Python+matplotlib实现计算两个信号的交叉谱密度实例


Posted in Python onJanuary 08, 2018

 计算两个信号的交叉谱密度

结果展示:

Python+matplotlib实现计算两个信号的交叉谱密度实例

完整代码:

import numpy as np
import matplotlib.pyplot as plt


fig, (ax1, ax2) = plt.subplots(2, 1)
# make a little extra space between the subplots
fig.subplots_adjust(hspace=0.5)

dt = 0.01
t = np.arange(0, 30, dt)

# Fixing random state for reproducibility
np.random.seed(19680801)


nse1 = np.random.randn(len(t))         # white noise 1
nse2 = np.random.randn(len(t))         # white noise 2
r = np.exp(-t / 0.05)

cnse1 = np.convolve(nse1, r, mode='same') * dt  # colored noise 1
cnse2 = np.convolve(nse2, r, mode='same') * dt  # colored noise 2

# two signals with a coherent part and a random part
s1 = 0.01 * np.sin(2 * np.pi * 10 * t) + cnse1
s2 = 0.01 * np.sin(2 * np.pi * 10 * t) + cnse2

ax1.plot(t, s1, t, s2)
ax1.set_xlim(0, 5)
ax1.set_xlabel('time')
ax1.set_ylabel('s1 and s2')
ax1.grid(True)

cxy, f = ax2.csd(s1, s2, 256, 1. / dt)
ax2.set_ylabel('CSD (db)')
plt.show()

总结

以上就是本文关于Python+matplotlib实现计算两个信号的交叉谱密度实例的全部内容,希望对大家有所帮助。感兴趣的朋友可以继续参阅本站其他相关专题,如有不足之处,欢迎留言指出。感谢朋友们对本站的支持!

Python 相关文章推荐
python通过定义一个类实例作为ftp回调方法
May 04 Python
Python OS模块常用函数说明
May 23 Python
Python正则表达式常用函数总结
Jun 24 Python
Python实现excel转sqlite的方法
Jul 17 Python
python判断字符串或者集合是否为空的实例
Jan 23 Python
python下载微信公众号相关文章
Feb 26 Python
解决python3中的requests解析中文页面出现乱码问题
Apr 19 Python
flask应用部署到服务器的方法
Jul 12 Python
Django多数据库的实现过程详解
Aug 01 Python
详解pycharm2020.1.1专业版安装指南(推荐)
Aug 07 Python
详解python对象之间的交互
Sep 29 Python
Python中的pprint模块
Nov 27 Python
python matplotlib 注释文本箭头简单代码示例
Jan 08 #Python
Python自定义简单图轴简单实例
Jan 08 #Python
[原创]python爬虫(入门教程、视频教程)
Jan 08 #Python
小米5s微信跳一跳小程序python源码
Jan 08 #Python
Python实现判断字符串中包含某个字符的判断函数示例
Jan 08 #Python
Python实现的字典值比较功能示例
Jan 08 #Python
python基础之包的导入和__init__.py的介绍
Jan 08 #Python
You might like
PHP中for循环语句的几种变型
2006/11/26 PHP
PHP Array交叉表实现代码
2010/08/05 PHP
php获取服务器端mac和客户端mac的地址支持WIN/LINUX
2014/05/15 PHP
PHP编程快速实现数组去重的方法详解
2017/07/22 PHP
PHP实现Redis单据锁以及防止并发重复写入
2018/04/10 PHP
PHP迭代器和迭代的实现与使用方法分析
2018/04/19 PHP
PHP设计模式之策略模式原理与用法实例分析
2019/04/04 PHP
php框架知识点的整理和补充
2021/03/01 PHP
capacityFixed 基于jquery的类似于新浪微博新消息提示的定位框
2011/05/24 Javascript
javaScript实现浮点数转十六进制字符
2013/10/29 Javascript
jQuery探测位置的提示弹窗(toolTip box)详细解析
2013/11/14 Javascript
js获取当前地址 JS获取当前URL的示例代码
2014/02/26 Javascript
js实现select跳转菜单新窗口效果代码分享(超简单)
2015/08/21 Javascript
微信JSSDK上传图片
2015/08/23 Javascript
JavaScript:Array类型全面解析
2016/05/19 Javascript
浅谈在fetch方法中添加header后遇到的预检请求问题
2017/08/31 Javascript
js 中rewrap-ajax.js插件实例代码
2017/10/20 Javascript
Vue 创建组件的两种方法小结(必看)
2018/02/23 Javascript
详解如何用babel转换es6的class语法
2018/04/03 Javascript
React事件处理的机制及原理
2018/12/03 Javascript
Vue+Express实现登录注销功能的实例代码
2019/05/05 Javascript
layui表格内放置图片,并点击放大的实例
2019/09/10 Javascript
Python查询阿里巴巴关键字排名的方法
2015/07/08 Python
django模型层(model)进行建表、查询与删除的基础教程
2017/11/21 Python
对python生成业务报表的实例详解
2019/02/03 Python
numpy.linalg.eig() 计算矩阵特征向量方式
2019/11/29 Python
浅谈Python 参数与变量
2020/06/20 Python
html5移动端价格输入键盘的实现
2019/09/16 HTML / CSS
html2canvas生成清晰的图片实现打印的示例代码
2019/09/30 HTML / CSS
西班牙拥有最佳品牌的动物商店:Animalear.com
2018/01/05 全球购物
旅游管理毕业生自荐信
2013/11/05 职场文书
宣传稿格式范文
2015/07/23 职场文书
2016年公司中秋节致辞
2015/11/26 职场文书
2016幼儿园新学期寄语
2015/12/03 职场文书
win server2012 r2服务器共享文件夹如何设置
2022/06/21 Servers
向Spring IOC 容器动态注册bean实现方式
2022/07/15 Java/Android