资源描述:
《数字信号实验报告new》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数字信号处理导论实验报告姓名:张伦裕学号:200905090205实验一信号、系统及系统响应当系统的输入输出差分方程为:Y(n)-0.8y(n-1)-0.5y(n-2)=0.7x(n)+0.3x(n-1)通过MATLAB编程实现并考虑如下问题:(1)当系统的输入为单位冲激函数时,分别利用filter函数和impz函数得到的系统单位冲激响应曲线。(2)当系统的输入为单位阶跃函数时,分别利用filter函数和impz函数得到的系统单位阶跃响应曲线。(3)对于不同的输入,系统的输出有什么变化,试讨论之。一、实验原理系统的零状态响应就是在系统初始状态为零的条件下微分方程的
2、解。y=filter(b,a,f)其中b,a分别为系统差分方程左右端系数向量,f为数输入向量。impz函数可直接求出离散系统的单位冲激响应,为:impz(b,a)。conv函数可用于计算离散序列的线性卷积。为:y=conv(x,h)。可用impz函数和conv函数计算对任意输入的离散系统的响应,即由impz求出单位冲激响应,再用单位冲激应与输入序列卷积。二、程序n=0:20;A=[1,-0.8,-0.5];B=[0.7,0.3];y1=filter(B,A,u1);[y2]=impz(B,A,21);Y1=filter(B,A,u2);[h]=impz(B,A,2
3、1);Y22=conv(h,u2);Y2=Y22(1:21);subplot(221);stem(n,y1);title('filterfunction');gridon;subplot(222);stem(n,y2);title('impzfunction');gridon;subplot(223);stem(n,Y1);title('filterfunction');gridon;subplot(224);stem(n,Y2);title('impzfunction');gridon;三、结果及分析:实验二使用FFT作频谱分析编程实现下列问题并讨论所得到的结果
4、。(1)使用FFT对MATLAB中的数据noissin.dat进行谱分析。(2)使用FFT对语音数据noisyspeech.wav进行谱分析。一、实验原理(1)离散傅里叶变换(DFT)公式为:X(k)=∑x(n)W^nk;x(n)=∑X(k)W^-nk;其中w=e^(2∏nk/N),N为离散序列的长度。(2)快速傅里叶变换(FFT)是利用w因子的取值特点,减少DFT的复数乘法的次数。其中一种是时间抽取基2算法,它将时间按奇偶逐级分开,直到两点的DFT。MATLAB提供了fft函数可用于计算FFT,器调用形式为;X=fft(x)或X=fft(x,N),N为2的整数次
5、幂,若x的长度小于N,则补零,若超过则舍去N后的数据。(3)函数形式[y,fs,bits]=wavread('Blip',[N1N2]);用于读取语音,采样值放在向量y中,fs表示采样频率(Hz),bits表示采样位数。[N1N2]表示读取从N1点到N2点的值(若只有一个N的点则表示读取前N点的采样值)。函数形式s=noissin(n1,n2)用于读取MALAB的噪声信号。二、程序xx=wavread('noisyspeech.wav');fs=100;N=128;x=xx(1:N);n=1:N;X=abs(fft(x,N));subplot(221);plot(
6、n,x);xlabel('n');ylabel('x(n)');gridon;subplot(222);plot(n,X);xlabel('k');ylabel('
7、X(k)
8、');gridon;loadnoissin;s=noissin(1:20);S=fft(s);subplot(223);plot(abs(s));xlabel('n');ylabel('
9、s(n)
10、');gridon;subplot(224);plot(abs(S));xlabel('k');ylabel('
11、S(k)
12、');gridon;三、结果与分析实验三使用双线性Z变换设计IIR滤波器
13、编程实现下列问题并讨论所得到的结果。(1)使用FFT对MATLAB中的数据noissin.dat进行谱分析。(2)使用FFT对语音数据noisyspeech.wav进行谱分析。一、实验原理(1)设计滤波器就是要设计一个系统是其能让一定频率的波段通过或滤去,对IIR滤波器,器转移函数是:H(Z)=(∑bz^(-r))/(1+∑az^(-k))。(2)设计的一般原则:若使滤波器拒绝某个频率,应在单位园上相应的频率处设置一个零点,反之则设置一个极点。(3)低通数字IIR滤波器设计步骤:给出数字低通滤波器的技术要求→映射为模拟低通的技术要求→归一化为模拟低通滤波器的技术要
14、求→设计出