正弦波及FFT变换py程序
import numpy as np
from scipy.fftpack import fft,ifft
from matplotlib.pylab import plt# Fs = 1000 #采样频率
# T = 1/Fs #采样周期,只相邻两数据点的时间间隔
# L = 1500 #信号长度Fs = 5000000000 #采样频率
T = 1/Fs #采样周期,只相邻两数据点的时间间隔
L = 10000#信号长度nnn = 500000t = np.arange(L)*T #间隔 T * 个数
npp = np.piX = np.sin(2*np.pi *nnn*t)
# X = S + np.random.rand(L)plt.plot(t[:nnn], X[:nnn])
plt.xlabel("Time(s)")
plt.ylabel("Amplitude")
plt.title("Signol with random noise")
plt.show()Y = fft(X)
p2 = np.abs(Y) # 双侧频谱
p1 = p2[:int(L/2)]
f = np.arange(int(L/2))*Fs/L;
# plt.plot(f,2*p1/L)
plt.plot(f, p1)
plt.title('Single-Sided Amplitude Spectrum of X(t)')
plt.xlabel('f (Hz)')
plt.ylabel('|P1(f)|')
plt.show()
正弦波及FFT变换py程序相关推荐
- matlab fft反变换,Matlab实现FFT变换
Matlab实现FFT变换 Matlab实现FFT变换(单边谱及双边谱) 以前对于Fourier Transform从来没有细究,不管在LabVIEW还是Matlab里都有现成的FFT(快速Fo ...
- 对固定长度序列进行不同长度下的FFT变换的区别
在matlab中函数fft负责对信号序列进行离散傅里叶变换,Y = fft(X)表示为对输入序列X进行离散傅里叶变化,Y为变化后的频域序列. Y = fft(X)中未对"采用几点fft&qu ...
- 深入浅出解释FFT(六)——深入理解fft变换
(如需交流,请关注公众号:神马观止) FFT(FastFourier Transform,快速傅立叶变换)是离散傅立叶变换的快速算法,也是我们在数字信号处理技术中经常会提到的一个概念.在大学的理工科课 ...
- 模数转换实验中断方式c语言,DSP实验报告--模拟信号的AD+FFT变换
DSP实验报告--模拟信号的AD+FFT变换,完整的实验报告,内含源代码(C语言). 一. 实验开发环境 1.通用 PC机一台,安装 Windows2000 或 WindowsXP 操作系统且已安 装 ...
- python读取wav文件并进行FFT变换
一.wav格式文件 WAV为微软公司(Microsoft)开发的一种声音文件格式,它符合RIFF(Resource Interchange File Format)文件规范,用于保存Windows平台 ...
- DSP28335通过FFT变换实现高频滤波
这次项目中,由于AD采集到的高频开关噪声大,对数据处理存在极大干扰,经过示波器测量,发现噪声频率主要集中在80kHz,寻求老师帮助,老师推荐使用FFT滤波,有同学使用的是FIR滤波,通过Matlab中 ...
- FFT变换的C语言实现
FFT变换的C语言实现 1.安装FFTW3库文件 FFTW是用来计算一维或者多维的离散傅里叶变换,输入可以为实数序列也可以为复数序列的C语言的子函数库,FFTW是免费软件,是作为fft函数库的各种应用 ...
- matlab中的fft变换后的频率的缩放系数问题。
本文的目的在于说明fft变换之后的幅值对应的真实频率到底是多少. 为了方便大家理解fft的基础知识,在下面附上<小波分析与傅里叶变换基础>的pdf,dft,fft相关在第三章. 链接:ht ...
- STM32F407 DSP+FPU进行FFT变换(2)
STM32F407 DSP+FPU进行FFT变换 接着上一篇继续,要用FFT运算的话,F4有FPU和DSP库,可以很方便让我们去对数据进行傅氏变换.首先得配置好DSP库和FPU. 配置DPS库和FPU ...
- fft 估计载波频率程序_基于FFT载频估计和Costas环的非协作通信载波同步系统的制作方法...
基于FFT载频估计和Costas环的非协作通信载波同步系统的制作方法 [技术领域] [0001] 本发明涉及本发明属于数字通信技术领域,特别是涉及一种通信载波同步系统. [背景技术] [0002] 在 ...
最新文章
- javascript:正则表达式、一个表单验证的例子
- Debian8 远程登录Permission Denied,please try again
- antd 动态添加表单_react Ant Design 动态生成表单,并带验证
- mysql补充(1)校对集utf8_unicode_ci与utf8_general_ci
- 迭代器 生成器 装饰器 匿名函数
- Acwing 135 最大子序和
- 【机器学习实战之一】:C++实现K-近邻算法KNN
- 超硬核的 OCR 开发利器推荐:从场景案例到实操演示,快速掌握OCR模型训练
- mysql client 升级_解决consider upgrading MySQL client问题
- Layer数据表格监听排序切换
- 专访 | 阿里前辈李智慧:如何培养架构思维,我把20年的架构经讲给你听
- c语言追赶法求方程组的解,MATLAB-追赶法求解三对角方程组的算法原理例题与程序...
- 时域信号处理——时域采样(抽样)
- 客运售票员_客运室优秀售票员事迹材料
- 使用turtle画随机颜色的圆圈
- 【Java】初识泛型(带你从初学者角度切入,通俗易懂,速进)
- vue3使用vuedraggable实现拖拽
- UVA 1603 Square Destroyer
- 财务管理软件QuickBooks 2020 for Mac
- ios 穿山甲广告联盟集成_今日头条发布穿山甲联盟 如何实现广告联盟模式进阶?...
热门文章
- matlab 读取shp面文件,在matlab中将处理结果输出为shp文件
- c语言使用的编程环境是,啊哈C语言编程-第3课-C语言编程环境
- Server 安装 caffee
- 新H5手机微商城运营源码 简约轻量版
- C++中的万能头文件
- 安装Selenium+Firefox+Firepath+Firebug
- 荒野行动android模拟,荒野行动用模拟器玩教程 荒野行动模拟器不支持机型解决方法...
- 进制转化(北理乐学编程题目)
- 软件项目管理 1.2.PMBOK与软件项目管理知识体系
- 【Unity游戏开发】Unity基础(一)面板介绍