资源描述：
《数字信号处理实验例题.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验一离散时间信号与系统例1-1用MATLAB计算序列{-201–13}和序列{120-1}的离散卷积。解MATLAB程序如下：a=[-201-13];b=[120-1];c=conv(a,b);M=length(c)-1;n=0:1:M;stem(n,c); xlabel('n');ylabel('幅度');图1.1给出了卷积结果的图形，求得的结果存放在数组c中为：{-2-413151-3}。 例1-2用MATLAB计算差分方程当输入序列为时的输出结果。解MATLAB程序如下：N=41;a
2、=[0.8-0.440.360.22];b=[10.7-0.45-0.6];x=[1zeros(1,N-1)];k=0:1:N-1;y=filter(a,b,x);stem(k,y)xlabel('n');ylabel('幅度') 图1.2给出了该差分方程的前41个样点的输出，即该系统的单位脉冲响应。 例1-3用MATLAB计算例1-2差分方程所对应的系统函数的DTFT。解例1-2差分方程所对应的系统函数为： 其DTFT为
3、 用MATLAB计算的程序如下：k=256;num=[0.8-0.440.360.02];den=[10.7-0.45-0.6];w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h));gridtitle('实部')xlabel('omega/pi');ylabel('幅度')subplot(2,2,2);plot(w/pi,imag(h));gridtitle('虚部')xl
4、abel('omega/pi');ylabel('Amplitude')subplot(2,2,3);plot(w/pi,abs(h));gridtitle('幅度谱')xlabel('omega/pi');ylabel('幅值')subplot(2,2,4);plot(w/pi,angle(h));gridtitle('相位谱')xlabel('omega/pi');ylabel('弧度') 实验二离散傅里叶变换及其快速算法 例2-1对连续的单一频率周期信号按采样频率
5、采样，截取长度N分别选N=20和N=16，观察其DFT结果的幅度谱。解此时离散序列，即k=8。用MATLAB计算并作图，函数fft用于计算离散傅里叶变换DFT，程序如下：k=8;n1=[0:1:19];xa1=sin(2*pi*n1/k);subplot(2,2,1)plot(n1,xa1)xlabel('t/T');ylabel('x(n)');xk1=fft(xa1);xk1=abs(xk1);subplot(2,2,2)stem(n1,xk1)xlabel('k');ylabel('X(k)');n2=[0:1
6、:15];xa2=sin(2*pi*n2/k);subplot(2,2,3)plot(n2,xa2)xlabel('t/T');ylabel('x(n)');xk2=fft(xa2);xk2=abs(xk2); subplot(2,2,4)stem(n2,xk2)xlabel('k');ylabel('X(k)'); 计算结果示于图2.1，(a)和(b)分别是N=20时的截取信号和DFT结果，由于截取了两个半周期，频谱出现泄漏；(c)和(d)分别是
7、N=16时的截取信号和DFT结果，由于截取了两个整周期，得到单一谱线的频谱。上述频谱的误差主要是由于时域中对信号的非整周期截断产生的频谱泄漏。例2-2用FFT计算两个序列的互相关函数。解用MATLAB计算程序如下：x=[13-112331];y=[21-1120-13];k=length(x);xk=fft(x,2*k); yk=fft(y,2*k);rm=real(ifft(conj(xk).*yk));rm=[rm(k+2:2*k)rm(1:k)];m=(-k+1):(k-1
8、);stem(m,rm)xlabel('m');ylabel('幅度'); 其计算结果如图2.2所示。 例2-3计算两个序列的的互相关函数，其中x(n)={23521–100123530–1–2012}y(n)=x(n-4)+e(n),e(n)为一随机噪声，在MATLAB中可以用随机函数ran
