版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

wpynwpu

木蟲 (正式寫手)

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

[求助] 如何在一個程序中實現(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;                %氮氣密度
b=1.873e-5;             %氮氣黏度
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,氮氣熱導率，單位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);

回復此樓

» 猜你喜歡

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

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

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

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

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

change0618

鐵桿木蟲 (著名寫手)

方丈大師

應(yīng)助: 44 (小學生)
金幣: 17724.5
紅花: 17
帖子: 2413
在線: 546.7小時
蟲號: 496517
注冊: 2008-01-19
專業(yè): 化學反應(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;                  %氮氣密度

b=1.873e-5;                %氮氣黏度

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,氮氣熱導率，單位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);

贊一下