数字信号处理MATLAB
第1页
习题数字信号处理MATLAB习题
M1-1 已知g1(t)?cos(6?t)，g2(t)?cos(14?t)，g3(t)?cos(26?t)，以抽样频率fsam?10Hz对上述三个信号进行抽样。在同一张图上画出g1(t)，
g2(t)和g3(t)及抽样点，对所得结果进行讨论。
解：
第2页
从以上两幅图中均可看出，三个余弦函数的周期虽然不同，但它们抽样后相应抽样点所对应的值都相同。那么这样还原回原先的函数就变成相同的，实际上是不一样的。这是抽样频率太小的原因，我们应该增大抽样频率才能真实还原。如下图：f=50Hz
第3页
程序代码
f=10;
t=-0.2:0.001:0.2; g1=cos(6.*pi.*t); g2=cos(14.*pi.*t); g3=cos(26.*pi.*t); k=-0.2:1/f:0.2; h1=cos(6.*pi.*k); h2=cos(14.*pi.*k); h3=cos(26.*pi.*k); % subplot(3,1,1);
% plot(k,h1,'r.',t,g1,'r'); % xlabel('t'); % ylabel('g1(t)'); % subplot(3,1,2);
% plot(k,h2,'g.',t,g2,'g'); % xlabel('t'); % ylabel('g2(t)'); % subplot(3,1,3);
% plot(k,h3,'b.',t,g3,'b'); % xlabel('t'); % ylabel('g3(t)');
plot(t,g1,'r',t,g2,'g',t,g3,'b',k,h1,'r.',k,h2,'g.',k,h3,'b.')
第4页
xlabel('t'); ylabel('g(t)');
legend('g1(t)','g2(t)','g3(t)');
M2-1 利用DFT的性质，编写一MATLAB程序，计算下列序列的循环卷积。
(1) g[k]={1,-3,4,2,0,-2,},h[k]={3,0,1,-1,2,1}； (2) x[k]=cos(?k/2),y[k]=3k,k=0,1,2,3,4,5。 解： (1)
循环卷积结果
6.0000 -3.0000 17.0000 -2.0000 7.0000 -13.0000
程序代码
第5页