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

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

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) 使用这样的方法求解,但从结果看好像你的方程有问题!

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

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

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

由于方程比较复杂,解析解不能用初等函数表示 只是要获得图像的话,用数值计算的办法可能更方便 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()函数,求出微分方程的解。 dsolve('Dp=0.2*p*(1-p/5)-0.21','p(0)=0') 其二,用plot()函数,绘出t~p(t)关系图 plot(t,p) 实现过程如下图:

看标题以为你要求微分方程呐,结果是画dR/dr vs. r % 画出图中的公式% 定义微分方程函数dRdr = @(r) 0.89 ./ r .* exp(-(log(r) + 0.84).^2 / 0.086);% 在(0, 10]上画图r = linspace(0.01, 3, 500);dR = dRdr(r);figure(10);plot(r, dR, '^b', '...

参考代码: m=10600; k=31986497; c=8.9E4;dy = @(t,y)[y(2); (k*y(1)-c*y(2)/abs(y(2))*y(2)^2)/m];ode45(dy,[0 5],[0 3.66]) 其中,蓝线为x,绿线为x'。

I use Matlab code to reply your question. You should run successfully by pasting it into your code file. In order to use built-in function of ODE45, I first off set a function file called 'myode.m' function [xp]=myode(t,x) xp=-...

matlab code: clear allclcf=@(t,y)([-y(2)+y(1)*(1-y(1)^2-y(2)^2); y(1)+y(2)*(1-y(1)^2-y(2)^2)]);[t y]=ode45(f,0:0.1:10,[5 5]);plot(t,y(:,1),t,y(:,2),y(:,2),y(:,1))legend('x-t','y-t','x-y')

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