数学建模社区-数学中国
标题:
大家好我想问如何在解集中保留实数解剔除复数解
[打印本页]
作者:
lhrman
时间:
2015-5-17 22:44
标题:
大家好我想问如何在解集中保留实数解剔除复数解
[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}；计算三元方程组的解
' a* K! x s+ L: b- L+ t
2 r V: `3 Z l) z6 f- {! p# Y( ^' [4 K
; a6 p& c0 f1 c# K# b, l
quen2 := fsolve(quen); 用fsolve只能得到
一个解
/ Y' f8 G7 K% B: X8 `: ^1 C
5 p3 i" H1 r) R0 g9 o& P
quen3 := allvalues(solve(quen));可以将所有的解求出，但是其中涵盖太
多的复数解
，我们只需要实数解（
实数解不只一个
）
4 ~* `; m4 p0 M. E% N' o% n5 \
% X: ~/ L1 T* s2 L! S! G
map(subs, [quen3], quen)；验证方程的解，quen3中所有的解都是方程的解。
" q: [2 C; }# S0 P: t2 p
3 r& C* n9 ~, Q1 J+ ` @8 Z* d
我想求出所有解之后把复数解都删掉
$ x C' P d" q# h
8 I: Z ]6 v& S1 v( \
作者:
haoyongle
时间:
2015-5-18 15:36
quen3(imag(quen3)~=0)=[];
+ Y0 r, [4 r/ o9 y2 l4 O* f+ j( G! ^
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5