 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 284 | 回復(fù): 1 | ||||
xuanshan新蟲 (初入文壇)
|
[求助]
程序問題 已有1人參與
|
|
function dy=rigid(t,y) a=0.06;%輪胎的接地印跡長(zhǎng)度 w=63;%輪胎的轉(zhuǎn)速 R=0.317;%輪胎半徑 u=1;%外界風(fēng)速 ls=1.293;%空氣密度 lg=2.5*10^3;%磨粒密度 d=10^(-5);%磨粒直徑 g=9.806;%重力加速度 m=1.79*10^(-6);%空氣動(dòng)力粘度 k1=3*ls/(4*lg*d); g1=(1-ls/lg)*g; dy=zeros(4,1);%一個(gè)列向量 V=sqrt(y(2)^2+y(4)^2); U=w*sqrt(R^2-(y(1)^2+y(3)^2)); Re1=ls*d*abs(u-V)/m; Re2=ls*d*abs(U-V)/m; Cd1=24*(1+0.0975*Re1-0.636*Re1^2)/Re1; Cd2=24*(1+0.0975*Re2-0.636*Re2^2)/Re2; k2=3*Cd2*ls/(4*lg*d); dy(1)=y(2); dy(3)=y(4); dy(2)=k1*(Cd1*(u-y(2))^2+Cd2*(w*(y(1)-a)/(sqrt(R^2-((y(1)-a)^2+(y(3)+sqrt(R^2-a^2))^2))+y(2))^2)); dy(4)=-g1+k2*(w*(y(3)+sqrt(R^2-a^2))/sqrt(R^2-((y(1)-a)^2+(y(3)+sqrt(R^2-a^2))^2)+y(4))^2); dy=[dy(1) dy(2) dy(3) dy(4)]'; end 這是我的程序代碼，但是在運(yùn)算的時(shí)候出現(xiàn)下列情況 tspan=linspace(0,1.5,1000); [T Y]=ode23s(@rigid,tspan,[0 19.6 0 3.87]); Warning: Failure at t=7.361808e-06. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (2.615440e-20) at time t. > In ode23s at 402 麻煩哪位大神給指導(dǎo)一下，謝謝！ |
回復(fù)此樓» 猜你喜歡
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有5人回復(fù)
-
297求調(diào)劑
已經(jīng)有5人回復(fù)
-
302求調(diào)劑
已經(jīng)有10人回復(fù)
-
266求調(diào)劑
已經(jīng)有8人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有17人回復(fù)
-
344求調(diào)劑
已經(jīng)有5人回復(fù)
-
312求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有9人回復(fù)
-
288求調(diào)劑，一志愿華南理工大學(xué)071005
已經(jīng)有4人回復(fù)
-
304求調(diào)劑
已經(jīng)有13人回復(fù)
512095518
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +52 - 仿真EPI: 1
- 應(yīng)助: 118 (高中生)
- 金幣: 7665.4
- 紅花: 78
- 帖子: 712
- 在線: 186.1小時(shí)
- 蟲號(hào): 3195576
- 注冊(cè): 2014-05-10
- 性別: GG
- 專業(yè): 大氣環(huán)境與全球氣候變化
- 管轄: 計(jì)算模擬
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
應(yīng)該是你的常微分方程組四個(gè)解寫得有問題吧，估計(jì)是太過于復(fù)雜才導(dǎo)致ode23求解不了。直接用數(shù)值會(huì)好一些�；蚴悄阍囅耾de45看看。具體你還需要自己去查下類似的方程該怎么寫。 下面是我用過的供你參考，其中a-n都是參數(shù)，都是提前算好的。 function dx = hanshu(t, x) global a b c d e f g h i j k m l n dx = zeros(4,1); dx(1) = a + b*x(3) + c*x(2) - d*x(1); dx(2) = e + f*x(1) + g*x(3) + h*x(4) - i*x(2); dx(3) = j + k*x(1) - l*x(3); dx(4) = m*x(2) - n*x(4); end |
2樓2015-10-14 15:31:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xuanshan 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 材料專碩306英一數(shù)二 +9 | z1z2z3879 2026-03-16 | 12/600 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|