| 5 | 1/1 | 返回列表 |
| 查看: 1056 | 回復: 4 | ||||
hailiang金蟲 (正式寫手)
|
[求助]
請教MATLAB中下面的提示是什么意思?
|
|
請教MATLAB中下面的提示是什么意思? In an assignment A(I) = B, the number of elements in B and I must be the same. Error in hlduffing (line 5) dy(2)=f*cos(0.8*t)-0.2*y(2)-y(1)-y(1)^3; Error in odearguments (line 88) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode45 (line 114) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in hltest2013 (line 7) [t,y]=ode45('hlduffing',t0,[1;1]) |
至尊木蟲 (著名寫手)
驃騎將軍

金蟲 (正式寫手)
|
麻煩幫我看看下面求解微分方程的程序問題出在哪? 不對f循環(huán)的時候沒有問題的,做循環(huán)就出問題了 function dy = hlduffing(t,y) global f dy=zeros(3,1); dy(1)=y(2); dy(2)=f*cos(0.8*t)-0.2*y(2)-y(1)-y(1)^3; end function hlduffing1_bifur_f_getmax % 最大值法求解分岔圖 global f clear all t0=[0 100];%積分時間 %bifurcation for f=linspace(10,32,1000); % f=10:32:1000; %r的變化精度 [t,y]=ode45('hlduffing',t0,[1;1]); [Xmax]=getmax(y(:,1)); plot(f,Xmax,'k','markersize',1) hold on clear Xmax end |
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗: +3 |
|
這里有兩個錯誤 : 1. clear all 的位置 它應該在程序的最上方,而在你的程序里,它在 global f 之后,這樣運行到 clear all 之后,全局變量 f 會被刪除,你的子程序就無法獲得 f 的值了 2. dy=zeros(3,1); 子程序里面,dy 的預設的錯誤,你的程序假設它是 (2, 1) : dy(1) / dy(2), 但是在程序一開始,卻執(zhí)行 dy=zeros(3,1); 這樣就會出現(xiàn) 長度錯誤,正確的應該是 dy=zeros(2,1); 不過因為我沒有你的 getmax 程序,無法運行,但是至少能解決你開始所提出的問題了 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +6 | 生物工程調(diào)劑 2026-03-17 | 10/500 |
|
|
[考研] 一志愿中南化學(0703)總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|