求助利用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)*
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/23 14:42:10
求助利用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)*
求助利用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)
求助利用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)*
这个数字精度应该是可以设置的,float型或者是double型都可以,这样就不会取整了
cuoleba