新）程佩青数字信号处理第三版笔记.docVIP

程佩青41页，两个余弦信号采样，采样值一样 % P41, sin sampling clc; clear; fs = 100; N = 512; n = 0:N-1; x1 = cos(2*pi*20*n/fs); x2 = cos(2*pi*120*n/fs); figure(1); subplot(2,1,1);stem(x1) subplot(2,1,2);stem(x2) F1 = fft(x1,N); F2 = fft(x2,N); figure(2); subplot(2,1,1);plot(abs(F1)); subplot(2,1,2);plot(abs(F2)); 原因，两个信号频率关于fs对称，高频混叠到低频处与低频正好重合 所以，10，110；30，130；40，140等等，都会有这个问题。 P43页习题14，对信号做fft % p43, exercise14 clc; clear; fs = 3000; N = 512; n = 0:N-1; x = ( 1+cos(2*pi*100*n/fs) ).*cos(2*pi*600*n/fs); F = fft(x,N); plot(abs(F)); ylabel(|X(k)|); xlabel(k); 三个峰值分别对应500，600，700Hz，对应k值分别为 86 103 12