用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')------------------------