ukman 发表于 2009-8-4 22:39

大家觉得solve解符号方程可靠吗

我尝试过解一个7元一次方程组,由于在各个方程中都有好几个变量,因此解出来的解很长很长,也不知道是对是错。手算又非常麻烦,大家凭经验觉得sovle可靠吗?

friendchj 发表于 2009-8-5 14:41

解出来的结果应该没有问题,你可以试着化简一下结果

ukman 发表于 2009-8-5 15:55

你是说用simplfy吗,解出来的结果已经不能再化简了

[ 本帖最后由 无水1324 于 2009-8-7 22:23 编辑 ]

friendchj 发表于 2009-8-5 16:38

回复 板凳 ukman 的帖子

能否给出你的例子,vpa也可用于化简

ukman 发表于 2009-8-5 17:30

好的,我整理下,谢谢先

zhoumath 发表于 2009-8-5 17:37

估计用maple更好一些

ukman 发表于 2009-8-5 19:46

方程定义如下,
syms A1 A2 A3 A4 A5 A6 A7 A8 A9 A10 A11 A12 A13 A14 A15 A16 A17 A18 B1 B2 B3 B4 B5 B6 B7 B8 B9 B10 B11 B12 B13 B14 B15 B16 B17 C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 D1 D2 D3 D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 D16 D17 D18 D19 E1 E2 E3 E4 E5 E6 E7 E8 E9 E10 E11 E12 E13 E14 E15 E16 E17 E18 E19 F1 F2 F3 F4 F5 F6 F7 F8 F9 F10 F11 F12 F13 F14 F15 F16 F17 F18 G1 G2 G3 G4 G5 G6 a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11 a12 a13 a14 X1 X2 X3 X4 X5 X6 X7 x0 x00 y0 y00;

=solve('A1*X1+A2*a2+A3*a1+A4*X3+A5*a6+A6*a5+A7*X4+A8*a8+A9*a7+A10*X5+A11*a10+A12*a9-A13*X6-A14*a10-A15*a11+A16*X7-A17*x00-A18*X0=0','A1*X2+B1*a4+B2*a3-B3*X3-B4*a6-B5*a5-B6*X4-B7*a8-B8*a7-B9*X5-B10*a10-B11*a9+B12*X6+B13*a12+B14*a11+B15*X7-B16*y00-B17*y0=0','C1*X1+C2*a2+C3*a1-C4*X2-C5*a4-C6*a3+C7*X3+C8*a6+C9*a5+C10*X4+C11*a8+C12*a7+C13*X5+C14*a10+C15*a9+C16*a12+C17*a11-C18*x00-C19*x0+C20*y00+C21*y0=0','D1*X1+D2*a2+D3*a1-D4*X2-D5*a4-D6*a3+D7*X3+D8*a6+D9*a5+D10*X4+D11*a8+D12*a7+D13*X5+D14*a10+D15*a9-D16*x00-D17*x0+D18*y00+D19*y0=0','E1*X1+E2*a2+E3*a1-E4*X2-E5*a4-E6*a3+E7*X3+E8*a6+E9*a5+E10*X4+E11*a8+E12*a7+E13*X5+E14*a10+E15*a9-E16*x00-E17*x0+E18*y00+E19*y0=0','F1*X1-F2*a2-F3*a1-F4*X2+F5*a4+F6*a3+F7*X6+F8*a12+F9*a11+F10*a6+F11*a5+F12*X7+F13*a14+F14*a13+F15*x00+F16*x0-F17*y00-F18*y0=0','G1*X1+G2*X2+G3*X6+G4*a12+G5*a11+G6*X7+G4*a14+G5*a13=0',X1, X2, X3, X4, X5, X6, X7)

我试了一下,结果 Output truncated.Text exceeds maximum line length of 25,000 characters for Command Window display.无奈啊!

[ 本帖最后由 ChaChing 于 2009-8-8 00:19 编辑 ]

friendchj 发表于 2009-8-6 14:23

其实就是求解:
A1*X1+……+A7*X7=C
可以把方程组化简一下,所有的常数记为一个符号,试试:
syms A1 A2 A3 A4 A5 A6 A7 A8 A9 A10 A11 A12 A13 A14 A15 A16 A17 A18 B1 B2...
   B3 B4 B5 B6 B7 B8 B9 B10 B11 B12 B13 B14 B15 B16 B17 C1 C2 C3 C4 C5...
   C6 C7 C8 C9 C10 C11 C12 C13 C14 C15 C16 C17 C18 C19 C20 C21 D1 D2 D3...
   D4 D5 D6 D7 D8 D9 D10 D11 D12 D13 D14 D15 D16 D17 D18 D19 E1 E2 E3...
   E4 E5 E6 E7 E8 E9 E10 E11 E12 E13 E14 E15 E16 E17 E18 E19 F1 F2 F3...
   F4 F5 F6 F7 F8 F9 F10 F11 F12 F13 F14 F15 F16 F17 F18 G1 G2 G3 G4 G5...
   G6X1 X2 X3 X4 X5 X6 X7 Z1 Z2 Z3 Z4 Z5 Z6 Z7

=solve('X1+A4*X3+A7*X4+A10*X5+A13*X6+A16*X7=Z1',...
                                          'X2+B3*X3+B6*X4+B9*X5+B12*X6+B15*X7=Z2',...
                                          'X1+C4*X2+C7*X3+C10*X4+C13*X5=Z3',...
                                          'X1+D4*X2+D7*X3+D10*X4+D13*X5=Z4',...
                                          'X1+E4*X2+E7*X3+E10*X4+E13*X5=Z5',...
                                          'X1+F4*X2+F7*X6+F12*X7=Z6',...
                                          'X1+G2*X2+G3*X6+G6*X7=Z7',...
                                          X1, X2, X3, X4, X5, X6, X7);

ukman 发表于 2009-8-7 16:52

恩,我后来也是这里处理的。方便多了。
页: [1]
查看完整版本: 大家觉得solve解符号方程可靠吗