MATLAB 常微分方程（ordinary differential equation）及常微分方程组的解析解

最新推荐文章于 2024-05-24 16:07:48 发布

我爱吃猪猪包

最新推荐文章于 2024-05-24 16:07:48 发布

阅读量853

点赞数 2

文章标签： matlab 开发语言线性代数

本文链接：https://blog.csdn.net/yuandebaba/article/details/126271100

版权

desolve函数求解常微分方程（组）

一.一阶微分方程

1.不带入初值

示例：y’ = 2x

// An MATLAB block
syms y(x)  %声明y是x的一个函数
dsolve(diff(y) == 2*x) %diff(y)表示对y求导

结果：在这里插入图片描述
可看出结果为一个不定积分

2.带入初值

示例：y’ = 2x，初值y(0)=1

// An MATLAB block
syms y(x)  %声明y是x的一个函数
dsolve(diff(y) == 2*x,y(0)==1) %diff(y)表示对y求导

结果：
在这里插入图片描述

二.二阶微分方程

示例：y’’ = uy’ + y, 初值 y(0）=1，y’(0)=0;

1.不带人初值

// An MATLAB block
syms y(x)  u
eqn = diff(y,2)==u*diff(y)+y;  %将等式赋值给变量
dsolve(eqn)

结果：
在这里插入图片描述
可看出结果为一个不定积分

2.带人初值

// An MATLAB block
syms y(x)  u
eqn = diff(y,2)==u*diff(y)+y;  %将等式赋值给变量
cond1 = y(0) == 1;
dy=diff(y);
cond2=dy(0)==0;
dsolve(eqn,cond1,cond2) %diff(y)表示对y求导

结果：
在这里插入图片描述
当方程较复杂时不能求出解析解，只能求出数值解

三.常微分方程组

示例：方程组（1）y’(x)=z(x) (2) z’(x)=-y(x) ;初值 y(0)=1 , z(0)=1

1.不带入初值

// An MATLAB block
syms y(x)  z(x)
eqn1 = diff(y)==z; 
eqn2 = diff(z)==-y; 
res=dsolve(eqn1,eqn2) ;
res.y
res.z

结果：
在这里插入图片描述
可看出结果为一个不定积分

2.带入初值

// An MATLAB block
syms y(x)  z(x)
eqn1 = diff(y)==z; 
eqn2 = diff(z)==-y; 
cond1 = y(0) == 1;
cond2 = z(0) == 1;
res=dsolve(eqn1,eqn2,cond1,cond2) ;
res.y
res.z

结果：
在这里插入图片描述
笔记学习视频：
https://www.bilibili.com/video/BV1fB4y1h7bC?p=38&vd_source=9a25f2f4c96b69c013330f038adef247