四川定制网站建设,开发网站服务器,网络推广外包怎么样,网站推广要点 优帮云function [x,y,Y]levydistrintheorytest(t0,lambda,T,a,b);% 计算理论的二维levy分布的图象 画出来的3维曲面图% t0 #xff1a;levy walk 趋于零的时间% lambda #xff1a;幂律值% T#xff1a;时间% a: 离散的大步长的最大值% b#xff1a;小步长的最大值e…function [x,y,Y]levydistrintheorytest(t0,lambda,T,a,b);% 计算理论的二维levy分布的图象 画出来的3维曲面图% t0 levy walk 趋于零的时间% lambda 幂律值% T时间% a: 离散的大步长的最大值% b小步长的最大值error0.125/180*pi; %假定平行于坐标轴的误差角度syms uc-gamma(2-lambda)*t0^(lambda-1)*cos(pi*lambda/2)/lambda;h120; %大的离散步长h20.1; %小的离散步长x1-a:h1:-b;x2-(b-h2):h2b-h2);x3b:h1:a;x[x1,x2,x3];disp(%%%%%%%%%%%计算的次数为%%%%%%%%%%%%%%%)Nlength(x)yx;Zzeros(N,N);for i1:1 %x坐标ifor j1:N % y坐标thetaangle(complex(x(i),y(j)));if thetaerror theta-error % x轴正半轴Z(j,i)double(int(cos(u.*x(i)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./picos(1/u.*x(i)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));elseif thetapi/2error thetapi/2-error % y轴正半轴Z(j,i)double(int(cos(u.*y(j)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./picos(1/u.*y(j)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));elseif thetapi-error || theta-pierror %x轴负半轴Z(j,i)double(int(cos(u.*x(i)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./picos(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/2error %y轴负半轴Z(j,i)double(int(cos(u.*y(j)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./picos(1/u.*y(j)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1));elseZ(j,i)double(int(cos(u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-(abs(u))^lambda)./picos(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)./picos(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)./picos(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;endendendYZ;end错误提示Error using mupadmexError in MuPAD command: DOUBLE cannot convert the input expression into adouble array.If the input expression contains a symbolic variable, use the VPA functioninstead.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)./picos(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(
