a1.rar_LFMCW_echo_lfmcw时频图_radar_锯齿LFMCW

clear all; %清除变量 A=1; F0=20e6; T=200e-6; B=8e6; K=B/T; Fs=100e6; N=round(T*Fs);%采样点个数 t=linspace(0,T,N);%横坐标范围设定 %St=A*exp(j*2*pi*(F0*t+(K*t.^2)/2));%单周期chirp信号,设A=1，f0=20Mhz,B=8MHz, T=200us, 可得调频斜率等于4.0x10^10。 St=cos(2*pi*(F0*t+(K*t.^2)/2)); figure(1) plot(linspace (0,T,N)*1e6,real(St)) xlabel('时间/us'); ylabel('幅度'); title('单周期线性调频信号时域图'); figure(2) freq=linspace(-Fs/2, Fs/2, N); plot (freq*1e-6, fftshift(abs(fft(St)))); axis([0, 50 0,600]); xlabel('频率/MHz' ); ylabel('幅度'); title('单周期线性调频信号频谱图'); Si=St; m=100;%扩展周期数 for k=1:1:(m-1)%周期延拓 Si=[Si, St]; end figure (3) %周期延拓后的时域图 plot(linspace (0, m*T, m*N)*1e6, real(Si)); xlabel('时间/us' ); ylabel('幅度'); title('多周期线性调频信号时域图'); figure (4) %FFT频谱图 freq=linspace( -Fs/2, Fs/2, m*N) ; plot (freq*1e-6, fftshift(abs(fft(Si)))); axis([0, 50 0,6000]); xlabel('频率MHz'); ylabel('幅度'); title('多周期线性调频信号频谱图'); grid on; n=length(Si);%取信号长度 M=16;%每次滑动保留重叠点数 L=512;%窗口长度 k=L-M; %实际每次滑动点数 t=fix (n/k) ;%滑动次数 length1=L; r= (rectwin (length1))' ;%矩形窗 for a=1: t%窗口滑动进行STFT n1=(L-M)*(a-1)+1; n2= (L-M)* (a- 1)+length1; y=Si (n1:n2) ;%截取一段信号 sf=fft(y.*conj(r), L);%信号片段与窗函数的复共轭相乘 Y(n1:n2)=sf(1:length1);%结果赋给Y [g h]=max (abs (sf));%幅值最大点 Q(a)=h/L*Fs;%幅值最大的频率点作为本窗口内的频率 end figure (5);%作出STFT频谱图 plot(abs(Y)); title('STFT频谱图'); xlabel('频率/Hz'); ylabel('幅度'); figure (6)%多周期chirp信号时频分析图 x=([0:t-1]*1e-3/(t/m)); plot(x*1e6, abs(Q),'-*') xlabel('时间/us'); ylabel ('频率/Hz'); title('多周期线性调频信号时频分析图');