清华杨顶辉数学实验十二作业

  • 格式:pdf
  • 大小:129.49 KB
  • 文档页数:5

第6题

(1)

h=[172169169171167171168165169168166168164

170165160175173172168155176172169176

173168169167170166161173175158170169173

164165167171166166172173178163169169178

177170167169173170160179172163173165176

162165172177182175170170169186174163172

176166167172177177169166182176172173174171

175165169168177184166171170];

w=[755564654762675262656265

5958645567746457576469585758

50527257554957765163636159

62536170635360645754666066

56545873586562504767586352

66596669756062637766505960

76635758677250636856596459

685665626470497159];

>>hist(h);得到身高直方图如下:

>>hist(w);得到体重直方图如下:

>>jbtest(h)

ans=

0

>>jbtest(w)

ans=

0

【结论】接受总体服从正态分布的假设

(2)>>[musigmamucisigmaci]=normfit(h)

mu=

170.2500

sigma=

5.4018

muci=

169.1782

171.3218

sigmaci=

4.7428

6.2751

>>[musigmamucisigmaci]=normfit(w)

mu=

61.2700

sigma=

6.8929

muci=

59.9023

62.6377

sigmaci=

6.0520

8.0073

【结论】身高/cm体重/kg

均值点估计170.250061.2700

均值区间估计[169.1782,171.3218][59.9023,62.6377]

标准差估计5.40186.8929标准差区间估计[4.7428,6.2751][6.0520,8.0073]

(3)>>ttest(h,167.5)

ans=

1

>>ttest(w,60.2)

ans=

0

所以身高拒绝原假设,体重接受原假设。

【结论】身高发生明显变化,体重没有发生明显变化。

【心得】hist命令可以用来画直方图

第8题

通过计算出x1的均值x1a,进行单个总体均值的检验,总体方差未知。原假设H0:x2所在

的总体的均值mu=x1a,备择假设H1:x2所在总体的均值mu≠x1a

>>x1=[9385799078768185886892738884907069838385];

>>x2=[8889868587887593887886868089857978888890];

>>x1a=mean(x1);

>>ttest(x2,x1a)

ans=

1

所以拒绝原假设,即x2所在总体的均值mu不等于x1a

【结论】两次测验难度有差异

第10题

简写中国为c,日本为j,男生为b,女生为g,均值为a,标准差为s

>>T=[7124.55.7122.55.4123.45.4121.85.4

8129.45.6128.15.5128.45.5127.65.7

9134.66.0133.45.4134.36.2133.56.3

10139.36.6138.95.9140.06.9140.26.6

11145.17.2144.96.7146.77.0146.76.7

12151.28.1152.07.8152.56.6151.96.2

13160.08.0159.67.6156.36.0155.15.4

14165.17.0165.16.8157.75.5156.75.2

15168.36.3168.56.2158.95.6157.45.0

16170.16.3170.05.9159.35.4157.95.3

17171.06.0170.86.0159.35.4158.15.0

18170.85.8171.15.9159.15.3158.25.1];

>>cba=T(:,2);cbs=T(:,3);jba=T(:,4);jbs=T(:,5);

>>cga=T(:,6);cgs=T(:,7);jga=T(:,8);jgs=T(:,9);

假设中国、日本每个年龄段的样本容量都是5000,不同样本容量结果可能不同。取5000是

为了便于和课后答案比对。原假设H0:中日学生在该年龄段没有差异。H1:中日学生在该年

龄段有差异。把这个问题看做两个总体均值的检验,且总体方差未知,假设2总体方差相等。

由于系统的ttest2检验需要知道原始数据,但是没有。所以自己写一个t检验的函数f10

functionh=f10(a1,a2,s1,s2,n1,n2)

temp=((n1-1).*s1.^2+(n2-1).*s2.^2)/(n1+n2-2);

t=(a1-a2)./sqrt(temp./n1+temp./n2);

alpha=0.05;

fori=1:1:12

ifabs(t(i))<=tinv(1-alpha/2,n1+n2-2);

h(i,1)=0;

else

h(i,1)=1;

end

end

>>b=f10(cba,jba,cbs,jbs,5000,5000);

>>g=f10(cga,jga,cgs,jgs,5000,5000);

>>[T(:,1),b,g]

ans=

711

811

911

1010

1100

1211

1311

1401

1501

1601

1701

1811

【结论】上表中第1列表示不同的年龄段,第2列表示男生比较结果,第3列表示女生比较

结果。其中0表示接受原假设,即对应年龄段学生身高没有差异,1表示否定原假设,即对

应年龄段学生身高有差异。

【心得】大概看了一下,好像中国男生、女生身高比日本高一点。本来以为日本营养好、孩

子运动更多,身高会高一些,看来不是这样啊。民族自豪感稍微有了一些。