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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

wpynwpu

木蟲 (正式寫手)

應(yīng)助: 4 (幼兒園)
金幣: 1469.1
紅花: 1
帖子: 550
在線: 262.5小時
蟲號: 438036
注冊: 2007-09-13
性別: GG
專業(yè): 金屬材料的凝固與結(jié)晶學(xué)

[求助] 如何在一個程序中實(shí)現(xiàn)兩個微分方程的求解（其中第一個方程的解要被第二方程用到）

function wpy531

options=odeset('outputfcn','odeplot');
[t,U]=ode45(@fly1,[0:0.0001:0.005],0.79,options);
U1=U(:,1)
plot(t,U1)
title('Droplet Velocity VS Time')
xlabel('Time/s')
ylabel('Droplet Velocity/m/s')
options=odeset('outputfcn','odeplot');

[t,TT]=ode45(@temperaturel,[0:0.0001:0.005],1173,options);
TT1=TT(:,1)
plot(t,TT1)
title('Droplet Velocity VS Time')
xlabel('Time/s')
ylabel('Droplet Velocity/m/s')

function y=fly1(t,U)
a=1.1616;                %氮?dú)饷芏?br /> b=1.873e-5;             %氮?dú)怵ざ?br /> d=1071e-6;                   %熔滴直徑
c=2385;                   %熔滴密度
g=9.8;                   %重力加速度
f=d*a*U/b;                %Re數(shù)
h=0.28+6/f.^0.5+21/f;    %拖曳系數(shù)
y=((1-a/c)*g-3*(0.28+6/f.^0.5+21/f)*a*U.^2/(4*d*c))/(1+0.5*a/c);

function x=temperaturel(t,TT)% 液相冷卻
rhod=2368;             % rho-d, 熔滴密度，單位為kg*m-3,-0.28*(T-933)
cpd=1086;             %cpd, 熔滴比熱容，單位為Jm-3K-1
zhj=1.102e-3;          %dd,熔滴直徑，單位m
d=0.6;                %熔滴黑度，
f=5.6032e-8;          %Stefan-Boltzmann常數(shù)，單位W*m-3*K-4
g=293;                %T0，環(huán)境溫度，單位K
ug=2.217e-5;             %氬氣粘度
cpg=520.3;                %氬氣熱容
kg=0.01712;             %kg,氮?dú)鉄釋?dǎo)率，單位W*m-1*K-1
rhog=1.784;             %rhog, 氣體密度，單位kg*m-3
%sudu=3.2+7.27*t;                %V，熔滴速度，單位m*s-1
h=kg/zhj*(2.0+0.6*(rhog*zhj*U./ug)^0.5*(ug*cpg/kg)^(1/3))  %h, 熔滴對流傳熱系數(shù)，W/(m2*K)
x=-6*[0.6*(5.6032e-8)*(T^4-7.3701e+009)+h*(T-293)]/(rhod*cpd*zhj);

回復(fù)此樓

» 猜你喜歡

085700資源與環(huán)境308求調(diào)劑已經(jīng)有4人回復(fù)
生物學(xué)調(diào)劑招人�。�！已經(jīng)有4人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有17人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有12人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有13人回復(fù)
288求調(diào)劑已經(jīng)有14人回復(fù)

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

一階微分方程求解，并優(yōu)化參數(shù)值已經(jīng)有14人回復(fù)
如何求解描述振動的二階微分方程已經(jīng)有21人回復(fù)
隱式歐拉法求解一階常微分方程已經(jīng)有7人回復(fù)
求解一個微分方程，謝謝。已經(jīng)有11人回復(fù)
請問如何求解二元一階微分方程組已經(jīng)有24人回復(fù)
一個偏微分方程組（包括邊界條件）如何求解得到自振頻率和振型？已經(jīng)有4人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】如何求解下列偏微分方程已經(jīng)有13人回復(fù)
【求助】一個二階常微分方程的解法（重獎金幣！）已經(jīng)有16人回復(fù)

1樓 2012-05-31 14:49:11

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

change0618

鐵桿木蟲 (著名寫手)

方丈大師

應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時
蟲號: 496517
注冊: 2008-01-19
專業(yè): 化學(xué)反應(yīng)工程

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
臭水溝: 金幣+3, 謝謝交流~~ 2012-05-31 18:58:27
wpynwpu: 金幣+50, ★★★★★最佳答案, 你對matlab很精通呀，羨慕，謝謝! 2012-05-31 23:15:49

CODE:

function wpy531

clear all

clc

[t,U]=ode45(@fly1,[0:0.0001:0.005],0.79);

figure(1)

plot(t,U)

title('Droplet Velocity VS Time')

xlabel('Time/s')

ylabel('Droplet Velocity/m/s')

[t,T]=ode45(@temperaturel,[0:0.0001:0.005],1173,[],U(end));

figure(2)

plot(t,T)

title('Droplet Velocity VS Time')

xlabel('Time/s')

ylabel('Droplet Velocity/m/s')

function y=fly1(t,U)

a=1.1616;                  %氮?dú)饷芏?br />
b=1.873e-5;                %氮?dú)怵ざ?br />
d=1071e-6;                    %熔滴直徑

c=2385;                    %熔滴密度

g=9.8;                     %重力加速度

f=d*a*U/b;                 %Re數(shù)

h=0.28+6/f.^0.5+21/f;      %拖曳系數(shù)

y=((1-a/c)*g-3*(0.28+6/f.^0.5+21/f)*a*U.^2/(4*d*c))/(1+0.5*a/c);

function x=temperaturel(t,T,U)% 液相冷卻

rhod=2368;                % rho-d, 熔滴密度，單位為kg*m-3,-0.28*(T-933)

cpd=1086;                %cpd, 熔滴比熱容，單位為Jm-3K-1

zhj=1.102e-3;            %dd,熔滴直徑，單位m

d=0.6;                  %熔滴黑度，

f=5.6032e-8;             %Stefan-Boltzmann常數(shù)，單位W*m-3*K-4

g=293;                  %T0，環(huán)境溫度，單位K

ug=2.217e-5;              %氬氣粘度

cpg=520.3;                 %氬氣熱容

kg=0.01712;                %kg,氮?dú)鉄釋?dǎo)率，單位W*m-1*K-1

rhog=1.784;               %rhog, 氣體密度，單位kg*m-3

sudu=3.2+7.27*t;                   %V，熔滴速度，單位m*s-1

h=kg/zhj*(2.0+0.6*(rhog*zhj*U/ug)^0.5*(ug*cpg/kg)^(1/3));  %h, 熔滴對流傳熱系數(shù)，W/(m2*K)

x=-6*(0.6*(5.6032e-8)*(T^4-7.3701e+009)+h*(T-293))/(rhod*cpd*zhj);

贊一下

回復(fù)此樓

2樓2012-05-31 18:27:08

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

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