【发布时间】:2015-02-27 09:08:21
【问题描述】:
我正在尝试求解 Maxima(即 wxMaxima)中的一阶 ODE。从ode2,我获得了一个包含不定积分的表达式,我想在我定义的限制下对其进行评估,即(0,t)。
pta: 'diff(q,t)+(1/taub)*q = 1/(c0+c1)*'diff(Q,t)+1/(c0*taub)*Q;
Qproblem: (Q1/t1)*t;
ptap: ev(pta,Q=Qproblem);
assume(t1>0, Q1>0, t>0);
sol: ode2(ptap,q,t);
grind(sol);
如果我在解决方案表达式上应用grind,我会得到以下结果
q = %e^-(t/taub)*((c0*t1*taub*'integrate(%e^(t/taub)*'diff(t*Q1/t1,t,1),t)
+c1*(t*taub-taub^2)*%e^(t/taub)*Q1
+c0*(t*taub-taub^2)*%e^(t/taub)*Q1)
/((c0*c1+c0^2)*t1*taub)
+%c)$
除非这些积分变成确定的,否则无法正确评估此表达式。因此,我想以某种方式转换它们。我在上面提到了所需的限制。
我应该使用哪种方法? ode2 是否有标志我不知道这对此有帮助吗?或者是否有一个函数可以从表达式中提取和替换部分,这样我就可以在'integrates 中添加两个参数,同时仍然保留第一个参数?
【问题讨论】:
标签: differential-equations maxima