山东大学第一二次数字信号处理实验报告.docxVIP

数字信号处理实验报告实验一、离散时间信号与系统1、实验目的（1）熟悉Matlab编程环境、编程方法、常用数学运算、函数和画图函数的应用。（2）用Matlab实现离散时间序列的基本操作。（3）区分conv和filter函数的异同。（4）掌握抽样定理。2、实验内容（1）请画出函数,其中取请画在一幅图上，以便于对比。实验程序：clc;clear all;close all;n=0:20;xn=0.5*sin(2*pi*0.01*n);yn=0.75*cos(2*pi*0.15*n);subplot(331);stem(n,xn)%画出x(n)散点图title(x(n))subplot(332);stem(n,yn)%画出y(n)散点图title(y(n))subplot(333);stem(n,xn+yn) %画出x(n)+y(n)散点图title(x(n)+y(n))subplot(334);stem(n,xn.*yn)%画出x(n)×y(n)散点图title(x(n)*y(n))subplot(335);stem(n,xn./yn)%画出x(n)/y(n)散点图title(x(n)/y(n))subplot(336);stem(n+3,xn)%画出x(n+3)散点图title(x(n+3))subplot(337);stem(n-5,xn)%画出x(n-5)散点图title(x(n-5))subplot(338);stem(-n,xn)%画出x(-n)散点图title(x(-n))实验结果：（2）产生复序列,画出其实部、虚部、模和相位函数。实验程序：clc;clear all;close all;n=-10:10;xn=exp(-0.2+0.3*n*i) %x(n)函数subplot(221)stem(n,real(xn))%real函数求实部title(实部)subplot(222)stem(n,imag(xn))%imag函数求虚部title(虚部)subplot(223)stem(n,abs(xn))%abs函数求绝对值，也就是模title(模)subplot(224)stem(n,phase(xn)) %phase函数求相位title(相位)实验结果：xn =1 至 10 列：-0.8105 - 0.1155i 、-0.7402 - 0.3499i 、 -0.6037 - 0.5530i、 -0.4133 - 0.7067i 、 -0.1860 - 0.7973i 0.0579 - 0.8167i 、 0.2967 - 0.7631i 、 0.5089 - 0.6413i 、 0.6757 - 0.4623i 、 0.7822 - 0.2420i11 至 20 列：0.8187 + 0.0000i 、 0.7822 + 0.2420i 、 0.6757 + 0.4623i 、 0.5089 + 0.6413i 、0.2967 + 0.7631i 、 0.0579 + 0.8167i、 -0.1860 + 0.7973i 、 -0.4133 + 0.7067i 、-0.6037 + 0.5530i 、 -0.7402 + 0.3499i21 列：-0.8105 + 0.1155i（3）设。分别用conv和filter实现，画图比较二者异同。实验程序：clc;clear all;close all;n=0:20;%n取0到20xn=0.5*sin(2*pi*0.01*n);hn=[1,0,0.5,0.25];subplot(221)stem(hn) %画出h(n)散点图title(h(n));subplot(222)stem(n,xn)%画出x(n)散点图title(x(n));subplot(223)y1n=conv(xn,hn)%用conv函数实现卷积stem(y1n)title(conv实现x(n)与h(n)卷积);subplot(224)y2n=filter(hn,1,xn);%用filter函数实现卷积stem(y2n)title(filter实现x(n)与h(n)卷积);实验结果：由实验结果的图可知用filter函数和conv函数实现卷积的结果并不相同，若要让filter显示完整的卷积结果，可以在x(n)后补上（h(n)长度 -1）个0后再用filter函数进行卷积。（4）已知,请选择合适的抽样频率对其进行抽样，所得离散时间信号不发生混叠失真。如果是周期性信号，应保持周期性，请画图表示。实验程序：clc;clear all;close all;dt=0.001;t=0:0.001:0.2;xt=cos(2*pi*50*t)+0.75*cos(2*pi*80*t);subplot(221)plot(t,xt)