Question

求助利用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)

Answer 1

这个数字精度应该是可以设置的，float型或者是double型都可以，这样就不会取整了

Answer 2

cuoleba