求助利用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/05/04 19:44:43
求助利用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)=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