资源描述：
《实验一 离散系统时域分析与Z 域分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验一离散系统时域分析与Z域分析一．实验目的(1)熟悉Matlab的主要操作命令；(2)掌握典型离散信号的特性(3)掌握离散系统的频率响应二．实验原理1．系统单位脉冲响应与系统频率响应系统单位脉冲响应：系统函数：单位圆上序列的Z变换为序列的傅里叶变换2．系统幅频响应与相频响应系统频率响应三．实验中调用的函数eig(x)：求解x矩阵的特征相量和特征值。fliplr(x)：是将x矩阵的列倒序。[h,w]=freqz(b,a)：系统函数的频率响应h，w表示0～π范围，其中b表示系统函数的分子系数，按降幂顺序排列，a表示系统函数的分母系数，按降幂顺序排列。20*log10(abs(h))：对
2、系统函数的的频率响应h取绝对值得到幅频响应h（分贝）angle(h)：是对系统函数的的频率响应h求取相频响应四．实验内容(1)数组的加、减、乘、除和乘方运算。输入A=[1234],B=[3456],求C=A+B,D=A-B,E=A.*B,F=A./B,G=A.^B,用stem语句画出A，B，C，D，E，F，G。程序如下：a=[1234]b=[3456]c=a+bd=a-be=a.*bf=a./bg=a.^bsubplot(2,4,1)-7-stem(a)title('a')subplot(2,4,2)stem(b)title('b')subplot(2,4,3)stem(c)titl
3、e('c')subplot(2,4,4)stem(d)title('d')subplot(2,4,5)stem(e)title('e')subplot(2,4,6)stem(f)title('f')subplot(2,4,7)stem(g)title('g')输出：a=1234b=3456c=46810d=-2-2-2-2e=381524f=0.33330.50000.60000.6667-7-g=1162434096(2)用Matlab实现下列序列：a)程序如下：N=16n=0:15xn=0.8.^nstem(xn)输出：b)程序如下：N=16n=0:15-7-xn=exp((0.
4、2+3*j)*n)stem(xn)输出：c)程序如下：N=16n=0:15xn=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)stem(xn)输出：d)将c)中的x(n)扩展为以16为周期的函数，绘出四个周期。程序如下：N=64n=0:N-1xn=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)fori=1:4forj=1:16yn((i-1)*16+j)=xn(j)endend-7-stem(yn)grid输出：e)将c)中的x(n)扩展为以10为周期的函数，绘出四个周期。程序如下：N
5、=40n=0:N-1xn=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)fori=1:4forj=1:10yn((i-1)*10+j)=xn(j)endendstem(yn)grid输出：(3)绘出下列时间函数的图形，对x轴、y轴以及图形上方须加适当的标注。a)程序如下：t=0:0.01:10-7-xt=sin(2*pi*t)plot(t,xt)xlabel('t/s')ylabel('Amplitude')legend('x(t)=sin(2*pi*t')输出：b)程序如下：t=0:0.01:4xt=cos(100*pi*t).*s
6、in(pi*t)plot(t,xt)xlabel('t/s')ylabel('Amplitude')legend('x(t)=cos(100*pi*t)*sin(pi*t)')输出：(4)一因果系统,求出并绘制H（Z）的幅频响应与相频响应。程序如下：b=[1,sqrt(2)]-7-a=[1,-0.67,0.9][h,w]=freqz(b,a)am=20*log(abs(h))subplot(1,2,1)plot(w,am)title('Amplitude')ph=angle(h)subplot(1,2,2)plot(w,ph)title('Phase')输出：五.实验思考题1.在MA
7、TLAB命令窗里运行demo指令,运行其中样例程序,体会MATLAB强大功能。编辑MATLAB的库函数要点是什么?答：-7-
