数字信号处理96105.docVIP

西安院 A sinusoid of frequency ω0=0.1πand duration of 300 samples, that is ,0≤n300, is input to a (causal) filter with transfer function ,where a=0.97。Adjust the scale factor b such that the filter’s gain at ω0 is unity。 Determine and plot the input x(n) and the output y(n) of the filter over the interval 0≤n450, by iterating the difference equation of the filter. At the same time, plot the filter’s magnitude response. 实验目的 给定一个系统的传输函数，输入一个正弦序列，经过系统的传输函数后求其差分方程，同时考虑零点，极点对幅频特性曲线的影响。 实验原理 一般时域离散系统或网络可以用差分方程，单位脉冲响应以及系统函数进行描述。 实验步骤简述 求系数b的值，利用其模值求出b的值。求出其差分方程y(n)，利用滤波函数将其传输函数进行逆变换。求出传输函数的幅频特性曲线，画出输入正弦函数的图形，画出差分方程的图形 实验程序框图 实验源程序 （1）close all; clear all; clc %tic a1=0.97; w=0.1*pi; b=abs(1-a1*exp(-j*w)); N=300; n=1:N; xn=sin(n.*w); x1=[xn zeros(1,150)]; a=[1 -a1]; yn=filter(b,a,x1); figure(1); subplot(211); plot(xn); title( x(n)); xlabel(n); ylabel(·ù?è); %t1=toc subplot(212); plot(yn); title( y(n)); xlabel(n); ylabel(·ù?è); figure(2); freqz(b,a); 实验结果： （2）close all; clear all; clc %tic a1=0.20; w=0.1*pi; b=abs(1-a1*exp(-j*w)); N=300; n=1:N; xn=sin(n.*w); x1=[xn zeros(1,150)]; a=[1 -a1]; yn=filter(b,a,x1); figure(1); subplot(211); plot(xn); title( x(n)); xlabel(n); ylabel(·ù?è); %t1=toc subplot(212); plot(yn); title( y(n)); xlabel(n); ylabel(·ù?è); figure(2); freqz(b,a); 实验结果： 6. 实验心得：通过实验理解了极点愈靠近圆，峰值愈高与尖锐，零点愈靠近单位园，谷值愈接近零 实验五：DTMF信号的生成是通过将两个有限长数字序列相加而实现；而对DTMF信号的检测是通过计算DTMF信号的DFT；然后测量在给定8个频率上的能量而实现。已知采样频率为8KHz，DTMF信号点数N＝205，对DTMF信号进行205点DFT。右表为DTMF数字。 1209Hz 1336Hz 1447Hz 1663Hz 1 697Hz 1 2 3 A 2 770Hz 4 5 6 B 3 852Hz 7 8 9 C 4 941Hz * 0 # D 要求：① 通过键盘任意输入16个键之一，生成DTMF信号。 ② 试实现对该输入信号的检测，并在屏幕上显示。③试从计算量角度考虑为什么采用DFT而不用FFT进行DTMF信号检测。④ 判断出每个频率对应的DFT的频率采样点。 提示：各键对应的ASCII码为0—9对应48—57，A—D对应65—68，#对应35，*对应42。 实验目的 实现两个数字序列信号的叠加，实现对DTMF信号进行DFT变换 实验原理 FFT是DFT快速算法，如果用算法执行这个DFT的计算,计算的次数是Nlog2N. 在这种情况下，如果仅想计算DFT中的M个点,Mlog2N,那么直接计算DFT更有效，利用DFT对输入信号进行检测。 实验步骤简述 生成两个信号，同时将两个进行叠加生成DTMF。将十六个键符转换成ASCII码值，将8个频率对应的采样点数算出。利用两重循环判断输入的键符是否为十六个键符之一，同时确定其对应的两个频率和采样点数。调