Question

怎样用matlab求解一个超越方程组

sin(x1)-sin(x2)=m1*0.601
sin(58/180*pi-2*x1-x2）-sin(x1)=m2*0.601
0<x1,x2<pi/2; m1,m2=0,+1,-1

Answer 1

求解这类超越方程组，可以通过循环语句，用fsolve（）函数求解其方程组的数值解。

求解结果

m1=-1  m2=-1

x1=0.22105     0.96186  x2=8.1657e-13  3.7059e-13

m1=-1  m2=0

x1=0.085333     0.75629  x2=6.4391e-10  4.9884e-13

m1=-1  m2=1

x1=-0.062527     0.56867  x2=2.4836e-13 -2.5296e-12

m1=0  m2=-1

x1=0.40647     0.40647  x2=8.8012e-08  8.6441e-08

m1=0  m2=0

x1=0.25307     0.25307  x2=7.2798e-10 -4.7684e-11

m1=0  m2=1

x1=0.085333    0.085333  x2=1.5973e-12 -1.7877e-12

根据0<x1,x2<pi/2;条件，可以自行确定其x1、 x2的值。

实施代码