数字信号处理实验材料-2010级.pdf

数字信号处理实验材料 （2010级数字媒体技术专业） Deng Z.H Qian P.J 整理 试验1 离散时间信号的MATLAB实现 第一部分 序列的生成 主要 第二部分 离散序列的基本运算与波形变换的MATLAB实现 内容 第三部分 线性时不变系统求解的MATLAB实现 第四部分 序列的傅里叶变换与信号重构MATLAB实现 试验 Matlab 平台 平台 试验 以第四部分序列的傅里叶变换与信号重构MATLAB实现为内容写一份试验 报告 报告 第一部分 序列的生成 1．正弦序列 离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函 π ⎛ ⎞ 数来画出序列的波形。下面就是正弦序列sin k 的MATLAB源程序。程序运行结果 ⎜ ⎟ 6 ⎝ ⎠ 如图1.1所示。 例1-1 %正弦序列实现程序 k=0:39; fk=sin(pi/6*k); stem(k,fk) 图1.1 正弦序列波形 2．指数序列 离散指数序列的一般形式为cak ，可用MATLAB中的数组幂运算（即点幂运算） c*a.^ k 来实现。下面为用MATLAB编写绘制离散时间实指数序列波形的函数。 实指数序列 function dszsu(c,a,k1,k2) %c：指数序列的幅度 %a：指数序列的底数 %k1：绘制序列的起始序号 %k2：绘制序列的终止序号 k=k1:k2; x=c*(a.^k); stem(k,x,filled) hold on plot([k1,k2],[0,0]) hold off 利用上述函数，实现实指数波形MATLAB程序如下（其中 a 值分别为 5 3 5 3 ， ， ， ）。 − − 4 4 4 4 例1-2 %离散时间实指数序列实现程序 subplot 221; dszsu(1,5/4,0,20); xlabel(k); title(f1[k]); subplot 222 dszsu(1,3/4,0,20); xlabel(k); title(f2[k]); subplot 223; dszsu(1,-5/4,0,20); xlabel(k); title(f3[k]); subplot 224;