代数方程式と微分方程式を含む方程式系を解こうとしています。これを象徴的に行うには、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 を、私が試した方程式またはその解と組み合わせても、有用な結果が得られません。何か案は?