SOGO論壇
  登入   註冊   找回密碼
SOGO論壇 › 論壇 › ≡ 電腦數位 ≡ › 電子電機暨資訊工程館 › 數位訊號處理入門---傅立葉分析
返回列表
查看: 2762|回覆: 3
列印 上一主題 下一主題

[心得分享] 數位訊號處理入門---傅立葉分析 [複製連結]

alphi

Rank: 5Rank: 5

數位軟體勳章

狀態︰ 離線
跳轉到指定樓層
1樓
發表於 2012-7-8 14:12:42 |只看該作者 |倒序瀏覽 | x 1
本帖最後由 alphi 於 2012-7-8 14:38 編輯

一、前言:
法國數學家傅立葉(Fourier)在1830年提出的傅立葉級數理論(Fourier Series),大意是講任意週期函數可以用正弦波(Sin)與餘弦波(Cos)級數和表示。
如今已經成為頻譜解析(Spectrum Analysis)的基礎理論在通訊與機率、統計、物理許多領域都會使用到該理論。那有何謂傅立葉轉換(Fourier Transformation)就是具 有週期性不規則波形(例如ECG:心電圖)利用傅立葉級數的原理(任意週期性波形是由該週期的基本頻率與整數倍頻率的高次正弦波與餘閒波所組成)可以求得波形的特徵 與資訊。
二、傅立葉級數
傅立葉級數(Fourier Series)就是[有一個隨時間而變化的訊號x(t)具有週期T之連續級數時，x(t)為直流成分與週期T的倒數1/T(Hz)頻率成分與其整數倍的頻率成分n/T(Hz)所組成]。

x(t)傅立葉級數展開=(直流成分)+(係數1*基本波)+(係數2*2次高諧波)+(係數3*3次高諧波)+....(係數n*n次高諧波)，其公式如下
  

  直流成分:表示輸入訊號的平均值
  高頻率成分:表示各高次諧波所占比率，由這些係數可以了解輸入波形的特徵。因此可以供各種波形分析之用。

三、傅立葉分析
  對被取樣後有限個數訊號做傅立葉轉換稱為離散傅立葉轉換(DFT:Discrete Fourier Transformation)。傅立葉轉換是求傅立葉級數之係數所做的處裡，針對傅立葉轉
  換要處裡的訊號可以用下列複數(Complex Number)表示
  
  
  
  因此做DFT(離散傅立葉轉換)之後轉換結果也為複數，如以下表示
  

四、程式碼
  由於此程式碼強調可讀性而非效能
  1.   #include <stdio.h>
  2.   #include <stdlib.h>
  3.   #include <math.h>
  4.   #define PI 3.14159265358
  5.   
  6.   typedef struct Complex
  7.   {
  8.      double r;
  9.      double i;
  10.   }Complex

  12.   void DFT(Complex *x,unsigned int N);
  13.   int main(int argc,char **argv)
  14.   {
  15.   }
  16.   void DFT(double *x,unsigned int N)
  17.   {
  18.       Complex X;

  20.       unsigned int n,k;
  21.       
  22.       for(k=0;k<N;k++)
  23.       {
  24.          x[k].i=0;
  25.       }
  26.       for(n=0;n<N;n++)
  27.       
  28.            X.r=0;
  29.            X.i=0;

  31.            for(k=0;k<N;k++)
  32.            {
  33.                X.r+=x[k].r*cos(2.0*PI*n*k/N)+x[k].i*sin(2.0*PI*n*k/N);
  34.                X.i+=x[k].r*cos(2.0*PI*n*k/N)-x[k].i*sin(2.0*PI*n*k/N);
  35.            }
  36.            printf("%f %f",X.r,X.i);
  37.       }
  38.   }
  39.   
複製代碼
頻譜分析
   頻譜分析(Spectrum Analysis)就是將訊號從時域轉換成頻域之後要分析的資料，如著眼在振幅就成為振幅頻譜(Amplitude Specturm),目前常見的三種頻譜
   為
   (a)功率頻譜(Power Spectrum)
      
   (b)振幅頻譜(Amplitude Spectrum)
        
   (c)相位頻譜(Phase Spectrum)

   功率頻譜的函式
  1.    void Power(Complex *x,const unsigned int N)
  2.    {
  3.             Complex X;
  4.             unsigned int k;
  5.             for(k=0;k<N/2;k++)
  6.             {
  7.                  X.r=pow(x[k].r,2.0);
  8.                  X.i=pow(x[k].i,2.0);

  10.                 printf("%f %f",X.r,X.i);         
  11.             }
  12.    }
  13.    
複製代碼
振幅頻譜的函式
  1.    void Amplitude(Complex *x,const unsigned int N)
  2.    {
  3.             Complex X;
  4.             unsigned int k;
  5.             for(k=0;k<N/2;k++)
  6.             {
  7.                  X.r=sqrt(pow(x[k].r,2.0));
  8.                  X.i=sqrt(pow(x[k].i,2.0));

  10.                 printf("%f %f",X.r,X.i);         
  11.             }
  12.    }
  13.    
複製代碼
相位頻譜的函式
  1.    void Amplitude(Complex *x,const unsigned int N)
  2.    {
  3.             Complex X;
  4.             unsigned int k;
  5.             for(k=0;k<N/2;k++)
  6.             {
  7.                  X.r=atan((x[k].r/x[k].i));
  8.                  X.i=atan(x[k].r/x[k].i));

  10.                 printf("%f %f",X.r,X.i);         
  11.             }
  12.    }
  13.    
複製代碼
已有 1 人評分威望 收起 理由
紅塵孤鳥 + 2 您發表的文章內容豐富，無私分享造福眾人，.

總評分: 威望 + 2   查看全部評分

喜歡嗎？分享這篇文章給親朋好友︰
               感謝作者     

失敗只有一種:那就是半途而廢

檢舉

mm117777

Rank: 11Rank: 11Rank: 11Rank: 11

熱心參予論壇活動及用心回覆主題勳章 數位硬體勳章

狀態︰ 離線
2樓
發表於 2012-7-8 14:22:48 |只看該作者
支持加1希望能多看到您的文章感恩合十

檢舉

didolin

Rank: 2

狀態︰ 離線
3樓
發表於 2012-9-9 01:55:24 |只看該作者
接下來就要玩調變了
傅立葉轉換真的非常重要
尤其是在通訊原理或訊號處理中
個人認為傅立葉是19世紀最偉大的數學家之一

檢舉

q2399

Rank: 3Rank: 3

狀態︰ 離線
4樓
發表於 2014-10-26 17:11:37 |只看該作者
感謝大大的分析加上程式，提供驗證的方法

檢舉

返回列表
請注意︰利用多帳號發表自問自答的業配文或置入性行銷廣告者，將直接禁訪或刪除帳號及全部文章！
高級模式
B Color Image Link Quote Code Smilies
您需要登錄後才可以回覆 登入 | 註冊


本論壇為非營利自由討論平台，所有個人言論不代表本站立場。文章內容如有涉及侵權，請通知管理人員，將立即刪除相關文章資料。侵權申訴或移除要求：abuse@oursogo.com

GMT+8, 2024-11-30 12:49

品味生活|休閒天地|興趣嗜好|學術藝文|女性頻道|電腦數位|遊戲天堂|影視音樂|宗教綜合|綜合圖片|綜合影片|成人園地
© 2004-2024 SOGO論壇 OURSOGO.COM
回頂部