线性调频信号模糊函数仿真，包含线性调频模糊函数的各类图形.zip

共1个文件

m：1个

版权申诉

5星 · 超过95%的资源 171 浏览量 2023-07-26 09:48:17 上传评论收藏 1KB ZIP 举报

线性调频（Linear Frequency Modulation，LFM）信号是一种广泛应用在雷达、通信和电子战中的信号类型。这种信号的特点是其频率随时间线性变化，因此也被称为 chirp 信号。模糊函数（Ambiguity Function，AF）是分析LFM信号特性的重要工具，它反映了信号的时频分辨率，对于理解和优化雷达系统的性能至关重要。模糊函数是通过计算LFM信号的自相关函数来定义的，用于衡量两个LFM信号在不同时间延迟和频率偏移下的相关程度。在雷达系统中，模糊函数的峰值对应于目标的距离和速度信息。对于LFM信号，由于其频率的变化，模糊函数具有良好的时频局部化特性，可以同时提供高的距离分辨力和多普勒分辨力。 MATLAB是一个强大的数值计算和可视化软件，常被用来进行LFM信号模糊函数的仿真。在这个压缩包文件中，"线性调频信号模糊函数仿真，包含线性调频模糊函数的各类图形.m" 是一个MATLAB脚本，它很可能包含了以下内容： 1. **LFM信号生成**：使用MATLAB的`awgn`或`frequencies`函数来创建线性调频信号，其中参数包括初始频率、终止频率、调频斜率以及信号长度等。 2. **自相关函数计算**：利用MATLAB的`xcorr`函数计算LFM信号的自相关函数，这一步骤是获取模糊函数的关键。 3. **模糊函数的计算**：通过自相关函数的二维傅里叶变换，我们可以得到模糊函数。MATLAB的`fft2`函数可以完成这一转换。 4. **可视化**：利用MATLAB的`imagesc`或`surf`函数将模糊函数的二维图像或者三维表面图展示出来，便于观察和分析。通常，模糊函数的中心会有一个强峰，表示零时延和零频率偏移，而其他峰值则表示可能的混淆情况。 5. **参数调整与实验**：脚本可能会包含不同参数设置，如调频斜率的改变，来研究它们如何影响模糊函数的形状和性能。 6. **性能评估**：通过分析模糊函数的形状，可以评估LFM信号在雷达探测中的性能，比如距离分辨力和速度分辨力。这个MATLAB脚本提供了一个全面理解LFM信号模糊函数的平台，可以帮助我们深入研究和优化线性调频信号在雷达系统中的应用。通过运行和修改这个脚本，你可以更直观地了解LFM信号的特性，并进行相关实验。

资源推荐

资源详情

资源评论

收起资源包目录

线性调频信号模糊函数仿真，包含线性调频模糊函数的各类图形.zip （1个子文件）

线性调频信号模糊函数仿真，包含线性调频模糊函数的各类图形.m 5KB

%%%%%%%% T_3.m %%%% clear all clc clf taup=1; %脉冲宽度 100us b=10; %带宽 up_down=-1; %up_down=-1正斜率， up_down=1负斜率 x=lfm_ambg(taup,b,up_down); %计算模糊函数 taux=-1.1*taup:.01:1.1*taup; fdy=-b:.01:b; figure(1) mesh(100*taux,fdy./10,x) %画模糊函数 xlabel('Delay - \mus') ylabel('Doppler - MHz') zlabel('| \chi ( \tau,fd) |') title('模糊函数') figure(2) contour(100.*taux,fdy./10,x) %画等高线 xlabel('Delay - \mus') ylabel('Doppler - MHz') title('模糊函数等高线') grid on N_fd_0=(length(fdy)+1)/2; % fd=0 的位置 x_tau=x(N_fd_0,:); % 时间模糊函数 figure(3) plot(100*taux,x_tau) axis([-110 110 0 1]) xlabel('Delay - \mus') ylabel('| \chi ( \tau,0) |') title(' 时间模糊函数') grid on N_tau_0=(length(taux)+1)/2; % tau=0 的位置 x_fd=x(:,N_tau_0); % 速度模糊函数 figure(4) plot(fdy./10,x_fd) xlabel('Doppler - MHz') ylabel('| \chi ( 0,fd) |') title(' 速度模糊函数') grid on x_db=20*log10(x+eps); [I,J]=find(abs(x_db+6)<0.09); %取6db点的位置 I=(I-b/.01)/(1/.01); %Doppler维坐标变换 J=(J-1.1*taup/.01)/(1/.01); %时间维坐标变换 figure(5) %6db 的等高线 plot(J*100,I/10,'.') axis([-110 110 -1 1]) xlabel('Delay - \mus') ylabel('Doppler - MHz') title('模糊函数 6db 的等高线') grid on %- - - - 模糊函数 - - - function x=lfm_ambg(taup,b,up_down) % taup 脉冲宽度; % b 带宽; %up_down=-1正斜率， up_down=1负斜率 eps=0.0000001; i=0; mu=up_down*b/2./taup; for tau=-1.1*taup:.01:1.1*taup i=i+1; j=0; for fd=-b:.01:b j=j+1; val1=1-abs(tau)/taup; val2=pi*taup*(1-abs(tau)/taup); val3=(fd+mu*tau); val=val2*val3+eps; x(j,i)=abs(val1*sin(val)/val); end end %%%%%%%%%%%%%%%%%% %%%%%%%% T_3.m %%%% clear all clc clf taup=1; %脉冲宽度 100us b=10; %带宽 up_down=-1; %up_down=-1正斜率， up_down=1负斜率 x=lfm_ambg(taup,b,up_down); %计算模糊函数 taux=-1.1*taup:.01:1.1*taup; fdy=-b:.01:b; figure(1) mesh(100*taux,fdy./10,x) %画模糊函数 xlabel('Delay - \mus') ylabel('Doppler - MHz') zlabel('| \chi ( \tau,fd) |') title('模糊函数') figure(2) contour(100.*taux,fdy./10,x) %画等高线 xlabel('Delay - \mus') ylabel('Doppler - MHz') title('模糊函数等高线') grid on N_fd_0=(length(fdy)+1)/2; % fd=0 的位置 x_tau=x(N_fd_0,:); % 时间模糊函数 figure(3) plot(100*taux,x_tau) axis([-110 110 0 1]) xlabel('Delay - \mus') ylabel('| \chi ( \tau,0) |') title(' 时间模糊函数') grid on N_tau_0=(length(taux)+1)/2; % tau=0 的位置 x_fd=x(:,N_tau_0); % 速度模糊函数 figure(4) plot(fdy./10,x_fd) xlabel('Doppler - MHz') ylabel('| \chi ( 0,fd) |') title(' 速度模糊函数') grid on x_db=20*log10(x+eps); [I,J]=find(abs(x_db+6)<0.09); %取6db点的位置 I=(I-b/.01)/(1/.01); %Doppler维坐标变换 J=(J-1.1*taup/.01)/(1/.01); %时间维坐标变换 figure(5) %6db 的等高线 plot(J*100,I/10,'.') axis([-110 110 -1 1]) xlabel('Delay - \mus') ylabel('Doppler - MHz') title('模糊函数 6db 的等高线') grid on %- - - - 模糊函数 - - - function x=lfm_ambg(taup,b,up_down) % taup 脉冲宽度; % b 带宽; %up_down=-1正斜率， up_down=1负斜率 eps=0.0000001; i=0; mu=up_down*b/2./taup; for tau=-1.1*taup:.01:1.1*taup i=i+1; j=0; for fd=-b:.01:b j=j+1; val1=1-abs(tau)/taup; val2=pi*taup*(1-abs(tau)/taup); val3=(fd+mu*tau); val=val2*val3+eps; x(j,i)=abs(val1*sin(val)/val); end end %%%%%%%%%%%%%%%%%% %%%%%%%% T_3.m %%%% clear all clc clf taup=1; %脉冲宽度 100us b=10; %带宽 up_down=-1; %up_down=-1正斜率， up_down=1负斜率 x=lfm_ambg(taup,b,up_down); %计算模糊函数 taux=-1.1*taup:.01:1.1*taup; fdy=-b:.01:b; figure(1) mesh(100*taux,fdy./10,x) %画模糊函数 xlabel('Delay - \mus') ylabel('Doppler - MHz') zlabel('| \chi ( \tau,fd) |') title('模糊函数') figure(2) contour(100.*taux,fdy./10,x) %画等高线 xlabel('Delay - \mus') ylabel('Doppler - MHz') title('模糊函数等高线') grid on N_fd_0=(length(fdy)+1)/2; % fd=0 的位置 x_tau=x(N_fd_0,:); % 时间模糊函数 figure(3) plot(100*taux,x_tau) axis([-110 110 0 1]) xlabel('Delay - \mus') ylabel('| \chi ( \tau,0) |') title(' 时间模糊函数') grid on N_tau_0=(length(taux)+1)/2; % tau=0 的位置 x_fd=x(:,N_tau_0); % 速度模糊函数 figure(4) plot(fdy./10,x_fd) xlabel('Doppler - MHz') ylabel('| \chi ( 0,fd) |') title(' 速度模糊函数') grid on x_db=20*log10(x+eps); [I,J]=find(abs(x_db+6)<0.09); %取6db点的位置 I=(I-b/.01)/(1/.01); %Doppler维坐标变换 J=(J-1.1*taup/.01)/(1/.01); %时间维坐标变换 figure(5) %6db 的等高线 plot(J*100,I/10,'.') axis([-110 110 -1 1]) xlabel('Delay - \mus') ylabel('Doppler - MHz') title('模糊函数 6db 的等高线') grid on %- - - - 模糊函数 - - - function x=lfm_ambg(taup,b,up_down) % taup 脉冲宽度; % b 带宽; %up_down=-1正斜率， up_down=1负斜率 eps=0.0000001; i=0; mu=up_down*b/2./taup; for tau=-1.1*taup:.01:1.1*taup i=i+1; j=0; for fd=-b:.01:b j=j+1; val1=1-abs(tau)/taup; val2=pi*taup*(1-abs(tau)/taup); val3=(fd+mu*tau); val=val2*val3+eps; x(j,i)=abs(val1*sin(val)/val); end end %%%%%%%%%%%%%%%%%%

评论收藏

内容反馈

版权申诉