目錄

• • DFT的計算公式
  • 代碼
  • gnuplot畫圖和matlab仿真結果

DFT的計算公式

公式評論說有錯

$X(k)={\sum }_{n=0}^{N?1}x(n)exp(?j2\Pi /N?n?k)$

根據公式以及歐拉定理

$e^{ix}=cos(x)+i?sin(x)$

所以有

$X(k)={\sum }_{n=0}^{N?1}x(n)?(cos(2?\Pi /N?n?k)+i?sin(2?\Pi /N?n?k))$

有

$z=a+i?b=\sqrt{a^2+b^2}$

代碼

#include <stdio.h>//標準輸入輸出頭文 #include<math.h>//聲明了常用的一些數學運算 #define pi 3.1415 #define N 32//32個點typedef struct {double real,imag; } complex;//復數 結構 complex dft_out[100];//單個點計算k complex dft_one[100];//單個點計算n double amp[N];//DFT后的結果 int main() {int n,k;double xn;for(k=0; k<N; k++)//k循環{for(n=0; n<N; n++)//n循環{xn=cos(n*pi/6);//要DFT的信號dft_one[n].real=xn*cos(2*pi/N*n*k);//實部信號dft_one[n].imag=xn*sin(2*pi/N*n*k);//虛部信號dft_out[k].real+=dft_one[n].real;dft_out[k].imag+=dft_one[n].imag;//DFT后的實部，虛部相加}amp[k]=sqrt(dft_out[k].real*dft_out[k].real+dft_out[k].imag*dft_out[k].imag);//歐拉公式printf("%d %f\n",k,amp[k]);} }

gnuplot畫圖和matlab仿真結果

C：

matlab：

總結

以上是生活随笔為你收集整理的DFT C 语言的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。