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 相关文章推荐
以windows service方式运行Python程序的方法
Jun 03 Python
举例讲解Django中数据模型访问外键值的方法
Jul 21 Python
wxPython的安装图文教程(Windows)
Dec 28 Python
tensorflow 加载部分变量的实例讲解
Jul 27 Python
浅析Python 实现一个自动化翻译和替换的工具
Apr 14 Python
详解python数据结构和算法
Apr 18 Python
解决pyecharts在jupyter notebook中使用报错问题
Apr 23 Python
python Tcp协议发送和接收信息的例子
Jul 22 Python
Python终端输出彩色字符方法详解
Feb 11 Python
解决TensorFlow程序无限制占用GPU的方法
Jun 30 Python
Python unittest装饰器实现原理及代码
Sep 08 Python
详解用 python-docx 创建浮动图片
Jan 24 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
解决phpmyadmin中文乱码问题。。。
2007/01/18 PHP
深入PHP运行环境配置的详解
2013/06/04 PHP
php判断ip黑名单程序代码实例
2014/02/24 PHP
Session的工作机制详解和安全性问题(PHP实例讲解)
2014/04/10 PHP
php分页函数完整实例代码
2014/09/22 PHP
php对象在内存中的存在形式分析
2015/02/03 PHP
laravel框架中视图的基本使用方法分析
2019/11/23 PHP
点击广告后才能获得下载地址
2006/10/26 Javascript
二级域名转向类
2006/11/09 Javascript
json-lib出现There is a cycle in the hierarchy解决办法
2010/02/24 Javascript
js操作输入框中选择内容兼容IE及其他主流浏览器
2014/04/22 Javascript
Node.js开发之访问Redis数据库教程
2015/01/14 Javascript
Javascript闭包(Closure)详解
2015/05/05 Javascript
JavaScript使用RegExp进行正则匹配的方法
2015/07/11 Javascript
JS、jQuery中select的用法详解
2016/04/21 Javascript
AngularJS中directive指令使用之事件绑定与指令交互用法示例
2016/11/22 Javascript
Parcel.js + Vue 2.x 极速零配置打包体验教程
2017/12/24 Javascript
Vue路由钩子之afterEach beforeEach的区别详解
2018/07/15 Javascript
React 组件中的 bind(this)示例代码
2018/09/16 Javascript
JS数组进阶示例【数组的几种函数用法】
2020/01/16 Javascript
ant-design-vue中的select选择器,对输入值的进行筛选操作
2020/10/24 Javascript
Python字符串替换实例分析
2015/05/11 Python
Python统计日志中每个IP出现次数的方法
2015/07/06 Python
python使用Plotly绘图工具绘制散点图、线形图
2019/04/02 Python
Python中的字符串切片(截取字符串)的详解
2019/05/15 Python
python GUI库图形界面开发之PyQt5动态(可拖动控件大小)布局控件QSplitter详细使用方法与实例
2020/03/06 Python
python实现udp传输图片功能
2020/03/20 Python
Python实现列表索引批量删除的5种方法
2020/11/16 Python
python 实现aes256加密
2020/11/27 Python
加拿大时尚床上用品零售商:QE Home | Quilts Etc
2018/01/22 全球购物
自我评价怎么写好呢?
2013/12/05 职场文书
家长给幼儿园的表扬信
2014/01/09 职场文书
争先创优心得体会
2014/09/12 职场文书
处级干部反四风个人对照检查材料思想汇报
2014/09/27 职场文书
保险公司增员口号
2015/12/25 职场文书
高端收音机+蓝牙音箱,JBL TUNER FM带收音蓝牙音箱评测
2021/04/24 无线电