的解,则令
tspa
t0t1t2tf(要求ti单调),
y0初始条件。
2例5:求解常微分方程y2y2x2x,0≤x≤05,y01的MATLAB
程序如下:fu
i
li
e2y2xx2x;xyode23fu
0051
f结果为:x00040000900014000190002400029000340003900044000490005000y100000924708434077540719906764064400622206105060840615406179d2ydy1y2y0y01y002dt例6:求解常微分方程dt的解,并画出解的图形。分析:这是一个二阶非线性方程,用现成的方法均不能求解,但我们可以通过下面的变换,将二阶方程化为一阶方程组,即可求解。
dyx1y,x2dt,7,则得到:令:
dx1x101dtx2dx271x2xxx001212dt接着,编写vdpm如下:fu
ctio
fyvdptxfyx271x12x2x1再编写m文件sy12_6m如下:y010txode45vdp040y0yx1dyx2plottytdy
f三、实验内容
dy3y8y1.利用MATLAB求常微分方程的初值问题dx,x02的解。
22.利用MATLAB求常微分方程的初值问题1xy2xy,yx01,yx03
的解。
3.利用MATLAB求常微分方程y
4
2yy0的解。xt0
32
dydxt2dt4xdtyedx3xy04.利用MATLAB求常微分方程组dt
yt00
的特解。
25.求解常微分方程y21yyy0,0≤x≤30,y01,y00的特解,
并做出解函数的曲线图。
6.完成实验报告。
fr