雷达回波脉冲压缩Python代码

雷达发射LFM信号
参数如下:
中心频率 f c = 10 G H z f_c=10GHz fc​=10GHz
脉冲宽度 T = 5 μ s T=5\mu s T=5μs
带宽 B W = 10 M H z BW=10MHz BW=10MHz
PRT是10倍的脉宽

点目标位于 R 0 = 2000 m R0=2000m R0=2000m处

雷达回波信号进行匹配滤波,实现脉冲压缩,在目标位置处出现了峰值

import numpy as np
import matplotlib.pyplot as plt# 雷达参数
fc = 10e9
T = 5e-6
BW = 10e6
K = BW/T
fs = 2*BW
PRT = 10*T
N = PRT*fs
t = np.linspace(0,PRT,N)
c = 3e8# 矩形脉冲函数
def rect_wave(x,c,c0):if x>=(c+c0):r = 0.0elif x<=c0:r = 0.0else:r = 1return r# 得到匹配滤波器的频域
aa = np.array([rect_wave(tt,T,0) for tt in t])
bb = np.exp(1j*np.pi*2*fc*t+1j*np.pi*K*np.power(t,2))
x = np.multiply(aa,bb) # 雷达发射信号
xf = np.fft.fftshift(np.fft.fft(x))
MF_f = np.conjugate(xf) # 匹配滤波器是发射信号频谱的共轭# 目标回波
R0 = 2000
tao = 2*R0/c
aa_tao = np.array([rect_wave(tt-tao,T,0) for tt in t])
bb_tao = np.exp(1j*np.pi*2*fc*(t-tao)+1j*np.pi*K*np.power((t-tao),2))
st = np.multiply(aa_tao,bb_tao)# 回波与匹配滤波器在频域相乘
sf = np.fft.fftshift(np.fft.fft(st))
Target_MF_f = np.multiply(sf,MF_f)
# 匹配滤波器的时域输出
Target_MF_t = np.fft.ifft(Target_MF_f)
plt.plot(c*t/2,np.abs(Target_MF_t))
plt.xlabel('Range(m)')
plt.ylabel('Matched Filter Output')
plt.title('Matched Filter Output in Time Domain')
plt.show()

得到的图形如下图所示

雷达回波脉冲压缩python代码相关推荐

  1. 数字图像处理作业: 包含质量因子的 JPEG压缩 python代码

    JPEG图片压缩的Python实现 注意:比较压缩率是 比较 原图与gpj文件的 大小 而非 解压后的bmp ,不管Q选择多少 解压后的bmp大小都是一样的. import numpy as np i ...

  2. 如何保护Python代码?

    我正在用Python开发一款软件,该软件将分发给我的雇主的客户. 我的雇主希望通过限时许可文件来限制软件的使用. 如果我们分发.py文件甚至.pyc文件,则很容易(反编译和)删除检查许可证文件的代码. ...

  3. python 解压js压缩的图片_20行Python代码,无损压缩千百张图片!Python有多强!就不用说了...

    每个设计师.摄影师或有图片处理需求小编,都会面临批量 高清大图的困扰. 因为高清大图放到网站上会严重拖慢加载速度,或是有的地方明确限制了图片大小,因此,为了完成工作,他们总是需要先把图片压缩,再上传. ...

  4. 60行Python代码搞定,王者荣耀雷达图!

    点击上方"菜鸟学Python",选择"星标"公众号 超级无敌干货第一时间推给你!!! 现有以下表格,我需要画成第二张图中的雷达图(类似fifa足球.王者荣耀里面 ...

  5. 涨见识了,在终端执行 Python 代码的 6 种方式

    作者 | BRETT CANNON 译者 | 豌豆花下猫@Python猫 为了我们推出的 VS Code 的 Python 插件[1],我写了一个简单的脚本来生成变更日志[2](类似于Towncrie ...

  6. DeepMind提图像生成的递归神经网络DRAW,158行Python代码复现

    作者 | Samuel Noriega 译者 | Freesia 编辑 | 夕颜 出品 | AI科技大本营(ID: rgznai100) [导读]最近,谷歌 DeepMInd 发表论文( DRAW: ...

  7. 158行Python代码复现:DeepMind提图像生成的递归神经网络DRAW

    授权自AI科技大本营(ID:rgznai100) 本文约5200字,建议阅读10+分钟. 本文作者基于代码实现系统的思路,详细阐述了 DRAW 的概念.架构和优势等. [ 导读 ]最近,谷歌 Deep ...

  8. 深度学习中的正则化技术(附Python代码)

    作者:SHUBHAM JAIN 翻译:和中华 校对:丁楠雅 本文约3500字,建议阅读20分钟. 本文简单介绍了什么是正则化以及在深度学习任务中可以采用哪些正则化技术,并以keras代码具体讲解了一个 ...

  9. 30秒内便能学会的30个超实用Python代码片段

    许多人在数据科学.机器学习.web开发.脚本编写和自动化等领域中都会使用Python,它是一种十分流行的语言. Python流行的部分原因在于简单易学. 本文将简要介绍30个简短的.且能在30秒内掌握 ...

最新文章

  1. 【 C 】在单链表中插入一个新节点的尝试(二)
  2. linux gcc make cmake 三工具的关系
  3. DIET:Dual Intent and Entity Transformer
  4. java performance tools / NetBeans Profiler / Sun BTrace / Eclipse MAT / IBM ISA
  5. Java小结(四)——折半查找、选择排序、冒泡排序
  6. 为什么说学人工智能一定要学Python?
  7. java与c/c++进行socket通信的一些问题
  8. C++ 读取单个字符
  9. C++ 用new 动态创建多维数组
  10. 又一国产数据库诞生!腾讯发布 TGDB,实时查询比 Neo4j 快 20-150 倍!
  11. Mac(不限于)中几个有内涵的工具
  12. 如何发布一个BT种子文件
  13. 把“友商”装进芯里威联通运行黑群晖最新DSM系统
  14. 红旗Linux系统安装教程
  15. 代理服务器没有响应解决方案
  16. (三)Position Rank代码解读(一)
  17. 浅谈AS2 传输/连接
  18. 【阿里云产品公测】以开发者角度看ACE服务『ACE应用构建指南』
  19. python 疯狂讲义 笔记(变量和基本类型)
  20. JXTA的安装和使用

热门文章

  1. 周末苦逼码代码,为css3的强大所颤抖了
  2. 浅谈H5业务场景下的人脸识别
  3. android4.4.2游戏,安卓模拟器4.4.2内核之上的穹顶之战
  4. Access-Control-Max-Age是什么
  5. “速课小龙”项目冲刺2
  6. 六月的天空,飘过大片云朵,抬头仰望,心中的那一抹温柔
  7. Linux上使用docker运行酷Q机器人
  8. CSS 实现input自定义样式--文本框
  9. 《CSS世界》读书笔记(九)
  10. launchScreen与launchImage 、APP图标问题