N=128;                         % 采樣點數
T=1;                                 % 采樣時間終點
t=linspace(0,T,N);                % 給出N個采樣時間ti(I=1:N)
x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t);  % 求各采樣點樣本值x
dt=t(2)-t(1);                   % 采樣周期
f=1/dt;                          % 采樣頻率(Hz)
X=fft(x);                        % 計算x的快速傅立葉變換X
F=X(1:N/2+1);                   % F(k)=X(k)(k=1:N/2+1)
f=f*(0:N/2)/N;                  % 使頻率軸f從零開始
plot(f,abs(F),'-*')            % 繪制振幅-頻率圖
xlabel('Frequency');
ylabel('|F(k)|')