| 1 | 1/1 | 返回列表 |
| 查看: 495 | 回復(fù): 0 | ||
lin007zhen鐵蟲 (初入文壇)
|
[求助]
求助,微分代數(shù)方程組
|
|
function dx = yuanpan4(t,x) m=10;r=0.5;k=100;g=9.8;u=0.4; dx=zeros(5,1); if x(4)>0 dx(1)=x(2); dx(2)=(-k*x(1)+x(5))/m; dx(3)= (-2*x(5))/(m*r); dx(4)=x(2)-r*x(3)-x(4); dx(5)= x(5)+m*g*u; elseif x(4)==0 dx(1)=x(2); dx(2)=(-k*x(1)+x(5))/m; dx(3)=(-2*x(5))/(m*r); dx(4)=x(2)-r*x(3)-x(4); dx(5)=0; else dx(1)=x(2); dx(2)=(-k*x(1)+x(5))/m; dx(3)=(-2*x(5))/(m*r); dx(4)=x(2)-r*x(3)-x(4); dx(5)=x(5)-m*g*u; end end 這是函數(shù)。 tspan=[0,16]; x0=[3;3;5;0.5;-39.2]; M=diag([1,1,1,0,0]); options=odeset('Mass',M,'RelTol',1e-6,'AbsTol',[1e-10,1e-10,1e-10,1e-10,1e-10]); [t,x]=ode15s(@yuanpan4,tspan,x0,options);plot(t,x(:,1),'-b',t,x(:,4),'-r',t,x(:,5),'-g'); grid on; 這是主程序。 Warning: Failure at t=1.192445e-02. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.775558e-17) at time t. 這是錯誤提示 |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +3 | 苑豆豆 2026-03-20 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機?/a> 2026-03-16 | 3/150 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|