2012数字信号处理实验讲义-实验一报告

《2012数字信号处理实验讲义-实验一报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验一离散时间系统时域、变换域分析班级：电子101班姓名：李怡敏学号：101045一、实验目的1.熟悉信号处理软件MATLAB的使用。2.离散信号的基本运算实现。3.了解基本序列及复杂序列的产生方法。4.运用卷积方法观察系统的吋域特性。5.熟悉信号的采样与恢复。6.掌握线性时不变系统的变换域表示方法。7.掌握序列傅里叶变换性质。二、实验内容(一)验证性实验内容1.熟悉扩展函数2.运行例题程序(二)提高综合运用能力部分内容1.编写利用扩展函数产牛•下列序列并画图的程序(a)兀(斤)=2*/(〃+2)—3(斤一

2、4)-5<=n<=5>>n=[-5:5];>>x=2*quyang(-2,-5,5)-quyang(4,-5,5);>>stem(n,x);»xlabel('n');»ylabel('x(n)');>>axis([-5,5,-2,5]);5r4-Ckee-1・eeeeeo'-5-4-3-2-1012345n(b)x(/t)=cos(0.04加)和y(〃)=cos(0.04如)+0.2w(〃)0<=n<=50w(n)为白噪声函数为w=randn(size(n))n=0:50;x=cos(0.04*n*pi);

3、y二x+0.2*randn(size(n));subplot(2,1,1);stem(n,x);xlabel('n');ylabel('x(n)');subplot(2,1,2);stem(n,y);xlabel('n');ylabel('y(n)');n21,雀0-1'205101520i'25n3035404550或者:»»»»»stem(n,y);>>xlabel('n');>>ylabel('y(n)');y=Columnsn=[0:50];x=cos(0.04*pi*n);w=randn(size

4、(n));y二sigadd(x,n,w,n)1.61451.49992.66101.52110.23271.1894-0.28010.6179Columns9through160.48760.4258-0.00881.2824-1.81120.36540.70830.4219Columns17through240.1521-0.49550.0397-0.1601-1.0647-1.2538-1.2257-2.4437Columns25through32-1.2261-0.8816-0.67730.4749-

5、1.2808-0.2531-0.01000.2119Columns33through40-1.6295-0.3238-0.1879-1.3168-0.92941.0195-0.06870.5773Columns41through480.3970-0.2097-0.02371.0811-0.22091.59021.44530.1081Columns49through510.7030-0.1957-1.2023through811.编写程序求系统输出y(n)的零状态响应并画图。设线性移不变系统的抽样响应为h(n

6、)=(0.9)"u(n),输入序列为x(n)=u(n)-u(n-10)»n=0:l:10;»h二0.9.八n;»x=jieyue(0,0,10)-jieyue(10,0,10);»y=juanji(h,x)实际答案:n二0:30;h=(0.9.”n);x=jieyue(0,0,20)-jieyue(10,0,20);[y,ny]=conv_m(x,n,h,n);h=[hzeros(1,(length(ny)-length(h)))];x=[xzeros(1,(length(ny)-length(x)))];

7、subplot(3,1,1);stem(ny,h);subplot(3,1,2);stem(ny,x);subplot(3,1,3);stem(ny,y);1.编写程序实现序列傅里叶变换频移性质、时域卷积后的频域特性(时域卷积定理的程序设计实现)。function[f,k]=juanji(fl,f2,kl,k2)%f=fl*f2%kl:fl对应序列号%k2:f2对应序列号w=-pi:2*pi/255:pi;f=conv(fl,f2)hl=freqz(xl,1,w);h2=freqz(x2,1,w);hp=h

8、l.*h2;h3=freqz(f,1,w);subplot(2,2,1)stem(kl,fl)title('fl(k)')xlabel('k')ylabel('fl(k)')subplot(2,2,2)plot(w/pi,abs(hp));grid;subplot(2,2,2);plot(w/pi,abs(h3));grid;subplot(2,2,3);plot(w/pi,angle(hp));grid;sub