用solve解方程组问题求助

zhizi123410

>> N=20;h=0.2;J=10;k=2;r=pi/2;
>> S=solve('sin(a-r)-N*sin(a+b)+N/2*(N-1)*(u+w)*cos(a+b)-N/2*k*((N-1)*u*cos(2*a)-sin(2*a))-h*N*(1/2*(N-1)*u*cos(r-a)+sin(r-a))=0,-sin(b+r)+N*sin(a+b)+N/2*(N-1)*(u+w)*cos(a+b)-N/2*k*((N-1)*w*cos(2*b)-sin(2*b))+h*N*(1/2*(N-1)*w*cos(r+b)-sin(r+b))=0,-N/3*(2*N-1)*(u+w)*cos(a+b)+N/2*sin(a+b)+2*J*u-N*k/4*(sin(a*2)-2/3*u*(2*N-1)*cos(2*a))+h*N*(1/3*(2*N-1)*u*cos(r-a)+1/2*sin(r-a))=0,-N/3*(2*N-1)*(u+w)*cos(a+b)+N/2*sin(a+b)+2*J*w-N*k/4*(sin(b*2)-2/3*w*(2*N-1)*cos(2*b))+h*N*(1/3*(2*N-1)*w*cos(r+b)+1/2*sin(r+b))=0','a,b,u,w')
这是我的方程组，参数是我自己随便选择的，不知道为什么一直解不出来呢？有人让我用数值解法求解，但是这样不是可以求出数值解了嘛？为什么很久都得不到解呢？？:@Q :@Q :@Q :@Q

无水1324

那你用fsolve试一下，详细的可以在matlab中help fsolve