怎样求解变系数的微分方程

gigicool333

例如下面：u(t)是一个关于t的函数，而不是常数
function q = bianxishu(t,y)  
q = zeros(2,1);               
q(1) =y(2);
q(2) = u(t) -2*y(1)-3*y(2);
我想问的是：
1，ode45是不是只能解这种常系数的微分方程？还是变系数的也可以
2 如果可以解变系数的话，怎么写程序？
  如果不可以的，又怎么用别的办法求解变系数的情况？

咕噜噜

ode45不能解变系数微分方程吧，变系数微分方程如果是单自由度的有好多书上都有介绍，但是你给的这个方程不算是变系数微分方程，也就相当于一个单自由度阻尼系统的强迫振动，不难
或者还是你方程表述错了

gigicool333

那样的话象我这个问题用ode45可以求解出来吗？

咕噜噜

可以，假如u(t)=sint，为下列程序
fun=inline('[y(2);sin(t)-2*y(1)-3*y(2)]','t','y')
ode45(fun,[0 1],[0 1])或者[t,y]=ode45(fun,[0 1],[0 1])
积分区间和初始条件是我自己随便给你，你自己调整，把函数u改变一下也可以，我没用m函数，内联函数一样的，你可以自己编写m函数