线性调频信号模糊函数-多普勒敏感性-仿真

clear all;close all; 
% 线性调频信号模糊函数
%定义参数 
B=10e6;       % B ——带宽 
Fc=10e9;      % Fc——载频 
Fs=10*B;      % Fs——采样率 
T2=10e-6;     % T2——脉冲宽度 
Tr=375e-6;    % Tr——探测周期 
C=2.997e8; 

%产生线性调频脉冲信号 
mu=B/T2; 
NT2=[Fs*T2]; 
t=[-NT2/2:NT2/2-1]*T2/NT2;  
f=Fc+mu*t.*t/2; 
sig=exp（2*pi*j*f）; 
figure（1）; 
plot（t*1e6real（sig））; 
title（‘近程模式下线性调频脉冲信号‘）; 
xlabel（‘延时（一个脉宽时间内）/us‘）; 
ylabel（‘归一化幅度‘）; 
grid on; 
%线性调频脉冲信号频谱 
sigfft=fftshift（fft（sig））; 
Nf=length（sigfft）; 
f=[-Nf/2:Nf/2-1]*Fs/Nf; 
figure（2）; 
plot（f*1e-6abs（sigfft）/max（abs（sigfft）））; 
% plot（fabs（sigfft）/max（abs（sigfft）））; 
title（‘近程模式下的信号频谱图‘）; 
xlabel（‘频率/MHz‘）; 
ylabel（‘归一化幅度‘）; 
grid on; 


%模糊函数 
Naw=80; %？？？？？？？？？？、
Fdmax=B; %%？？？？？？？？？？？？？？？？？？？？
i=0; 
for t=[-Naw:Naw-1]*T2/Naw 
    i=i+1; 
    j=0; 
    for Fd=[-Na