实验十三.数理统计实验

最新推荐文章于 2024-07-15 16:27:13 发布

智多星736

最新推荐文章于 2024-07-15 16:27:13 发布

阅读量546

点赞数 1

分类专栏： matlab实验报告文章标签： matlab

本文链接：https://blog.csdn.net/qq_67421354/article/details/127684725

版权

matlab实验报告专栏收录该内容

15 篇文章 9 订阅

订阅专栏

实验内容

1.教材（MATLAB与数学实验，主编：刘二根、王广超、朱旭生。）P140,6.5练习1-6小题。要求：每一题均需写出完整的MATLAB程序

2.自动包装机每额定质量是50kg,开工后随机抽查12袋,质量为:50.2 49.5 49.9 50.1 50.2 50.0 49.2 49.8 50.1 49.6 49.7 50.2。若每袋质量服从正态分布,问包装机工作是否正常(α=0.05)?

3.随机调查15名女生和18名男生的某科成绩,如下,

女生:91 77 79 77 57 82 68 93 71 84 81 70 53 95 69

男生:80 78 99 80 55 76 48 68 67 68 62 90 31 77 86 83 88 69

问是否可以认为(1)男女生成绩相当?(2)女生成绩优于男生(α=0.05)?

4.某产品质量指标如下:4177 4069 3244 3123 3204 2025 3266 3618 3815 4944 3149 4124 4299 4441 2508 3712 3738 2492 2758 4582。问该产品的质量指标是否服从正态分布?

5.随机从A批导线中抽取4根、从B批中抽取5根,测得电阻如下:A：0.143 0.142 0.143 0.137 B：0.140 0.142 0.136 0.138 0.140; 设A、B批导线的电阻服从同方差的正态分布,且两样本相互独立,则求μ₁-μ₂的95%的置信区间。

6.对两个总体抽样,A：111 92 94 108 87 96 85 88 108 112 100 107 108 99 85 108 109 97 91 107;

B：40 125 106 160 56 130 95 137 146 86 99 119 230 153 69 42 124 132 137 95。问A、B两个样本对应的总体分布类型是否一致?

实验步骤、过程

第一题

在命令行输入

>> y1=pdf('exp',3,2)  %第(1)题
>> y2=pdf('f',1.58,11,7)  %第(2)题

第二题

在命令行输入

>> a1=cdf('norm',6,5,3);
>> a2=cdf('norm',-6,5,3);
>> y1=(1-a1)+a2;
>> y1

第三题

在命令行输入

>> y1=norminv(0.95)
>> y2=tinv(0.025,17)
>> y3=chi2inv(0.975,15)
>> y4=finv(0.95,7,11)

第四题

由X~N(μ,σ2)知E(X)=μ，D(X)=σ2。从而

α1=E(X)=μ，α2=E(X2)=D(X)+E2(X)=σ2+μ2

从而得μ=A1，σ2+μ2=A2

使用MATLAB求解上述方程组：

>> syms x y A1 A2
>> f1=x-A1;
>> f2=x^2+y^2-A2;
>> [p,q]=solve(f1,f2,'x','y')

第五题

>> syms x a b n
>> fx=x*a*exp(-n*x)*b;
>> lnf=log(fx);
>> dlnf=diff(lnf,x);
>> solve(dlnf,x)

第六题

>> x1=unifrnd(2,7,[1,100]);
>> [a,b]=unifit(x1,0.1)

>> x2=poissrnd(2,[1,100]);
>> [a,b]=poissfit(x2)

>> x3=normrnd(2,5,[1,100]);
>> [a,b,c,d]=normfit(x3)

>> x4=binornd(100,0.55,[1,100]);
>> [a,b]=binofit(x4,100,0.01)

>> x=[50.2 49.5 49.9 50.1 50.2 50.0 49.2 49.8 50.1 49.6 49.7 50.2];
>> n=length(x)
>> mx=mean(x);
>> sx=std(x);
>> (mx-50)/(sx/sqrt(n))
>> talfa=tinv(0.95,12)

>> x=[91 77 79 77 57 82 68 93 71 84 81 70 53 95 69];
>> y=[80 78 99 80 55 76 48 68 67 68 62 90 31 77 86 83 88 69];
>> n1=length(x);
>> n2=length(y);
>> mx=mean(x);
>> my=mean(y);
>> z=(mx=my)/(3*sqrt(1/n1+1/n2))
>> zalfa=norminv(0.95)

>> x=[4177 4069 3244 3123 3204 2025 3266 3618 3815 4944 3149 4124 4299 4441 2508 3712 3738 2492 2758 4582];
>> [h,p,ci]=ttest(x,3665,0.05,1)

>> A=[0.143 0.142 0.143 0.137];
>> B=[0.140 0.142 0.136 0.138 0.140];
>> m1=mean(A);
>> m2=mean(B);
>> lc=(m1-m2)*(1/finv(0.95,3,4))
>> uc=(m1-m2)*(1/finv(0.05,3,4))

>> A=[111 92 94 108 87 96 85 88 108 112 100 107 108 99 85 108 109 97 91 107];
>> B=[40 125 106 160 56 130 95 137 146 86 99 119 230 153 69 42 124 132 137 95];
>> [p,h,stats]=ranksum(A,B,0.05)

答案供参考,如有错误请自行改正，本人只是随意分享之前的作业，代码均来源于网络+个人稍微修改，认为有用可参考哈～