MATLAB学习之微分方程边值求解（十二）

微分方程边值求解

微分方程在给定的边值条件下求解。

例 求微分方程$x{y}^{\prime \prime }-5y^{\prime }+x^3=0$中$y(1)=0,y(5)=0$的解。

（1）首先利用dsolve函数求解方程的解

y = dsolve('x*D2y-5*Dy=-x^3','y(1)=0,y(5)=0','x')

y =

-(13*x^6)/2604 + x^4/8 - 625/5208

（2）利用MATLAB绘图绘制曲线

xn = -1:6;%绘制x轴的坐标变化范围从-1到6，步长为1
yn = subs(y,'x',xn)%将符号表达式y中的'x'替换为xn

yn =

[ 0, -625/5208, 0, 387/248, 592/93, 2835/248, 0, -1705/24]

ezplot(y,[-1,6])
hold on

plot([1,5],[0,0],'.r','MarkerSize',20)
text(1,1,'y(1) = 0')
text(4,1,'y(5) = 0')
title(['x*D2y - 5*Dy = -x^3',', y(1) = 0,y(5) = 0'])
hold off
yn

yn =

[ 0, -625/5208, 0, 387/248, 592/93, 2835/248, 0, -1705/24]

绘制曲线结果如下图所示