matlab求cos-1,求助利用matlab计算含三角函数的方程组的数值解,如何使数值解更为精确?13576000*(cos(x(1))+cos(a)-(1+b/2)^2*(1-b)*(cos(a)...

求助利用matlab计算含三角函数的方程组的数值解,如何使数值解更为精确?13576000*(cos(x(1))+cos(a)-(1+b/2)^2*(1-b)*(cos(a)+cos(x(2))))-x(3)*sin(c)=013576000*(sin(x(1))+(1+b/2)^2*(1-b)*sin(x(2)))-x(3)*cos(c)=0x(1)-(1+b/2)*(1-b)*

2015-11-13

求助利用matlab计算含三角函数的方程组的数值解,如何使数值解更为精确?

13576000*(cos(x(1))+cos(a)-(1+b/2)^2*(1-b)*(cos(a)+cos(x(2))))-x(3)*sin(c)=0

13576000*(sin(x(1))+(1+b/2)^2*(1-b)*sin(x(2)))-x(3)*cos(c)=0

x(1)-(1+b/2)*(1-b)*(pi-x(2))=0

利用[x,fv.ef.out.jac]=fsolve(@NSC,x0,[],pi/12,0.1,88*pi/180)

初值设为

0 ,3 ,70000000

得到的数值解不够精确,理论上x(1)应该略大于0,x(2)接近3.14,可惜得不到,

我是matlab菜鸟,如何进行次方程组的数值计算!

M文件：

function y=NSC(x,a,b,c)

y(1)=13576000*(cos(x(1))+cos(a)-(1+b/2)^2*(1-b)*(cos(a)+cos(x(2))))-x(3)*sin(c);

y(2)=13576000*(sin(x(1))+(1+b/2)^2*(1-b)*sin(x(2)))-x(3)*cos(c);

y(3)=x(1)-(1+b/2)*(1-b)*(pi-x(2));

matlab命令：

x0=[0,3,7e7]

[x,fv.ef.out.jac]=fsolve(@NSC,x0,[],pi/12,0.1,88*pi/180)