版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

zhaoqian59

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1514
散金: 8
紅花: 5
帖子: 194
在線: 57.4小時
蟲號: 2656296
注冊: 2013-09-15
性別: MM
專業(yè): 爆炸與沖擊動力學(xué)

[求助] 這個程序是不是哪里有問題？初始值跟實驗值是一樣的，但是模擬結(jié)果不對已有1人參與

f(x)
pl=1000;
phyd=110000;
P0=29.39E9;
a0=0.0341;
v=1.5;

x1=x(1);
x2=x(2);
pg=P0*(a0/x1)^(3*v);
y=[x2 (1/pl*(pg-phyd)-3/2*x2^2)/x1]';

Ts = 0.001; % 仿真步長
Tn = 1;    % 仿真終止時間
t  = 0:Ts:Tn; % 仿真時間范圍
N  = length(t);
x=[0.034 0]'; % 系統(tǒng)初值：x=[a a的導(dǎo)數(shù)]
%% Variable Declaration and Initialization
a=0.034*ones(1,N);
a_dt=zeros(1,N);
%x = [a a_d]'; % intial states
%% Simulation Iteration
for n = 1:N
a(n)=x(1);
a_dt(n)=x(2);
k1 = fx(x);
k2 = fx(x+Ts/2*k1);
k3 = fx(x+Ts/2*k2);
k4 = fx( x+Ts *k3);
x = x + Ts/6 * ( k1 + 2*k2 + 2*k3 + k4 );
end
%% Plot Response
plot(t,a)
xlabel('t/(s)')
ylabel('a')
是一個二階非線性常微分方程

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

matlab編程繪圖

» 猜你喜歡

0856調(diào)劑，是學(xué)校就去已經(jīng)有5人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有13人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有11人回復(fù)
復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有11人回復(fù)
321求調(diào)劑已經(jīng)有10人回復(fù)
東華理工大學(xué)化材專業(yè)26屆碩士博士申請已經(jīng)有8人回復(fù)
288求調(diào)劑已經(jīng)有11人回復(fù)
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有5人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有12人回復(fù)

爆炸力學(xué)

1樓 2015-08-11 09:29:52

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

樓主自己編寫的四階龍格庫塔來求解常微分方程？
能給出方程具體形式么，初值，以及各常數(shù)？

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2015-08-11 10:37:11

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

zhaoqian59

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1514
散金: 8
紅花: 5
帖子: 194
在線: 57.4小時
蟲號: 2656296
注冊: 2013-09-15
性別: MM
專業(yè): 爆炸與沖擊動力學(xué)

引用回帖:

2樓: Originally posted by 月只藍 at 2015-08-11 10:37:11
樓主自己編寫的四階龍格庫塔來求解常微分方程？
能給出方程具體形式么，初值，以及各常數(shù)？

aa''+1.5a'=1/1000*(P-Pwater)
P=P0*(a0/a)^v
a是關(guān)于t的函數(shù)，Pwater=1.1*100000，P0=29E9，a0=0.0341，v=3.謝謝你啦~

贊一下

回復(fù)此樓

爆炸力學(xué)

3樓2015-08-11 15:03:26

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

3樓: Originally posted by zhaoqian59 at 2015-08-11 15:03:26
aa''+1.5a'=1/1000*(P-Pwater)
P=P0*(a0/a)^v
a是關(guān)于t的函數(shù)，Pwater=1.1*100000，P0=29E9，a0=0.0341，v=3.謝謝你啦~...

t=0時，對應(yīng)a a'的初值呢？

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

4樓2015-08-11 15:08:12

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

zhaoqian59

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1514
散金: 8
紅花: 5
帖子: 194
在線: 57.4小時
蟲號: 2656296
注冊: 2013-09-15
性別: MM
專業(yè): 爆炸與沖擊動力學(xué)

引用回帖:

4樓: Originally posted by 月只藍 at 2015-08-11 15:08:12
t=0時，對應(yīng)a a'的初值呢？...

t=0時，a(t)=0.0341,a'=0

贊一下

回復(fù)此樓

爆炸力學(xué)

5樓2015-08-11 15:14:44

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

引用回帖:

記 a=u1，a'=u2；
原二階方程可將階為方程組：
u1'=u2;
u2'=1/u1* ( (P-Pwater)/1000 - 1.5*u2 );
其中，P=P0*(a0/u1)^v;

MATLAB代碼如下：

CODE:

function solve_odes_two_oder

clear all;clc

u0=[0.0341 0];

tspan=linspace(0,0.1,100);

[t u]=ode45(@odefun,tspan,u0);

figure(1)

plot(t,u(:,1),'bo--',t,u(:,2),'r-*'),legend('a',' da/dt ')

function f=odefun(t,u)

Pwater=1.1*100000;

P0=29E9;

a0=0.0341;

v=3;

P=P0*(a0/u(1))^v;

f(1)=u(2);

f(2)=1/u(1)*   (  (P-Pwater)/1000 - 1.5*u(2)     );

f=f';

結(jié)果中，a和a'迅速增大，即便在0.1的時間之內(nèi)，見附圖1。

附圖1.png

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

6樓2015-08-11 15:47:29

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

zhaoqian59

捐助貴賓 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1514
散金: 8
紅花: 5
帖子: 194
在線: 57.4小時
蟲號: 2656296
注冊: 2013-09-15
性別: MM
專業(yè): 爆炸與沖擊動力學(xué)

引用回帖:

6樓: Originally posted by 月只藍 at 2015-08-11 15:47:29
記 a=u1，a'=u2；
原二階方程可將階為方程組：
u1'=u2;
u2'=1/u1* ( (P-Pwater)/1000 - 1.5*u2 );
其中，P=P0*(a0/u1)^v;

MATLAB代碼如下：

function solve_odes_two_oder
clear all;clc
u0=;
...

復(fù)制過去之后運行不出來.....

贊一下

回復(fù)此樓

爆炸力學(xué)

7樓2015-08-12 09:40:23

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhaoqian59: 金幣+50 2015-08-12 12:40:35

引用回帖:

7樓: Originally posted by zhaoqian59 at 2015-08-12 09:40:23
復(fù)制過去之后運行不出來........

代碼完整復(fù)制到一個新建的m文件，運行即可，不要在主程序窗口運行。

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

8樓2015-08-12 10:17:04

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 zhaoqian59 的主題更新

返回列表

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

[考研] 271材料工程求調(diào)劑 +6	.6lL 2026-03-18	6/300	2026-03-19 15:41 by haoshis
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7	步川酷紫123 2026-03-13	7/350	2026-03-18 17:12 by 盡舜堯1
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 293求調(diào)劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、

24小時熱門版塊排行榜

[求助] 這個程序是不是哪里有問題？初始值跟實驗值是一樣的，但是模擬結(jié)果不對 已有1人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] 這個程序是不是哪里有問題？初始值跟實驗值是一樣的，但是模擬結(jié)果不對已有1人參與