matlab求解三角函数方程组

上机实跑，随用随Ctrl+v

举例：

cosx+cos(x-y)＝1.96
sinx+sin(x-y)＝0.5
求x和y分别为多少度

代码如下：

function solve_equations()
    % 定义方程组函数,如果有系数记得加上*号
    fun = @(x) [cos(x(1)) + cos(x(1) - x(2)) - 1.96; ...
               sin(x(1)) + sin(x(1) - x(2)) - 0.5];
    
    % 初始猜测值（可能需要调整以找到解）
    x0 = [0, 0]; % 初始猜测值，这里设为0度，但可能需要更合理的猜测
    
    % 使用fsolve求解方程组
    options = optimoptions('fsolve','Algorithm','trust-region-reflective'); % 设置算法选项
    sol = fsolve(fun, x0, options);
    
    % 将弧度转换为度数
    sol_deg = rad2deg(sol);
    
    % 显示结果
    fprintf('近似解为：x = %f度, y = %f度\n', sol_deg(1), sol_deg(2));
end