数字信号实验一
步骤 1：
X1=[1,1,1,1];X2=[1:4];X3=[4:-1:1];X4=[X2,X3];X5=[X3,X2];X1K8=fft(X1,8);X1K16=fft(X1,16);X4K8=fft(X4,8);X4K16=fft(X4,16);X5K8=fft(X5,8);X5K16=fft(X5,16);figure(1);subplot(2,3,1);stem(abs(X1K8));axis([0,8,0,1.2*max(abs(X1K8))]);title("x1 的 8 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,4);stem(abs(X1K16));axis([0,16,0,1.2*max(abs(X1K16))]);title("x1 的 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,2);stem(abs(X4K8));axis([0,8,0,1.2*max(abs(X4K8))]);title("x2 的 8 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,5);stem(abs(X4K16));axis([0,16,0,1.2*max(abs(X4K16))]);title("x2 的 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,3);stem(abs(X5K8));axis([0,8,0,1.2*max(abs(X5K8))]);title("x1 的 8 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,6);stem(abs(X5K16));axis([0,16,0,1.2*max(abs(X5K16))]);title("x1 的 8 点 FFT");xlabel("数字角频率 w");
ylabel("幅值 A");0 501234x1的 8点 FFT数 字 角 频 率 w幅值A0 1001234x1的 16点 FFT数 字 角 频 率 w幅值A0 505101520x2的 8点 FFT数 字 角 频 率 w幅值A0 1005101520x2的 16点 FFT数 字 角 频 率 w幅值A0 505101520x1的 8点 FFT数 字 角 频 率 w幅值A0 1005101520x1的 8点 FFT数 字 角 频 率 w幅值A
步骤 2：
X1=[0:7];X2=[0:15];Y18=cos(pi/4*X1);Y116=cos(pi/4*X2);Y28=cos(pi/4*X1)+cos(pi/8*X1);Y216=cos(pi/4*X2)+cos(pi/8*X2);Y18K8=fft(Y18,8);Y18K16=fft(Y18,16);Y116K16=fft(Y116,16);Y28K8=fft(Y28,8);Y28K16=fft(Y28,16);Y216K16=fft(Y216,16);figure(2);subplot(2,3,1);stem(abs(Y18K8));axis([0,8,0,1.2*max(abs(Y18K8))]);title("x4 的 8 点抽样 8 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,2);stem(abs(Y18K16));axis([0,16,0,1.2*max(abs(Y18K16))]);title("x4 的 8 点抽样 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,3);stem(abs(Y116K16));axis([0,16,0,1.2*max(abs(Y116K16))]);title("x4 的 16 点抽样 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,4);stem(abs(Y28K8));axis([0,8,0,1.2*max(abs(Y28K8))]);title("x5 的 8 点抽样 8 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,5);stem(abs(Y28K16));axis([0,16,0,1.2*max(abs(Y28K16))]);title("x5 的 8 点抽样 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(2,3,6);stem(abs(Y216K16));axis([0,16,0,1.2*max(abs(Y216K16))]);
title("x5 的 16 点抽样 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");0 501234x4的 8点 抽 样 8点 FFT数 字 角 频 率 w幅值A0 1001234x4的 8点 抽 样 16点 FFT数 字 角 频 率 w幅值A0 1002468x4的 16点 抽 样 16点 FFT数 字 角 频 率 w幅值A0 50246x5的 8点 抽 样 8点 FFT数 字 角 频 率 w幅值A0 100246x5的 8点 抽 样 16点 FFT数 字 角 频 率 w幅值A0 1002468x5的 16点 抽 样 16点 FFT数 字 角 频 率 w幅值A
步骤 3：
t16=[0:1/64:15/64];t32=[0:1/64:31/64];t64=[0:1/64:63/64];Y16=cos(8*pi*t16)+cos(16*pi*t16)+cos(20*pi*t16);Y32=cos(8*pi*t32)+cos(16*pi*t32)+cos(20*pi*t32);Y64=cos(8*pi*t64)+cos(16*pi*t64)+cos(20*pi*t64);YK16=fft(Y16,16);YK32=fft(Y32,32);YK64=fft(Y64,64);figure(3);subplot(1,3,1);stem(abs(YK16));axis([0,16,0,1.2*max(abs(YK16))]);title("x8 的 16 点抽样 16 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(1,3,2);stem(abs(YK32));axis([0,32,0,1.2*max(abs(YK32))]);title("x8 的 32 点抽样 32 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");subplot(1,3,3);stem(abs(YK64));axis([0,64,0,1.2*max(abs(YK64))]);title("x8 的 64 点抽样 64 点 FFT");xlabel("数字角频率 w");ylabel("幅值 A");
