多元方程的求解问题（符号求解）

TDF

我想用matlab里的矩阵符号运算功能解出这组四元一次方程的解。但是解不出来。求各位大侠 指点一下。谢谢啦。程序如下：
A=sym('-e^2*a^2-1,2*c*e*a+2*f*d*e*a^2,-f*a^2,-2*f*d*e*a^2;-2*c*e+2*f*d*e*a^2,e^2*a^2-f*a^2+1,2*d*f*e*a^2,-f*a^2;-a^2,-a^3*e*d,-e^2*a^2+a^2+t/pi*e^2*a^2*(j+k),2*d*e*a^2+t/pi*(h+i)*e^2*a^2;2*d*e*a^2,-a^2,-2*d*e*a^2-t/pi*e^2*a^2*(h+i),-e^2*a^2+a^2+t/pi*e^2*a^2*(m+n)');
B=sym('0;b;-t/pi*e^2*a^2*l;t/pi*e^2*a^2*o');
D=A\B

运行后，提示出错：
??? Error using ==> maple
at offset 9, missing operator or `;`

Error in ==> sym.sym>char2sym at 394
         if str2num(char(maple('nops',['{' S '}']))) > 1

Error in ==> sym.sym at 92
   S = char2sym(x);

是B矩阵有问题吗？该如何解决呢？求高人指点啊

TDF

怎么没人回答啊:'(

ChaChing

我使用R2009a试试, 报错讯息并不一样!

??? Error using ==> sym.mldivide at 28
First dimensions must agree.
A是4*4sym, 而B是1*1sym, 怎能D=A\B?

引起个人好奇的是, LZ的A,B写法类似, 为何A可以顺利转换, 而B不行!?
所以特地以debug模式切入看个究竟?
发现其中逻辑竟然是以逗号(,)判断, 所以B的写法无逗号当然无法正常转换成4*1 !  使用逗号再转置即可正常转换成4*1

1. clc; clear
   
2. A=sym('-e^2*a^2-1,2*c*e*a+2*f*d*e*a^2,-f*a^2,-2*f*d*e*a^2;-2*c*e+2*f*d*e*a^2,e^2*a^2-f*a^2+1,2*d*f*e*a^2,-f*a^2;-a^2,-a^3*e*d,-e^2*a^2+a^2+t/pi*e^2*a^2*(j+k),2*d*e*a^2+t/pi*(h+i)*e^2*a^2;2*d*e*a^2,-a^2,-2*d*e*a^2-t/pi*e^2*a^2*(h+i),-e^2*a^2+a^2+t/pi*e^2*a^2*(m+n)');
   
3. B=sym('0,b,-t/pi*e^2*a^2*l,t/pi*e^2*a^2*o'); B=B';

复制代码

其实个人比较习惯使用syms

1. clc; clear; syms a b c d e f h j k l m n o t
   
2. A=[-e^2*a^2-1,2*c*e*a+2*f*d*e*a^2,-f*a^2,-2*f*d*e*a^2;-2*c*e+2*f*d*e*a^2,e^2*a^2-f*a^2+1,2*d*f*e*a^2,-f*a^2;-a^2,-a^3*e*d,-e^2*a^2+a^2+t/pi*e^2*a^2*(j+k),2*d*e*a^2+t/pi*(h+i)*e^2*a^2;2*d*e*a^2,-a^2,-2*d*e*a^2-t/pi*e^2*a^2*(h+i),-e^2*a^2+a^2+t/pi*e^2*a^2*(m+n)];
   
3. B=[0;b;-t/pi*e^2*a^2*l;t/pi*e^2*a^2*o];

复制代码

但无论那一种方式, 进入D=A\B;时就一直busy了, 其实个人并不意外, 那麼多参数那麼复杂式子!
个人以为可以数值就儘量数值!

还有标题请明确些!

TDF

好的，多谢前辈指点啊。其实后来我自己又运算了下。结果是出来了，就是有号几十页。一看到就很头疼啊。