求解方程组

cboboc

这里是要解一个四个未知数六个方程的方程组
我写的是syms a1 a3 b1 b3;
[a1 a3 b1 b3]=solve('-sigma*a1+k*h*b1+(k/2)*(a2*b1+b2*a1+a3*b2+a2*b3)=0','-sigma*a2+k*h*b2+(k/2)*(a1*b1+a3*b1+a1*b3)=0','-sigma*a3+k*h*b3+(h/2)*(a2*b1+b2*a1)=0','-sigma*b1+(k/4)*(b1*b1)+(k/2)*(b3*b1)+g*k*a2+(B-(1/3))*(h^2)*4*(k^2)*sigma*b2-B*g*(h^3)*4*(k^3)*a2=0','-sigma*b1+(k/2)*(b2*b1+b3*b2)+g*k*a1+(B-(1/3))*(h^2)*(k^2)*sigma*b1-B*g*(h^3)*(k^3)*a1=0','-sigma*b3+(k/2)*b2*b1+g*k*a3+(B-(1/3))*(h^2)*9*(k^2)*sigma*b3-B*g*(h^3)*9*(k^3)*a3=0','a1','a3','b1','b3');
但是提示我：Warning: Explicit solution could not be found，应该怎样改呢，原题目是

[ 本帖最后由 ChaChing 于 2010-3-7 12:44 编辑 ]

elitecn

你是不是题目看错了？6个方程解4个未知数？题目中应该意思是a1--a3，b1--b3是未知数，不然怎么会用到6个方程？

cboboc

谢谢你这么快回复我，我是想求出a1，a3，b1，b3，不就能求出来他们的差了吗 ？

ChaChing

a1-a3, b1-b3是未知数
应该是这意思
a1~a3, b1~b3是未知数
所以有6个未知数!? 并非要求出来他们的差!?

cboboc

嗯是的。因为是帮给别人解的，看错了题目。楼上说的是对的，谢谢你

siyaoming

是啊，我也在解这方面的方程，苦于一直找不到类似的程序，离自己主动编程的距离有差的太远，等待中。。。。。。。

dingd

如果是数值解，可用1stOpt试试，求解方程比Matlab强大方便的多！