多元分析实验项目

1、 多元统计分析实验报告实验序号：01班级： 姓名： 学号：实验名称正态分布随机向量的生成、可视化多元正态分布函数实验目的：1、通过实践，掌握仿真数据的生成实验内容：1、 通过matlab产生正态分布随机数，利用产生的随机数，进一步产生正态分布的随机向量。 2、计算产生的正态分布随机向量的均值及协方差阵，作出相应的概率密度函数图像 实验要求：1、独立完成上述实验内容。2、要求有计算程序和结果。实验示例%multi_variableanalysis_1.m x=normrnd(0,1,1,5000);y=normrnd(0,1,1,5000);x=x;y;figure(1)plot(x(1,:),

2、x(2,:),'r.'); %d(ax)=a'*d(x)*a;%a=eye(2);a(2,2)=2;xx=a*x;figure(2)plot(xx(1,:),xx(2,:),'m.'); a=cos(0.5),sin(0.5);-sin(0.5),cos(0.5);a=a*1,0;0,3;%a(2,2)=5;xx=a*x;figure(3)plot(xx(1,:),xx(2,:),'.');s,v=eig(a*a');%- covmatrix=inv(a*a');x0=-15:0.5:15;y0=15:-0.5:-15;x

3、0,y0=meshgrid(x0,y0);z=zeros(size(x0)for i=1:size(x0,1) for j=1:size(y0,2); currv=x0(j,i),y0(j,i)' scale=-0.5*currv'*covmatrix*currv; z(i,j)=exp(scale); endendfigure(4)imagesc(z)数学实验报告实验序号：02班级： 姓名： 学号：实验名称样本均值、样本协差阵、样本相关阵的计算实验目的：1、通过实践，掌握多元统计的基本概念。实验内容： 1、仿真方法产生多元正态分布的样本，计算其样本均值、样本协差阵、样本相关阵

4、。实验要求：1、独立完成上述实验内容。2、要求有计算程序和结果。实验示例%multi_variableanalysis_2.m % simulated dataset x=normrnd(0,1,1,5000);y=normrnd(0,1,1,5000);x=x;y; a=cos(0.5),sin(0.5);-sin(0.5),cos(0.5)a=a*1,0;0,2;%a(2,2)=5;xx=a*x;covmatrix=a*a'meanv=zeros(2,1); meanvs=mean(xx')'covmatrixs=zeros(2,2);for i=1:size(xx

5、,2) covmatrixs=covmatrixs+(xx(:,i)-meanvs)*(xx(:,i)-meanvs)'endcovmatrixs0= 1/size(xx,2)*covmatrixs;covmatrixs1= 1/(size(xx,2)-1)*covmatrixs;corrmatrixs=covmatrixs0;for i=1:size(corrmatrixs,1) corrmatrixs(i,:)=corrmatrixs(i,:)./sqrt(covmatrixs0(i,i)end for j=1:size(corrmatrixs,2) corrmatrixs(:,j

6、)=corrmatrixs(:,j)./sqrt(covmatrixs0(j,j);end corrmatrixs数学实验报告实验序号：03班级： 姓名： 学号：实验名称距离判别分析实验目的：1、通过实践，掌握距离判别分析实验内容： 1、仿真方法产生数据，并采用距离判别分析生成的数据。 2、完成书本习题4.10实验要求：1、独立完成上述实验内容。2、要求有计算程序和结果。实验示例数学实验报告实验序号：04班级： 姓名： 学号：实验名称费希尔判别分析实验目的：1、通过实践，掌握费希尔判别分析方法实验内容： 1、完成课本习题4.9, 4.10实验要求：1、独立完成上述实验内容。2、要求有计算程序和

7、结果。实验示例数学实验报告实验序号：05班级： 姓名： 学号：实验名称k-均值聚类分析实验目的：1、通过实践，掌握k-均值聚类分析方法实验内容： 1、仿真方法产生数据，验证k-均值聚类分析； 实验要求：1、独立完成上述实验内容。2、要求有计算程序和结果。实验示例%k_means.m%-simulated dataxx1=5*(rand(3,50)-0.5);xx2=5*(rand(3,50)+0.5);xx3=5*(rand(3,50)+1.5);plot3(xx1(1,:),xx1(2,:),xx1(3,:),'.')hold onplot3(xx2(1,:),xx2(2,:

8、),xx2(3,:),'.')plot3(xx3(1,:),xx3(2,:),xx3(3,:),'.')%-k=3;xx=xx1,xx2,xx3;length_xx=size(xx,2);aa=zeros(k,150);centerp=zeros(3,k);idexx=round(length_xx*rand(1,k)+0.5);idexx=min(idexx,length_xx);centerp=xx(:,idexx);for jjtk=1:10 aa=zeros(k,length_xx);for i=1:length_xx aa=xx(:,i); tt=ze

9、ros(1,k); for j=1:k tt(j)=norm(aa-centerp(:,j); end idexv,idextt=min(tt); aa(idextt,i)=1;endcenterp1=centerp;for ii=1:k idexaa=find(aa(ii,:)=1); centerp(:,ii)=mean(xx(:,idexaa)')'endaxxs=max(max(abs(centerp-centerp1);if axxs=0 breakendend%computing q_valueq_value=0;for jj=1:kidexx=find(aa(jj

10、,:)=1);ss=xx(:,idexx);aaq=0;for tt=1:size(ss,2) aaq=aaq+norm(ss(:,tt)-centerp(:,jj)2;endq_value=q_value+aaq/size(ss,2);endq_value figure() idexx=find(aa(1,:)=1); plot3(xx(1,idexx),xx(2,idexx),xx(3,idexx),'ro') hold on idexx=find(aa(2,:)=1); plot3(xx(1,idexx),xx(2,idexx),xx(3,idexx),'k*&#

11、39;) idexx=find(aa(3,:)=1); plot3(xx(1,idexx),xx(2,idexx),xx(3,idexx),'+') 数学实验报告实验序号：06班级： 姓名： 学号：实验名称主成分分析实验目的：1、通过实践，掌握主成分分析实验内容： 1、完成课本习题6.8实验要求：1、独立完成上述实验内容。2、要求有计算程序和结果。实验示例%pca.mdata=load('datap.txt');nx,ny=size(data);data_s=zeros(nx,ny);means=zeros(1,ny);stds=zeros(1,ny);for

12、i=1:ny means(i)=mean(data(:,i); stds(i)=std(data(:,i);end for i=1:ny ll=(data(:,i)-means(i)./stds(i); data_s(:,i)=ll;end r=zeros(ny,ny);for i=1:nx r=r+data_s(i,:)'*data_s(i,:);endr=r./(nx-1); s,v,d=svd(r);for i=1:ny if sum(s(:,i)<0 s(:,i)=-s(:,i); endend y=data_s*s;lamda=zeros(ny,1);for i=1:n

13、y lamda(i)=v(i,i);endlamda=lamda./sum(lamda);y*lamda数学实验报告实验序号：07班级： 姓名： 学号：实验名称综合实验分析1实验目的：1、通过实践，综合运用多元统计方法实验内容：两种蠓虫af和apf巳由生物学家wl grogan和w.w.wirth(1981年)根据它们的触角长度和翼长加以区分。现测得6只apf和9只af蠓虫的触角长度和翼长识别的依据是蠓虫的触角和翅膀的长度的数据（见附录一），现在要根据数据,制定一种方法,正确区分两类蠓虫；并且将建立的模型将用于已知触角长和翼长三个待定的样本进行识别。蠓虫的触角和翅膀的长度的数据：apf(g1)af(g2)序号触角长度翅膀长度序号触角长度翅膀长度11.141.7871.241.7221.181.9681.361.7431.201.8691.381.6441.262.00101.381.8251.282.00111.381.9061.301.96121.401