版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

bluethky

木蟲(chóng) (著名寫(xiě)手)

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

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

function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b);
%  計(jì)算理論的二維levy分布的圖象畫(huà)出來(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ù)此樓

» 猜你喜歡

工科材料085601 279求調(diào)劑已經(jīng)有8人回復(fù)
317求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有5人回復(fù)
274求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有9人回復(fù)
招收調(diào)劑碩士已經(jīng)有11人回復(fù)
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有4人回復(fù)
【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有14人回復(fù)
材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦已經(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)于分類(lèi)討論的問(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畫(huà)三維圖程序編寫(xiě)的問(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的回帖

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

wurongjun

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

專(zhuān)家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時(shí)
蟲(chóng)號(hào): 1368482
注冊(cè): 2011-08-14
性別: GG
專(zhuān)業(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

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時(shí)
蟲(chóng)號(hào): 2609601
注冊(cè): 2013-08-21
專(zhuān)業(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是 ...

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

贊一下(1人)

回復(fù)此樓

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

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

bluethky

木蟲(chóng) (著名寫(xiě)手)

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

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

贊一下(1人)

回復(fù)此樓

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

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

wurongjun

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

專(zhuān)家經(jīng)驗(yàn): +831
數(shù)學(xué)EPI: 9
應(yīng)助: 791 (博后)
貴賓: 0.308
金幣: 24609
散金: 310
紅花: 75
帖子: 3004
在線: 881.4小時(shí)
蟲(chóng)號(hào): 1368482
注冊(cè): 2011-08-14
性別: GG
專(zhuān)業(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的回帖

yogort

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 25 (小學(xué)生)
金幣: 5127
紅花: 20
帖子: 520
在線: 302小時(shí)
蟲(chóng)號(hào): 1060529
注冊(cè): 2010-07-19
專(zhuān)業(yè): 偏微分方程

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

贊一下(1人)

回復(fù)此樓

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

bluethky

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

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

普通回帖

bluethky

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 4693.8
散金: 849
紅花: 11
帖子: 1397
在線: 227.7小時(shí)
蟲(chóng)號(hào): 2609601
注冊(cè): 2013-08-21
專(zhuān)業(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的回帖

bluethky

木蟲(chóng) (著名寫(xiě)手)

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

引用回帖:

由于程序里是關(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的回帖

bluethky

木蟲(chóng) (著名寫(xiě)手)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

bluethky

木蟲(chóng) (著名寫(xiě)手)

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

送紅花一朵

引用回帖:

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

本來(lái)是一個(gè)0到無(wú)窮的積分，我在1處截?cái)嗔耍賹?后面的積分轉(zhuǎn)換到0到1，用數(shù)值的方法算出來(lái)了，謝謝大俠！

贊一下

回復(fù)此樓

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主逃課的小屁孩的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 0856調(diào)劑，是學(xué)校就去 +7	sllhht 2026-03-19	8/400	2026-03-20 09:30 by 每天只擺一小會(huì)
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 294求調(diào)劑材料與化工專(zhuān)碩 +14	陌の森林 2026-03-18	14/700	2026-03-19 22:38 by 學(xué)員8dgXkO
[考研] 梁成偉老師課題組歡迎你的加入 +9	一鴨鴨喲 2026-03-14	11/550	2026-03-19 17:22 by ！本暗一次！
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 311求調(diào)劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專(zhuān)利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 一志愿蘇州大學(xué)材料工程（085601）專(zhuān)碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專(zhuān)利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 材料與化工專(zhuān)碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么

24小時(shí)熱門(mén)版塊排行榜

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

» 猜你喜歡

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

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

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

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

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

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