lzth.net
当前位置:首页 >> 如何用mAtlAB求解微分方程并画图 >>

如何用mAtlAB求解微分方程并画图

function zd0412 function dy=odefun(t,y) dy=[y(2); 3*sin(2*t)+exp(t)*y(1)-t*y(2)]; end tspan=[0,2]; y0=[1;-1] [t,y]=ode45(@odefun,tspan,y0); plot(t,y(:,1)) end

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

syms t v=dsolve('Dv=(190.708-90.64*v^2)/47.27','v(0)=0','t'); t=0:0.00001:0.002; v=eval(v); plot(t,v) 使用这样的方法求解,但从结果看好像你的方程有问题!

说明 求解常微分方程一般使用ode系列函数(其中最常用的是ode45),基本做法就是把微分方程用函数表示出来,然后调用ode函数求解即可。对于本题而言,使用匿名函数描述微分方程更为方便一些。 参考代码 k0=8; k1=1; K=1; k2=5; n=2; % 常数定义d...

如何用matlab求解微分方程并画图,可以先用dsolve()或ode()求出其微分方程(组)的解析解或数值解,然后用plot()绘制其图形。 例如:解微分方程 y'=y-2t/y,y(0)=1,0

function zd0412 function dy=odefun(t,y) dy=[y(2); 3*sin(2*t)+exp(t)*y(1)-t*y(2)]; end tspan=[0,2]; y0=[1;-1] [t,y]=ode45(@odefun,tspan,y0); plot(t,y(:,1)) end

由于方程比较复杂,解析解不能用初等函数表示 只是要获得图像的话,用数值计算的办法可能更方便 fun=@(t,y) 1.44*(10^9)*(1-y).*exp(-109170./(8.314*t)); [T,Y] = ode23t(fun,500:600,0); DY=fun(T,Y); AX=plotyy(T,Y,T,DY); set(get(AX(1),'Xl...

求解微分方程的函数是dsolve()。如要画微分方程的解函数,可以用plot()函数。例如: 求解dy/dx=1+y^2,并且当y(0)=1的解,画出y=f(x)的图形。 y=dsolve('Dy=1+y^2','y(0)=1','x') x=-2*pi:pi/10:2*pi; y1=eval(y); plot(x,y1,'k-') xlabel('...

说明:微分方程dy/dx-2y/x+1=(x+1)^5/2 缺 初始条件。 如此微分方程的初始条件为 y(1)=0,则你问题可以这样来解决。 1、用dsolve()函数,求出微分方程的解析解。 y=dsolve(Dy-2*y/x+1==(x+1)^5/2,y(1)==0) 2、用plot()函数,绘出y(x)的图形。 ...

只需要提取x和y的前n个数字再plot,但必须保证x和y维数一致,下面是一个小例子 x=1:0.1:4; y=4:0.1:7; x1=x(1,1:10); y1=y(1,1:10); subplot(2,1,1) plot(x,y) subplot(2,1,2) plot(x1,y1) plot(x(1,1:10),y(1,1:10))

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