数学建模社区-数学中国
标题:
大家好我想问如何在解集中保留实数解剔除复数解
[打印本页]
作者:
lhrman
时间:
2015-5-17 23:36
标题:
大家好我想问如何在解集中保留实数解剔除复数解
[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}；计算三元方程组的解
* V/ `7 C5 c# X% `( @% h
3 ~4 a( V0 l! g0 k' j; c
/ P, b) J5 I, `& a3 J- M
quen2 := fsolve(quen); 用fsolve只能得到
一个解
# F: W& T! A* L. H9 ]* l8 _
* i5 [! f/ A# [* l
quen3 := allvalues(solve(quen));可以将所有的解求出，但是其中涵盖太
多的复数解
，我们只需要实数解（
实数解不只一个
）
+ p5 ]3 z$ o# r4 c( ^5 ]! t4 e& I
R, `, _- |5 I1 H1 z! |, K
map(subs, [quen3], quen)；验证方程的解，quen3中所有的解都是方程的解。
* ~/ P! S: p- D6 |( _7 Z! i/ ~
^8 j/ v7 Z% ~" S/ p
我想求出所有解之后把复数解
) ]4 D. L$ B$ p' Q0 g6 u
j0 J( e: H7 K4 O
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5