figure(2) subplot(2,1,1) plot(t, modulated) title('DSB调幅信号') xlim([0 0.02]) xlabel('时间') ylabel('幅度') subplot(2,1,2) NFFT = 2^nextpow2(length(filtered)); Y = fft(filtered, NFFT)/length(filtered); f = fs*linspace(0,1,NFFT+1);

subplot(2,1,2) 绘制的是对调幅信号进行滤波后的频谱图像。NFFT = 2^nextpow2(length(filtered)) 是为了对滤波后的信号进行FFT变换时使用的变换点数，nextpow2() 函数得到大于等于length(filtered)的最小2的整数次幂。Y = fft(filtered, NFFT)/length(filtered) 是对滤波后的信号进行FFT变换并归一化，得到频域中的能量值。f = fs*linspace(0,1,NFFT+1) 是生成频率向量，长度为NFFT+1，用于绘制频谱图。最后，subplot(2,1,2) 中使用 plot(f(1:NFFT/2),20*log10(abs(Y(1:NFFT/2)))) 绘制信号的单边频谱。其中，f(1:NFFT/2) 取频率向量的前一半，因为后一半是对称的。20*log10(abs(Y(1:NFFT/2))) 是对能量值取对数，并转换为分贝单位，可以更清晰地显示信号频域特征。