求解 模擬

Tc.normal

function test=tc normal(c,gamma,alpha,lambda)
%c=0.1;
%alpha=10;
%lamda=6;

n=1;
%gama=1;
k=1;
tc=0;
mu=0;
f1=1000;
f2=1000;
rtc_new=0;
while (tc==0);
n=n+1;
if (1-(gamma(alpha+n*k*0.5+gama)/gamma(alpha+n*k*0.5))*exp(gama)*((alpha+n*k*0.5)/(alpha+n*k*0.5+gama))^(alpha+n*k*0.5)/(alpha+n*k*0.5+gama)^gama<=c*n);
tc=n;
end;
end;


for t1=1:f1;
theta=gamrnd(alpha,1/lamda);
for t2=1:f2;
xn=0;   
xn=sum(((normrnd(mu,1/sqrt(theta),1,tc)-mu).^2)/2);
d(t1,t2)=(alpha+tc*k/2)/(lamda+xn);

end;
rtc_new=rtc_new+mean(1-(theta.*(exp(1-theta./d(t1,:)))./d(t1,:)).^gama+c*tc);
end;
RTC=rtc_new/f1;
RTC1=(RTC-2*sqrt(alpha)*sqrt(c))/(2*sqrt(alpha)*sqrt(c));
test(1)=tc;
test(2)=RTC;
test(3)=RTC1;


哪裡有錯 請指教