newjifen.m

来自「生成验证平稳分布用matlab」· M 代码 · 共 11 行

function sum=newjifen(tal)
sum=0;
pha=2*pi*rand(1);
 a=rand(1)+4.5;
 t=-50:0.01:50;
 for i=1:10001
xt(i)=a*cos(2*pi*t(i)+pha)*a*cos(2*pi*(t(i)+tal)+pha);
end;
 for z=1:10001
sum=(xt(z)*0.01)/100+sum;
end