wcnmaz

應(yīng)助: 0 (幼兒園)
金幣: 212
帖子: 131
在線: 34.7小時(shí)
蟲(chóng)號(hào): 2492252
注冊(cè): 2013-06-02
專業(yè): 機(jī)械測(cè)試?yán)碚撆c技術(shù)

[交流] MATLAB7.0出錯(cuò)了？已有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程序運(yùn)行到上面一行出錯(cuò)了？感覺(jué)也沒(méi)什么錯(cuò)誤啊
怎么改，謝謝了

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)];

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

MS7.0安裝問(wèn)題已經(jīng)有5人回復(fù)
matlab 運(yùn)行ode45出錯(cuò)，不知道什么原因已經(jīng)有4人回復(fù)
Matlab 7.0 基礎(chǔ)教程清華大學(xué) 已經(jīng)有29人回復(fù)
matlab一運(yùn)行就出現(xiàn)的問(wèn)題已經(jīng)有7人回復(fù)
運(yùn)行matlab出現(xiàn)錯(cuò)誤已經(jīng)有4人回復(fù)
matlab程序出錯(cuò)調(diào)試已經(jīng)有7人回復(fù)
Quantum-ESPRESSO4.3安裝問(wèn)題，求助！已經(jīng)有4人回復(fù)
問(wèn)下這個(gè)問(wèn)題怎么解決.....真心啊已經(jīng)有4人回復(fù)
matlab程序問(wèn)題，總是出錯(cuò)，麻煩大家看看問(wèn)題在那里？已經(jīng)有6人回復(fù)
xcysden-1.5.53.tar.gz安裝問(wèn)題求助已經(jīng)有5人回復(fù)
求助 matlab交點(diǎn)問(wèn)題已經(jīng)有19人回復(fù)
一筆畫(huà)問(wèn)題matlab實(shí)現(xiàn) 已經(jīng)有5人回復(fù)
matlab出錯(cuò) 求助高手解答已經(jīng)有5人回復(fù)
采用matlab中nlinfit擬合，出現(xiàn)問(wèn)題了，不知道是哪錯(cuò)了，該怎么解決，求高手指點(diǎn)一下已經(jīng)有9人回復(fù)
【求助】matlab 中，幾個(gè)m文件調(diào)用已經(jīng)有14人回復(fù)
MATLAB安裝后出現(xiàn)一堆java的問(wèn)題已經(jīng)有12人回復(fù)

1樓 2015-04-05 15:00:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

amoroso_xmc

新蟲(chóng) (初入文壇)

應(yīng)助: 16 (小學(xué)生)
金幣: 342.5
散金: 5
帖子: 45
在線: 5.8小時(shí)
蟲(chóng)號(hào): 2482315
注冊(cè): 2013-05-26
專業(yè): 流體力學(xué)

★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖

2.047和t之間漏了乘號(hào)*。
或者t不是一個(gè)標(biāo)量？

贊一下

回復(fù)此樓

2樓2015-04-05 20:41:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wyjjf

新蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 294.1
散金: 46
紅花: 3
帖子: 742
在線: 160.5小時(shí)
蟲(chóng)號(hào): 2588385
注冊(cè): 2013-08-09
專業(yè): 海洋工程

★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

2樓: Originally posted by amoroso_xmc at 2015-04-05 20:41:14
2.047和t之間漏了乘號(hào)*。
或者t不是一個(gè)標(biāo)量？

標(biāo)量？

回復(fù)此樓

3樓2015-04-06 07:06:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wcnmaz

新蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 212
帖子: 131
在線: 34.7小時(shí)
蟲(chóng)號(hào): 2492252
注冊(cè): 2013-06-02
專業(yè): 機(jī)械測(cè)試?yán)碚撆c技術(shù)

引用回帖:

2樓: Originally posted by amoroso_xmc at 2015-04-05 20:41:14
2.047和t之間漏了乘號(hào)*。
或者t不是一個(gè)標(biāo)量？

謝謝回答了
2.047*t了，這一行還是出錯(cuò)啊
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); 沒(méi)有錯(cuò)誤吧

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怎么設(shè)置大一些？可以犧牲一些準(zhǔn)確度

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)];

贊一下

回復(fù)此樓

4樓2015-04-06 14:43:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wcnmaz

新蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 212
帖子: 131
在線: 34.7小時(shí)
蟲(chóng)號(hào): 2492252
注冊(cè): 2013-06-02
專業(yè): 機(jī)械測(cè)試?yán)碚撆c技術(shù)

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);
這一行應(yīng)該沒(méi)有錯(cuò)誤了吧？？？

贊一下

回復(fù)此樓

5樓2015-04-06 19:16:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	10/500	2026-03-18 20:41 by Wangjingyue
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz