نقل قول:
'admin' pid='8524' dateline='1418740112'
لطفا کد جدید را قرار بدهید تا ببینم تغییراتی که دادید، چی بوده.
این کد جدیده.
من باید انترگرال تی رو حتما حساب کنم .
کد پیاچپی:
clc%calculate switching probability+Tact with T-distribution
clear all;
close all;
Nr=100;%intial value
sigma=Nr;
mu=0;
D=trnd(5,1,1);%generate random numbers in T distribution
b=0;
c=0;
q=0;
g=0;
z=0;
s=0;
for alfa=1:8
for landa=1:8
if D<-landa%dar vaziyat foroosh
%syms x
f=tpdf(-5:5,5)
b=int(f,-Inf,-landa)
elseif D>landa%dar vaziyat kharid
%syms x
j=tpdf(-5:5,5)
c=int(j,landa,Inf)
end
inact(alfa,landa)=b+c
end
for landa=1:8
if -landa<D<0
%syms x
o=tpdf(-5:5,5)
z=1/2*int(o,-landa-alfa,landa-alfa)
elseif 0<D<landa
%syms x
u=tpdf(-5:5,5)
s=1/2*int(u,-landa+alfa,landa+alfa)
end
act(alfa,landa)=z+s
end
end
for alfa=1:8
for landa=1:8
h(alfa,landa)=act(alfa,landa)+inact(alfa,landa);
Tact(alfa,landa)=inact(alfa,landa)/h(alfa,landa);
Tinact(alfa,landa)=act(alfa,landa)/h(alfa,landa);
k(alfa,landa)=((alfa^4+6*alfa^2+3)*Tact(alfa,landa)+3*Tinact(alfa,landa))/((alfa^2+1)*Tact(alfa,landa)+Tinact(alfa,landa))^2;
end
end
[x,y]=meshgrid(1:8)
contourf(x,y,Tact,15)
colorbar
xlabel('alfa');
ylabel('landa');
title('Tact')