matlab求三元方程组带三角函数

Question

代码如下 我想求a,b,t b是角度

>> syms a b t
>> eq1 = a*cos(b)/1020-t

eq1 =

(a*cos(b))/1020 - t

>> eq2 = a*sin(b)/1200-t

eq2 =

(a*sin(b))/1200 - t

>> eq3 = a*sin(b)/[a*cos(b)+1.89]-1200/1020

eq3 =

(a*sin(b))/(a*cos(b) + 189/100) - 20/17

Answer 1

用matlab求带三角函数的三元方程组，可以用fsolve()函数求出其近似值。其求解方法

[x,fval] = fsolve(solve_fun,x0)

其中：solve_fun——三元方程组自定义函数，x0——初值，x——变量值，fval——等式值

>> x0=[0.4447    0.6154    0.7919];

>> solve_fun=@(x)[x(1)*cos(x(2)*pi/180)/1020-x(3);x(1)*sin(x(2)*pi/180)/1200-x(3);x(1)*sin(x(2)*pi/180)/(x(1)*cos(x(2)*pi/180)+1.89)-1200/1020];

>> [x,fval] = fsolve(solve_fun,x0)

运行结果（注意，计算时用将角度转换成弧度）