微分方程数值解法

dsolve('x^2*D2y+4*x*Dy+2*y=0','y(1)=2','Dy(1)=-3','x') ans 

结果：代码：clearallclcf=@(x,y)([y(2);   0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]

请参考ode45函数的用法将所求值设为终值即可

……一定要用matlab吗,只会mathematica的路过.而且啊,你一面说要求方程的数值解,一面又只给一个全是字母系数的方程,你是要怎样……再问：呃，abcdef都是常量，可以随便假定的，就赋1、

当然不是了,首先解齐议程对应的特征方程r^2-r+1=0r=(1±√3i)/2所以齐次通解是y=e^(1/2x)(C1cos√3x+C2sin√3x)特解可能观察得得y=a因此非齐次通解为y=e^(1

你这个属于边值问题.得用打靶算法.如果想简单,你令一个变量为x=y-2001,就变成正常的微分方程.

matlab里面常使用龙格库塔方法求解常微分方程组,命令是ode45,还有其他一些函数,但是最常用的是ode45,lz可以help一下,很简单的,另外给你一个文档,讲的还是比较详细,希望可以帮到你ht

这个方程可以求符号解啊,虽然答案里用到了特殊函数,但是稍微带入下具体数值求个数值解也是可以的：eqn=r'[t]^2==cr[t]^2+2k^2r[t]-h^2/r[t];DSolve[eqn,r[t

不知道LZ的教材是哪一本,是不是高教出的王高雄,周之铭版的,常微的教材很多,你可以多参考几本：吉林大学出版社,王怀柔,伍卓群版的不错或者东北师大版的也不错但个人认为包络和奇解这个方向好像没什么深入的必

首先建立vdp1.m文件functiondx=vdp1(t,x)dx=zeros(2,1);dx(1)=x(2);dx(2)=(10000*x+(96.321*exp(-t/378.47108)+28

引进数值解法其实是迫不得已.如所说,并不是所有的微分方程都可以求出具体的解.我们所学的微分方程之所以都可以写出解是因为这些都是非常经典性质很好的微分方程.但对于工程实际上的微分方程就远没那么好了,各种

用幂级数法：设y=c0+c1x+c2x^2+...+cnx^n+...则y'=c1+2c2x+3c3x^2+...+ncnx^(n-1)y"=2c2+6c3x+12c4x^2+...+n(n-1)cn

这个方程最简单的方法是用shootingmethod.可wiki之.主要就假设一个y(a1),然后看,y'(a2)是否等于b2.可以用matlab写两个小函数,一个用在ode45（fun1）,然后y(

设u=积分(0~x)ydx,那么原来方程就是u''=-0.04u'^2-sin(u)+0,44u''(0)=0,u'=3,u(0)=0(积分上下界都为0)再问：����������룬��ʤ�м����

改动：functionS=history(t)S=ones(3,1);改为S=ones(4,1);OVER!PS.多翻翻help!

因为你x=0时2/x是无穷大呀,然后y'又是0,然后(2/x)y'就是nan了,所以后面算的全是nan了.

见图

dsolve('D2y-y*sin(x)=0,y(0)=1,Dy(0)=0','x')Warning:Explicitsolutioncouldnotbefound.>Indsolveat101ans

dy/dt+y+y²=2dy/dt=-y²-y+2dy/(-y²-y+2)=dt两边积分得1/3*㏑(y+2)/(y-1)=t+C

下载地址：http://www.khdaw.com/bbs/thread-41865-1-1.html自己注册号,下载啊.采纳哦