LFM回波信号仿真,加汉明窗,可用的matlab代码

运行结果(文尾有程序)

在数字信号处理过程中,每次FFT变换只能对有限长度的时域数据进行变换,因此,需要对时域信号进行信号截断。
即使是周期信号,如果截断的时间长度不是周期的整数倍(周期截断),那么,截取后的信号将会存在泄漏。为了将这个泄漏误差减少到最小程度(注意我说是的减少,而不是消除),我们需要使用加权函数,也叫窗函数。

理想频响是分段恒定的,在边界频率处有突变点,所以,这样得到的理想滤波器单位脉冲响应hd(n)是无限长序列,而且是非因果的。但FIR的h(n)是有限长的,怎样用一个有限长的序列去近似无限长的hd(n)?最简单的办法是直接截取一段 hd(n) 代替 h(n) 。这种截取可以形象地想像为h(n)是通过一个“窗口”所看到的一段hd(n),因此,h(n)也可表达为hd(n)乘以一个“窗函数”

加窗主要是为了使时域信号似乎更好地满足FFT处理的周期性要求,减少泄漏。这里面,有很多窗函数可以使用,不过,每种窗函数的使用场景是不同的。由于直接对信号(加矩形窗)截断会产生频率泄露,为了改善频率泄露的情况,加非矩形窗,一般都是加汉明窗,因为汉明窗的幅频特性是旁瓣衰减较大,主瓣峰值与第一个旁瓣峰值衰减可达40db。



%% Pulse Compress
clear all; close all; clc;% LFM parameter
B       =   4e+6;       % 4MHz带宽
Tao     =   200e-6;     % 200us时宽
T       =   2e-3;       % 2ms脉冲重复周期
fs      =   8e+6;       %采样频率
SNR     =   20;         % 信噪比20dB
dis     =   T*fs/2;     % 将目标设置在回波中间处% Generate LFM
t = -round(Tao*fs/2):1:round(Tao*fs/2)-1; % 脉冲采样点
lfm = (10^(SNR/20))*exp(1i*pi*B/Tao*(t/fs).^2);figure;
subplot(2,1,1); plot(real(lfm),'b'); title('线性调频信号实部');
subplot(2,1,2); plot(imag(lfm),'r'); title('线性调频信号虚部');% Generate echo
echo  = zeros(1,T*fs);
echo(dis:1:dis+Tao*fs-1) = lfm;
noise = normrnd(0,1,1,T*fs) + 1i*normrnd(0,1,1,T*fs);
echo = echo + noise;figure;
subplot(2,1,1); plot(real(echo),'b'); title('回波信号实部');
subplot(2,1,2); plot(imag(echo),'r'); title('回波信号虚部');% Generate filter coeff
coeff = conj(fliplr(lfm)).* hamming(Tao*fs).'; %翻转共轭/加汉明窗figure; freqz(coeff);% pulse compress
fft_n = 2^(floor(log2(T*fs)) + 1);
pc_res = ifft(fft(echo,fft_n).*fft(coeff,fft_n)); % 未截取不完全滤波点figure;
plot(db(abs(pc_res)/max(abs(pc_res))),'r');  title('脉压结果分贝图');

LFM回波信号仿真,加汉明窗,可用的matlab代码相关推荐

  1. 复杂海洋环境下运动目标回波信号建模及仿真

    [摘  要]面向现代海战场智能化.无人化.网络化的发展趋势,UUV承担着重要任务和使命.由于UUV工作环境复杂.干扰因素较多,传统回波信号仿真中自由场假设将不再适用.针对UUV等高动态目标实时探测问题 ...

  2. 基于matlab的回波,基于MATLAB回波信号产生与消除.doc

    基于MATLAB回波信号产生与消除 摘 要 MATLAB可以进行矩阵运算.绘制函数和数据.实现算法.创建用户界面.连接其他编程语言的程序等,主要应用于工程计算.控制设计.信号处理与通讯.图像处理.信号 ...

  3. wcdma matlab程序,一种基于MATLAB的WCDMA信号仿真方法.doc

    目 录 1 引言 ----------------------- 6 2 WCDMA系统 ------------------- 7 2.1 WCDMA通信系统的模型与结构 ----------- 7 ...

  4. python信号处理加汉明窗_SciPyTutorial-非整周期信号加窗

    17. Scipy Tutorial- 非整周期信号加窗 17.1 什么是加窗? 使用FFT分析信号的频率成分时,分析的是有限的数据集合. FFT认为波形是一组有限数据的集合,一个连续的波形是由若干段 ...

  5. SAR成像系列:【3】合成孔径雷达(SAR)的二维回波信号与简单距离多普勒(RD)算法 (附matlab代码)

    合成孔径雷达发射信号以线性调频信号(LFM)为基础,目前大部分合成孔径雷达都是LFM体制,为了减轻雷达重量也采用线性调频连续波(FMCW)体制:为了获得大带宽亦采用线性调频步进频(FMSF)体制. ( ...

  6. 雷达辐射源调制信号仿真

    雷达辐射源调制信号仿真 说明:通过Matlab进行单载频(CW),线性调频(LFM).非线性调频(NLFM).二相编码(BPSK).四相编码(QPSK).二频编码(BFSK).四频编码(QFSK)七种 ...

  7. 基于线性调频(LFM)信号分数阶傅里叶变换的多径时延估计

    仿真信号 仿真信号为线性调频信号,信号采样率为2kHz. 线性调频(LFM)信号参数为:初始频率为300Hz,带宽400Hz,脉宽1s,调频斜率为 k ^ = − F s Δ t c o t ( p ...

  8. [解疑][TI]TI毫米波雷达系列(三):调频连续波雷达回波信号3DFFT处理原理(测距、测速、测角)

    本文是经过参考多个文章并整理的,相关程序已经经过验证其可行性.在此感谢原文作者(文末有相关链接)的无私分享. 1.测距.测速 毫米波雷达测距主要是通过检测回波时延来计算目标距离:测速是通过检测目标运动 ...

  9. 信号归一化功率_线性调频(LFM)信号仿真分析

    线性调频(LFM)信号是指瞬时频率随时间成线性变化的信号.(设振幅归一化,初始相位为零).线性调频信号也称为鸟声(Chirp)信号,因为其频谱带宽落于可听范围,听着像鸟声,所以又称Chirp扩展频谱( ...

  10. 开源夏令营《基于HackRF开发GPS信号仿真模拟器》工作总结(三)

    2014.07.21-2014.07.27,项目已经开展三周了,接下来就要面对中期检查了.实验室已经放假了,现在只剩下为数不多的几个人了,其中包括我. 这周主要解决的问题 这周主要解决的问题分为三个方 ...

最新文章

  1. dell 服务器系统6,dell C6100服务器安装windows2008R2系统
  2. 网络故障排除工具NeoTraceProTrial325
  3. IT培训课、视频教程和书本之PK
  4. Android 混淆打包
  5. php odbc 分页,用php实现odbc数据分页显示一例_php技巧
  6. 【转载】ASP.NET自定义404和500错误页面
  7. 做了极度危险的事情各种奔忙
  8. Visio中插入公式
  9. 毕业设计 - 题目:基于深度学习的中文汉字识别 - 深度学习 卷积神经网络 机器视觉 OCR
  10. mysql with,mysql+with
  11. IM即时通讯需要解决的问题
  12. 电脑常用快捷键及组合键
  13. 查看电脑系统是否永久激活
  14. strchr、strstr函数
  15. 什么是Blazor?
  16. 录音文件转换成文字的方法
  17. 深入浅出使用python编程_深入浅出Python元编程
  18. 编码(1)学点编码知识又不会死:Unicode的流言终结者和编码大揭秘
  19. 汽车理论第六章matlab,matlab汽车理论第一章
  20. 杰理科技AD15N系列玩具类音频开发

热门文章

  1. 点分十进制IP和网络字节序IP互转
  2. js刻度尺插件_JS鼠标拖动刻度标尺(尺子)效果
  3. RK3568-ANDROID11-4G-EC20-驱动篇(移远模块)
  4. 用 Javascript 生成二维码
  5. 8uftp使用教程 图解,8uftp使用教程 图解
  6. php代挂程序什么原理,小新云全套代挂加速 来源于网络资源简介:01、本源码调用官方接口 联合开发网 - pudn.com...
  7. 用C语言写的爬虫项目
  8. 做微信小程序商城需要什么材料?
  9. python导入第三方库失败_史上最详细 Python第三方库添加方法 and 错误解决方法
  10. photoshop--抽出滤镜