DSP实验讲稿.docVIP

DSP实验讲稿 ?第1章? 离散时间系统和信号 ?第2章? 离散傅里叶变换及其快速算法 第3章? 无限长单位脉冲响应(IIR)滤波器的设计方法 第4章? 有限长单位脉冲响应(FIR)滤波器的设计方法 ?第5章? 数字信号处理系统的实现 ? 第章? 多采样率信号处理第１章? 离散时间信号与系统 例用MATLAB计算序列{-2? 0? 1? –1? 3}和序列{1? 2? 0? -1}的离散卷积。 解 MATLAB程序如下： ??????????????????????? a=[-2 0 1 -1 3]; ??????????????????????? b=[1 2 0 -1]; ??????????????????????? c=conv(a,b); ??????????????????????? M=length(c)-1; ??????????????????????? n=0:1:M; ??????????????????????? stem(n,c); ?????????????? xlabel(n); ylabel(幅度); ?????????????????????? ????????? 图1.1给出了卷积结果的图形，求得的结果存放在数组c中为：{-2? -4? 1? 3? 1? 5? 1? -3}。 ? 例 用MATLAB计算差分方程 ?????????? 当输入序列为 时的输出结果 。 解 MATLAB程序如下： ?????????????????? N=41; ?????????????????? a=[0.8 -0.44 0.36 0.2]; ?????????????????? b=[1 0.7 -0.45 -0.6]; ?????????????????? x=[1 zeros(1,N-1)]; ?????????????????? k=0:1:N-1; ?????????????????? y=filter(a,b,x); ?????????????????? stem(k,y) ?????????????????? xlabel(n);ylabel(幅度)?????????? 图 1.2 给出了该差分方程的前41个样点的输出，即该系统的单位脉冲响应。 ??????????????????? ? 　 　 用MATLAB计算例差分方程 ?????????????? 所对应的系统函数的DTFT。 解差分方程所对应的系统函数为： ?????????????????????????????? ??????? 其DTFT为 ?????????????????????????????? ????? 用MATLAB计算的程序如下：?? ?????????????????????????? k=256; ?????????????????????????? num=[0.8 -0.44 0.36 0.02]; ?????????????????????????? den=[1 0.7 -0.45 -0.6]; ?????????????????????????? w=0:pi/k:pi;%改为w=-pi:pi/k:pi,比较结果 ?????????????????????????? h=freqz(num,den,w); ?????????????????????????? subplot(2,2,1); ?????????????????????????? plot(w/pi,real(h));grid ?????????????????????????? title(实部) ?????????????????????????? xlabel(\omega/\pi);ylabel(幅度) ?????????????????????????? subplot(2,2,2); ?????????????????????????? plot(w/pi,imag(h));grid ?????????????????????????? title(虚部) ?????????????????????????? xlabel(\omega/\pi);ylabel(Amplitude) ?????????????????????????? subplot(2,2,3); ?????????????????????????? plot(w/pi,abs(h));grid ?????????????????????????? title(幅度谱) ?????????????????????????? xlabel(\omega/\pi);ylabel(幅值) ??????????