| 5 | 1/1 | 返回列表 |
| 查看: 1043 | 回復(fù): 4 | ||||
hailiang金蟲 (正式寫手)
|
[求助]
請(qǐng)教MATLAB中下面的提示是什么意思?
|
|
請(qǐng)教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]) |
至尊木蟲 (著名寫手)
驃騎將軍

金蟲 (正式寫手)
|
麻煩幫我看看下面求解微分方程的程序問題出在哪? 不對(duì)f循環(huán)的時(shí)候沒有問題的,做循環(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];%積分時(shí)間 %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)驗(yàn): +3 |
|
這里有兩個(gè)錯(cuò)誤 : 1. clear all 的位置 它應(yīng)該在程序的最上方,而在你的程序里,它在 global f 之后,這樣運(yùn)行到 clear all 之后,全局變量 f 會(huì)被刪除,你的子程序就無法獲得 f 的值了 2. dy=zeros(3,1); 子程序里面,dy 的預(yù)設(shè)的錯(cuò)誤,你的程序假設(shè)它是 (2, 1) : dy(1) / dy(2), 但是在程序一開始,卻執(zhí)行 dy=zeros(3,1); 這樣就會(huì)出現(xiàn) 長(zhǎng)度錯(cuò)誤,正確的應(yīng)該是 dy=zeros(2,1); 不過因?yàn)槲覜]有你的 getmax 程序,無法運(yùn)行,但是至少能解決你開始所提出的問題了 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +6 | Liwangman 2026-03-15 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] (081700)化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|