两样本T检验matlab程序-均值差置信区间matlab程序

作业一1.用 matlab 编写样本均值和方差的程序X=input（’请输入样本值： \n’）;％以矩阵形式输入%a=mean（X） %求均值%b=var(X) ％求方差％2.编写均值差的置信区间程序%两正态总体方差未知，但方差相等%X=input(’请输入正态总体 X 的样本值: \n'）； Y=input('请输入正态总体 Y 的样本值: \n’）；S1=var（X); ％求样本 X 的方差％S2=var（Y)； ％求样本 Y 的方差％n11=size（X）； n22=size（Y)；n1=n11（：,2）; n2=n22(:，2)；%计算样本 X，Y 的个数%Sw=sqrt（（(n1-1)＊S1+(n2-1)＊S2)/(n1+n2—2））； %求 Sw％a=mean（X）； %求样本 X 的均值％b=mean（Y）; %求样本 Y 的均值%L=input（’请输入置信度： \n’);alpha=1—L；t=tinv（1—alpha/2,n1+n2-2)； %置信度为 L 的 T 值％d1=(a—b）-t＊Sw.＊sqrt(1。/n1+1./n2); %置信区间下界％d2=(a—b)+t*Sw。*sqrt(1./n1+1。/n2)； %置信区间上界％c=［d1，d2］ %均值差置信区间%％两正态总体方差已知,为 sigma1，sigma2％X=input（’请输入正态总体 X 的样本值: \n’)； Y=input('请输入正态总体 Y 的样本值: \n’);a=mean(X); %求样本 X 的均值%b=mean（Y）; %求样本 Y 的均值%n11=size(X）； n22=size(Y)；n1=n11(：,2)； n2=n22(:，2）; %计算样本 X，Y 的个数％L=input('请输入置信度： \n’)；alpha=1-L;u= norminv(1—alpha/2）； ％置信区间为 1—alpha 的 U 值％sigma1=input(’请输入总体 X 的方差平方根：\n'）;sigma2=input('请输入总体 Y 的方差平方根：\n');d1=a—b-u＊sqrt(sigma1.^2/n1+sigma2.^2/n2）； %置信区间下界%d2=a-b+u*sqrt(sigma1。^2/n1+sigma2。^2/n2)； ％置信区间上界%c=[d1,d2] ％均值差置信区间%3.编写两样本 T 检验程序%两样本 t 检验法,总体方差未知但相等％X=input(’请输入正态总体 X 的样本值: \n'）; Y=input（’请输入正态总体 Y 的样本值: \n');S1=var(X）； %求样本 X 的方差%S2=var(Y); %求样本 Y 的方差％n11=size（X）; n22=size（Y)；n1=n11(：,2）; n2=n22（：，2);％计算样本 X，Y 的个数％Sw=sqrt（（(n1-1）*S1+（n2—1）＊S2)/（n1+n2-2））； %求 Sw%a=mean(X）; %求样本 X 的均值%b=mean（Y）； ％求样本 Y 的均值％L=input(’请输入置信度： \n’）；alpha=1-L;t=tinv（1—alpha/2,n1+n2-2）； ％置信度为 L 的 T 值% %% if (abs（a—b)〉Sw。*t*sqrt（1/n1+1/n2）) ％双侧检验，u1=u2% %％ if (（a—b）〉Sw。*t*sqrt(1/n1+1/n2)) ％单侧检验，u1—u2〈=0％ %％ if （（a—b)<—（Sw.*t＊sqrt（1/n1+1/n2））) ％单侧检验，u1—u2〉=0％ disp('拒绝假设’）else disp（'假设成立')end4.证明样本方差是总体方差的无偏估量量(详见概率课本 P159）