% 简单脉冲信号的模糊函数
clc
close all;
clear all;
T = 1;
t = -T:5e-2:T;
fd = -10*T:2e-1:10*T;
fd = fd + eps;
[xx,yy] = meshgrid(t,fd);
AF = abs((sin(pi*yy.*(T-abs(xx))))./(T*pi*yy));
figure(1)
mesh(xx,yy,AF);
title('脉宽为T,具有单位能量简单脉冲的模糊函数');
xlabel('时延');
ylabel('多普勒频移');
zlabel('模糊函数AF');
figure(2)
v = [10^(-0.5/20),10^(-3/20),10^(-10/20),10^(-20/20)];%%%%%%%%
contour(xx,yy,AF,v);
text(0.5,-8,{'-0.5,-3,-10和','-20dB等高线图'},'Edgecolor',[0 0 0]);%%%%%%
title('简单脉冲模糊函数的等高线图');
xlabel('时延');
ylabel('多普勒频移');
t = -T:1e-6:T;
At = 1-abs(t);
figure(3)
plot(t,At);
title('简单脉冲AF的零多普勒截线');
xlabel('时延');
ylabel('幅值');
fd = -10*T:1e-2:10*T;
fd = fd + eps;
Af = abs((sin(pi*fd*T))./(T*pi*fd));
figure(4)
plot(fd,Af);
title('零延迟截线');
xlabel('多普勒频移');
ylabel('幅值');
% t = -T:1e-6:T;
figure(5)
fd = 0.31;
Atf = abs((sin(pi*fd*(T-abs(t))))./(T*pi*fd));
plot(t,Atf,'r')
hold on
fd = 1.72;
Atf = abs((sin(pi*fd*(T-abs(t))))./(T*pi*fd));
plot(t,Atf,'g')
hold on
fd = 7.5;
Atf = abs((sin(pi*fd*(T-abs(t))))./(T*pi*fd));
plot(t,Atf,'b')
title('时延剖面叠加');
xlabel('时延');
ylabel('幅值');