用matlab求解常微分方程

我有以下建议：

b'(t) + k16*b(t) = k15*a(t)

其中k15和k16是常数

你知道怎么解决吗

谢谢！

---

Amit

这是一首一阶颂歌。有一个解析解（只需使用积分因子）。不需要集成。

但是，如果要在MATLAB中求解它：

>> k15 = 0.2; k16 = 0.3; % type your constants here
>> a = @(t) t^2; % type your expression for a here
>> dbdt = @(t,b) -k16*b + k15*a(t);
>> tf = 10; % final time of integration
>> b0 = 1; % initial value of b
>> [t,y] = ode45(@dbdt,[0 tf],b0)
>> plot(t,y) % display solution.

---

这是一首一阶颂歌。有一个解析解（只需使用积分因子）。不需要集成。

但是，如果要在MATLAB中求解它：

>> k15 = 0.2; k16 = 0.3; % type your constants here
>> a = @(t) t^2; % type your expression for a here
>> dbdt = @(t,b) -k16*b + k15*a(t);
>> tf = 10; % final time of integration
>> b0 = 1; % initial value of b
>> [t,y] = ode45(@dbdt,[0 tf],b0)
>> plot(t,y) % display solution.