雷达模糊函数

admin

[mw_shl_code=applescript,true]clear all
clc
clf
taop = 1;   %脉冲宽度 100us
B = 10;        %带宽 10MHz
up_down = -1;    %正调频
x = lfm_ambg(taop,B,up_down);%计算模糊函数
taux = -1.1*taop:0.01:1.1*taop;
fdy = -B:0.01:B;
figure(1)
mesh(100*taux,fdy./10,x);
xlabel('Delay - \mus');
ylabel('Doppler - MHz');
zlabel('|\chi(\tau,fd)|');
title('模糊函数');
figure(2)
contour(100*taux,fdy./10,x)
xlabel('Delay - \mus');
ylabel('Doppler - MHz');
title('模糊函数等高线');
grid on
N_fd_0 = (length(fdy) + 1)/2;%fd = 0位置
x_tau = x(N_fd_0,;%时间模糊函数
figure(3)
plot(100*taux,x_tau)
axis([-110 110 0 1]);
xlabel('Delay - \mus')
ylabel('|\chi(\tau,0)|')
title('时间模糊函数')
grid on
N_tau_0 = (length(taux) + 1)/2;% tau = 0位置
x_fd = x(:,N_tau_0);%速度模糊函数
figure(4)
plot(fdy./10,x_fd);
xlabel('Doppler - MHz');
ylabel('|\chi(0,fd)|')
title('速度模糊函数')
grid on
x_db = 20*log10(x+eps);
[I,J] = find(abs(x_db+6)<0.09);
I = (I-B/0.01)/(1/0.01);
J = (J-1.1*taop/0.01)/(1/0.01);
figure(5)
plot(J*100,I/10,'.');
xlabel('Delay - \mus');
ylabel('Doppler - MHz');
axis([-110 110 -1 1]);
title('模糊函数6dB的等高线')
grid on[/mw_shl_code]