xwbm.m

来自「各种雷达脉冲压缩信号的时域」· M 代码 · 共 26 行

close all;clear all;clc;
T1=5;
P=13;

x=-5:0.08:5;
y=-5:0.08:5;

x1=x(1):((x(end)-x(1))/P):x(end);
y1=y(1):((y(end)-y(1))/P):y(end);

[X,Y]=meshgrid(x1,y1);

for m=0:(P-1)
{   a(m+1)=exp(j*pi*Y(m+1)*(T1-(X(m+1)-m*T1)));
    b(m+1)=pi*Y(m+1)*(T1-abs(X(m+1)-m*T1))+eps;
    c(m+1)=(sin(b)./b).*(1-abs(X(m+1)-m*T1)/T1)/sqrt(P);

    d=[1,1,1,1,1,0,0,1,1,0,1,0,1];
    for k=0:(P-1-m)
    e=d(k+1)*d(k+1+m)*exp(j*2*pi*T1);
    Z(k+1)=a(k+1)*c(k+1)*e
    end
}
end

surf(X,Y,Z),title('线性调频信号模糊图'),xlabel('时间'),ylabel('频率')