我想寻求y
特别是ODE y'' - y' - 2y = 4x^2
我制作了以下剧本:
syms x A0 A1 A2
ypa = A2*x^2+A1*x+A0; % y_p assume
cyp = diff(ypa,2) - diff(ypa) - 2*ypa % according to ODE
P1 = 4*x^2; P2 = cyp ; % Equating P1 and P2
C = coeffs(P1 - P2,x);
A0 = solve(C(1),A0)
A1 = solve(C(2),A1)
A2 = solve(C(3),A2)
我得到A2 = -2
的正确答案。但我没有得到A0
(应该是-3
)和A1
(应该是2
)。如何自动获取它们?
P.S我正在使用MATLAB R2013a。
答案 0 :(得分:2)
不是在C
的每个等式上调用solve 3次,而是应该在整个方程组上调用一次,以便进行适当的替换以获得数值结果每个变量:
>> [A0, A1, A2] = solve(C)
A0 =
-3
A1 =
2
A2 =
-2