function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b);

%??計(jì)算理論的二維levy分布的圖象? ?畫出來的3維曲面圖

%? ?t0 ：levy walk??趨于零的時(shí)間

%??lambda ：冪律值

%? ?T：時(shí)間

%? ?a: 離散的大步長的最大值

%? ?b：小步長的最大值

error=0.125/180*pi;??%假定平行于坐標(biāo)軸的誤差角度

syms u

c=-gamma(2-lambda)*t0^(lambda-1)*cos(pi*lambda/2)/lambda;

h1=20;??%大的離散步長

h2=0.1; %小的離散步長

x1=-a:h1:-b;

x2=-(b-h2):h2b-h2);

x3=b:h1:a;

x=[x1,x2,x3];

disp('%%%%%%%%%%%計(jì)算的次數(shù)為%%%%%%%%%%%%%%%')

N=length(x)

y=x;

Z=zeros(N,N);

for i=1:1? ?%x坐標(biāo)

i

for j=1:N??% y坐標(biāo)

theta=angle(complex(x(i),y(j)));

if theta<error && theta>-error??% x軸正半軸

Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));

elseif? ???theta<pi/2+error && theta>pi/2-error??% y軸正半軸

Z(j,i)=double(int(cos(u.*y(j)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));

elseif? ?theta>pi-error || theta<-pi+error??%x軸負(fù)半軸

Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));

elseif? ?theta>-pi/2-error && theta<-pi/2+error? ?%y軸負(fù)半軸

Z(j,i)=double(int(cos(u.*y(j)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));

else

Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1))*(double(int(cos(u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1)));

%zy=(double(int(cos(u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1)));

%Z(j,i)=zx*zy;

end

end

end

Y=Z;

end

錯(cuò)誤提示；

Error using mupadmex

Error in MuPAD command: DOUBLE cannot convert the input expression into a

double array.

If the input expression contains a symbolic variable, use the VPA function

instead.

Error in sym/double (line 710)

Xstr = mupadmex('symobj::double', S.s, 0);

Error in levydistrintheorytest (line 37)

Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1))*(double(int(cos(u.*y(

總結(jié)

以上是生活随笔為你收集整理的matlab 打开软件报错,matlab程序错误，提示如下【弄了半夜还是不行】 - 数学 - 小木虫 - 学术 科研 互动社区...的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。