zhangleiii

版主

應(yīng)助: 0 (幼兒園)
金幣: 7.5
帖子: 19
在線: 8.2小時
蟲號: 1468594
注冊: 2011-10-31
專業(yè): 光學(xué)

[求助] matlab出錯求助高手解答

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

運行出錯：
??? 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論文潤色/翻譯怎么收費? 已經(jīng)有122人回復(fù)
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會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)價值貼推薦，對您同樣有幫助:

【求助】拜求高手解釋下matlab的報錯-----“out of memory” 已經(jīng)有13人回復(fù)
【求助】請各位高手告訴我：這種想法能不能用matlab實現(xiàn)？200金幣求助已經(jīng)有16人回復(fù)
【整理】今天整理硬盤上的資料，偶然發(fā)現(xiàn)的好玩的小東西就直接發(fā)上來吧已經(jīng)有105人回復(fù)
【求助】CASTEP計算出錯，運行失敗，請高手幫忙。已經(jīng)有11人回復(fù)
【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進來指點下）已經(jīng)有5人回復(fù)
【求助】硅化反應(yīng)，求高手解答已經(jīng)有16人回復(fù)
【求助】大家看看，matlab安裝過程的截圖，找出錯已經(jīng)有14人回復(fù)
【求助】MS4.4在REDHAT5.0集群上安裝出錯，附步驟，望高手指正已經(jīng)有16人回復(fù)

1樓 2011-11-23 10:16:04

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

zhangleiii

版主

應(yīng)助: 0 (幼兒園)
金幣: 7.5
帖子: 19
在線: 8.2小時
蟲號: 1468594
注冊: 2011-10-31
專業(yè): 光學(xué)

謝謝你~~上個問題解決了。我又重新調(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
運行出錯
??? 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

我實在弄不懂了，謝謝幫忙啊~~

贊一下

回復(fù)此樓

4樓2011-11-23 21:17:27

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

查看全部 6 個回答

zhangleiii

版主

應(yīng)助: 0 (幼兒園)
金幣: 7.5
帖子: 19
在線: 8.2小時
蟲號: 1468594
注冊: 2011-10-31
專業(yè): 光學(xué)

怎么沒人回復(fù)呀~先謝謝各位了！

贊一下

回復(fù)此樓

2樓2011-11-23 16:27:10

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

dbb627

新蟲

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

zhangleiii(金幣+1): 謝謝！ 2011-11-24 18:46:58

Roberg at 12
你的Roberg積分程序有問題
這里沒有你的Roberg程序，看不出來
估計下錯誤應(yīng)該是
subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

3樓2011-11-23 18:06:04

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

dbb627

無蟲

專家經(jīng)驗: +4
計算強帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計算模擬

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

zhangleiii(金幣+2): 謝謝你的意見啦，問題解決了~~ 2011-11-25 09:12:55

你的程序一定使用了遞歸調(diào)用，MATLAB限制了遞歸調(diào)用的次數(shù)，默認(rèn)是500次

所謂遞歸就是自己調(diào)用自己，因為這樣比較占用資源，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

贊一下(1人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

5樓2011-11-24 19:40:24

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

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

[考研] 考研調(diào)劑 +11	Amber00 2026-03-31	11/550	2026-04-01 11:32 by wangjy2002
[考研] 275求調(diào)劑 +7	waltzh 2026-04-01	7/350	2026-04-01 11:23 by 無際的草原
[考研] 材料0856 英一數(shù)二 323 求調(diào)劑 +8	袁sy 2026-04-01	8/400	2026-04-01 10:11 by 唐沐兒
[考研] 求0861交通運輸專碩or材料專碩調(diào)劑 +4	勒布朗@ 2026-03-31	4/200	2026-04-01 09:54 by 一只好果子?
[考研] 0805求調(diào)劑 +6	是水分 2026-03-31	6/300	2026-04-01 09:05 by oooqiao
[考研] 一志愿中國科學(xué)院大學(xué)265求調(diào)劑 +7	恬淡ye 2026-03-31	8/400	2026-04-01 08:38 by JourneyLucky
[考研] 土木304求調(diào)劑 +5	頂級擦擦 2026-03-31	5/250	2026-04-01 08:15 by fdcxdystjk￥
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +5	Ripcord上岸 2026-03-27	5/250	2026-03-31 19:52 by mg1014
[考研] 生醫(yī)工0831調(diào)劑求推薦 +4	小熊睿睿_s 2026-03-27	6/300	2026-03-31 17:12 by 記事本2026
[考研] 材料專碩調(diào)劑 +13	椰椰。 2026-03-29	13/650	2026-03-31 16:37 by hypershenger
[考研] 289求調(diào)劑 +3	Acesczlo 2026-03-29	4/200	2026-03-31 14:48 by 熱情沙漠
[考研] 334求調(diào)劑 +7	Trying] 2026-03-31	7/350	2026-03-31 12:33 by 無際的草原
[考研] 274求調(diào)劑 +6	xiao愛同學(xué) 2026-03-30	6/300	2026-03-31 10:04 by cal0306
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[考研] 328求調(diào)劑 +8	嗯滴的基本都 2026-03-27	8/400	2026-03-30 17:20 by Wang200018
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8	AZMK 2026-03-25	12/600	2026-03-28 18:15 by AZMK
[考研] 壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 復(fù)試調(diào)劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落�？伤�
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-27	8/400	2026-03-28 12:42 by 唐沐兒
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀