大家好我想问如何在解集中保留实数解剔除复数解

lhrman

[img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\289623144\QQ\WinTemp\RichOle\N`H1F`X@4FAYE0A41M_[O~3.jpg[/img][img=0,1]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\289623144\QQ\WinTemp\RichOle\N`H1F`X@4FAYE0A41M_[O~3.jpg[/img]quen := {s0*s2 = 3/8, s0+s1+s2 = 1, 2*s0*s1+2*s1*s2 = -9/2}；计算三元方程组的解
. G& Q2 {: W5 F6 }! ]
0 j( H& G# x) t7 z9 d) Q9 N. P
! F( y2 H. {$ N: d* K, t6 `$ E9 m. Gquen2 := fsolve(quen);  用fsolve只能得到一个解

- i1 H- ~6 U% i  oquen3 := allvalues(solve(quen));可以将所有的解求出，但是其中涵盖太多的复数解，我们只需要实数解（实数解不只一个）

2 H- T+ N) N9 L2 W" b4 ]map(subs, [quen3], quen)；验证方程的解，quen3中所有的解都是方程的解。

& o) i! e9 O9 y( K0 P: J/ {8 N我想求出所有解之后把复数解都删掉
4 w. D0 J2 j% }: A* y0 n
' p; a9 F+ R) S4 _

haoyongle

quen3(imag(quen3)~=0)=[];