clzm.net
当前位置:首页 >> mAtlAB解微分方程 >>

mAtlAB解微分方程

dsolve('D3y+D2y/x-Dy/(x^2)=x','y(500)=0,Dy(500)=0','x'),用这个能够得到y的一个含参数的解,然后根据x=0时,y为有限值来确定参数的值,之后就能解出相应的y’’ 和y’ ,后面那个貌似用MATLAB不太好实现,最好是自己来分析。

这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数 就一个例子(选自matlab帮助文档) 解入下微分方程组 1,建立方程组函数 function dy = rigid(t,y)dy = zeros(3,1); % a column vectordy(1) = y(2) * y(3);dy(2) = -y(1) * y(3);d...

对于相对简单的微分方程或微分方程组,可以dsolve()函数得到解析值(精确解)。 例①: >>syms x(t) a >>dsolve(diff(x) == -a*x) ans = C2*exp(-a*t) >>dsolve(diff(x) == -a*x, x(0) == 1) ans = exp(-a*t) 例②: syms y(t) a >> Dy = diff(y...

你这个程序明显不对。没看懂是什么意思。跟你说下错在哪吧: 你看下diff()的帮助,返回值是比输入的向量少一个元素的,你输入的是个标量(可以理解为长度为1),返回值当然是长度为0(即空元素【 】)。matlab中任何数与空【】计算得到空【】。...

function dx=differential(t,x) dx = zeros(2,1);; dx(1)=x(1)*(1-x(1))*30; dx(2)=x(2)*(1-x(2))*(160*x(1)-20); end clc options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-5]); [T,Y]=ode45(@differential,[0 5],[0.5 0.5],options); figur...

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

(1)用matlab求微分方程初值问题的符号解: syms y(x) Dy=diff(y,1);D2y=diff(y,2); y=dsolve(D2y+4*Dy+29*y==0,y(0)==0,Dy(0)==15) (2)用matlab求微分方程初值问题的数值解: x0=[0 15]; [x,y]=ode45(@func,[0 20],x0) func——微分方程自定义...

方程里面使用了中文符号 >> dsolve('Dp^2=15*(7-x)/((7-x)^2+4)^0.5/2','p(0)=0','x') ans = 1/2/(53-14*x+x^2)^(1/2)*(-30*(53-14*x+x^2)^(1/2)*(-7+x))^(1/2)*t -1/2/(53-14*x+x^2)^(1/2)*(-30*(53-14*x+x^2)^(1/2)*(-7+x))^(1/2)*t

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

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

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