 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 506 | 回復(fù): 1 | ||
xuanshan新蟲 (初入文壇)
|
[求助]
程序問題 已有1人參與
|
|
function dy=rigid4(t,y) w=78;%輪胎的轉(zhuǎn)速 R=0.317;%輪胎半徑 u=10;%外界風(fēng)速 ls=1.225;%空氣密度 lg=2.5*10^3;%磨粒密度 d=10^(-5);%磨粒直徑 g=9.806;%重力加速度 m=1.79*10^(-5);%空氣動(dòng)力粘度 n=m/ls;%運(yùn)動(dòng)粘度 k=3*ls/(4*lg*d); g1=(1-ls/lg)*g; dy=zeros(4,1);%一個(gè)列向量 V=sqrt(y(2)^2+y(4)^2); u0=w*R^2/(R+y(1)); U1=(0.32*y(1)*u^1.5*n*y(3)-0.005*y(1)^3*u^2.5)/(n^1.5*y(3)^1.5)+u0; U2=(0.08*n*y(1)^2*u^1.5*y(3)-0.002*y(1)^4*u^2.5)/(n^1.5*y(3)^2.5); U=sqrt(U1^2+U2^2); Re=ls*d*abs(U-V)/m; Cd=19.5/Re^0.633; dy(1)=y(2); dy(3)=y(4); dy(2)=k*Cd*abs(U2-y(2))*(U2-y(2)); dy(4)=-g1+k*Cd*abs(U1-y(4))*(U1-y(4)); dy=[dy(1) dy(2) dy(3) dy(4)]'; end 上面是主程序，下面是命令窗口中的程序 tspan=linspace(0,1,1000); [t,y]=ode45(@rigid4,tspan,[0 2.3 0 2.1]); plot(y(:,1),y(:,3)) axis([0 0.003 0 1]) 得出的結(jié)果只顯示坐標(biāo)軸不顯示圖像，麻煩哪位大神給指導(dǎo)指導(dǎo)，在這里謝過！ |
回復(fù)此樓» 猜你喜歡
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有4人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有4人回復(fù)
-
086000生物與醫(yī)藥調(diào)劑
已經(jīng)有7人回復(fù)
-
305求調(diào)劑
已經(jīng)有7人回復(fù)
-
348求調(diào)劑
已經(jīng)有3人回復(fù)
-
332求調(diào)劑
已經(jīng)有6人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有6人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有10人回復(fù)
-
356求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
已經(jīng)有7人回復(fù)
文刀星
木蟲 (小有名氣)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 3651.2
- 散金: 500
- 紅花: 5
- 帖子: 60
- 在線: 100.9小時(shí)
- 蟲號(hào): 3693767
- 注冊(cè): 2015-02-22
- 性別: GG
- 專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
function dy=rigid4(t,y) w=78;%輪胎的轉(zhuǎn)速 R=0.317;%輪胎半徑 u=10;%外界風(fēng)速 ls=1.225;%空氣密度 lg=2.5*10^3;%磨粒密度 d=10^(-5);%磨粒直徑 g=9.806;%重力加速度 m=1.79*10^(-5);%空氣動(dòng)力粘度 n=m/ls;%運(yùn)動(dòng)粘度 k=3*ls/(4*lg*d); g1=(1-ls/lg)*g; dy=zeros(4,1);%一個(gè)列向量 V=sqrt(y(2)^2+y(4)^2); u0=w*R^2/(R+y(1)); U1=(0.32*y(1)*u^1.5*n*y(3)-0.005*y(1)^3*u^2.5)/(n^1.5*y(3)^1.5)+u0;%%%% U1=NAN 有問題，問題在于y(3)=0,做除數(shù)了 U2=(0.08*n*y(1)^2*u^1.5*y(3)-0.002*y(1)^4*u^2.5)/(n^1.5*y(3)^2.5);%%% 同樣，U2=NAN，也由于y(3)=0,在分母上 U=sqrt(U1^2+U2^2); % NAN Re=ls*d*abs(U-V)/m;% NAN Cd=19.5/Re^0.633;% NAN dy(1)=y(2); dy(3)=y(4); dy(2)=k*Cd*abs(U2-y(2))*(U2-y(2)); %NAN dy(4)=-g1+k*Cd*abs(U1-y(4))*(U1-y(4));% NAN end 你求出的數(shù)值解y都是ＮＡＮ，當(dāng)然畫不出圖來(lái)啦！出現(xiàn)NANd的的情況我已經(jīng)在上面的代碼中注明了。還是好好檢查檢查的的微分方程組是否有錯(cuò)。把初值改為非零，但是也還積不出來(lái)。感覺可能是你的微分方程有問題！ |
2樓2015-12-16 14:09:03
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xuanshan 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +3 | 小懶蟲不懶了 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 311（085601）求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +11 | zhn03 2026-03-25 | 12/600 |
|
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，各位老師看看我�。。� +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|