版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

[求助] 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('%%%%%%%%%%%計算的次數(shù)為%%%%%%%%%%%%%%%')
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(

回復此樓

» 猜你喜歡

288求調劑已經(jīng)有5人回復
281求調劑（0805）已經(jīng)有19人回復
材料考研調劑已經(jīng)有3人回復
346求調劑[0856] 已經(jīng)有6人回復
294求調劑材料與化工專碩已經(jīng)有9人回復
328求調劑，英語六級551，有科研經(jīng)歷已經(jīng)有12人回復
一志愿西安交通大學材料工程專業(yè) 282分求調劑已經(jīng)有5人回復
085601材料工程專碩求調劑已經(jīng)有7人回復
321求調劑已經(jīng)有6人回復
一志愿中海洋材料工程專碩330分求調劑已經(jīng)有7人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

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

回帖支持 ( 顯示支持度最高的前 50 名 )

wurongjun

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

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

【答案】應助回帖

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

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

贊一下(2人)

回復此樓

善惡到頭終有報,人間正道是滄桑.

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

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

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

引用回帖:

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

我把它設置成 int（f1,'u',0,1）這種類型的，系統(tǒng)提示對Undefined function or variable 'u'。。。。。很詭異

贊一下(1人)

回復此樓

5樓2014-05-10 10:14:42

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

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

我在每個int前面加了一個simplify后可以算出結果，但是速度不是一般的慢，你還有其他好方法沒有？

贊一下(1人)

回復此樓

6樓2014-05-10 12:27:28

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

wurongjun

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

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

【答案】應助回帖

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

贊一下(1人)

回復此樓

善惡到頭終有報,人間正道是滄桑.

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

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

yogort

鐵桿木蟲 (正式寫手)

應助: 25 (小學生)
金幣: 5127
紅花: 20
帖子: 520
在線: 302小時
蟲號: 1060529
注冊: 2010-07-19
專業(yè): 偏微分方程

這是不是就是分段函數(shù)要畫圖啊？

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

bluethky

9樓2014-05-11 18:46:44

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

普通回帖

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

程序本身在命令窗口可以計算，，，但是在m文件中一直報錯。。。。我也沒找到錯哪了，，，求大俠給點意見，問度娘也沒有找出原因。。。。。
太困了。。。。先謝謝了！
早安！

贊一下

回復此樓

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

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

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

引用回帖:

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

贊一下

回復此樓

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

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

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

引用回帖:

7樓: Originally posted by wurongjun at 2014-05-10 17:31:40
符號計算與數(shù)值計算混合使用是Matlab程序的大忌!
喪失了快速計算得優(yōu)點!
你可以把(0,1)區(qū)間劃分為更小的區(qū)間,然后考慮數(shù)值積分!

確實速度老慢了。。。。區(qū)間不能再少了，，，我現(xiàn)在考慮用inline 然后用quad，但是由于其中有其他變量，提示錯誤

贊一下

回復此樓

8樓2014-05-10 19:31:24

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

bluethky

木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時
蟲號: 2609601
注冊: 2013-08-21
專業(yè): 植物學研究的新技術、新方

送紅花一朵

引用回帖:

9樓: Originally posted by yogort at 2014-05-11 18:46:44
這是不是就是分段函數(shù)要畫圖啊？

本來是一個0到無窮的積分，我在1處截斷了，再將1后面的積分轉換到0到1，用數(shù)值的方法算出來了，謝謝大俠！

贊一下

回復此樓

10樓2014-05-12 12:50:59

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

相關版塊跳轉我要訂閱樓主逃課的小屁孩的主題更新

返回列表

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

[考研] 281求調劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 能源材料化學課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	11/550	2026-03-19 08:11 by 脫穎而出
[考研] 0703化學 305求調劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 344求調劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 一志愿中國海洋大學，生物學，301分，求調劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 085601專碩，總分342求調劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考博] 環(huán)境領域全國重點實驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 工科材料085601 279求調劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 301求調劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 考研調劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 277材料科學與工程080500求調劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 0703化學調劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 0856求調劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 求老師收留調劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 中科大材料與化工319求調劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 266求調劑 +4	學員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 材料工程調劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 307求調劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 328化工專碩求調劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky

24小時熱門版塊排行榜

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

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

【答案】應助回帖

【答案】應助回帖

» 本帖已獲得的紅花（最新10朵）

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

» 本主題相關價值貼推薦，對您同樣有幫助: