ntjm.net
当前位置:首页 >> mAtlAB解微分方程画图 >>

mAtlAB解微分方程画图

用matlab求解微分方程组并画图,其实现过程:1、首先应根据已知微分方程组,编写其函数 dy(1)=-0.10982*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3)); dy(2)=1441.5074*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3));2、用ode45函数求解其微分方程组的数值解 [t,y]=ode45(@ode_fun,tspan,y0)3、然后用plot绘图函数,绘出t-Cg和t-θ关系曲线图4、编程后运行结果

最低0.27元开通文库会员,查看完整内容> 原发布者:加菲鼠go 用MATLAB求解微分方程1.微分方程的解析解求微分方程(组)的解析解命令:dsolve('方程1','方程2',…'方程n','初始条件','自变量')记号:在表达微分方程时,用

第一个是线性方程,可以直接求解析解 用dsolve求数值解:%%%%%%%%%%%%%%%%M-file 以下要单存一个文件在工作路径下%%%%%%%%%%%%%function dotx=myfun(t,y)qin=50;q=50;Cain=1.3;V=1000;K=0.05;dotx=(qin*Cain-q*y)/V-K

能不能把你的方程组重新整理一下,我看不明白. 你给出的是微分方程吗,我还是看不懂! v(1)=300 v(2)=100 v(3)=0 以上上个都是常数 则可以知道: dv(1)/dt=0 dv(2)/dt=0 dv(3)/dt=0 怎么又会有后面的: dv(2)/dt=100v(3)x dv(3)/dt=-100v(2)x 呢

初始条件y(t)=0是不是有问题啊,如果是这样,直接解出y=0,是不是y(0)=0啊? clc;clear y=dsolve('d2y=g-k*s*v/m','dy(0)=0,y(0)=0') y1=diff(y,'t') v=solve(y1-'v','v') y=subs(y) 结果: y = 1/2*(g-k*s*v/m)*t^2 y1 = (g-k*s*v/m)*t v = g*t/(k*s*t+m)*m y = 1/2*(g-k*s*g*t/(k*s*t+m))*t^2

该二阶微分方程可以ode()函数求得,其方程的数值解,然后用plot()函数绘制其xy、xdy/dx函数图形,以及ydy/dx函数的相平面图形.主要执行代码

syms x y t;s=dsolve('dx=-10*x+10*y','dy=28*x-y-x*z','dz=-8/3*z+x*y','x(0)=0','y(0)=0','z(0)=1','t');x=s.xy=s.yz=s.zhold on;ezplot(x);ezplot(y);ezplot(z);%祝楼主学有所成 嘿嘿

如何用matlab求解微分方程并画图 dsolve()求解微分方程后,使用ezplot()作图,参考代码:12>> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间

||把T换成抄x,把alpha换成袭y>> syms x y>> z=dsolve('Dy=1.44*(10^知9)*(1-y)*exp(-109170/(8.314*x))')z =C2*exp(-1440000000*t*exp(-131308636035602598027423622804907/(10000000000000000000000000000*x))) + 1>> pretty(z) / /

function dz=caonim(t,z) %保存为 caonim.m 文件dz=[-z(1)^3-z(2);z(1)-z(2)^3];------------------------------------运行[T,Z]=ode45('caonim',[0 30],[1,0.5]);plot(T,Z)legend('x-t','y-t')------------------------

网站首页 | 网站地图
All rights reserved Powered by www.ntjm.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com