（全新）程佩青数字信号处理第三版笔记.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 120 计算方法： 采样间隔：Ts = 1/fs; 信号持续时间：T = Ts*N 频域分辨率为：F = 1/T; 每个k代表的频率值， K1 = f1/F= K2 = f2/F= K3 = f3/F= Ps：在使用data cursor时候，我们会发现，每次只能显示一个点的数据，当需要显示多个点的时候就会觉得很不方便。告诉大家一个小窍门，当你选择data cursor工具时，按住alt键，点左键选择曲线上的点。那你就可以想要多少个显示就要多少个了。/viewthread.php?tid=59001sid=6enFyn DFT末尾补零效果 原来x(n)的DFT为X(k)，在末尾补零rN个后成为y(n)，则 Y(k) = X(k/r)，在中间相当于插值，在频域插值。 然后，末尾补零x(n)的Z变换和DTFT是不会改变的仍为X(z), X(exp(jw)) % padding zero in the end, check effect to Xk clc; clear; f1=2.67; f2=3.75; f3=6.75; fs=20; N = 16; n = 0:N-1; x = sin(2*pi*f1*n/fs) + sin(2*pi*f2*n/fs) + sin(2*pi*f3*n/fs); X = fft(x, N); X = abs(X); k = fs/N * (0:N-1); %fs/N是频域分辨率 stem(k, X(1:N), o, r); hold on % padding with rN zeros r = 2; x(N+1:r*N) = 0; X = fft(x, r*N); X = abs(X); k = fs/(r*N) * (0:r*N-1); stem(k, X(1:r*N), *); hold on % padding with rN zeros r = 4; x(N+1:r*N) = 0; X = fft(x, r*N); X = abs(X); k = fs/(r*N) * (0:r*N-1); stem(k, X(1:r*N), , y); r = 2, 即补零N个时，X(k/2)，在频域每两个k中间插入r-1=1个值 r = 3, 即补零N个时，X(k/3)，在频域每两个k中间插入r-1=2个值 红色圆圈为不补零时结果 蓝色*为补零rN = 2N个时，在频域每两个k中间插入r-1=1个值 黄色为补零rN = 4N个时，在频域每两个k中间插入r-1=3个值 由重合的可以看出，原来的数值仍在原来的地方，只是重合。 序列插值之后变化 % Interpolation 0, DFT effect clc; clear; f1=2.67; f2=3.75; fs=20; N = 16; n = 0:N-1; x = sin(2*pi*f1*n/fs) + sin(2*pi*f2*n/fs); X = fft(x, N); X = abs(X); k = fs/N * (0:N-1); %fs/N是频域分辨率 subplot(2,1,1); stem(k, X(1:N), o, r); % 内插