作者hsuanboy (压力)
看板MATLAB
标题[讨论] 解联立方程式 小问题
时间Tue Jun 25 10:37:20 2013
以下是解二元二次方程式的小问题
第一个程式是简单的联立方程式 在解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
1F：推 kurt28:变数宣告的问题 06/25 23:56
2F：推 kurt28:前面syms Q W a b ;中间W Q不给 最後面subs(S,[Q W],[3 2]) 06/26 00:01
3F：推 kurt28:或者是像你这样 最後面直接加 subs(S,{'Q','W'},[Q W]) 06/26 00:06