happyup 发表于 2009-5-5 11:31

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

朋友、版主好;
请邦我运行一下,下面的这个方程组,只需把它拷进matlab中即可。我运行后出现了提示,我不明白是怎么回事。
       谢谢!!
=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 发表于 2009-5-5 11:42

回复 楼主 happyup 的帖子

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

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

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

happyup 发表于 2009-5-5 11:55

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

以上的是出错提示



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

friendchj 发表于 2009-5-5 13:21

回复 板凳 happyup 的帖子

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

happyup 发表于 2009-5-5 13:30

原帖由 friendchj 于 2009-5-5 13:21 发表 http://www.chinavib.com/forum/images/common/back.gif
Explicit solution 是 显式解 的意思;
说明方程没有显式解。

谢谢!
是无解吗?那怎么办呢?邦我一下,谢谢!

ChaChing 发表于 2009-5-5 13:44

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

happyup 发表于 2009-5-5 13:49

原帖由 ChaChing 于 2009-5-5 13:44 发表 http://www.chinavib.com/forum/images/common/back.gif
符号方面不熟! 报错好像告知无法求解了!
瞄了下LZ的方程, 真有够伤视力!
感觉好像若A为定数, 其实便是简单的AX=b问题!?
建议LZ自行想办法简化, matlab仅是一工具而已, 又非万能!
待高人路过!

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

friendchj 发表于 2009-5-5 13:50

应该可以求数值解

ChaChing 发表于 2009-5-5 13:56

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

happyup 发表于 2009-5-5 13:59

原帖由 ChaChing 于 2009-5-5 13:56 发表 http://www.chinavib.com/forum/images/common/back.gif
有点好奇! 原始工程问题是什麽?

是一些节点受力和位移。

happyup 发表于 2009-5-5 14:08

原帖由 friendchj 于 2009-5-5 13:50 发表 http://www.chinavib.com/forum/images/common/back.gif
应该可以求数值解

怎么做呢,请明示。
            谢谢!!

ChaChing 发表于 2009-5-5 14:12

回复 10楼 happyup 的帖子

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

happyup 发表于 2009-5-5 14:17

原帖由 ChaChing 于 2009-5-5 14:12 发表 http://www.chinavib.com/forum/images/common/back.gif
不解, 节点受力和位移怎会有那些非线性项?

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

friendchj 发表于 2009-5-5 14:19

回复 11楼 happyup 的帖子

主程序:
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 发表于 2009-5-5 14:30

谢谢friendchj朋友,这里的解是近似解吧?
其实这个方程组中,A的值是最量要的,能不能让A准确些呢?
            谢谢!!
页: [1] 2
查看完整版本: 求助,请邦我用matlab运行一下这个方程组