假设检验 MATLAB编程

假设检验MATLAB编程

function Pow = Example61a(n)

Mu = 119:.05:121;

M = 1e4;

Pow = zeros(length(Mu),1);

for s = 1:length(Mu)

mu = Mu(s);

Flag = zeros(M,1);

for m = 1:M

X =normrnd(mu,ones(n,1));

Tx = (mean(X) - 120)/(std(X)/sqrt(n)); Flag(m) = 1 - (Tx>= -1.96)*(Tx<= 1.96); end

Pow(s) = mean(Flag);

end

%plot(Mu,Pow);