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

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}；计算三元方程组的解
" L, x; r3 S* C- V. K! ^
0 o: B, j  ]# d! B! P6 g- C: I, a
! |. A5 ^* W( @; [7 W: j8 aquen2 := fsolve(quen);  用fsolve只能得到一个解

9 z7 n3 I8 x6 V7 M- x6 r2 I* Bquen3 := allvalues(solve(quen));可以将所有的解求出，但是其中涵盖太多的复数解，我们只需要实数解（实数解不只一个）
" D9 K2 a2 x) T+ \4 O" Y
* ~- ~- ]7 }, m, Qmap(subs, [quen3], quen)；验证方程的解，quen3中所有的解都是方程的解。
  r) v5 E/ s9 q# |
我想求出所有解之后把复数解都删掉

0 n2 T. B: Q4 c, z. A

haoyongle

quen3(imag(quen3)~=0)=[];
' d* j4 m6 P$ x( x; o( c& K