4

代数方程式と微分方程式を含む方程式系を解こうとしています。これを象徴的に行うには、dsolve と solve を組み合わせる必要があります (そうですか?)。

次の例を考えてみましょう: 3 つの基本方程式があります。

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,a) を解くと、

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

solve および/または dsolve を、私が試した方程式またはその解と組み合わせても、有用な結果が得られません。何か案は?

4

1 に答える 1