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中replace方法实例分析
Aug 20 Python
Python使用chardet判断字符编码
May 09 Python
python获取指定路径下所有指定后缀文件的方法
May 26 Python
python实现的DES加密算法和3DES加密算法实例
Jun 03 Python
python中偏函数partial用法实例分析
Jul 08 Python
视觉直观感受若干常用排序算法
Apr 13 Python
Python安装lz4-0.10.1遇到的坑
May 20 Python
python全局变量引用与修改过程解析
Jan 07 Python
什么是Python中的匿名函数
Jun 02 Python
使用bandit对目标python代码进行安全函数扫描的案例分析
Jan 27 Python
Django中session进行权限管理的使用
Jul 09 Python
Pandas 稀疏数据结构的实现
Jul 25 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和ACCESS写聊天室(七)
2006/10/09 PHP
PHP下10件你也许并不了解的事情
2008/09/11 PHP
fetchAll()与mysql_fetch_array()的区别详解
2013/06/05 PHP
微信API接口大全
2015/04/15 PHP
ThinkPHP中order()使用方法详解
2016/04/19 PHP
Yii2实现ajax上传图片插件用法
2016/04/28 PHP
php支持断点续传、分块下载的类
2016/05/02 PHP
Smarty模板引擎缓存机制详解
2016/05/23 PHP
PHPStrom 新建FTP项目以及在线操作教程
2016/10/16 PHP
Laravel 错误提示本地化的实现
2019/10/22 PHP
Google韩国首页图标动画效果
2007/08/26 Javascript
用Javascript获取页面元素的具体位置
2013/12/09 Javascript
jQuery中position()方法用法实例
2015/01/16 Javascript
Javascript函数的参数
2015/07/16 Javascript
全面解析Bootstrap图片轮播效果
2015/12/03 Javascript
jQuery实现的自定义滚动条实例详解
2016/09/20 Javascript
jquery实时获取时间的简单实例
2017/01/26 Javascript
JavaScript实现替换字符串中最后一个字符的方法
2017/03/07 Javascript
Jquery中attr与prop的区别详解
2017/05/27 jQuery
Angularjs自定义指令实现分页插件(DEMO)
2017/09/16 Javascript
使用Anaconda3建立虚拟独立的python2.7环境方法
2018/06/11 Python
tensorflow 输出权重到csv或txt的实例
2018/06/14 Python
python框架中flask知识点总结
2018/08/17 Python
使用PM2+nginx部署python项目的方法示例
2018/11/07 Python
python实现抽奖小程序
2020/04/15 Python
python使用Pandas库提升项目的运行速度过程详解
2019/07/12 Python
Pyqt5自适应布局实例
2019/12/13 Python
基于Python生成个性二维码过程详解
2020/03/05 Python
详解Html5 Canvas画线有毛边解决方法
2018/03/01 HTML / CSS
台湾网友喜爱的综合型网路购物商城:Yahoo! 奇摩购物中心
2018/03/10 全球购物
美国婴儿服装购物网站:Gerber Childrenswear
2020/05/06 全球购物
一道SQL存储过程面试题
2016/10/07 面试题
会计助理的岗位职责
2013/11/29 职场文书
医院护士的求职信
2014/01/03 职场文书
党员年终个人总结
2015/02/14 职场文书
Python 数据可视化神器Pyecharts绘制图像练习
2022/02/28 Python