原创力文档- 1、本文档共3页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数学实验第八章习题答案
%ex8_1
clear;
data=[459, 362, 624, 542, 509, 584, 433, 748, 815, 505, 612, 452, 434, 982, ...
640, 782, 742, 565, 706, 593, 680, 926, 653, 164, 487, 734, 608, 428, 1153, ...
593, 844, 527, 552, 513, 781, 474, 388, 824, 538, 862, 659, 775, 859, 755, ...
649, 697, 515, 628, 954, 771, 609, 2, 960, 885, 610, 292, 837, 473, 677, 358, ...
638, 699, 634, 555, 570, 84, 416, 606, 1062, 484, 120, 447, 654, 564, 339, ...
280, 246, 687, 539, 790, 581, 621, 724, 531, 512, 577, 496, 468, 499, 544, ...
645, 764, 558, 378, 765, 666, 763, 217, 715, 310, 851];
s=sort(data);
s([1:3,end-2:end])=[];%删除上下各3个值
hist(s,10) %直方图基本符合正态分布
mean(s),std(s)
%ex8_2
clear;
x=[480 500 380 1100 1100 490 230 250 300 510 1300];
y=[180 150 170 350 460 240 60 90 110 250 200];
plot(x,y,o) %散点图可看出是
corrcoef(x,y)%系数0.7373
%ex8_3
%c, f, b, e, a, d
%ex8_4
%用1表示正面,2表示反面
N=100000;R=unidrnd(2,1,N);P=sum(R==1)/N
%ex8_5
function s=ex8_5(n,p)
r=rand(1,n);s=0;
for i=1:n
if r(1,i)p
s=s+1;
end
end
%例如:用蒙特卡罗法求B(10,0.2)的数学期望和方差
N=100000;R=zeros(1,N);
for i=1:N
R(i)=ex8_5(10,0.2);
end
E=mean(R)
V=std(R)^2
%可见E接近10*0.2, V接近10*0.2*0.8
%ex8_6
N=100000;n=10;p=0.3;
R=binornd(n,p,1,N);
R=(R-n*p)/sqrt(n*p*(1-p));
x=linspace(-4,4,20);
c=zeros(1,20);
for i=1:20
c(i)=sum(R=x(i))/N;
end
plot(x,c,x,normcdf(x,0,1),r:)
legend(二项分布函数,正态分布函数)
%ex8_7
%(1)
clear;N=10000;
x=unifrnd(0,1,1,N);y=unifrnd(0,1,1,N);
Nc=sum(yexp(-x.^2/2));
I=Nc/N/sqrt(2*pi)
%(2)
clear;N=100000;
x=unifrnd(0,2*pi,1,N);y=unifrnd(0,exp(pi),1,N);
Nc=sum(y(exp(x/2).*sin(x).^2));
I=Nc/N*2*pi*exp(pi)
%(3)
clear;N=10000;
x=unifrnd(0,pi,N,1);y=unifrnd(0,1,N,1);z=unifrnd(0,1,N,1);
c1=ysin(x);c2=zexp(-x.^2-y.^2);
Nc=sum(c1c2);I=Nc/N*pi
%用int验算
syms x y;I=int(exp(-x^2-y^2),y,0,sin(x));
I2=int(I,x,0,pi);
vpa(I2)
%ex8_8
clear;N=1000;
x=unifrnd(-1.5,1.5,N,1);y=unifrnd(-1.5,1.5,N,1);
[x,y]=meshgrid(x,y);
f=(x.^2+2*y.^2+x.*y).*exp(-x.^2-y.^2);
mesh(x,y,f)
[fmax,id]=max(f(:))
xmax=x(id),ymax=y(id)
%用fminsearch验算
fun=inline(-(x(1)^2+2*x(2)^2+x(
文档评论（0）