我正在尝试求解包含代数和微分方程的方程系统。要象征性地做到这一点,我需要结合dsolve和解决(我呢?)。将求解和dsolve结合起来求解具有微分和代数方程的方程系统
考虑下面的例子: 我们有三个基方程
a == b + c; % algebraic equation
diff(b,1) == 1/C1*y(t); % differential equation 1
diff(c,1) == 1/C2*y(t); % differential equation 2
求解两个微分方程,消除INT(Y,0..t),然后求解C = F(C1,C2,一个)收益率
C1*b == C2*c or C1*(a-c) == C2*c
c = C1/(C1+C2) * a
我该如何说服Matlab给我那个结果?这是我的尝试:
syms a b c y C1 C2;
Eq1 = a == b + c; % algebraic equation
dEq1 = 'Db == 1/C1*y(t)'; % differential equation 1
dEq2 = 'Dc == 1/C2*y(t)'; % differential equation 2
[sol_dEq1, sol_dEq2]=dsolve(dEq1,dEq2,'b(0)==0','c(0)==0'); % this works, but no inclusion of algebraic equation
%[sol_dEq1, sol_dEq2]=dsolve(dEq1,dEq2,Eq1,'c'); % does not work
%solve(Eq1,dEq1,dEq2,'c') % does not work
%solve(Eq1,sol_dEq_C1,sol_dEq_C2,'c') % does not work
解决方案和/或dsolve与方程或他们的解决方案没有组合我尝试给了我一个有用的结果。有任何想法吗?