matlab|dsolve解决常微分初值与讲解(含实例使用)

90 阅读 0 评论 60 点赞

我是靠谱客的博主自然跳跳糖，这篇文章主要介绍matlab|dsolve解决常微分初值与讲解(含实例使用)，现在分享给大家，希望可以做个参考。

本博文源于matlab中dsolve命令，这个命令主要解决常微分方程。涉及求解

复制代码

1
2
dsolve('eqn1','eqn2',...,t)

例子：求解x’(t)+4x(t)=sin2t

复制代码

1
2
3
4
5
6
7
8
>> syms x
>> dsolve('Dx+4*x=sin(2*t)')
 
ans =
 
-1/10*cos(2*t)+1/5*sin(2*t)+exp(-4*t)*C1

例子：对于方程降阶x’(t)=t

复制代码

1
2
3
4
5
6
7
8
9
10
 
>> dsolve('Dx=t')
 
ans =
 
1/2*t^2+C1
 
 
>>

例子：求解x’’(t)+2x’(t)-3x(t)=e^t

复制代码

1
2
3
4
5
6
7
8
9
>> dsolve('D2x+2*Dx-3*x=exp(t)')
 
ans =
 
exp(-3*t)*C2+exp(t)*C1+1/4*t*exp(t)
 
 
>>

例子：求解常微分方程组：x’(t)=-3y(t) y’(t)=4x(t)

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
>> syms t;
>> [x,y] =dsolve('Dx=-3*y','Dy=4*x')
 
x =
 
C1*sin(2*3^(1/2)*t)+C2*cos(2*3^(1/2)*t)
 
 
 
y =
 
-2/3*3^(1/2)*(cos(2*3^(1/2)*t)*C1-C2*sin(2*3^(1/2)*t))
 
 
>>

例子：求解常微分方程初值问题

在这里插入图片描述

复制代码

1
2
3
4
5
6
7
8
9
10
>> syms x 
>> dsolve('D2y+3*Dy+2*y=sin(2*x)','y(0)=1','Dy(0)=1')
 
ans =
 
1/2*sin(2*x)-exp(-2*t)*(-1/2*sin(2*x)+2)+exp(-t)*(-sin(2*x)+3)
 
 
>>