lzth.net
当前位置:首页 >> 用MATLAB求微分方程Dy/Dx%2y/x+1=(x+1)∧5/2,并... >>

用MATLAB求微分方程Dy/Dx%2y/x+1=(x+1)∧5/2,并...

说明:微分方程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)的图形。 ...

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

微分方程的话可以用dsolve函数,格式可以总结为dsolve('equation','condition',自变量)的形式。自变量默认为t。 这种方法解微分方程还好。如果是微分方程组的话建议用ode45函数(该方法也用于微分方程) 首先建立M文件写方程组的表达形式 functi...

此题有误,两个方程的左侧应是dx/dt,dy/dt [x y]=dsolve('Dx=4*x-2*y', 'Dy=2*x-y','t') x = 2*C2*exp(3*t)+1/2*C1 y = C1+C2*exp(3*t) [x y]=dsolve('Dy=4*x-2*y', 'Dx=2*x-y','t') x = 1/4*C1+1/2*C1*t+1/2*C2 y = C1*t+C2

>> clear >> syms a b c d e; >> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e') y = (2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(...

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可以解这类延迟微分方程,但是是数值解法;所以需要之到一个初始条件 x(0)的值; 你能给出x(0)的值我可以帮你解 首先编写关于延迟函数的M文件; function dx=yanchi(t,x,z) tau=z;%定义延迟时间 dx=x*(1-tau);%延...

1、对于解析值,你可以用dsolve()函数求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解 >> syms y(x) >>Dy = diff(y);D2y = diff(y, 2); >>dsolve(x*D2y+x*(Dy)^2-Dy==0,'x') 运行结果 2、对于数值解,你可以查看二阶常微分方程这个例子,http...

先变形为dy/dx=y/x-x,再用dsolve求通解或ode45求数值解。如: syms y(x) y=dsolve(diff(y)==y/x-x) 结果是: y = - x^2 + C1*x

对于一般的微分方程用matlab的dsolve()函数命令求得其解析值(即函数表达式)。例如 1、求的解析值。 syms y(t) a Dy = diff(y); D2y = diff(y,2); dsolve(D2y == -a^2*y, y(0) == 1, Dy(pi/a) == 0) 运行结果 y=exp(-a*t*1i)/2 + exp(a*t*1i)...

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