當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

bluethky

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 4698.8
散金: 849
紅花: 11
帖子: 1398
在線: 227.7小時(shí)
蟲號(hào): 2609601
注冊(cè): 2013-08-21
專業(yè): 植物學(xué)研究的新技術(shù)、新方

[求助] matlab程序錯(cuò)誤，提示如下【弄了半夜還是不行】已有2人參與

function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b);
%  計(jì)算理論的二維levy分布的圖象畫出來(lái)的3維曲面圖
% t0 ：levy walk  趨于零的時(shí)間
%  lambda ：冪律值
% T：時(shí)間
% a: 離散的大步長(zhǎng)的最大值
% b：小步長(zhǎng)的最大值

error=0.125/180*pi;  %假定平行于坐標(biāo)軸的誤差角度
syms u
c=-gamma(2-lambda)*t0^(lambda-1)*cos(pi*lambda/2)/lambda;
h1=20;  %大的離散步長(zhǎng)
h2=0.1; %小的離散步長(zhǎng)
x1=-a:h1:-b;
x2=-(b-h2):h2

b-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(

回復(fù)此樓

» 猜你喜歡

考研調(diào)劑已經(jīng)有4人回復(fù)
281求調(diào)劑已經(jīng)有4人回復(fù)
0805 316求調(diào)劑已經(jīng)有6人回復(fù)
085601求調(diào)劑總分293英一數(shù)二已經(jīng)有3人回復(fù)
08工學(xué)調(diào)劑已經(jīng)有17人回復(fù)
340求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
食品專碩一志愿雙一流 328 已經(jīng)有4人回復(fù)
材料調(diào)劑已經(jīng)有6人回復(fù)
300分，材料，求調(diào)劑，英一數(shù)二已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab 運(yùn)行ode45出錯(cuò)，不知道什么原因已經(jīng)有4人回復(fù)
MATLAB解方程提示錯(cuò)誤‘is not a valid expression equation’ 已經(jīng)有10人回復(fù)
matlab為啥出現(xiàn)以下這種錯(cuò)誤已經(jīng)有6人回復(fù)
以下matlab代碼出現(xiàn)問(wèn)題已經(jīng)有10人回復(fù)
請(qǐng)大俠指教matlab的這個(gè)錯(cuò)誤提示到底是個(gè)神馬意思？已經(jīng)有12人回復(fù)
用matlab求解非線性方程組說(shuō)無(wú)解，一定是方程組本身無(wú)解，還是有可能程序有問(wèn)題呢？已經(jīng)有11人回復(fù)
matlab的m文件關(guān)聯(lián)問(wèn)題已經(jīng)有17人回復(fù)
請(qǐng)matlab高手看一下這個(gè)問(wèn)題怎樣編程才好，附程序代碼，謝謝！已經(jīng)有9人回復(fù)
matlab-出錯(cuò)提示Undefined method for input arguments type 'double'. 已經(jīng)有5人回復(fù)
matlab語(yǔ)句問(wèn)題已經(jīng)有5人回復(fù)
求助一個(gè)matlab關(guān)于分類討論的問(wèn)題已經(jīng)有9人回復(fù)
大家愛(ài)幫我看看這段MATLAB程序?yàn)樯冻鲥e(cuò)了已經(jīng)有5人回復(fù)
問(wèn)下這個(gè)問(wèn)題怎么解決.....真心啊已經(jīng)有4人回復(fù)
Matlab兩個(gè)小問(wèn)題已經(jīng)有10人回復(fù)
同樣的MATLAB程序，為何在別人的電腦上運(yùn)行不了已經(jīng)有21人回復(fù)
VS2005下matlab調(diào)用已經(jīng)有8人回復(fù)
【求助】用MATLAB畫三維圖程序編寫的問(wèn)題已經(jīng)有6人回復(fù)
【求助】程序出錯(cuò) 已經(jīng)有4人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wurongjun

專家顧問(wèn) (職業(yè)作家)

專家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時(shí)
蟲號(hào): 1368482
注冊(cè): 2011-08-14
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
管轄: 數(shù)學(xué)

【答案】應(yīng)助回帖

符號(hào)計(jì)算與數(shù)值計(jì)算混合使用是Matlab程序的大忌!
喪失了快速計(jì)算得優(yōu)點(diǎn)!
你可以把(0,1)區(qū)間劃分為更小的區(qū)間,然后考慮數(shù)值積分!

贊一下(1人)

回復(fù)此樓

善惡到頭終有報(bào),人間正道是滄桑.

7樓2014-05-10 17:31:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個(gè)回答

bluethky

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 4698.8
散金: 849
紅花: 11
帖子: 1398
在線: 227.7小時(shí)
蟲號(hào): 2609601
注冊(cè): 2013-08-21
專業(yè): 植物學(xué)研究的新技術(shù)、新方

程序本身在命令窗口可以計(jì)算，，，但是在m文件中一直報(bào)錯(cuò)。。。。我也沒(méi)找到錯(cuò)哪了，，，求大俠給點(diǎn)意見(jiàn)，問(wèn)度娘也沒(méi)有找出原因。。。。。
太困了。。。。先謝謝了！
早安！

贊一下

回復(fù)此樓

2樓2014-05-10 01:35:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wurongjun

專家顧問(wèn) (職業(yè)作家)

專家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時(shí)
蟲號(hào): 1368482
注冊(cè): 2011-08-14
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
管轄: 數(shù)學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

你的那個(gè)u值有問(wèn)題!
你定義的它是一個(gè)符號(hào)變量:syms u
程序里又把它當(dāng)數(shù)值變量使用,所以出錯(cuò)!
至于你說(shuō)的在命令行可以運(yùn)行,原因可能是,你的work space里面已經(jīng)給u賦值了,你可以清空space再看一下!
縱觀程序 u是需要賦值的!

贊一下(2人)

回復(fù)此樓

善惡到頭終有報(bào),人間正道是滄桑.

3樓2014-05-10 06:53:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

bluethky

木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 4698.8
散金: 849
紅花: 11
帖子: 1398
在線: 227.7小時(shí)
蟲號(hào): 2609601
注冊(cè): 2013-08-21
專業(yè): 植物學(xué)研究的新技術(shù)、新方

引用回帖:

3樓: Originally posted by wurongjun at 2014-05-10 06:53:34
你的那個(gè)u值有問(wèn)題!
你定義的它是一個(gè)符號(hào)變量:syms u
程序里又把它當(dāng)數(shù)值變量使用,所以出錯(cuò)!
至于你說(shuō)的在命令行可以運(yùn)行,原因可能是,你的work space里面已經(jīng)給u賦值了,你可以清空space再看一下!
縱觀程序 u是 ...

由于程序里是關(guān)于u在（0,1）的積分，我們不是都是對(duì)其先設(shè)置成，符號(hào)變量在進(jìn)行積分嗎？我不是太明白。。不知我應(yīng)該怎樣將其設(shè)置成數(shù)值變量？求老師指點(diǎn)一下。。。
不過(guò)你的解答應(yīng)該是正確的

贊一下

回復(fù)此樓

4樓2014-05-10 10:03:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑323材料與化工 +3	1124361 2026-03-24	3/150	2026-03-24 20:24 by peike
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6	小羊睡著了? 2026-03-23	9/450	2026-03-24 20:21 by 小羊睡著了?
[考研] 085602 289分求調(diào)劑 +5	WWW西西弗斯 2026-03-24	5/250	2026-03-24 18:51 by jhhcooi
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +7	holy J 2026-03-21	7/350	2026-03-24 17:24 by xiaohai104
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 336求調(diào)劑 +4	收到VS 2026-03-20	4/200	2026-03-23 19:02 by macy2011
[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風(fēng)里
[考研] 315分，誠(chéng)求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 286求調(diào)劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 329求調(diào)劑 +9	想上學(xué)吖吖 2026-03-19	9/450	2026-03-20 22:01 by luoyongfeng
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO