模糊函数仿真代码

clear all
close all
clc

fs=20e6;%采样率
t=10e-6;%脉宽
tt=0:1/fs:t-1/fs;%时间刻度
b=10e6;%带宽
Nz=round(b*t);
u=b/(2*t);%调频斜率
y=exp(1i*2*pi*u*tt.^2);%线性调频信号
len_y=length(y);
x=[zeros(1,2*len_y),y,zeros(1,2*len_y)];
x1=[zeros(1,len_y),y,zeros(1,len_y)];
N1=length(x1);
val1=sum(abs(y).^2);
tau1=-t*fs:1:t*fs;
fd1=-b*t:1:b*t;
k=0;
for tau=tau1
    k=k+1;
    l=0;
   for fd=fd1
       l=l+1;
        for m=1:N1
             val2(m)=conj(x(round(-tau1(1)+m+round(tau))));%时间因子
             val3(m)=exp(j*2*pi*fd*(m-1)/len_y);%频率因子
             val(m)=x1(m)*val2(m)*val3(m);
        end
    MF(l,k)=abs(sum(1/val1.*val));
   end
end 
figure
surfc(tau1/fs*1e6,fd1/t/1e6,MF,'edgecolor','none')  
xlabel('时延(us)')
ylabel('多普勒频移(MHz)')
zlabel('模糊函数值')
title('线性调频模糊函数','FontSize',14)