 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 817 | 回復: 4 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
wcnmaz新蟲 (小有名氣)
|
[交流]
MATLAB7.0出錯了？ 已有2人參與
|
||
|
y(i,6)=-178075801.6*y(i,4)*y(i,5)+490910.1296*y(i,3)*y(i,2)+1128974.05+7481.104004*sin(2.047t); 為什么MATLAB7.0程序運行到上面一行出錯了？感覺也沒什么錯誤啊 怎么改，謝謝了  ODE45_main文件: clc;clear tspan=[0,1]; y0=[0; 0;0;0;0]; [t,y]=ode45('ODE45_fun',tspan,y0); [m,n]=size(y); for i=1:m y(i,6)=-178075801.6*y(i,4)*y(i,5)+490910.1296*y(i,3)*y(i,2)+1128974.05+7481.104004*sin(2.047t); end data=[t,y]; save ODE45_data.txt data -ascii subplot(2,3,1),plot(t,y(:,1)),title('y(1)') Xlabel('t');Ylabel('y'); subplot(2,3,2),plot(t,y(:,2)),title('y(2)') Xlabel('t');Ylabel('y'); subplot(2,3,3),plot(t,y(:,3)),title('y(3)') Xlabel('t');Ylabel('y'); subplot(2,3,4),plot(t,y(:,4)),title('y(4)') Xlabel('t');Ylabel('y'); subplot(2,3,5),plot(t,y(:,5)),title('y(5)') Xlabel('t');Ylabel('y'); subplot(2,3,6),plot(t,y(:,6)),title('y(6)') Xlabel('t');Ylabel('y'); grid on ODE45_fun文件： function dy=ODE45_fun(t,y) dy(1)=-1.918298553*y(3)*y(4)-121.6697369*y(5)*y(2)+0.006472085*y(2)*y(2)+15.25250926*y(5)*y(5)-0.518363603*sin(2.047t)+0.001124759; dy(2)=0.007229182*y(5)*y(1)-0.013867729*y(3)-0.005151943*sin(2.047t)+33.43424564*y(4)*y(5)-0.092169794*y(3)*y(2)-0.698266828; dy(3)=72.10986245*y(2)+0.52129529*y(4)*y(1)+0.025471074*y(3)+921.886526*y(4)/y(1)-0.47870471*y(4)+0.025471074*y(1)*t-0.38220722*cos(2.047t)-4.62279911; dy(4)=-57.37263009*y(5)+0.001053501*y(3)*y(1)+20.28825016/y(1)-0.001053501*y(3)+0.064741605*y(4)+20.28825016*t-0.006201915*sin(2.047t)-3651.904822; dy(5)=0.017284293*y(4)+0.00278644*y(1)*y(2)-0.551218454*y(2)*y(5)*y(5)+0.010839281*y(2)*y(2)*y(5)+0.020353114*cos(2.047t)+0.110594984; dy=[dy(1);dy(2);dy(3);dy(4);dy(5)]; |
回復此樓» 猜你喜歡
材料與化工調(diào)劑一志愿大連海事085600，349
已經(jīng)有4人回復
-
甲基亞磺磺酸鈉和甲基磺酸酯反應機理
已經(jīng)有4人回復
-
理學07化學 303求調(diào)劑
已經(jīng)有3人回復
-
一志愿大連理工大學材料求調(diào)劑
已經(jīng)有4人回復
-
085600材料與化工調(diào)劑
已經(jīng)有11人回復
-
生物考研337分求調(diào)劑
已經(jīng)有4人回復
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有9人回復
-
哈爾濱工業(yè)大學材料與化工專碩378求調(diào)劑
已經(jīng)有3人回復
-
求調(diào)劑，一志愿 南京航空航天大學 ，080500材料科學與工程學碩，總分289分
已經(jīng)有9人回復
-
材料科學與工程求調(diào)劑
已經(jīng)有8人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- MS7.0安裝問題
已經(jīng)有5人回復
- matlab 運行ode45出錯，不知道什么原因
已經(jīng)有4人回復
- Matlab 7.0 基礎教程 清華大學
已經(jīng)有29人回復
- matlab一運行就出現(xiàn)的問題
已經(jīng)有7人回復
- 運行matlab出現(xiàn)錯誤
已經(jīng)有4人回復
- matlab程序出錯調(diào)試
已經(jīng)有7人回復
- Quantum-ESPRESSO4.3安裝問題，求助！
已經(jīng)有4人回復
- 問下這個問題怎么解決.....真心啊
已經(jīng)有4人回復
- matlab程序問題，總是出錯，麻煩大家看看問題在那里？
已經(jīng)有6人回復
- xcysden-1.5.53.tar.gz安裝問題求助
已經(jīng)有5人回復
- 求助 matlab交點問題
已經(jīng)有19人回復
- 一筆畫問題matlab實現(xiàn)
已經(jīng)有5人回復
- matlab出錯 求助高手解答
已經(jīng)有5人回復
- 采用matlab中nlinfit擬合，出現(xiàn)問題了，不知道是哪錯了，該怎么解決，求高手指點一下
已經(jīng)有9人回復
- 【求助】matlab 中，幾個m文件調(diào)用
已經(jīng)有14人回復
- MATLAB安裝后出現(xiàn)一堆java的問題
已經(jīng)有12人回復
wcnmaz
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 212
- 帖子: 131
- 在線: 34.7小時
- 蟲號: 2492252
- 注冊: 2013-06-02
- 專業(yè): 機械測試理論與技術(shù)
|
謝謝回答了 2.047*t了，這一行還是出錯啊 y(i,6)=-178075801.6*y(i,4)*y(i,5)+490910.1296*y(i,3)*y(i,2)+1128974.05+7481.104004*sin(2.047*t); 還有：[t,y]=ode45('ODE45_fun',tspan,y0); 沒有錯誤吧 Warning: Failure at t=5.501826e-001【【【【y0=換不同的初始值，F(xiàn)ailure at t不同】】】】. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.776357e-015) at time t. the smallest value allowed怎么設置大一些？可以犧牲一些準確度 ODE45_main文件: clc;clear tspan=[0,180]; y0=[1e-5;1e-5;1e5;1e-5;1e-5]; [t,y]=ode45('ODE45_fun',tspan,y0); [m,n]=size(y); for i=1:m y(i,6)=-178075801.6*y(i,4)*y(i,5)+490910.1296*y(i,3)*y (i,2)+1128974.05+7481.104004*sin(2.047*t); end data=[t,y]; save ODE45_data.txt data -ascii subplot(2,3,1),plot(t,y(:,1)),title('y(1)') Xlabel('t');Ylabel('y'); subplot(2,3,2),plot(t,y(:,2)),title('y(2)') Xlabel('t');Ylabel('y'); subplot(2,3,3),plot(t,y(:,3)),title('y(3)') Xlabel('t');Ylabel('y'); subplot(2,3,4),plot(t,y(:,4)),title('y(4)') Xlabel('t');Ylabel('y'); subplot(2,3,5),plot(t,y(:,5)),title('y(5)') Xlabel('t');Ylabel('y'); subplot(2,3,6),plot(t,y(:,6)),title('y(6)') Xlabel('t');Ylabel('y'); % plot(t,y(:,1),'bo',t,y(:,2),'rx',t,y(:,3),'gv',t,y(:,4),'r-'); grid on ODE45_fun文件： function dy=ODE45_fun(t,y) dy(1)=-1.918298553*y(3)*y(4)-121.6697369*y(5)*y(2)+0.006472085*y(2)*y(2)+15.25250926*y(5)*y(5)-0.518363603*sin(2.047*t)+0.001124759; dy(2)=0.007229182*y(5)*y(1)-0.013867729*y(3)-0.005151943*sin(2.047*t)+33.43424564*y(4)*y(5)-0.092169794*y(3)*y(2)-0.698266828; dy(3)=72.10986245*y(2)+0.52129529*y(4)*y(1)+0.025471074*y(3)+921.886526*y(4)/y(1)-0.47870471*y(4)+0.025471074*y(1)*t-0.38220722*cos(2.047*t)-4.62279911; dy(4)=-57.37263009*y(5)+0.001053501*y(3)*y(1)+20.28825016/y(1)-0.001053501*y(3)+0.064741605*y(4)+20.28825016*t-0.006201915*sin(2.047*t)-3651.904822; dy(5)=0.017284293*y(4)+0.00278644*y(1)*y(2)-0.551218454*y(2)*y(5)*y(5)+0.010839281*y(2)*y(2)*y(5)+0.020353114*cos(2.047*t)+0.110594984; dy=[dy(1);dy(2);dy(3);dy(4);dy(5)]; |
4樓2015-04-06 14:43:44
amoroso_xmc
新蟲 (初入文壇)
- 應助: 16 (小學生)
- 金幣: 342.5
- 散金: 5
- 帖子: 45
- 在線: 5.8小時
- 蟲號: 2482315
- 注冊: 2013-05-26
- 專業(yè): 流體力學
2樓2015-04-05 20:41:14
3樓2015-04-06 07:06:14
wcnmaz
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 212
- 帖子: 131
- 在線: 34.7小時
- 蟲號: 2492252
- 注冊: 2013-06-02
- 專業(yè): 機械測試理論與技術(shù)
5樓2015-04-06 19:16:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科求調(diào)劑286 +4 | tgs_001 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 11408總分309，一志愿東南大學求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 286分調(diào)劑 +10 | Faune 2026-03-30 | 11/550 |
|
|
|
[考研]
 105500藥學求調(diào)劑，一志愿山東大學藥學，348分
+3
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
|
[考研] 26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 材料專碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13 | 丹青奶蓋 2026-03-26 | 15/750 |
|
|
|
[考研]
一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+11
|
nnnnnnn5 2026-03-25 | 11/550 |
|
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
|
[考研] 312，生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 一志愿天津大學339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|