[討論] 解聯立方程式 小問題

[hsuanboy (壓力)]

以下是解二元二次方程式的小問題

第一個程式是簡單的聯立方程式 在解a和b的值 可以容易被解出來
syms  a b
S=solve('3*a^2+2*a+2*b^2-2*b=1','3*a+2*b=3');
S=[S.a S.b]


S =

[ 1/3,   1]
[ 1/5, 6/5]

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

第二個程式是比較複雜的聯立方程式 在解a和b的值 沒辦法容易的被解出來
a和b的值會呈現下方的樣子 a和b的值會連帶Q和W 為什麼它不會直接算出跟
第一個程式一樣的值呢 我不太曉得原因 想請教高手 有沒有什麼方式
可以利用第二個程式來解出a和b的值 或是有其他的指令可以使用的呢?
syms  a b
Q=3;
W=2;
S=solve('Q*a^2 + 2*a + W*b^2 - 2*b = 1','Q*a+ W*b = 3');
S=[S.a S.b]


S =

[ -((W*(Q + 4*W + (Q^2*W + Q^2 + Q*W^2 - 7*Q*W + W^2)^(1/2)))/(W^2 + Q*W) -
3)/Q, (Q + 4*W + (Q^2*W + Q^2 + Q*W^2 - 7*Q*W + W^2)^(1/2))/(W^2 + Q*W)]
[ -((W*(Q + 4*W - (Q^2*W + Q^2 + Q*W^2 - 7*Q*W + W^2)^(1/2)))/(W^2 + Q*W) -
3)/Q, (Q + 4*W - (Q^2*W + Q^2 + Q*W^2 - 7*Q*W + W^2)^(1/2))/(W^2 + Q*W)]

希望有好心的大大能幫我解答 謝謝阿!~

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 120.126.11.120
推 kurt28:變數宣告的問題
推 kurt28:前面syms Q W a b ;中間W Q不給 最後面subs(S,[Q W],[3 2])
推 kurt28:或者是像你這樣  最後面直接加 subs(S,{'Q','W'},[Q W])