信号与系统实验指导书(带源程序)解析.docVIP

实验一 离散时间系统与MATLAB 一. 实验目的 1. 进一步加深对离散时间系统的理解。 2. 学习在MATLAB中怎样表示离散时间信号。 3. 熟悉离散时间信号的作图。 二. 实验步骤 1. 复习离散时间系统的有关内容。 2. 复习MATLAB的基本语法。 3. 按实验内容熟悉 stem。 4. 编写程序。 5. 输出结果，总结结论，按要求写出实验报告。 三. 实验内容 1．掌握stem函数 STEM(Y) plots the data sequence Y as stems from the x axis terminated with circles for the data value. STEM(X,Y) plots the data sequence Y at the values specified in X. 例：t=[0:0.1:2]; x=cos(pi*t+0.6); stem(t,x); xn=[4,2,2,3,6,7]; stem(xn); 思考：STEM(Y)与STEM(X,Y)有什么不同？STEM与PLOT函数有什么不同？ 2．掌握subplot函数 H = SUBPLOT(m,n,p), or SUBPLOT(mnp), breaks the Figure window into an m-by-n matrix of small axes, selects the p-th axes for the current plot, and returns the axis handle. The axes are counted along the top row of the Figure window, then the second row, etc. 例： n1=0:3;x1=[1,1,1,1];subplot(221);stem(n1,x1);title(x1序列); n2=0:7;x2=[1,2,3,4,4,3,2,1];subplot(222);stem(n2,x2);title(x2序列); n3=0:7;x3=[4,3,2,1,1,2,3,4];subplot(223);stem(n3,x3);title(x3序列); n4=0:7;x41=cos((pi/4)*n4);subplot(224);stem(n4,x41);title(x4序列); 思考：subplot是怎样分配各个作图分区的顺序号的？ 3．信号的运算 ，，请作出，的图形。 思考：假如与长度不同，序列的求和和乘积运算能否执行，结果怎样？ 程序：x1=[1,0.7,0.4,0.1,0] x2=[0.1,0.3,0.5,0.7,0.9] subplot(121);stem(x1+x2);title(x1+x2); subplot(122);stem(x1.*x2);title(x1*x2); 4．掌握freqs与freqz函数 freqs用于s域的频率响应的计算，f reqz用于z域的频率响应的计算。请参看两者的help文件。 例：建立2个函数文件，分别输入为系统函数的分子，分母，和角频率的最大值，输出为系统频率响应的db，幅度，相位，和角频率。参看下面的程序。 function[db,mag,pha,Omega]=freqs_m(b,a,Omega_Max) % s域频率响应的计算 Omega=[0:1:500]*Omega_Max/500; H=freqs(b,a,Omega); mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); function[db,mag,pha,w]=freqz_m(b,a,w_Max); % z域频率响应的计算 w=[0:1:499]*w_Max/500; H=freqz(b,a,w); mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); （1）若，请画出该系统幅频和相频特性。 程序： b=1 a=[1 5 6] Omega_Max=2*pi [db,mag,pha,Omega]=freqs_m(b,a,Omega_Max) subplot(221);plot(Omega,mag);title(AF的幅度响应); subplot(222);plot(Omega,db);title(AF的幅度响应db); subplot(223);plot(Omega,pha);title(AF的相位响应); （2）若，请画出该的幅频和相频特性。 程序： b=[1 1] a=[1 0.8 0.64] w_Max=2*pi [db,mag,p