gigicool333 发表于 2007-5-13 10:04

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

例如下面：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 如果可以解变系数的话，怎么写程序？
如果不可以的，又怎么用别的办法求解变系数的情况？

咕噜噜 发表于 2007-5-13 10:31

回复 #1 gigicool333 的帖子

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

gigicool333 发表于 2007-5-13 10:45

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

咕噜噜 发表于 2007-5-13 11:19

回复 #3 gigicool333 的帖子

可以，假如u(t)=sint，为下列程序
fun=inline('','t','y')
ode45(fun,,)或者=ode45(fun,,)
积分区间和初始条件是我自己随便给你，你自己调整，把函数u改变一下也可以，我没用m函数，内联函数一样的，你可以自己编写m函数
页: [1]
查看完整版本: 怎样求解变系数的微分方程