 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1038 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
zhangleiii管理員
|
[求助]
matlab出錯(cuò) 求助高手解答
|
|||
|
function fun=weifen11 global r x ; A ,B, C ,D, F (r) ,G(r) ,H(r) 賦值 x=zeros(2,1); x0=[1 100]; tspan=[0 1]; fun=@(t,x)[I*(x(2)+M)*jifen1-J*x(2)*jifen2-0.1*x(2) -K*x(1)*jifen3-0.01*x(1)]; function a=jifen1 a=Roberg(f1,0,1.25*10^-6) function y1= f1 y1=F.*2.*pi.*r.*A.*(D.*H.*x(1,1)+C.*G.*x(2,1))/(1+C.*G.*x(2,1)+B.*F.*x(2,1)+D.*H.*x(1,1)); %n2 end end function b=jifen2 b=Roberg(f2,0,1.25*10^-6) function y2=f2 y2=G.*2.*pi.*r.*A.*(1+B.*F.*x(2))/(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); end end function c=jifen3 c=Roberg(f3,0,1.25*10^-6) function y3=f3 y3=H*2*pi.*r.*A.*(1+B.*F.*x(2))/(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %n1 end end options=odeset('reltol',1e-6,'abstol',1e-8); [t,x]=ode45(fun,tspan,x0,options); plot(t,x) end 運(yùn)行出錯(cuò)： ??? Subscripted assignment dimension mismatch. Error in ==> Roberg at 12 T(1,1)=(h/2)*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b)); Error in ==> weifen11>jifen1 at 22 a=Roberg(f1,0,1.25*10^-6) Error in ==> weifen11>@(t,x)[I*(x(2)+M)*jifen1-J*x(2)*jifen2-0.1*x(2);-K*x(1)*jifen3-0.01*x(1)] at 18 fun=@(t,x)[I*(x(2)+M)*jifen1-J*x(2)*jifen2-0.1*x(2) Error in ==> funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ==> ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in ==> weifen11 at 40 [t,x]=ode45(fun,tspan,x0,options); |
回復(fù)此樓» 猜你喜歡
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有57人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有22人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】拜求高手解釋下matlab的報(bào)錯(cuò)-----“out of memory”
已經(jīng)有13人回復(fù)
- 【求助】請(qǐng)各位高手告訴我：這種想法能不能用matlab實(shí)現(xiàn)？200金幣求助
已經(jīng)有16人回復(fù)
- 【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧
已經(jīng)有105人回復(fù)
- 【求助】CASTEP計(jì)算出錯(cuò)，運(yùn)行失敗，請(qǐng)高手幫忙。
已經(jīng)有11人回復(fù)
- 【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進(jìn)來指點(diǎn)下）
已經(jīng)有5人回復(fù)
- 【求助】硅化反應(yīng)，求高手解答
已經(jīng)有16人回復(fù)
- 【求助】大家看看，matlab安裝過程的截圖，找出錯(cuò)
已經(jīng)有14人回復(fù)
- 【求助】MS4.4在REDHAT5.0集群上安裝出錯(cuò)，附步驟，望高手指正
已經(jīng)有16人回復(fù)
dbb627
新蟲
-
專家經(jīng)驗(yàn): +4 - 計(jì)算強(qiáng)帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時(shí)
- 蟲號(hào): 149791
- 注冊(cè): 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計(jì)算模擬
3樓2011-11-23 18:06:04
zhangleiii
超級(jí)版主
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.5
- 帖子: 19
- 在線: 8.2小時(shí)
- 蟲號(hào): 1468594
- 注冊(cè): 2011-10-31
- 專業(yè): 光學(xué)
2樓2011-11-23 16:27:10
zhangleiii
主管區(qū)長
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.5
- 帖子: 19
- 在線: 8.2小時(shí)
- 蟲號(hào): 1468594
- 注冊(cè): 2011-10-31
- 專業(yè): 光學(xué)
|
謝謝你~~上個(gè)問題解決了。我又重新調(diào)試了一下 如下： function fun=weifen11 global r x ; A B C D F(r) G(r) H(r) 賦值 fun=@(t,x)[-K.*x(1).*jifen3-0.01.*x(1) I.*(x(2)+M).*jifen1-J.*x(2).*jifen2-0.1.*x(2)]; function a=jifen1 a=Roberg(f1,0,1.25*10^-6) function y1= f1 y1=F.*2.*pi.*r.*A.*(D.*H.*x(1)+C.*G.*x(2))./(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %isen2 end end function b=jifen2 b=Roberg(f2,0,1.25*10^-6) function y2=f2 y2=G.*2.*pi.*r.*A.*(1+B.*F.*x(2))./(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %isan1 end end function c=jifen3 c=Roberg(f3,0,1.25*10^-6) function y3=f3 y3=H.*2.*pi.*r.*A.*(1+B.*F.*x(2))./(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %ipn1 end end x0=[10 10]; tspan=[0 1]; options=odeset('reltol',1e-6,'abstol',1e-8); [t,x]=ode45(weifen11,tspan,x0,options) plot(t,x) end 運(yùn)行出錯(cuò) ??? Maximum recursion limit of 500 reached. Use set(0,'RecursionLimit',N) to change the limit. Be aware that exceeding your available stack space can crash MATLAB and/or your computer. Error in ==> besschk at 50 end 我實(shí)在弄不懂了，謝謝幫忙啊~~ |
4樓2011-11-23 21:17:27
dbb627
新蟲
-
專家經(jīng)驗(yàn): +4 - 計(jì)算強(qiáng)帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時(shí)
- 蟲號(hào): 149791
- 注冊(cè): 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計(jì)算模擬
【答案】應(yīng)助回帖
zhangleiii(金幣+2): 謝謝你的意見啦，問題解決了~~ 2011-11-25 09:12:55
|
你的程序一定使用了遞歸調(diào)用，MATLAB限制了遞歸調(diào)用的次數(shù)，默認(rèn)是500次 所謂遞歸就是自己調(diào)用自己，因?yàn)檫@樣比較占用資源，MATLAB為了防止死循環(huán)，所以限制了遞歸次數(shù) 你可以使用 set(0,'RecursionLimit',1000)將遞歸次數(shù)修改為1000次 我覺得你的函數(shù)修改下比較好 function fun=weifen11(t,x) global r x ; A B C D F(r) G(r) H(r) 賦值 a=Roberg(f1,0,1.25*10^-6) b=Roberg(f2,0,1.25*10^-6) c=Roberg(f3,0,1.25*10^-6) fun=[-K.*x(1).*c-0.01.*x(1) I.*(x(2)+M).*a-J.*x(2).*b-0.1.*x(2)]; function y1= f1(x) y1=F.*2.*pi.*r.*A.*(D.*H.*x(1)+C.*G.*x(2))./(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %isen2 function y2=f2(x) y2=G.*2.*pi.*r.*A.*(1+B.*F.*x(2))./(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %isan1 function y3=f3(x) y3=H.*2.*pi.*r.*A.*(1+B.*F.*x(2))./(1+C.*G.*x(2)+B.*F.*x(2)+D.*H.*x(1)); %ipn1 |
5樓2011-11-24 19:40:24
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +22 | lllllos 2026-03-30 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4 | 1段時(shí)間 2026-03-31 | 5/250 |
|
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
 130+4
|
Yishengeryi 2026-03-30 | 4/200 |
|
|
|
[考研] 315求調(diào)劑 +6 | akie... 2026-03-28 | 7/350 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-30 | 6/300 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 085601材料工程找調(diào)劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|