python傅里叶变换例子

本质：

就是对信号进行基变换，只不过新的基是正弦/余弦波罢了

我有一个东西，长得特别难看，数学性质特别差，傅立叶变换就是用性质很好的三角函数一部分一部分的去逼近它，最后得出一个看起来很复杂实际上数学性质很好的东西。

代码：

  

   

    

     

    
    

     

      # -*- coding: utf-8 -*-
     
    
   

    

     

    
    

     

      import os
     
    
   

    

     

    
    

     

      import numpy as np
     
    
   

    

     

    
    

     

      import cv2
     
    
   

    

     

    
    

     
 
     
    
   

    

     

    
    

     

      import numpy as np
     
    
   

    

     

    
    

     

      from scipy.fftpack import fft, ifft
     
    
   

    

     

    
    

     

      import matplotlib.pyplot as plt
     
    
   

    

     

    
    

     

      from matplotlib.pylab import mpl
     
    
   

    

     

    
    

     

      from mpl_toolkits.mplot3d import Axes3D
     
    
   

    

     

    
    

     

      import time
     
    
   

    

     

    
    

     
 
     
    
   

    

     

    
    

     

      mpl.rcParams['font.sans-serif'] = ['SimHei']  # 显示中文
     
    
   

    

     

    
    

     

      mpl.rcParams['axes.unicode_minus'] = False  # 显示负号
     
    
   

    

     

    
    

     
 
     
    
   

    

     

    
    

     

      # 采样点选择1400个，因为设置的信号频率分量最高为600赫兹，根据采样定理知采样频率要大于信号频率2倍
     
    
   

    

     

    
    

     

      a = 0
     
    
   

    

     

    
    

     

      while a < 100:
     
    
   

    

     

    
    

     

       a += 1
     
    
   

    

     

    
    

     

       x = np.linspace(0, 1, 1024)
     
    
   

    

     

    
    

     
 
     
    
   

    

     

    
    

     
 # 设置需要采样的信号，频率分量有200，400和600
     
    
   

    

     

    
    

     

       i = np.random.randint(10, size=3)
     
    
   

    

     

    
    

     

       fr = np.random.randint(512, size=3)
     
    
   

    

     

    
    

     
 
     
    
   

    

     

    
    

     

       y = i[0] * np.sin(2 * np.pi * fr[0] * x) + i[1] * np.sin(2 * np.pi * fr[1] * x) &
     
    
  
 

文章来源: blog.csdn.net，作者：网奇，版权归原作者所有，如需转载，请联系作者。

原文链接：blog.csdn.net/jacke121/article/details/111881350