3楼: Originally posted by 月只蓝 at 2014-12-09 21:56:02
令y1=y y2=y'
原二阶常微分初值问题 化为 一阶常微分方程组初值问题
化成方程组：
y1'=y2
y2'=-(2*y2+y1)

初值条件：
y1(0)=4, y2(0)=2

解析解见附图1。

数值解MATLAB代码：

function solve_odes
...

function solve_odes
clear all;clc
[t y]=ode45(@odes,[0:0.05:20],[4 2]);
figure(33)
plot(t,y(:,1),'b-s',t,y(:,2),'r-d'),axis([-0.5 20 -2 4.5]),legend('y1','y2')


function f=odes(t,y)
f(1)=y(2);
f(2)=-(2*y(2)+y(1));
f=f';

4楼: Originally posted by 依旧是则徐 at 2014-12-10 10:09:25
谢谢，能再请问一个问题吗，如果是V``(x)+(a+bcos2x）v=0这样的方程编写思路是怎样的，其中a，b 未知，可以变化。我会追加金币的。非常感谢！！！...