hqulhh

新蟲 (初入文壇)

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

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

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

回復(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ù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)
330求調(diào)劑已經(jīng)有3人回復(fù)

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

用matlab畫一個(gè)函數(shù)，能幫我看看程序哪里有問題嗎已經(jīng)有8人回復(fù)
求助正態(tài)分布的函數(shù)是什么分布已經(jīng)有8人回復(fù)
matlab求函數(shù)積分出錯(cuò) 已經(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這個(gè)函數(shù)是什么意思？已經(jīng)有8人回復(fù)
如何在matlab里寫.m函數(shù)，老是出錯(cuò) 有誰知道具體格式不已經(jīng)有4人回復(fù)
matlab快速切換圖已經(jīng)有6人回復(fù)
求助一個(gè)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ù)
問個(gè)關(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)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

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

看這個(gè)方程就差不多，其中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)驗(yàn): +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

★
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小時(shí)
蟲號: 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)代替就可以。

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

贊一下

回復(fù)此樓

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

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

csgt0

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

彩色掛圖

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

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

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

你沒發(fā)現(xiàn)你的式子里沒有P么，大概這個(gè)情況，你再改改吧。需要注意不能出現(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小時(shí)
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

5樓: Originally posted by csgt0 at 2012-08-02 17:21:37
你沒發(fā)現(xiàn)你的式子里沒有P么，大概這個(gè)情況，你再改改吧。需要注意不能出現(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 ...

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

贊一下

回復(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小時(shí)
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

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

贊一下

回復(fù)此樓

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

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

csgt0

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

彩色掛圖

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

引用回帖:

7樓: Originally posted by hqulhh at 2012-08-03 16:23:55
我這個(gè)計(jì)算過程中，單位要怎么統(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小時(shí)
蟲號: 1263377
注冊: 2011-04-12
專業(yè): 零件成形制造

引用回帖:

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

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

贊一下

回復(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小時(shí)
蟲號: 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ā)表

[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[教師之家] 焦慮 +8	水冰月月野兔 2026-03-13	12/600	2026-03-18 15:27 by 咪嗚喵嗚
[考研] 化工學(xué)碩306求調(diào)劑 +10	42838695 2026-03-12	10/500	2026-03-18 14:42 by haxia
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +4	pupcoco 2026-03-17	7/350	2026-03-18 12:14 by djl2006
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[基金申請] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 考研調(diào)劑 +4	芬達(dá)46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky