数字信号处理教程——MATLAB释义与实现 (第3版)

数字信号处理教程——MATLAB释义与实现 (第3版)

本教程系统地讲解了数字信号处理的基本理论和方法，注重物理概念的阐述，尽量把MATLAB的多媒体功能用于图形、动画和声音等形象的演示，是一本自始至终使用MATLAB来阐述问题和进行计算的数字信号处理教材。本教程把《数字信号处理》课程定位为讲述“用数字方法处理（任何）信号”的技术，因此，对模拟信号的计算机处理给予了充分的篇幅，而这正是目前许多数字信号处理教材的弱点。

本书适合作为电子信息专业的大学本科学生教材，对于学过数字信号处理课程，但没有用MATLAB来解决有关问题的研究生和工程技术人员，也有很好的参考价值。

Retrieve Companion Software

1. %《数字信号处理教程——MATLAB释义与实现》第一章例1.1.1程序hc111
   
2. % 连续信号、离散信号和数字信号
   
3. % 电子工业出版社出版  陈怀琛编著 2004年9月
   
4. %
   
5. dt=0.001;t=0:dt:6;                % 建立连续自变量向量
   
6. xa=sqrt(t)+cos(t);                % 原始的连续时间信号xa(t)
   
7. T=0.5;n=0:6/T;                                % T为采样周期，建立离散自变量向量
   
8. x=sqrt(n*T)+cos(n*T);                 % 采样周期为T的离散时间信号x(n)
   
9. deltax=0.5;                                        % deltax为x的量化单位
   
10. xq=round(x/deltax)* deltax;                % 舍入量化后的数字信号xq
    
11. subplot(1,2,1),
    
12. plot(t,xa,':'),hold on,grid on                % 画出连续时间信号曲线
    
13. plot(n*T,x,'o')                                                % 画出离散时间信号曲线
    
14. stem(n*T,xq,'*')                                        % 画出画出数字信号曲线
    
15. grid on
    
16. legend('连续信号xa','离散时间信号x','数字信号xq')  % 画出图例标注
    
17. subplot(1,2,2)
    
18. stairs(n*T,xq),        grid on        % 画出数字信号采样保持后恢复的连续信号曲线
    
19. legend('将数字信号采样保持','恢复后连续信号曲线')    % 画出图例标注
    
20. set(gcf,'color','w')     % 将图的背景色置为白色

复制代码

1. %《数字信号处理教程——MATLAB释义与实现》第二章例2.2.1程序hc221
   
2. % 离散信号的基本序列
   
3. % 电子工业出版社出版  陈怀琛编著 2004年9月
   
4. %
   
5. clear,n0=0;nf=10;ns=3;n03=-2;
   
6. % n1=n0:nf;x1=[zeros(1,ns-n0),1,zeros(1,nf-ns)];  %单位脉冲序列的产生
   
7. n1 = n0:nf; x1=[(n1-ns)==0];         % 显然,用逻辑式是比较高明的方法
   
8. n2=n0:nf;x2=[zeros(1,ns-n0),ones(1,nf-ns+1)];  %单位阶跃序列的产生
   
9. % 也有类似的用逻辑比较语句的方法,留给读者思考
   
10. n3 = n03:nf; x3=exp((-0.2+0.5j)*n3);        % 复数指数序列
    
11. subplot(2,2,1),stem(n1,x1);title('单位脉冲序列')
    
12. axis([0,10,0,1.1])
    
13. subplot(2,2,3),stem(n2,x2,'.');title('单位阶跃序列') % 用小圆点画序列
    
14. axis([0,10,0,1.1])
    
15. subplot(2,2,2),stem(n3,real(x3),'x');line([-5,10],[0,0]) % 画横坐标
    
16. title('复指数序列'),ylabel('实部')
    
17. subplot(2,2,4),stem(n3,imag(x3),'filled');  % 用实心圆点画序列
    
18. line([-5,10],[0,0]),ylabel('虚部')
    
19. set(gcf,'color','w')                        % 置图形背景色为白

复制代码

1. %《数字信号处理教程——MATLAB释义与实现》第二章例2.2.1程序hc221
   
2. % 离散信号的基本序列
   
3. % 电子工业出版社出版  陈怀琛编著 2004年9月
   
4. %
   
5. clear,n0=0;nf=10;ns=3;n03=-2;
   
6. % n1=n0:nf;x1=[zeros(1,ns-n0),1,zeros(1,nf-ns)];  %单位脉冲序列的产生
   
7. n1 = n0:nf; x1=[(n1-ns)==0];         % 显然,用逻辑式是比较高明的方法
   
8. n2=n0:nf;x2=[zeros(1,ns-n0),ones(1,nf-ns+1)];  %单位阶跃序列的产生
   
9. % 也有类似的用逻辑比较语句的方法,留给读者思考
   
10. n3 = n03:nf; x3=exp((-0.2+0.5j)*n3);        % 复数指数序列
    
11. subplot(2,2,1),stem(n1,x1);title('单位脉冲序列')
    
12. axis([0,10,0,1.1])
    
13. subplot(2,2,3),stem(n2,x2,'.');title('单位阶跃序列') % 用小圆点画序列
    
14. axis([0,10,0,1.1])
    
15. subplot(2,2,2),stem(n3,real(x3),'x');line([-5,10],[0,0]) % 画横坐标
    
16. title('复指数序列'),ylabel('实部')
    
17. subplot(2,2,4),stem(n3,imag(x3),'filled');  % 用实心圆点画序列
    
18. line([-5,10],[0,0]),ylabel('虚部')
    
19. set(gcf,'color','w')                        % 置图形背景色为白

复制代码

1. %《数字信号处理教程——MATLAB释义与实现》第二章例2.2.3程序hc223
   
2. % 离散信号序列的移位、相加和相乘
   
3. % 电子工业出版社出版  陈怀琛编著 2004年9月
   
4. %
   
5. x1=[0,1,2,3,4,3,2,1,0];ns1=-2;                % 给定x1及ns1
   
6. x2=[2,2,0,0,0,-2,-2]; ns2=2;                % 给定x2及ns2
   
7. nf1=ns1+length(x1)-1; nf2=ns2+length(x2)-1;
   
8. ny= min(ns1,ns2):max(nf1,nf2);                 % y(n)的位置向量
   
9. y1 = zeros(1,length(ny)); y2 = y1;       % 延拓序列初始化
   
10. y1(find((ny>=ns1)&(ny<=nf1)==1))=x1; % 给y1赋值x1
    
11. y2(find((ny>=ns2)&(ny<=nf2)==1))=x2; %  给y2 赋值x2
    
12. ya = y1 + y2;                                     % 序列相加
    
13. yp = y1.* y2;                                     % 序列相乘
    
14. subplot(4,2,1), stem(ns1:nf1,x1,'.')                % 绘图
    
15. xlabel('nx1'),ylabel('x1'),axis([-5,10,0,4])
    
16. subplot(4,2,3), stem(ns2:nf2,x2,'.'),axis([-5,10,-2,2])
    
17. xlabel('nx2'),ylabel('x2')
    
18. subplot(4,2,2), stem(ny,y1,'.')                % 绘图
    
19. xlabel('ny'),ylabel('y1')
    
20. subplot(4,2,4), stem(ny,y2,'.')
    
21. xlabel('ny'),ylabel('y2')
    
22. line([ny(1),ny(end)],[0,0])                % 画x轴
    
23. subplot(4,2,6), stem(ny,ya,'.')
    
24. xlabel('ny'),ylabel('ya')
    
25. line([ny(1),ny(end)],[0,0])                % 画x轴
    
26. subplot(4,2,8), stem(ny,yp,'.')
    
27. xlabel('ny'),ylabel('yp')
    
28. line([ny(1),ny(end)],[0,0])                % 画x轴
    
29. set(gcf,'color','w')            % 置图形背景色为白

复制代码