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的迭代器与生成器实例详解
Jul 16 Python
Python返回真假值(True or False)小技巧
Apr 10 Python
python图像处理之反色实现方法
May 30 Python
简单谈谈Python中的反转字符串问题
Oct 24 Python
python3制作捧腹网段子页爬虫
Feb 12 Python
如何用itertools解决无序排列组合的问题
May 18 Python
Python实现数据库并行读取和写入实例
Jun 09 Python
flask框架使用orm连接数据库的方法示例
Jul 16 Python
Python获取Redis所有Key以及内容的方法
Feb 19 Python
在tensorflow中实现去除不足一个batch的数据
Jan 20 Python
Python3和PyCharm安装与环境配置【图文教程】
Feb 14 Python
Python关键字及可变参数*args,**kw原理解析
Apr 04 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
十大催泪虐心动漫电影,有几部你还没看
2020/03/04 日漫
php笔记之:有规律大文件的读取与写入的分析
2013/04/26 PHP
UTF-8正则表达式如何匹配汉字
2015/08/03 PHP
PHP实现补齐关闭的HTML标签
2016/03/22 PHP
PHP下的浮点运算不准的解决方法
2016/10/27 PHP
jQuery 表单验证扩展(四)
2010/10/20 Javascript
jQuery平滑旋转幻灯片特效代码分享
2015/09/07 Javascript
JavaScript 身份证号有效验证详解及实例代码
2016/10/20 Javascript
微信小程序开发入门基础教程
2017/04/19 Javascript
Spring Boot/VUE中路由传递参数的实现代码
2018/03/02 Javascript
iview同时验证多个表单问题总结
2018/09/29 Javascript
在 Vue-CLI 中引入 simple-mock实现简易的 API Mock 接口数据模拟
2018/11/28 Javascript
vue.js仿hover效果的实现方法示例
2019/01/28 Javascript
利用d3.js实现蜂巢图表带动画效果
2019/09/03 Javascript
基于JS正则表达式实现模板数据动态渲染(实现思路详解)
2020/03/07 Javascript
Node.js API详解之 module模块用法实例分析
2020/05/13 Javascript
vue实现图书管理系统
2020/12/29 Vue.js
手写Vue源码之数据劫持示例详解
2021/01/04 Vue.js
python简单的函数定义和用法实例
2015/05/07 Python
python进行TCP端口扫描的实现
2018/12/21 Python
元组列表字典(莫烦python基础)
2019/04/03 Python
详解用python实现基本的学生管理系统(文件存储版)(python3)
2019/04/25 Python
Python提取转移文件夹内所有.jpg文件并查看每一帧的方法
2019/06/27 Python
深入解析神经网络从原理到实现
2019/07/26 Python
python实现音乐播放和下载小程序功能
2020/04/26 Python
scrapy与selenium结合爬取数据(爬取动态网站)的示例代码
2020/09/28 Python
佳能加拿大网上商店:Canon eStore Canada
2018/04/04 全球购物
英国123鲜花网站:123 Flowers
2019/07/07 全球购物
教师实习的自我鉴定
2013/10/26 职场文书
校运会入场式解说词
2014/02/10 职场文书
合作经营协议书
2014/04/17 职场文书
竞聘演讲稿精彩开头和结尾
2014/05/14 职场文书
《草船借箭》教学反思
2016/02/23 职场文书
初一语文教学反思
2016/03/03 职场文书
Python Pandas模块实现数据的统计分析的方法
2021/06/24 Python
未发现nvidia显卡怎么办?Win11系统中未检测到nvidia显卡解决教程
2022/04/08 数码科技