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程序相关推荐

  1. matlab fft反变换,Matlab实现FFT变换

    Matlab实现FFT变换 Matlab实现FFT变换(单边谱及双边谱)­ ­ 以前对于Fourier Transform从来没有细究,不管在LabVIEW还是Matlab里都有现成的FFT(快速Fo ...

  2. 对固定长度序列进行不同长度下的FFT变换的区别

    在matlab中函数fft负责对信号序列进行离散傅里叶变换,Y = fft(X)表示为对输入序列X进行离散傅里叶变化,Y为变化后的频域序列. Y = fft(X)中未对"采用几点fft&qu ...

  3. 深入浅出解释FFT(六)——深入理解fft变换

    (如需交流,请关注公众号:神马观止) FFT(FastFourier Transform,快速傅立叶变换)是离散傅立叶变换的快速算法,也是我们在数字信号处理技术中经常会提到的一个概念.在大学的理工科课 ...

  4. 模数转换实验中断方式c语言,DSP实验报告--模拟信号的AD+FFT变换

    DSP实验报告--模拟信号的AD+FFT变换,完整的实验报告,内含源代码(C语言). 一. 实验开发环境 1.通用 PC机一台,安装 Windows2000 或 WindowsXP 操作系统且已安 装 ...

  5. python读取wav文件并进行FFT变换

    一.wav格式文件 WAV为微软公司(Microsoft)开发的一种声音文件格式,它符合RIFF(Resource Interchange File Format)文件规范,用于保存Windows平台 ...

  6. DSP28335通过FFT变换实现高频滤波

    这次项目中,由于AD采集到的高频开关噪声大,对数据处理存在极大干扰,经过示波器测量,发现噪声频率主要集中在80kHz,寻求老师帮助,老师推荐使用FFT滤波,有同学使用的是FIR滤波,通过Matlab中 ...

  7. FFT变换的C语言实现

    FFT变换的C语言实现 1.安装FFTW3库文件 FFTW是用来计算一维或者多维的离散傅里叶变换,输入可以为实数序列也可以为复数序列的C语言的子函数库,FFTW是免费软件,是作为fft函数库的各种应用 ...

  8. matlab中的fft变换后的频率的缩放系数问题。

    本文的目的在于说明fft变换之后的幅值对应的真实频率到底是多少. 为了方便大家理解fft的基础知识,在下面附上<小波分析与傅里叶变换基础>的pdf,dft,fft相关在第三章. 链接:ht ...

  9. STM32F407 DSP+FPU进行FFT变换(2)

    STM32F407 DSP+FPU进行FFT变换 接着上一篇继续,要用FFT运算的话,F4有FPU和DSP库,可以很方便让我们去对数据进行傅氏变换.首先得配置好DSP库和FPU. 配置DPS库和FPU ...

  10. fft 估计载波频率程序_基于FFT载频估计和Costas环的非协作通信载波同步系统的制作方法...

    基于FFT载频估计和Costas环的非协作通信载波同步系统的制作方法 [技术领域] [0001] 本发明涉及本发明属于数字通信技术领域,特别是涉及一种通信载波同步系统. [背景技术] [0002] 在 ...

最新文章

  1. javascript:正则表达式、一个表单验证的例子
  2. Debian8 远程登录Permission Denied,please try again
  3. antd 动态添加表单_react Ant Design 动态生成表单,并带验证
  4. mysql补充(1)校对集utf8_unicode_ci与utf8_general_ci
  5. 迭代器 生成器 装饰器 匿名函数
  6. Acwing 135 最大子序和
  7. 【机器学习实战之一】:C++实现K-近邻算法KNN
  8. 超硬核的 OCR 开发利器推荐:从场景案例到实操演示,快速掌握OCR模型训练
  9. mysql client 升级_解决consider upgrading MySQL client问题
  10. Layer数据表格监听排序切换
  11. 专访 | 阿里前辈李智慧:如何培养架构思维,我把20年的架构经讲给你听
  12. c语言追赶法求方程组的解,MATLAB-追赶法求解三对角方程组的算法原理例题与程序...
  13. 时域信号处理——时域采样(抽样)
  14. 客运售票员_客运室优秀售票员事迹材料
  15. 使用turtle画随机颜色的圆圈
  16. 【Java】初识泛型(带你从初学者角度切入,通俗易懂,速进)
  17. vue3使用vuedraggable实现拖拽
  18. UVA 1603 Square Destroyer
  19. 财务管理软件QuickBooks 2020 for Mac
  20. ios 穿山甲广告联盟集成_今日头条发布穿山甲联盟 如何实现广告联盟模式进阶?...

热门文章

  1. matlab 读取shp面文件,在matlab中将处理结果输出为shp文件
  2. c语言使用的编程环境是,啊哈C语言编程-第3课-C语言编程环境
  3. Server 安装 caffee
  4. 新H5手机微商城运营源码 简约轻量版
  5. C++中的万能头文件
  6. 安装Selenium+Firefox+Firepath+Firebug
  7. 荒野行动android模拟,荒野行动用模拟器玩教程 荒野行动模拟器不支持机型解决方法...
  8. 进制转化(北理乐学编程题目)
  9. 软件项目管理 1.2.PMBOK与软件项目管理知识体系
  10. 【Unity游戏开发】Unity基础(一)面板介绍