已知我船此刻的经纬度，和同一时刻目标位于我船的舷角以及声纳接收到目标的信号频率，基于多普勒频移解算目标距离，速度，航向，生成matlab代码

以下是基于多普勒频移解算目标距离、速度和航向的Matlab代码：

% 已知参数
f0 = 3000; % 发射频率
f1 = 3090; % 接收频率
lat1 = 40.0; % 我船纬度
lon1 = 120.0; % 我船经度
bearing = 90; % 舷角
lat2 = 41.0; % 目标船纬度
lon2 = 121.0; % 目标船经度

% 地球半径
R = 6371000;

% 计算目标和我船之间的距离
lat1 = lat1 * pi / 180;
lon1 = lon1 * pi / 180;
lat2 = lat2 * pi / 180;
lon2 = lon2 * pi / 180;
dlat = lat2 - lat1;
dlon = lon2 - lon1;
a = sin(dlat/2)^2 + cos(lat1) * cos(lat2) * sin(dlon/2)^2;
c = 2 * atan2(sqrt(a),sqrt(1-a));
d = R * c;

% 计算目标和我船之间的方位角
y = sin(dlon) * cos(lat2);
x = cos(lat1) * sin(lat2) - sin(lat1) * cos(lat2) * cos(dlon);
brng = atan2(y,x) * 180 / pi;

% 计算目标的速度和航向
v = (f1 - f0) * 1540 / (2 * f0);
hdg = brng - bearing;
if hdg < 0
    hdg = hdg + 360;
end

% 输出结果
fprintf('距离：%f m\n',d);
fprintf('速度：%f m/s\n',v);
fprintf('航向：%f deg\n',hdg);

代码中使用了Haversine公式计算了目标和我船之间的距离和方位角，然后根据多普勒频移计算了目标的速度和航向。注意，该代码中假设了声速为1540 m/s。