解方程组

rogen

我在解方程的时候出现了这样的解，请各位专家学者给解释一下,谢谢！

clear all;clc
w1=0.2713;f1=1;f2=0.4654;
[x1,x2,x3,x4,x5,x6,x7,y1,y2,y3,y4,y5,y6]=solve(...
    'y1*(f2*x3*x6+f2*(x4+x5))+y2*(f1*x3*x6-w1*f1*x6-f2*x5)+y3*(-f2*x5)+y4*(w1*f2*x7+f2*x5)+y5*x6*f1+y6*(f1*(x4+x5)-f2*x5)=0',...
    'y1*(w1*f1*x6-(x4+x5)*f1)+y2*(x5*f1)+y3*(w1*f2*x6-f2*x3*x7+x5*f1)+y4*(f1*x3*x7-x5*f1)+y5*f2*x7+y6*(f2*x5-f1*x5)=0',...
    'x4+y1*(-f2*x1*x6)+y2*(f1*x1*x6)+y3*(-f2*x2*x7)+y4*(f1*x2*x7)=0',...
    'x3-w1*x7+y1*(-(f1*x2-f2*x1))+y6*(f1*x1)=0',...
    'y1*(f2*x1-f1*x2)+y2*(f1*x2-f2*x1)+y3*(f1*x2-f2*x1)-y4*(f1*x2-f2*x1)+y6*(f1*x1-f2*x1+f2*x2)=0',...
    'w1^2*x7+y1*(w1*f1*x2-f2*x1*x3)+y2*(f1*x1*x3-w1*f1*x1)+y3*(w1*f2*x2)+y5*f1*x1=0',...
    '-w1-w1*x4+w1^2*x6+y3*(-f2*x2*x3)+y4*(-w1*f2*x1+f1*x2*x3)+y5*f2*x2=0',...
    'w1*f1*x2*x6-f2*x1*x3*x6-(x4+x5)*(f1*x2-f2*x1)=0',...
    'f1*x1*x3*x6-w1*f1*x6*x1+x5*(f1*x2-f2*x1)=0',...
    'w1*f2*x6*x2-f2*x2*x3*x7+x5*(f1*x2-f2*x1)=0',...
    '-w1*f2*x1*x7+f1*x2*x3*x7-x5*(f1*x2-f2*x1)=0',...
    'f1*x1*x6+f2*x2*x7=0',...
    '(f1*(x4+x5)-f2*x5)*x1+(f2*x5-f1*x5)*x2=0')

这是我的程序

[ 本帖最后由 ChaChing 于 2010-3-15 23:56 编辑 ]

rogen

得到的结果如下：

x1 =

         0
         0
         0
-f2*x2/f1
        x1
        x1
        x1
        x1
  f1*x2/f2
        x1



x2 =

         0
         0
        x2
        x2
         0
        x2
-f1*x1/f2
-f1*x1/f2
        x2
-f1*x1/f2



x3 =

                                 0
                                 0
                          y1*f1*x2
          f2*w1*(f2+f1)/(f2-f1)/f1
     -1-y2*f2*x1-y3*f2*x1-y6*f2*x1
        y1*f1*x2-y1*f2*x1-y6*f1*x1
                                w1
-x1*(y1*f1^2+f2^2*y1+f2*y6*f1)/f2
                                 0
          f2*w1*(f2+f1)/(f2-f1)/f1



x4 =

                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
1/2*(f2^6*w1+f2^5*f1*w1-f1^2*f2^4+2*f2^4*f1^2*w1+2*f1^3*f2^3+f2^3*f1^3*w1-f1^4*f2^2-f1^6*w1)/f2^2/f1^2*(f2+f1)/(-f2*f1+f1^2+f2^2)/(f2-f1)
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0



x5 =

                                                                                                                                        x5
                                                                                                                                         0
                                                                                                                                         0
1/2*(f2^6*w1+f2^5*f1*w1-f1^2*f2^4+2*f2^4*f1^2*w1+2*f1^3*f2^3+f2^3*f1^3*w1-f1^4*f2^2-f1^6*w1)/(f2^2-2*f2*f1+f1^2)/f1/f2/(-f2*f1+f1^2+f2^2)
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0
                                                                                                                                         0



x6 =

                                                                                                                               1/w1
                                                                                                                               1/w1
                                                                                                                                  0
1/2*(f2^6*w1+f2^5*f1*w1-f1^2*f2^4+2*f2^4*f1^2*w1+2*f1^3*f2^3+f2^3*f1^3*w1-f1^4*f2^2-f1^6*w1)/w1/(-f2*f1+f1^2+f2^2)/f2^2/f1/(f2-f1)
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0



x7 =

                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
1/2*(f2^6*w1+f2^5*f1*w1-f1^2*f2^4+2*f2^4*f1^2*w1+2*f1^3*f2^3+f2^3*f1^3*w1-f1^4*f2^2-f1^6*w1)/w1/(-f2*f1+f1^2+f2^2)/f2^2/f1/(f2-f1)
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0
                                                                                                                                  0



y1 =

                                                                                                                   y1
                                                                                                            -f2*y3/f1
                                                                                                                   y1
-1/2*(3*f2^5*w1+f2^3*f1^2-f1^3*f2^2+f2^2*f1^3*w1-f1^4*w1*f2+f1^5*w1)/x2/f2^2/(f2^4-f2^3*f1+2*f2^2*f1^2-f1^3*f2+f1^4)
                                                                        (y2*f2*x1+y3*f2*x1+1-y6*f1*x1+y6*f2*x1)/f2/x1
                                                                                                                   y1
                                                                           -(y4*f1^2*x1+f2+y4*f2^2*x1)/x1/(f1^2+f2^2)
                                                                                                                   y1
                                                                                                                   y1
                               -f2*(f2^4*w1+f2^3*f1*w1-f2^2*f1^2+2*f1^3*f2-f1^4)/f1^3/x1/(-f2^2*f1+f2*f1^2-f1^3+f2^3)

[ 本帖最后由 ChaChing 于 2010-3-15 23:55 编辑 ]

rogen

结果还有一些，由于版面的问题，我就不在列举了，希望大家帮我看看！

列出的是不是方程的多个解啊？

问题解决了，希望楼主把这帖子删了？

[ 本帖最后由 ChaChing 于 2010-3-15 23:57 编辑 ]

ChaChing

建议与大家分享你的成果, 做个结束!