求助，请邦我用matlab运行一下这个方程组

happyup

朋友、版主好；
  请邦我运行一下，下面的这个方程组，只需把它拷进matlab中即可。我运行后出现了提示，我不明白是怎么回事。
       谢谢！！
[X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13,X14,X15,X16,X17,X18,X19,X20,X21,X22,A]=solve('A*X1+A*X2+A*X3+2*X4+A*X5+X6+X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+X14+A*X15+X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=0','X1+A*X2+A*X3+X4+A*X5+2*X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+X14+A*X15+A*X16+A*X17+X18+A*X19+A*X20+ X21+A*X22=2','A*X1+A*X2+A*X3+X4+X5+X6+A*X7+X8+A*X9+A*X10+A*X11+X12+A*X13+X14+A*X15+A*X16+A*X17+A*X18+A*X19+X20+X21+A*X22=6','2*X1+A*X2+A*X3+X4+A*X5+X6+A*X7+A*X8+A*X9+A*X10+X11+X12+A*X13+X14+A*X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=3','A*X1+A*X2+X3+X4+X5+X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+X13+X14+A*X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=5','A*X1+A*X2+A*X3+X4+X5+X6+A*X7+A*X8+X9+A*X10+A*X11+X12+A*X13+X14+X15+A*X16+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=7','A*X1+A*X2+A*X3+X4+X5+X6+A*X7+X8+A*X9+A*X10+A*X11+X12+A*X13+X14+A*X15+A*X16+ A*X17+A*X18+A*X19+A*X20+X21+A*X22=4','A*X1+A*X2+A*X3+X4+A*X5+X6+A*X7+A*X8+X9+X10+A*X11+X12+A*X13+X14+A*X15+A*X16+ A*X17+A*X18+A*X19+X20+X21+A*X22=3','A*X1+X2+X3+X4+A*X5+X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+X13+X14+A*X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+X22=5','A*X1+A*X2+A*X3+X4+X5+X6+A*X7+A*X8+X9+A*X10+A*X11+X12+X13+X14+A*X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=6','A*X1+A*X2+X3+A*X4+X5+X6+X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+2+X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=9','A*X1+A*X2+A*X3+A*X4+X5+X6+X7+A*X8+X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+X17+A*X18+A*X19+A*X20+X21+A*X22=4','A*X1+X2+X3+A*X4+X5+X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+A*X19+X20+X21+A*X22=1','A*X1+A*X2+A*X3+X4+X5+X6+X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+X22=1','X1+A*X2+A*X3+A*X4+X5+2*X6+A*X7+A*X8+A*X9+A*X10+A*X11+2*X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=3','A*X1+A*X2+A*X3+A*X4+X5+X6+X7+A*X8+X9+A*X10+A*X11+X12+A*X13+A*X14+2*X15+A*X16+A*X17+A*X18+A*X19+A*X20+X20+A*X22=3','X1+A*X2+A*X3+A*X4+X5+X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+X17+A*X18+A*X19+A*X20+X21+A*X22=8','A*X1+A*X2+X3+A*X4+2*X5+X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+X19+A*X20+X21+A*X22=8','X1+A*X2+A*X3+A*X4+X5+2*X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+X22=3','A*X1+A*X2+A*X3+A*X4+2*X5+X6+A*X7+X8+A*X9+A*X10+A*X11+2*X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=3','A*X1+A*X2+A*X3+A*X4+X5+X6+A*X7+A*X8+X9+X10+A*X11+X12+A*X13+X14+X15+A*X16+A*X17+A*X18+A*X19+A*X20+X21+A*X22=8','A*X1+A*X2+X3+A*X4+2*X5+X6+A*X7+A*X8+A*X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+X17+A*X18+A*X19+A*X20+X21+A*X22=4','A*X1+A*X2+A*X3+A*X4+2*X5+X6+A*X7+A*X8+X9+A*X10+A*X11+X12+A*X13+A*X14+X15+A*X16+A*X17+A*X18+X19+A*X20+X21+A*X22=4')

运行后，没有解出来，有提示，我不明白。

[ 本帖最后由 ChaChing 于 2009-5-5 11:44 编辑 ]

ChaChing

给齐完整格式：出错代码和出错提示!

这麽复杂又长, 我看又没人有耐心看了!:@)

[ 本帖最后由 ChaChing 于 2009-5-5 11:48 编辑 ]

happyup

Warning: Explicit solution could not be found.
> In solve at 140

以上的是出错提示



我写的方程组是完整格式的，拷入程序后就可以运算。

friendchj

Explicit solution 是 显式解 的意思；
说明方程没有显式解。

happyup

原帖由 friendchj 于 2009-5-5 13:21 发表

                               
登录/注册后可看大图

Explicit solution 是 显式解 的意思；
说明方程没有显式解。

谢谢！
  是无解吗？那怎么办呢？邦我一下，谢谢！

ChaChing

符号方面不熟! 报错好像告知无法求解了!
瞄了下LZ的方程, 真有够伤视力!
感觉好像若A为定数, 其实便是简单的AX=b问题!?
建议LZ自行想办法简化, matlab仅是一工具而已, 又非万能!
待高人路过!

happyup

原帖由 ChaChing 于 2009-5-5 13:44 发表

                               
登录/注册后可看大图

符号方面不熟! 报错好像告知无法求解了!
瞄了下LZ的方程, 真有够伤视力!
感觉好像若A为定数, 其实便是简单的AX=b问题!?
建议LZ自行想办法简化, matlab仅是一工具而已, 又非万能!
待高人路过!

说的及是，但是A也是变量，如是定数则太好解了。

friendchj

应该可以求数值解

ChaChing

有点好奇! 原始工程问题是什麽?

happyup

原帖由 ChaChing 于 2009-5-5 13:56 发表

                               
登录/注册后可看大图

有点好奇! 原始工程问题是什麽?

是一些节点受力和位移。

happyup

原帖由 friendchj 于 2009-5-5 13:50 发表

                               
登录/注册后可看大图

应该可以求数值解

怎么做呢，请明示。
            谢谢！！

ChaChing

不解, 节点受力和位移怎会有那些非线性项?

happyup

原帖由 ChaChing 于 2009-5-5 14:12 发表

                               
登录/注册后可看大图

不解, 节点受力和位移怎会有那些非线性项?

不是通常的单元体。显然你是这方面的高手。

friendchj

主程序：
clear
clc
x0=zeros(23,1);
F=fsolve(@myfun,x0);
disp(F.')
子程序：
function F=myfun(x)
F=[x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+2*x(4)+x(23)*x(5)+x(6)+x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(23)*x(15)+x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22);
   x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(23)*x(5)+2*x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+x(23)*x(17)+x(18)+x(23)*x(19)+x(23)*x(20)+ x(21)+x(23)*x(22)-2;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(5)+x(6)+x(23)*x(7)+x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(20)+x(21)+x(23)*x(22)-6;
   2*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(23)*x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(11)+x(12)+x(23)*x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-3;
   x(23)*x(1)+x(23)*x(2)+x(3)+x(4)+x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-5;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(15)+x(23)*x(16)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-7;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(5)+x(6)+x(23)*x(7)+x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+ x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-4;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(23)*x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(9)+x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+ x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(20)+x(21)+x(23)*x(22)-3;
   x(23)*x(1)+x(2)+x(3)+x(4)+x(23)*x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(22)-5;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(13)+x(14)+x(23)*x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-6;
   x(23)*x(1)+x(23)*x(2)+x(3)+x(23)*x(4)+x(5)+x(6)+x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+2+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-9;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+x(5)+x(6)+x(7)+x(23)*x(8)+x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-4;
   x(23)*x(1)+x(2)+x(3)+x(23)*x(4)+x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(20)+x(21)+x(23)*x(22)-1;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(4)+x(5)+x(6)+x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(22)-1;
   x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+x(5)+2*x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+2*x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-3;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+x(5)+x(6)+x(7)+x(23)*x(8)+x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+2*x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(20)+x(23)*x(22)-3;
   x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-8;
   x(23)*x(1)+x(23)*x(2)+x(3)+x(23)*x(4)+2*x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-8;
   x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+x(5)+2*x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(22)-3;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+2*x(5)+x(6)+x(23)*x(7)+x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+2*x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-3;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(9)+x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-8;
   x(23)*x(1)+x(23)*x(2)+x(3)+x(23)*x(4)+2*x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(23)*x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(17)+x(23)*x(18)+x(23)*x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-4;
   x(23)*x(1)+x(23)*x(2)+x(23)*x(3)+x(23)*x(4)+2*x(5)+x(6)+x(23)*x(7)+x(23)*x(8)+x(9)+x(23)*x(10)+x(23)*x(11)+x(12)+x(23)*x(13)+x(23)*x(14)+x(15)+x(23)*x(16)+x(23)*x(17)+x(23)*x(18)+x(19)+x(23)*x(20)+x(21)+x(23)*x(22)-4];
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
其中X1相当于x(1)，……，X22相当于x(22)，A相当于x(23)
运行结果(初值随便取的，可能不是最优解)：
Maximum number of function evaluations reached:
increase options.MaxFunEvals.
  Columns 1 through 12
    1.8513   -1.2056    1.2049   -1.1751   -0.0800   -0.4339    1.2981    1.0365   -0.0541   -0.0713   -2.8257   -0.2233
  Columns 13 through 23
    0.9443    3.6688    1.5892   -2.2378    0.7666   -2.2905    1.9104    0.1621    2.5969   -0.3174   -0.8440

[ 本帖最后由 friendchj 于 2009-5-5 14:22 编辑 ]

happyup

谢谢friendchj朋友，这里的解是近似解吧？
  其实这个方程组中，A的值是最量要的，能不能让A准确些呢？
              谢谢！！