其實是早應(yīng)該想到的，今天寫程序的時候想到了。然后，跟同學說起來，說，原來國外的教材很多都是矩陣的形式來表示離散傅里葉變換的。但是國內(nèi)的教材沒有這么表達。

一個是，自己看的東西還是少了，一個是，matlab本來應(yīng)該是混飯吃的東西，用的也不熟，理解還欠很多，另外一個是，對于線性代數(shù)，矩陣，的理解還真是淺薄。這么有用的東西，自己還是應(yīng)該抽空看看的。

對于x[n]的離散時間傅里葉變換

離散時間傅里葉變換

對于N點序列,離散傅里葉變換

離散傅里葉變換

可以看到，先是乘法，然后是求和。很容易想到，對于矩陣的乘法，就是先一行乘以一列，再把各個值相加，為一項的值，再去做下一步。

到最后需要得到的是一個1xN(或1xW)的矩陣，所以矩陣乘法中，左邊的應(yīng)該是一個一行的矩陣，為x[n]，exp中，需要得到一個NxN(或者NxW)的矩陣，由，nk得到，形式為n’*w。

clc;

fs=100; %采樣頻率

N=128; %數(shù)據(jù)點數(shù)

n=0:N-1;

t=n/fs; %時間序列

x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); %信號

%此處x[n]，就是需要的輸入的序列

%注意，matlab標號從一開始，所以需要有n+1

w=[0:999]

%w為一千個點，在這里相當于，w是連續(xù)的，若w的數(shù)量等于N，則就是N點的DFT。

y=x(n+1)*exp(-j*2*pi/N*n'*w)

%y[w]即傅里葉變換系數(shù)

總結(jié)

以上是生活随笔為你收集整理的matlab dtft的函数,DTFT的Matlab矩阵计算的理解的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。