传统基于逗留相位原理的窗函数反求法-非线性调频（NLFM）信号设计，脉冲压缩

%传统基于逗留相位原理的窗函数反求法-NLFM信号 clc;close all;clear all; fs=100e6;%100Mhz fc=10e6;%10Mhz T=200e-6; B=10e6; N=T*fs; f=(-N/2:N/2-1)*B/N; W=0.08+0.92*cos(pi.*f/B).^2; %k=0.08,n=2时，为海明窗 % W=cos(pi*f/B).^2; %k=0,n=2时，为余弦二次方窗，即汉宁窗 % W=cos(pi*f/B).^3; %k=0,n=3时，为余弦三次方窗 % W=cos(pi*f/B).^4; %k=0,n=4时，为余弦四次方窗 另一版本w4=cos(0.4*2*pi.*t/T).^4 % W=0.333+0.667*cos(pi*f/B).^2; %k=0.333,n=2时为3:1锥比加权函数 % W=exp(-12.*(f/B).^2); %取k=1或12，为高斯窗 % W=0.42+0.5.*cos(2.*pi.*f/B)+0.08*cos(4*pi.*f/B) %布莱克曼（Blackman）窗 % W=1+0.8565*cos(2*pi*f/B)+0.0004*cos(4*pi*f/B)+0.00854*cos(6*pi.*f/B)-0.00386*cos(8*pi.*f/B)+0.00148*cos(10*pi.*f/B); %5阶泰勒（Taylor）窗 sumi=0; for i=1:N W(i)=W(i)*(B/N)+sumi; sumi=W(i); end K1=T/sumi;%常数 Tf=K1.*W; figure(1) subplot(211) plot(f,Tf);%群时延函数 grid on ; title('群时延函数Tf/f') xlabel('frequency/hz') ylabel('time/s') subplot(212); % f_t=interp1(Tf,f,'linear');%三次样条插值法 [p,~,mu]=polyfit(Tf,f,20); %多项式拟合 t=0:1/fs:T-1/fs; f_t=polyval(p,t,[],mu); plot(t,f_t); grid on ; title('时频函数f/Td') xlabel('time/s') ylabel('frequency/hz') sumj=0; %数值积分 for i=1:N f_t(i)=f_t(i)*(1/fs)+sumj; sumj=f_t(i); end phi=2*pi*f_t;%相位 s=exp(1j*phi);%S形NLFM信号 figure(2) subplot(211) plot(t,real(s)) title('NLFM基带信号实部'); subplot(212); plot(t,imag(s)); title('NLFM基带信号虚部'); figure(3) subplot(211); plot(t,phi); title('NLFM相位') S=fftshift(abs(fft(s))); subplot(212) plot(f,S);%NLFM频谱 title('NLFM信号频谱') %脉冲压缩 figure(4) ht=conj(s);%匹配滤波器 Y=fft(s).*fft(ht);%计算循环卷积 时域卷积，频域相乘 y=abs(fftshift(ifft(Y))); subplot(311) plot(t,y); title('脉压输出信号'); subplot(312); plot(t,20*log10(y/max(y)));%经扩展 dB title('脉压输出信号/dB'); grid on; subplot(313) plot(f,abs(fftshift(Y)));%脉压输出信号 title('脉压输出信号频谱');