假设平面上有两个基站，发射站为A位于点（xA，yA），接收站为B位于点（xB，yB），接收站方位角为目标点的坐标为θ，写一个求解目标坐标（x，y）的matlab函数

假设目标点到基站A的距离为d1，到基站B的距离为d2，则有以下公式：

d1 = sqrt((x-xA)^2 + (y-yA)^2)
d2 = sqrt((x-xB)^2 + (y-yB)^2)

同时，根据方位角θ，可以得到以下方程：

tan(θ) = (y-yB)/(x-xB)

化简可得：

y = tan(θ) * x - tan(θ) * xB + yB

将y代入d2的公式中，可得到关于x的一元二次方程：

(x-xB)^2 + (tan(θ)*x - tan(θ)*xB + yB - y)^2 = d2^2

化简可得：

(x^2 - 2*xB*x + xB^2) + (tan(θ)^2 * x^2 - 2*tan(θ)^2 * x * xB + tan(θ)^2 * xB^2 - 2*tan(θ)*(y-yB)*x + 2*tan(θ)*(y-yB)*xB + (y-yB)^2 - d2^2) = 0

将其表示为ax^2 + bx + c = 0的形式，可得到：

a = 1 + tan(θ)^2
b = -2*xB - 2*tan(θ)*(y-yB)
c = xB^2 + (y-yB)^2 - d2^2 - (yB - tan(θ)*xB)^2

解一元二次方程即可得到目标点的坐标。

下面是对应的Matlab代码：

function [x, y] = calculate_target(xA, yA, xB, yB, theta, d1, d2)
    a = 1 + tan(theta)^2;
    b = -2*xB - 2*tan(theta)*(yA-yB);
    c = xB^2 + (yA-yB)^2 - d2^2 - (yB - tan(theta)*xB)^2;

    delta = b^2 - 4*a*c;

    if delta < 0
        x = NaN;
        y = NaN;
    else
        x1 = (-b + sqrt(delta)) / (2*a);
        x2 = (-b - sqrt(delta)) / (2*a);

        y1 = tan(theta) * x1 - tan(theta) * xB + yB;
        y2 = tan(theta) * x2 - tan(theta) * xB + yB;

        if abs(d1 - sqrt((x1-xA)^2 + (y1-yA)^2)) < abs(d1 - sqrt((x2-xA)^2 + (y2-yA)^2))
            x = x1;
            y = y1;
        else
            x = x2;
            y = y2;
        end
    end
end

其中，输入参数为xA、yA、xB、yB、theta、d1、d2，分别表示基站A、基站B的坐标，接收站方位角，目标到基站A、B的距离。输出为目标点的坐标x、y。如果无解，则返回NaN。