数理统计应用概要1.pptVIP

数理统计应用 1.解题思路： 运用回归系数的最小二乘估计。 2.解题过程： （1）.首先，将两种蠓虫的触角长度与翅膀长度以（xi,yi）为坐标，在平面直角坐标系中画出散点图。 散点图程序： %散点图 clear clc Y1=[1.14 1.18 1.20 1.26 1.28 1.30];%apt触角长度 Y2=[1.78 1.96 1.86 2.00 1.96 2.00];%apt翅膀长度 axis([1 2.1 1.6 2.1]); scatter(Y1,Y2,k); hold on; Y3=[1.30 1.24 1.36 1.38 1.38 1.38 1.40 1.49 1.54];%af触角长度 Y4=[1.72 1.74 1.64 1.82 1.90 1.70 1.82 1.82 2.08];%af翅膀长度 axis([1 2.1 1.6 2.1]); scatter(Y3,Y4,r); （2）.设出线性回归方程并计算出其中的参数。 线性回归方程为：y1=a1+b1x1; 计算apt的触角长度Y1与翅膀长度Y2的均值。 ①.E(Y1)=(1/n1)∑(Y1i)=1.2267;(n1=6) ②.E(Y2)=(1/n2)∑(Y2i)=1.9267;(n2=6) 估计a1与b1的点值： ①.Z1=∑(Y1i-E(Y1))^2=0.0197; ②.V1=∑(Y1i-E(Y1))*(Y2i-E(Y2))=0.0217; ③.b1=Z1/V1=1.1014; ④.a1=E(Y1)-E(Y2)*b1=0.5757; (3).设出线性回归方程并计算出其中的参数. 线性回归方程为：y2=a2+b2x2; 计算af的触角长度Y3与翅膀长度Y4的均值。 ①.E(Y3)=(1/n3)∑(Y3i)=1.3856;(n3=10) ②.E(Y4)=(1/n4)∑(Y4i)=1.8044;(n4=10) 估计a2与b2的点值： ①.Z2=∑(Y3i-E(Y3))^2=0.0642; ②.V2=∑(Y3i-E(Y3))*(Y4i-E(Y4))=0.0652; ③.b2=Z2/V2=1.0149; ④.a2=E(Y3)-E(Y4)*b2=0.3983; (4).线性回归方程的图像： (5).对(1.24,1.80)(1.28,1.84)(1.40,2.04)三个样本进行辨别。 偏离线性回归方程的误差： e1=(y1-x1*a1-b1)^2;(apt) e2=(y2-x2*a2-b2)^2;(af) 将上面数据分别带入，并比较e1与e2大小，数值越小越接近回归方程。 通过判断得出结论： apt,af,apt. (6).程序： %回归系数的最小二乘估计 clear clc Y1=[1.14 1.18 1.20 1.26 1.28]; Y2=[1.78 1.96 1.86 2.00 1.96]; c1=mean(Y1); c2=mean(Y2); n1=5; for i=1:5 z1(i)=(Y1(i)-c1)^2; v1(i)=(Y1(i)-c1)*(Y2(i)-c2); end Z1=sum(z1); V1=sum(v1); b1=V1/Z1; a1=c2-b1*c1; disp(a1); disp(b1); x1=0:0.01:2.1; y1=b1*x1+a1; plot(Y1,Y2,o); hold on; plot(x1,y1); hold on; %2 Y3=[1.30 1.24 1.36 1.38 1.38 1.38 1.40 1.49 1.54 1.56]; Y4=[1.72 1.74 1.64 1.82 1.90 1.70 1.82 1.82 2.08 1.96]; c3=mean(Y3); c4=mean(Y4); n2=10; for i=1:10 z2(i)=(Y3(i)-c3)^2; v2(i)=(Y3(i)-c3)*(Y4(i)-c4); end Z2=sum(z2); V2=sum(v2); b2=V2/Z2; a2=c4-b2*c3; disp(a2); disp(b2); x2=0:0.01:2.1; y2=b2*x2+a2; plot(Y3,Y4,*); plot(x2,y2); %辨别 Y5=[1.24 1.28 1.40]; Y6=[1.80 1.84 2.04]; for j=1:3 K1(j)=Y6(j)-b1*Y5(j)-a1; K2(j)=Y6(j)-b2*Y5(j)-a2; if(K1(j)^2K2(j)^2) O(j)=1;%af else O(j)=0;%apt end end disp(