function [Xk, w] = DFT(x, M, str)%離散傅里葉變換及x(n),FT[x(n)]的繪圖函數    n = 0 : M - 1;    %單位圓離散譜點數    N = 200;    %畫2個周期    k = -N : N - 1;    %離散頻率橫坐標    w = (2 * pi / N) * k;    %x(e^jw)=∑x(n)e^(jwn)，使用矩陣乘法實現    Xk = x * exp( -j * (2 * pi / N) ) .^ ( n' * k );    magX = abs(Xk);        subplot(2, 2, 1);    t = max(x);    n = 0 : M - 1;    stem(n, x, '.');    xlabel('n');    ylabel('x(n)');    title('時域信號波形');    text( (0.3 * M), (0.8 * t), str );        subplot(2, 2, 3)    plot(w/pi, magX, 'b');    xlabel('ω/π');    ylabel('|X(e^j^w)|');    title('信號的幅度譜');%     %使用FFT算法%     subplot(2, 2, 4)%     Xk = fft(x, N);%     Xk = [Xk, Xk];%補足2個周期%     magX = abs(Xk);%     stem(w/pi, magX, 'r');%     xlabel('ω/π');%     ylabel('|X(e^j^w)|');%     title('信號的FFT幅度譜');