Matlab随笔之判别分析

从概率论角度,判别分析是根据所给样本数据,对所给的未分类数据进行分类。

如下表,已知有t个样本数据,每个数据关于n个量化特征有一个值,又已知该样本数据的分类,据此,求s个未分类数据的分类情况class。

Matlab 的统计工具箱提供了判别函数

[class,err] = classify(sample,training,group, type)

其中,err 给出了分类误判率的估计值,type为分类方法,缺省值为‘linear‘,即线性分类,type 还可取值‘quadratic‘, ‘mahalanobis‘( mahalanobis 距离)。

eg:

程序:

training=[13.54,14.36,87.46,566.3,0.09779
13.08,15.71,85.63,520,0.1075
9.504,12.44,60.34,273.9,0.1024
17.99,10.38,122.8,1001,0.1184
20.57,17.77,132.9,1326,0.08474
19.69,21.25,130,1203,0.1096
11.42,20.38,77.58,386.1,0.1425
20.29,14.34,135.1,1297,0.1003];
sample=[16.6,28.08,108.3,858.1,0.08455
20.6,29.33,140.1,1265,0.1178
7.76,24.54,47.92,181,0.05263];
group=[zeros(3,1);ones(5,1)];%0表示良性肿瘤,1表示恶性肿瘤
[class,err]=classify(sample,training,group)

运行结果:

class =

     0
     0
     0

err =

     0

  

所以,三个样本的分类都为良性肿瘤。

时间: 2024-10-13 00:11:00

Matlab随笔之判别分析的相关文章

Matlab随笔之线性规划

原文:Matlab随笔之线性规划   LP(Linear programming,线性规划)是一种优化方法,在优化问题中目标函数和约束函数均为向量变量的线性函数,LP问题可描述为:min xs.t. A·x b Aeq·x=beq vlb x vub其中 ,b,beq均为向量,A,Aeq为矩阵,x为向量变量.矩阵A和向量b是线性不等式约束条件的系数, Aeq和beq是等式约束条件的系数. 在MATLAB中,用于LP的求解函数为linprog.其调用格式为:[x,fval,lambda]=linp

Matlab随笔之模拟退火算法

问题描述: 我方有一个基地,经度和纬度为( 70,40).假设我方飞机的速度为 1000 公里/小时. 我方派一架飞机从基地出发,侦察完敌方所有目标,再返回原来的基地.在敌方每一目 标点的侦察时间不计,求该架飞机所花费的时间(假设我方飞机巡航时间可以充分长). 这是一个旅行商问题.我们依次给基地编号为 1,敌方目标依次编号为 2, 3,…, 101, 最后我方基地再重复编号为 102(这样便于程序中计算). 距离矩阵 D = ( dij )102×102 , 其中 dij 表示表示 i, j 两

Matlab随笔之矩阵入门知识

直接输入法创建矩阵 – 矩阵的所有元素必须放在方括号“[ ]”内: – 矩阵列元素之间必须用逗号“,”或空格隔开,每行必须用“;”隔开 – 矩阵元素可以是任何不含未定义变量的表达式.可以是实数,或者是复数. – 例a=[1,2;3,4] 或 a=[2 1+3j;sqrt(4) 5] 创建基本矩阵的函数 – 空阵 [ ] — matlab允许输入空阵,当一项操作无结果时,返回空阵 – ones(N,M) —全部元素都为1的矩阵 – zeros(N,M) —全部元素都为0的矩阵 – rand(N,M

Matlab随笔(3)之指派问题的整数规划

注:除了指派问题外,一般的整数规划问题无法直接利用Matlab函数,必须Matlab编程实现分支定界法和割平面解法. 常用Lingo等专用软件求解整数规划问题. %指派矩阵为 % 3, 8, 2,10, 3 % 8, 7, 2, 9, 7 % 6, 4, 2, 7, 5 % 8, 4, 2, 3, 5 % 9,10, 6, 9,10 % %解题过程 c=[3 8 2 10 3; 8 7 2 9 7; 6 4 2 7 5; 8 4 2 3 5; 9 10 6 9 10];%指派矩阵 c=c(:);

Matlab随笔(1)之画图函数总结

MATLAB函数画图 MATLAB不但擅长於矩阵相关的数值运算,也适合用在各种科学目视表示(Scientific visualization).本节将介绍MATLAB基本xy平面及xyz空间的各项绘图命令,包含一维曲线及二维曲面的绘制.列印及存档. 基本XY平面绘图命令 Plot 是绘制一维曲线的基本函数,但在使用此函数之前,我们需先定义曲线上每一点的x及y座标.下例可画出一条正弦曲线: close all; %linspace(5,100,20)和5:5:100的区别在于前者已知元素总个数而不

Matlab随笔之插值与拟合(上)

1.拉格朗日插值 新建如下函数: function y=lagrange(x0,y0,x) %拉格朗日插值函数 %n 个节点数据以数组 x0, y0 输入(注意 Matlat 的数组下标从1开始), %m 个插值点以数组 x 输入,输出数组 y 为 m 个插值 n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end

Matlab随笔之三维图形绘制

1.三维曲线 用到了plot3函数.plot3(x,y,z)用来绘制3维曲线图,而不能绘制曲面图!就是把所有的(x,y,z)点连接在一起. t=linspace(0,10*pi,2000); x=sin(t); y=cos(t); z=t; plot3(x,y,z) 2.网格曲面 meshgrid是MATLAB中用于生成网格采样点的函数: mesh(X,Y,Z)函数功能是绘制由线条框构成的曲面. x=-3:0.1:3;y=-5:0.1:5; x1=ones(size(y'))*x;y1=y'*o

matlab随笔1

1.查看.mat文件:load('文件名')即可,其中有这样的语句: load('ex3data1.mat'); % training data stored in arrays X, y m = size(X, 1); 这其中在文件中已经定义了X,所以在程序的内容中X已经有了数据 2.关于~isempty(XXX)的使用: matlab中 ~表示逻辑运算符"非",也就是取反:~isempty(msg) 表示如果msg是空元素(空元素表示未赋值的元素,0并不是空元素),结果为0(fal

Matlab随笔(4)之分段线性函数化为线性规划

eg:         10x,            0<=x<=500 c(x)=1000+8x,    500<=x<=1000         3000+6x,    1000<=x<=1500   解法一: 可引入0-1变量,令z1=1,z2=1,z3=1分别表示0<=x<=500,500<=x<=1000,1000<=x<=1500,则 500z2<=x1<=500z1, 500z3<=x2<=50