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");
数字信号处理实验(实验一)数字信号实验报告vfp实验一必修一实验