其实是早应该想到的，今天写程序的时候想到了。然后，跟同学说起来，说，原来国外的教材很多都是矩阵的形式来表示离散傅里叶变换的。但是国内的教材没有这么表达。
一个是，自己看的东西还是少了，一个是，matlab本来应该是混饭吃的东西，用的也不熟，理解还欠很多，另外一个是，对于线性代数，矩阵，的理解还真是浅薄。这么有用的东西，自己还是应该抽空看看的。
对于x[n]的离散时间傅里叶变换
离散时间傅里叶变换
对于N点序列,离散傅里叶变换
离散傅里叶变换
可以看到，先是乘法，然后是求和。很容易想到，对于矩阵的乘法，就是先一行乘以一列，再把各个值相加，为一项的值，再去做下一步。
到最后需要得到的是一个1xN（或1xW）的矩阵，所以矩阵乘法中，左边的应该是一个一行的矩阵，为x[n]，exp中，需要得到一个NxN（或者NxW）的矩阵，由，nk得到，形式为n’*w。
clc;
fs=100; %采样频率
N=128; %数据点数
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是连续的，若w的数量等于N，则就是N点的DFT。
y=x(n+1)*exp(-j*2*pi/N*n'*w)
%y[w]即傅里叶变换系数
