cwlr.net
当前位置:首页 >> 求MATLAB编程实现微分画图 >>

求MATLAB编程实现微分画图

把这部分保存为M函数,注意文件名一定是:Lemaitrechaboche function dy=Lemaitrechaboche(t,y)r=13.19;A=1209;x=620;a0=13.2478;a1=0.00007865;a2=-0.0001924;z=733.25;k=a0+a1*(x-z)+a2*(x-z)^2;dy=[0;(x/A)^r*(1-y(2))^(-k)];end ————————分...

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

M=[m1 0 0 0 0;0 m2 0 0 0;0 0 m3 0 0;0 0 0 m4 0;0 0 0 0 m5] 这里多了个0

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

由于方程比较复杂,解析解不能用初等函数表示 只是要获得图像的话,用数值计算的办法可能更方便 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', '...

此微分方程组可以用数值的方法求得其数值解。用ode45()函数,得到微分方程组解的图形。 执行代码和微分方程组定义函数func(t,y)

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

1定义函数: function y=fun(t,x) y=zeros(3,1); x1=x(1);x2=x(2);x3=x(3); y(1)=x1*(1-x1/150000-0.5*x2/30000); y(2)=x2*(-1-x2/30000+0.5*x1/150000-2*x3/2200); y(3)=x3*(-1-x3/2200+2*x2/30000); 2解方程并画图 [t,x]=ode23(@fun,[0,1],[12...

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