Fisher最优求解算法Matlab代码.docVIP

Fisher最优求解算法Matlab代码 代码: %===========dim.m======================%函数:计算类的直径function y =dim(i,j,q)y=0;qb=0;for l=i:j? ? qb=qb+q(l,:);endqb=qb/(j-i+1);for l=i:j? ? y=y+(q(l,:)-qb)*(q(l,:)-qb);end%================主程序main.m======================function main()%读入文件q=load(data1.txt);%得到行列数[m,n]=size(q); 以下内容需要回复才能看到 复制内容到剪贴板 代码: fprintf(请输入一个0-%d之间的,m);K=input(分类数量:);if K1 || Km? ? fprintf(你的输入有误,超过了分类范围.\n);? ? return;end%得到直径序列for i=1:m? ? for j=i:m? ?? ???e(i,j)=dim(i,j,q);? ? endend%生成归类,K是类的数量,Kmclassify=ones(K);classify=classify(1,:);for k=1:K-1? ? %生成误差序列? ? for n=1:m? ?? ???err(n)=inf;? ? end? ? for i=1:m? ? %假设i是归类点,不重复? ? %将i放入分类,排序.? ?? ???if isempty(find(classify==i)) ? ?? ?? ???temp=classify;? ?? ?? ???temp(2)=i;? ?? ?? ???temp=sort(temp);? ?? ?? ???%得到该分类误差? ?? ?? ?? ?err(i)=0;? ?? ?? ?? ?for j=1:K-1? ?? ?? ?? ?? ? if temp(j)temp(j+1)? ?? ?? ?? ?? ???err(i)=err(i)+e(temp(j),temp(j+1)-1);? ?? ?? ?? ?? ? end? ?? ?? ?? ?end? ?? ?? ?? ?err(i)=err(i)+e(temp(K),m);? ?? ???end ? ? end? ? %找出最小的误差类? ? %注意,有可能出两个最小值的情况.我们取第一个最小的? ? tem=find(err==min(err));? ? classify(2)=tem(1);? ? classify=sort(classify);end%显示结果classify(K+1)=m+1;fprintf(分%d类结果:,K);for i=1:K? ? fprintf(();? ? for j=classify(i):classify(i+1)-1? ?? ???fprintf( %d ,j);? ? end? ? fprintf() );endfprintf(\n);