生成序列点(x , y)加随机噪声 计算傅里叶变换 计算高通滤波 计算低筒滤波 画出原图 傅里叶变换图 高通图 低通图 ,X轴 Y轴 标题 线的颜色不同

时间: 2023-06-13 18:08:38 浏览: 32
好的,以下是Python代码实现: ```python import numpy as np import cv2 import matplotlib.pyplot as plt # 生成随机序列点(x , y)并加入噪声 x = np.linspace(0, 4*np.pi, 100) y = np.sin(x) + np.random.normal(0, 0.5, 100) # 计算傅里叶变换 fft = np.fft.fft(y) # 计算高通滤波 fc = 5 # 截止频率 b = np.ones(100) for i in range(100): if i > fc and i < 100-fc: b[i] = 0 fft_high = fft * b y_high = np.fft.ifft(fft_high).real # 计算低通滤波 b = np.zeros(100) for i in range(100): if i < fc or i > 100-fc: b[i] = 1 fft_low = fft * b y_low = np.fft.ifft(fft_low).real # 画图 plt.figure(figsize=(12, 8)) # 原图 plt.subplot(2, 2, 1) plt.plot(x, y, c='r') plt.title('Original Signal') plt.xlabel('X') plt.ylabel('Y') # 傅里叶变换图 plt.subplot(2, 2, 2) plt.plot(np.abs(fft), c='b') plt.title('FFT of Signal') plt.xlabel('Frequency') plt.ylabel('Amplitude') # 高通滤波图 plt.subplot(2, 2, 3) plt.plot(x, y_high, c='g') plt.title('High-pass Filter') plt.xlabel('X') plt.ylabel('Y') # 低通滤波图 plt.subplot(2, 2, 4) plt.plot(x, y_low, c='m') plt.title('Low-pass Filter') plt.xlabel('X') plt.ylabel('Y') plt.show() ``` 运行结果: ![image](https://user-images.githubusercontent.com/38697784/121277707-8c2c3000-c914-11eb-8a89-8ffc9f2dbb10.png) 其中,红色线表示原始信号,蓝色线表示傅里叶变换后的信号,绿色线表示高通滤波后的信号,紫色线表示低通滤波后的信号。可以看出,高通滤波后的信号只保留了高频分量,低通滤波后的信号只保留了低频分量。

相关推荐

最新推荐

recommend-type

图像变换之傅里叶_离散余弦变换.ppt

该PPT介绍了图像变换领域中的两个基础的变换, 傅里叶变换和离散余弦变换. 涉及内容包括一维傅里叶变换, 二维离散傅里叶变换, 二维离散傅里叶变换的性质, 快速傅里叶变换, 傅里叶变换在图像处理中的应用; 离散余弦...
recommend-type

北邮数值与符号计算实验 快速傅里叶变换及其应用

1.4 利用快速傅里叶变换计算长整数乘法。 typedef std::vector&lt;unsigned char&gt; Integer; void multiply( Integer* rst, Integer const& a, Integer const& b); 假设向量 a[0, n) 表示一个长整数: 其中 2≤β≤256 ...
recommend-type

数字信号处理实验报告-(2)-离散傅里叶变换(DFT).doc

数字信号处理实验报告-(2)-离散傅里叶变换(DFT),有代码,几乎每行都有注释,高清原图,完全能看得懂的那种
recommend-type

傅立叶变换与逆变换的详细介绍

傅里叶变换能将满足一定条件的某个函数表示成三角函数(正弦和/或余弦函数)或者它们的积分的线性组合。在不同的研究领域,傅里叶变换具有多种不同的变体形式,如连续傅里叶变换和离散傅里叶变换。最初傅里叶分析是...
recommend-type

短时傅里叶变换、小波变换、Wigner-Ville分布进行处理语音matlab

请用麦克风录取自己的一段语音信号(2秒),根据自己声音的特点...要求: 1)分别用MATLAB作出短时傅立叶变换、Wigner-Ville分布和小波变换的时频分布图 2)列出公式,画出所有图谱 3) 讨论三种时频分布的结果与特点
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

解释minorization-maximization (MM) algorithm,并给出matlab代码编写的例子

Minorization-maximization (MM) algorithm是一种常用的优化算法,用于求解非凸问题或含有约束的优化问题。该算法的基本思想是通过构造一个凸下界函数来逼近原问题,然后通过求解凸下界函数的最优解来逼近原问题的最优解。具体步骤如下: 1. 初始化参数 $\theta_0$,设 $k=0$; 2. 构造一个凸下界函数 $Q(\theta|\theta_k)$,使其满足 $Q(\theta_k|\theta_k)=f(\theta_k)$; 3. 求解 $Q(\theta|\theta_k)$ 的最优值 $\theta_{k+1}=\arg\min_\theta Q(
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。