untitled2.m

来自「产生正态变量的随机数,进行未知参数的区间估计 产生两个正态总体的随机数, 进行」· M 代码 · 共 48 行

本题两问，第一问用SPSS实现，第二问用MATLAB实现,以下是第二问，产生了两个正态分布的
1000个随机数，均值差用正态分布进行区间估计（由于均值差服从正态分布），方差比用f分布进行区间估计，
其中u=icdf('f',0.975,9999,9999);表示f分布的0.975下分位数。
在指令窗口运行如下程序
y=random('normal',10,2,1,10000);
z=random('normal',6,3,1,10000);
x=0;s=0;m=0;n=0;
for i=1:10000
  x=x+y(i);
end
x=x/10000;
for i=1:10000
    s=s+z(i);
end
s=s/10000;
for i=1:10000
    m=m+(y(i)-x)^2;
end
m=m/9999;
for i=1:10000
    n=n+(z(i)-s)^2;
end
n=n/9999;
u=icdf('f',0.975,9999,9999);
p=icdf('f',0.025,9999,9999);
t=icdf('normal',0.975,4,0.0361);
lower=x-s-t*sqrt(m/9999+n/9999)（表示均值差区间估计的下界）
uper=x-s+t*sqrt(m/9999+n/9999)（表示均值差区间估计的上界）
low=m/(n*u)（表示方差比区间估计的下界）
up=m/(n*p)（表示方差比区间估计的上界）
lower =

    3.8384


uper =

    4.1331


low =

    0.4307


up =

    0.4658