matlab中lambertw,MATLAB解常微分方程

在MATLAB中，由函数dsolve()解决常微分方程(组)的求解问题，其具体格式如下：

X=dsolve(‘eqn1’,’eqn2’,…)

函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。例1：求解常微分方程

的MATLAB程序为：dsolve('Dy=1/(x+y)','x')，注意，系统缺省的自变量为t，因此这里要把自变量写明。

结果为：-lambertw(-C1*exp(-x-1))-x-1 其中：Y=lambertw(X)表示函数关系Y*exp(Y)=X。

例2：求解常微分方程的MATLAB程序为：Y2=dsolve('y*D2y-Dy^2=0’,’x’) 结果为： Y2 =[

exp((x+C2)/C1)] [ C2] 我们看到有两个解，其中一个是常数。例3：求常微分方程组通解的MATLAB程序为：

[X,Y]=dsolve('Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t)','t') 例4：求常微分方程组

通解的MATLAB程序为：

[X,Y]=dsolve('Dx+2*x-Dy=10*cos(t),Dx+Dy+2*y=4*exp(-2*t)','x(0)=2','y(0)=0')

以上这些都是常微分方程的精确解法，也称为常微分方程的符号解。但是，我们知道，有大量的常微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面，MATLAB具有丰富的函数，我们将其统称为solver，其一般格式为：

[T,Y]=solver(odefun,tspan,y0)

该函数表示在区间tspan=[t0,tf]上，用初始条件y0求解显式常微分方程。

solver为命令ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb之一，这些命令各有特点。

例5：求解常微分方程，，

的MATLAB程序如下：fun=inline('-2*y+2*x*x+2*x')；[x,y]=ode23(fun,[0,0.5],1)

结果为： x =

0,0.0400,0.0900,0.1400,0.1900,0.2400,0.2900,0.3400,0.3900,0.4400,0.4900,0.5000

y =

1.0000,0.9247,0.8434,0.7754,0.7199,0.6764,0.6440,0.6222,0.6105,0.6084,0.6154,0.6179

例6：求解常微分方程的解，并画出解的图形。

分析：这是一个二阶非线性方程，用现成的方法均不能求解，但我们可以通过下面的变换，将二阶方程化为一阶方程组，即可求解。令：，，

，则得到：接着，编写vdp.m如下： function fy=vdp(t,x)

fy=[x(2);7*(1-x(1)^2)*x(2)-x(1)]; 再编写m文件sy12_6.m如下： y0=[1;0]

[t,x]=ode45(@vdp,[0,40],y0); y=x(:,1);dy=x(:,2);

plot(t,y,t,dy)

matlab中lambertw,MATLAB解常微分方程相关推荐

matlab中linprog函数解线性规划方程
线性规划问题,建立线性规划问题模型的基本步骤: 一.明确目标二.确定问题的决策变量三.构建目标函数四.确定约束条件建立好线性规划方程之后,紧接着而来的就是如何解方程,我们使用 matlab中的 ...
【Matlab】一、解常微分方程ODE
文章目录求解常微分方程 ODE (1)求解解析解 (2)求解数值解求解常微分方程 ODE 在matlab中,我们可以求解常微分方程的解析解,和数值解,一般使用dsolve来求解常微分方程的解析 ...
matlab中facedetector,Matlab人脸检测方法（Face Parts Detection）详解
今天同学让我帮忙制作一个人脸表情识别的样本库,其中主要是对人脸进行裁剪,这里用到了一个相对较新的Matlab人脸检测方法Face Parts Detection,网上百度了一下发现关于Matlab人脸 ...
matlab中desolve,MATLAB考试题的.doc
MATLAB考试题的一.填空题 1.在MATLAB中,清除所有变量的命令为_clear___________. 2. 在MATLAB中,保存和载入数据分别使用命令___save_________. ...
matlab中dblquad,MATLAB积分.ppt
MATLAB积分 (2) 用int命令求解: 先将重积分化成累次积分: MATLAB代码为 clear; syms x y z; int(int(int('y*sin(x)+z*cos(x)',z,- ...
matlab中normfit,MATLAB中如何得到一组统计数据的分布特征
MATLAB中如何得到一组统计数据的分布特征我想通过MATLAB知道一组统计数据的分布特征,如属于哪种分布,均值方差等,请问如何实现?谢谢! function f=p_judge(A,alpha) ...
matlab中fdyn,Matlab的用法总结
1. 对序列进行洗牌 randperm() randperm()产生随机的序列 %if filepaths 是一个5*1的结构体,then cshuffle = randperm(length(fil ...
matlab中rsa,matlab - Matlab中的RSA代码 - 堆栈内存溢出
我想加密一条消息,例如" HELO1234" ,然后解密以获取原始消息.我已经在matlab中编写了无法正常工作的RSA代码. 参数计算 temp=1; range=1:10; k ...
在matlab中function,Matlab中function函数使用操作方法
亲们或许不知道Matlab中function函数如何使用,那么今天小编就讲解Matlab中function函数使用操作方法哦,希望能够帮助到大家呢. Matlab中function函数使用操作方法打 ...

matlab中lambertw,MATLAB解常微分方程

matlab中lambertw,MATLAB解常微分方程相关推荐

最新文章

热门文章