hqulhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 220.3
紅花: 1
帖子: 26
在線: 40.9小時
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

[求助] MATLAB函數(shù)轉(zhuǎn)化問題。請教

請教下面公式如何在MATLAB實(shí)現(xiàn)的程序

回復(fù)此樓

» 猜你喜歡

材料專碩英一數(shù)二306 已經(jīng)有5人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有9人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有12人回復(fù)
0703化學(xué) 305求調(diào)劑已經(jīng)有4人回復(fù)
0703化學(xué)調(diào)劑，求各位老師收留已經(jīng)有10人回復(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ù)

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

用matlab畫一個函數(shù)，能幫我看看程序哪里有問題嗎已經(jīng)有8人回復(fù)
求助正態(tài)分布的函數(shù)是什么分布已經(jīng)有8人回復(fù)
matlab求函數(shù)積分出錯已經(jīng)有12人回復(fù)
使用lsqnonlin函數(shù)優(yōu)化動力學(xué)參數(shù)，總是得不到合理的結(jié)果已經(jīng)有6人回復(fù)
請教大家關(guān)于matlab初值的問題已經(jīng)有4人回復(fù)
COMSOL4.2與MATLAB2011b的連接問題？求助已經(jīng)有9人回復(fù)
請教matlab用數(shù)組給參數(shù)矩陣賦值的問題已經(jīng)有11人回復(fù)
請教MATLAB中的LMI求解問題，急�。�！已經(jīng)有7人回復(fù)
Matlab這個函數(shù)是什么意思？已經(jīng)有8人回復(fù)
如何在matlab里寫.m函數(shù)，老是出錯有誰知道具體格式不已經(jīng)有4人回復(fù)
matlab快速切換圖已經(jīng)有6人回復(fù)
求助一個matlab關(guān)于分類討論的問題已經(jīng)有9人回復(fù)
matlab 數(shù)組求導(dǎo) 已經(jīng)有3人回復(fù)
scanf()函數(shù)中讀入字符串問題已經(jīng)有4人回復(fù)
matlab里面怎么實(shí)現(xiàn)結(jié)構(gòu)體數(shù)組的可視化？？？？？？已經(jīng)有13人回復(fù)
matlab的函數(shù)調(diào)用問題已經(jīng)有3人回復(fù)
問個關(guān)于matlab里find函數(shù)的問題！已經(jīng)有4人回復(fù)
matlab生成隨機(jī)行走網(wǎng)格模擬材料生長已經(jīng)有26人回復(fù)
【求助】matlab生成旋轉(zhuǎn)體？已經(jīng)有9人回復(fù)
MATLAB安裝后出現(xiàn)一堆java的問題已經(jīng)有12人回復(fù)
matlab做函數(shù)圖象的問題已經(jīng)有6人回復(fù)
【討論】WMSN能量仿真NS2還是MATLAB 已經(jīng)有10人回復(fù)

1樓 2012-08-02 15:29:04

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

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

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
hqulhh: 金幣+2, ★★★很有幫助 2012-08-02 16:54:13
jjdg: 金幣+1, 感謝參與 2012-08-02 18:25:43

看這個方程就差不多，其中y(1)對應(yīng)你的L，y(2)對應(yīng)G

CODE:

[t,l]=ode45(@fun,[0 1],ones(2,1))

function dy=fun(t,y)

dy=zeros(2,1);

dy(1)=y(1)*t*y(2)^2;

dy(2)=5/y(2);

end

贊一下

回復(fù)此樓

showmethemoney

2樓2012-08-02 16:29:29

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

★
jjdg: 金幣+1, 感謝參與 2012-08-02 18:25:33

其他的參數(shù)都有，所以方程中所有L都用y(1), dL/dt用dy(1), G都用y(2),dG/dt用dy(2)代替就可以。

贊一下

回復(fù)此樓

showmethemoney

3樓2012-08-02 16:32:27

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

hqulhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 220.3
紅花: 1
帖子: 26
在線: 40.9小時
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

3樓: Originally posted by csgt0 at 2012-08-02 16:32:27
其他的參數(shù)都有，所以方程中所有L都用y(1), dL/dt用dy(1), G都用y(2),dG/dt用dy(2)代替就可以。

還有個Ρ是時間t的函數(shù) Ρ=kt

贊一下

回復(fù)此樓

4樓2012-08-02 16:53:28

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

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

★ ★ ★ ★ ★ ★ ★ ★
hqulhh: 金幣+5, ★★★很有幫助 2012-08-02 17:31:08
jjdg: 金幣+3, 感謝參與 2012-08-02 18:25:23

你沒發(fā)現(xiàn)你的式子里沒有P么，大概這個情況，你再改改吧。需要注意不能出現(xiàn)讓K=0

CODE:

L0=1;

[t,y]=ode45(@fun,0:100,[L0;1]);

L=y(:,1);

G=y(:,2);

delta=(L-L0)/L0

end

function dy=fun(t,y)

L0=1;

alpha=0.01;

deltaT=10;

C=1;

r=1;

A=1;

B=1;

v=1;

D=1;

Q=1;

R=8.314;

rou=0.5+0.005*t;

sigma=C*rou^2/r;

T=273+10*t;

dy=zeros(2,1);

dy(2)=D*exp(-Q/(R*T))/y(2);

yita=A*T*y(2)^3*rou^2*exp(B/T);

K=yita/3/(1-2*v);

dy(1)=y(1)*(alpha*deltaT-sigma/(3*K));

end

贊一下

回復(fù)此樓

showmethemoney

5樓2012-08-02 17:21:37

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

hqulhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 220.3
紅花: 1
帖子: 26
在線: 40.9小時
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

5樓: Originally posted by csgt0 at 2012-08-02 17:21:37
你沒發(fā)現(xiàn)你的式子里沒有P么，大概這個情況，你再改改吧。需要注意不能出現(xiàn)讓K=0

L0=1;
=ode45(@fun,0:100,);
L=y(:,1);
G=y(:,2);
delta=(L-L0)/L0
end

function dy=fun(t,y)
L0=1;
alpha=0.01;
delt ...

非常感謝，我講的那個P就是那個相對密度

贊一下

回復(fù)此樓

6樓2012-08-02 17:30:35

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

hqulhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 220.3
紅花: 1
帖子: 26
在線: 40.9小時
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

我這個計算過程中，單位要怎么統(tǒng)一，在MATLAB中？謝啦。

贊一下

回復(fù)此樓

7樓2012-08-03 16:23:55

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

引用回帖:

7樓: Originally posted by hqulhh at 2012-08-03 16:23:55
我這個計算過程中，單位要怎么統(tǒng)一，在MATLAB中？謝啦。...

matlab不管單位，自己先換算好

贊一下

回復(fù)此樓

showmethemoney

8樓2012-08-03 17:02:49

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

hqulhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 220.3
紅花: 1
帖子: 26
在線: 40.9小時
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

8樓: Originally posted by csgt0 at 2012-08-03 17:02:49
matlab不管單位，自己先換算好...

G要不要賦給初值，是不是這個式中[t,y]=ode45(@fun,0:100,[L0;1]); 1代表的那個

贊一下

回復(fù)此樓

9樓2012-08-03 20:19:05

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

hqulhh

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 220.3
紅花: 1
帖子: 26
在線: 40.9小時
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

用plot畫出G值不變，假設(shè)G初值為10，怎么回事？
假設(shè)各初值L0=100;alpha=14.6;deltaT=10;C=2;r=16;A=8;B=2;
v=0.28;D=1e-12;Q=315;R=8.31;
幫我看下。謝啦

贊一下

回復(fù)此樓

10樓2012-08-03 21:12:30

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

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

[考研] 271材料工程求調(diào)劑 +5	.6lL 2026-03-18	5/250	2026-03-19 03:07 by 無懈可擊111
[考研] 一志愿天大材料與化工（085600）總分338 +4	蔡大美女 2026-03-13	4/200	2026-03-18 22:47 by walc
[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 307求調(diào)劑 +3	冷笙123 2026-03-17	3/150	2026-03-18 09:55 by macy2011
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky