24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bluethky

木虫 (著名写手)

应助: 1 (幼儿园)
金币: 4460.8
散金: 849
红花: 11
帖子: 1354
在线: 227.3小时
虫号: 2609601
注册: 2013-08-21
专业: 物理学II

[求助] matlab程序错误，提示如下【弄了半夜还是不行】已有2人参与

function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b);
%  计算理论的二维levy分布的图象画出来的3维曲面图
% t0 ：levy walk  趋于零的时间
%  lambda ：幂律值
% T：时间
% a: 离散的大步长的最大值
% b：小步长的最大值

error=0.125/180*pi;  %假定平行于坐标轴的误差角度
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):h2

b-h2);
x3=b:h1:a;
x=[x1,x2,x3];
disp('%%%%%%%%%%%计算的次数为%%%%%%%%%%%%%%%')
N=length(x)
y=x;
Z=zeros(N,N);
for i=1:1 %x坐标
i
for j=1:N  % y坐标
      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轴负半轴
         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));
      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

错误提示；

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(

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有23人回复
导师想让我从独立一作变成了共一第一已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab 运行ode45出错，不知道什么原因已经有4人回复
MATLAB解方程提示错误‘is not a valid expression equation’ 已经有10人回复
matlab为啥出现以下这种错误已经有6人回复
以下matlab代码出现问题已经有10人回复
请大侠指教matlab的这个错误提示到底是个神马意思？已经有12人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab的m文件关联问题已经有17人回复
请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！已经有9人回复
matlab-出错提示Undefined method for input arguments type 'double'. 已经有5人回复
matlab语句问题已经有5人回复
求助一个matlab关于分类讨论的问题已经有9人回复
大家爱帮我看看这段MATLAB程序为啥出错了已经有5人回复
问下这个问题怎么解决.....真心啊已经有4人回复
Matlab两个小问题已经有10人回复
同样的MATLAB程序，为何在别人的电脑上运行不了已经有21人回复
VS2005下matlab调用已经有8人回复
【求助】用MATLAB画三维图程序编写的问题已经有6人回复
【求助】程序出错已经有4人回复

1楼2014-05-10 01:32:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

bluethky

木虫 (著名写手)

应助: 1 (幼儿园)
金币: 4460.8
散金: 849
红花: 11
帖子: 1354
在线: 227.3小时
虫号: 2609601
注册: 2013-08-21
专业: 物理学II

引用回帖:

3楼: Originally posted by wurongjun at 2014-05-10 06:53:34
你的那个u值有问题!
你定义的它是一个符号变量:syms u
程序里又把它当数值变量使用,所以出错!
至于你说的在命令行可以运行,原因可能是,你的work space里面已经给u赋值了,你可以清空space再看一下!
纵观程序 u是 ...

我把它设置成 int（f1,'u',0,1）这种类型的，系统提示对Undefined function or variable 'u'。。。。。很诡异