求助，请邦我用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 编辑 ]
回复 楼主 happyup 的帖子
给齐完整格式：出错代码和出错提示!这麽复杂又长, 我看又没人有耐心看了!:@)
[ 本帖最后由 ChaChing 于 2009-5-5 11:48 编辑 ] Warning: Explicit solution could not be found.
> In solve at 140
以上的是出错提示
我写的方程组是完整格式的，拷入程序后就可以运算。
回复 板凳 happyup 的帖子
Explicit solution 是 显式解 的意思；说明方程没有显式解。 原帖由 friendchj 于 2009-5-5 13:21 发表 http://www.chinavib.com/forum/images/common/back.gif
Explicit solution 是 显式解 的意思；
说明方程没有显式解。
谢谢！
是无解吗？那怎么办呢？邦我一下，谢谢！ 符号方面不熟! 报错好像告知无法求解了!
瞄了下LZ的方程, 真有够伤视力!
感觉好像若A为定数, 其实便是简单的AX=b问题!?
建议LZ自行想办法简化, matlab仅是一工具而已, 又非万能!
待高人路过! 原帖由 ChaChing 于 2009-5-5 13:44 发表 http://www.chinavib.com/forum/images/common/back.gif
符号方面不熟! 报错好像告知无法求解了!
瞄了下LZ的方程, 真有够伤视力!
感觉好像若A为定数, 其实便是简单的AX=b问题!?
建议LZ自行想办法简化, matlab仅是一工具而已, 又非万能!
待高人路过!
说的及是，但是A也是变量，如是定数则太好解了。 应该可以求数值解 有点好奇! 原始工程问题是什麽? 原帖由 ChaChing 于 2009-5-5 13:56 发表 http://www.chinavib.com/forum/images/common/back.gif
有点好奇! 原始工程问题是什麽?
是一些节点受力和位移。 原帖由 friendchj 于 2009-5-5 13:50 发表 http://www.chinavib.com/forum/images/common/back.gif
应该可以求数值解
怎么做呢，请明示。
谢谢！！
回复 10楼 happyup 的帖子
不解, 节点受力和位移怎会有那些非线性项? 原帖由 ChaChing 于 2009-5-5 14:12 发表 http://www.chinavib.com/forum/images/common/back.gif不解, 节点受力和位移怎会有那些非线性项?
不是通常的单元体。显然你是这方面的高手。
回复 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 编辑 ] 谢谢friendchj朋友，这里的解是近似解吧？
其实这个方程组中，A的值是最量要的，能不能让A准确些呢？
谢谢！！
页:
[1]
2