LFM.zip_LFM 模糊函数_LFM 模糊图_LFM模糊函数_lfm模糊度_模糊图等高线

close all;clear;clc; %% 线性调频信号 f0=10e3; %中心频率 B=10e3; T=10e-3; K=B/T; fs=8*f0; N=T*fs; dt=1/fs; t=-T/2:dt:T/2-dt; u=1/sqrt(T)*rectpuls(t/T).*exp(1j*2*pi*(f0*t+K/2*t.^2/2)); %发射信号 y=fftshift((fft(u,N))); figure(1); subplot(211),plot(t,real(u)); title('信号'),xlabel('时间'),ylabel('信号幅度'); subplot(212),plot(-fs/2:fs/N:fs/2-fs/N,abs(y)); title(['采样频率为',num2str(fs),'Hz']),xlabel('频率-HZ'),ylabel('信号幅度'); %% 模糊函数 tm=1; %脉冲宽度 b=10; %多普勒带宽 m=-1; %m=-1正斜率， m=1负斜率 tao=-tm:.01:tm; fd=-b/2:.01:b/2; M=m*b/2/tm; for ii=1:length(tao) for jj=1:length(fd) t=1-abs(tao(ii))/tm; phi=pi*(fd(jj)+M*tao(ii))*tm*(1-abs(tao(ii))/tm); X(jj,ii)=abs(t*sin(phi)/phi); end end %% 画模糊函数 figure(2),mesh(100*tao,fd./10,X) xlabel('时延\tau '); ylabel('归一化多普勒fd'); zlabel('| \chi ( \tau,fd) |'); title('模糊函数'); %% 画等高线 figure(3) contour(100.*tao,fd./10,X) xlabel('时延\tau '); ylabel('归一化多普勒fd'); title('模糊函数等高线') grid on %% 时间模糊函数 N_fd_0=(length(fd)+1)/2; % fd=0 的位置 x_tau=X(N_fd_0,:); figure(4) plot(100*tao,x_tau) axis([-110 110 0 1]) xlabel('时延\tau '); ylabel('| \chi ( \tau,0) |') title(' 时间模糊函数') grid on %% 速度模糊函数 N_tau_0=(length(tao)+1)/2; % tao=0 的位置 X_fd=X(:,N_tau_0); figure(5) plot(fd./10,X_fd) xlabel('归一化多普勒fd'); ylabel('| \chi ( 0,fd) |') title(' 速度模糊函数') grid on %% 6db 的等高线 X_db=20*log10(X+eps); [I,J]=find(abs(X_db+3)<0.09); %取6db点的位置 I=(I-b/2/.01)/(1/.01); %Doppler维 坐标变换 J=(J-1.1*tm/.01)/(1/.01); %时间维 坐标变换 figure(6) plot(J*100,I/10,'.') axis([-110 110 -1 1]) xlabel('时延\tau '); ylabel('归一化多普勒fd'); title('模糊函数 3db 的等高线') grid on