matlab解微分方程dydt=(cot(alfa)+tan(glma)*R)/(y*(R-y*tan(alfa)))*(K_ab*K*((P_a^(n-1))/(P_t^n))*Vc*y*sigma+K_ad*exp(a*T)*Vc*sigma+K_di*sqrt(Vc*y)*exp(-K_Q/(T+273)));alfa=15*pi/180;glma=15*pi/180;R=0.8;K_ab=0.0295;K_ad=1.4761E-14;a=9.0313E-4;K_di=
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/18 23:28:25
![matlab解微分方程dydt=(cot(alfa)+tan(glma)*R)/(y*(R-y*tan(alfa)))*(K_ab*K*((P_a^(n-1))/(P_t^n))*Vc*y*sigma+K_ad*exp(a*T)*Vc*sigma+K_di*sqrt(Vc*y)*exp(-K_Q/(T+273)));alfa=15*pi/180;glma=15*pi/180;R=0.8;K_ab=0.0295;K_ad=1.4761E-14;a=9.0313E-4;K_di=](/uploads/image/z/14897557-37-7.jpg?t=matlab%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8Bdydt%3D%28cot%28alfa%29%2Btan%28glma%29%2AR%29%2F%28y%2A%28R-y%2Atan%28alfa%29%29%29%2A%28K_ab%2AK%2A%28%28P_a%5E%28n-1%29%29%2F%28P_t%5En%29%29%2AVc%2Ay%2Asigma%2BK_ad%2Aexp%28a%2AT%29%2AVc%2Asigma%2BK_di%2Asqrt%28Vc%2Ay%29%2Aexp%28-K_Q%2F%28T%2B273%29%29%29%3Balfa%3D15%2Api%2F180%3Bglma%3D15%2Api%2F180%3BR%3D0.8%3BK_ab%3D0.0295%3BK_ad%3D1.4761E-14%3Ba%3D9.0313E-4%3BK_di%3D)
matlab解微分方程dydt=(cot(alfa)+tan(glma)*R)/(y*(R-y*tan(alfa)))*(K_ab*K*((P_a^(n-1))/(P_t^n))*Vc*y*sigma+K_ad*exp(a*T)*Vc*sigma+K_di*sqrt(Vc*y)*exp(-K_Q/(T+273)));alfa=15*pi/180;glma=15*pi/180;R=0.8;K_ab=0.0295;K_ad=1.4761E-14;a=9.0313E-4;K_di=
matlab解微分方程
dydt=(cot(alfa)+tan(glma)*R)/(y*(R-y*tan(alfa)))*(K_ab*K*((P_a^(n-1))/(P_t^n))*Vc*y*sigma+K_ad*exp(a*T)*Vc*sigma+K_di*sqrt(Vc*y)*exp(-K_Q/(T+273)));
alfa=15*pi/180;
glma=15*pi/180;
R=0.8;
K_ab=0.0295;
K_ad=1.4761E-14;
a=9.0313E-4;
K_di=5.7204E6;
T=800;
sigma=3E8;
K=0.416;
n=7;
P_a=1800;
P_t=1600;
Vc=50/60;
K_Q=20460;
左边为dy/dt
结果能表示成y关于t的函数最好了,或者画成(y-t)图像
matlab解微分方程dydt=(cot(alfa)+tan(glma)*R)/(y*(R-y*tan(alfa)))*(K_ab*K*((P_a^(n-1))/(P_t^n))*Vc*y*sigma+K_ad*exp(a*T)*Vc*sigma+K_di*sqrt(Vc*y)*exp(-K_Q/(T+273)));alfa=15*pi/180;glma=15*pi/180;R=0.8;K_ab=0.0295;K_ad=1.4761E-14;a=9.0313E-4;K_di=
这个方程是解不出解析解的.
我用ODE45给您算了一下.
先赋个初值y(0)=10000
编写m文件
function f=exam2fun (t,y)
f=3.9464/(y*(0.8-y*0.2679))*(3.8873e3*y+7.6004e-006+0.0299*sqrt(0.8333*y));
f=f(:);
再在主窗口输入
[t,y]=ode45('exam2fun',[0:1:1000],10000)
plot(t,y)
可以求得0到1000内的y(t)
最后有一部分为虚数,若增大初值,虚数分界线变大
得到的图像如图,最后一部分虚数没有显示.
您可以选择自己的初值,然后选择范围求解并画图.