统计计算作业new

《统计计算作业new》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、统计计算作业一课上作业1.EM算法问题：2.bootstrap算法问题：利用bootstrap算法估计正态分布样本方差的方差。算法：（1）产生500个标准正态分布样本；（2）从已知样本分布中有放回的抽取B组大小为n的样本；（3）用生成的新样本估计样本方差的方差。程序：%%%%用bootstrap方法估计样本方差的估计值n=500;x=normrnd(0,1,1,n);%产生正太随机数fori=1:1000;%B=1000y=length(x);z1=ceil(n*rand(1,n));d=x(z1);%在产生的样本中重复抽样，抽取50

2、个样本s1=(1/n)*(d*d.');s2=(1/n)*sum(d);s(i)=s1-(s2)^2;%每一组用新产生的样本进行估计endm1=0.001*(s*s.');n1=(0.001*sum(s))^2;m1-n1%对1000组进行估计最后的s平方3.JICKNIFE算法问题：利用JICKNIFE算法估计正态分布样本方差的方差。算法：(1)程序：%%%%用JICKNIFE方法估计样本方差的估计值n=500;x=normrnd(0,1,1,n);%产生正太随机数fori=1:n;j=i;a=x*x.';b=x(j)^2;f(j)

3、=(1/(n-1))*(a-b);g(j)=(1/(n-1))*(sum(x)-x(j));t0=(1/n)*a-(sum(x)/n)^2;%s平方T(j)=f(j)-(g(j))^2;%JACKNIFE方法中去掉一个样本的统计量T1(j)=n*t0-(n-1)*T(j);%Tj星公式endJ=(1/n)*sum(T1);%JICKNIFEs平方的估计值h=T1*T1.';z1=(1/(n-1))*J^2;V=(1/(n*(n-1)))*h-z1%s平方方差的估计值T;V;4.遗传算法问题：利用遗传算法找出方程的最优解，其中，。算法：

4、程序：%%%遗传算法求方程的最优解clear;clc;x=1:7;z=ceil(7*rand(1,8));d=x(z);%m1=[d(1)d(2)];m2=[d(3)d(4)];m3=[d(5)d(6)];m4=[d(7)d(8)];n=1000;count=1;while(count

5、=f2/fsum;p3=f3/fsum;p4=f4/fsum;a1=dec2bin(d(1),3);%二进制三位编码a2=dec2bin(d(2),3);a3=dec2bin(d(3),3);a4=dec2bin(d(4),3);a5=dec2bin(d(5),3);a6=dec2bin(d(6),3);a7=dec2bin(d(7),3);a8=dec2bin(d(8),3);s1=strcat(a1,a2);%连接编码s2=strcat(a3,a4);s3=strcat(a5,a6);s4=strcat(a7,a8);pp1=p1

6、;%0~1区间概率pp2=pp1+p2;pp3=pp2+p3;pp4=pp3+p4;u=rand(1,4);%随机选择4个,代表Sifori=1:4;ifu(i)

7、ews4=s4;if(crossover<=0.8)%交叉1位点的前半部分fori=1:r(1)%前两个交叉%第一个串前半部分if(select(1)==1)news1(i)=s1(i);elseif(select(1)==2)news1(i)=s2(i);elseif(select(1)==3)news1(i)=s3(i);elsenews1(i)=s4(i);end%第二个串前半部分if(select(2)==1)news2(i)=s1(i);elseif(select(2)==2)news2(i)=s2(i);elseif(se

8、lect(2)==3)news2(i)=s3(i);elsenews2(i)=s4(i);endendfori=r(1)+1:6%第一个串前后部分if(select(2)==1)news1(i)=s1(i);elseif