%单基地方向余弦-多普勒谱
%-------------------------------------------------------------------------%
% Function : 单基地方向余弦-多普勒谱
% Date : 2018.05.29
% Author : ls
% Status : ok
%-------------------------------------------------------------------------%
clc
clear all
close all
%%
%参数定义（对于舰载高频地波雷达擦地角一般认为是0°）
f0=10e6; % 雷达工作频率
c=3e8; % 光速
g=9.8; % 重力加速度
lamda=c/f0; % 雷达工作波长
L=lamda/2; % 海浪波长
vb=sqrt(g*L/2/pi); % 海浪相速度
fMB=sqrt(g/pi/lamda); % 单基地Bragg散射频率
vr=30*1.852/36; % 船速度******************
theta_s = 0*pi/180; % 平台运动方向************************
theta_p = 0*pi/180; %偏航角
theta = (0:5:180)*pi/180; % 舰船运动速度与海杂波散射体之间的夹角
theta_r=theta-theta_s+theta_p; %实际波束指向
%%
%画图
figure;
fd1=fMB+2*vr*cos(theta_r)/lamda; %舰载一阶海杂波正方向多普勒频率
plot(fd1,cos(theta) ,'LineWidth',1.5);
hold on;
grid on;
fd2=-fMB+2*vr*cos(theta_r)/lamda; %舰载一阶海杂波负方向多普勒频率
plot(fd2,cos(theta), 'LineWidth',1.5);
hold on;
xlabel('归一化多普勒频率（Hz）');
ylabel('方向余弦');
title('单基地方向余弦-多普勒谱');
axis([-1 1 -1 1]);
