Python+matplotlib实现计算两个信号的交叉谱密度实例
作者：mengwei 发布时间：2023-05-12 23:29:55
标签：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()
来源：https://matplotlib.org/index.html
0
投稿猜你喜欢
在学习redis过程中提到一个缓存击穿的问题， 书中参考的解决方案之一是使用布隆过滤器， 那么就有必要来了解一下什么是布隆过滤器。在参考了许- 简介在生产环境上，为了避免数据的丢失，通常情况下都会定时的对数据库进行备份。而Linux的crontab指令则可以帮助我们实现对数据库定时进
- pandas中对DataFrame筛选数据的方法有很多的，以后会后续进行补充，这里只整理遇到错误的情况。1.使用布尔型DataFrame对数
- 默认情况下，SQL Server 保存 7 个 ErrorLog 文件，名为： ErrorLog ErrorLog.1 ErrorLog.2
- hashlibhashlib主要提供字符加密功能，将md5和sha模块整合到了一起，支持md5,sha1, sha224, sha256,
- 本文使用案例是基于 python2.7 实现以下内容均为个人使用 peewee 的经验和遇到的坑，不会涉及过多的基本操作。所以，没有使用过
- 你是否苦恼于网上无法下载的“小说在线阅读”内容？或是某些文章的内容让你很有收藏的冲动，却找不到一个下载的链接？是不是有种自己写个程序把全部搞
- Prometheus是什么Prometheus是一套开源监控系统和告警为一体，由go语言(golang)开发，是监控+报警+时间序列数据库的
- 很早之前就在PJ的blog上看到可以用VS2005调试ASP程序，但是没有写出具体的步骤，后来一次偶尔也让我找到了方法，但是一直没把它写出来
- 因为工作原因，需要定期清理某个文件夹下面创建时间超过1年的所有文件，所以今天集中学习了一下Python对于本地文件及文件夹的操作。网上 这篇
- 1、要点击链接，然后点击里面的上传tab,不熟悉的人可能找不到这个上传功能 2、插入的就是1个链接，我希望插入链接的同时插入1个图片代表文件
- 我就废话不多说了，大家还是直接看代码吧~import tensorflow as tfn1 = tf.constant(2)n2 = tf.
- 本文实例讲述了关于php中SimpleXML 函数的用法,此函数是允许您把 XML 转换为对象,分享给大家供大家参考。具体分析如下：Simp
- 第一个问题是重命名数据库问题：在企业管理器中是无法直接对数据库重命名的，只能在查询分析器中操作create proc killspid (@
- 【人工智能项目】Python Flask搭建yolov3目标检测系统后端代码from flask import Flask, request
- 简介滚动条小部件用于向下滚顶其他小部件的内容,如列表框,文本和画布,但是,我们也可以为Entry小部件创建水平滚动条,常常被用于实现文本,画
- 1. 规范简介本规范主要规定ASP源程序在书写过程中所应遵循的规则及注意事项。编写该规范的目的是使项目开发人员的源代码书写习惯保持一致。这样
- 一、对列表（list）进行排序推荐的排序方式是使用内建的sort()方法，速度最快而且属于稳定排序>>> a = [1,9
- 1.首先在pycharm上使用pip安装pip install html-table pip install jira2.初始化发件人邮箱，
- 在使用python做爬虫的过程中，有些页面的的部分数据是通过js异步加载的，js调用接口的请求中有时还带有些加密的参数很难破解无法使用req
