数字信号处理实验4报告.docVIP

广东海洋大学学生实验报告书（学生用表） 实验名称 课程名称 课程号 学院(系) 专业 班级 学生姓名 学号 实验地点 实验日期 11.4离散时间信号的频域分析 实验目的 学会运用MATLAB求离散时间的信号的z变换和z反变换。 学会运用MATLAB求信号的离散时间傅立叶变换。 实验内容 z变换和z反变换 (1)用ztrans函数求函数x(n)=a^ncos(pi*n)u(n)的z变换。 x=sym(a^n*cos(pi*n)); z1=ztrans(x); z=simplify(z1) 实验结果：z =z/(z+a) (2)用iztrans函数求函数X(z)=z(2z^2-11z+12)/[(z-1)(z-2)]的z变换。 z=sym(z*(2*z^2-11*z+12)/(z-1)/(z-2)^3); x=iztrans(z); simplify(x) 实验结果：ans = -3+3*2^n-1/4*2^n*n-1/4*2^n*n^2 (3)用MATLAB命令对函数X(z)=18/[18+3z^(-1)-4z^(-2)-z^(-3)]进行部分分式展开，并求其z反变换。 B=[18]; A=[18,3,-4,-1]; [R,P,K]=residuez(B,A) 实验结果：R = 0.3600 0.2400 0.4000 P = 0.5000 -0.3333 -0.3333 K = [] 序列的离散时间傅立叶变换 n=-10:10; x=(-0.8).^n; k=-200:200; w=(pi/100)*k; X=x*(exp(-j*pi/100)).^(n*k); subplot(4,1,1) plot(w/pi,real(X));grid title(X(e^{j\omega})实部) xlabel(\omega/\pi); ylabel(振幅); subplot(4,1,2) plot(w/pi,imag(X));grid title(X(e^{j\omega})虚部) xlabel(\omega/\pi); ylabel(振幅); subplot(4,1,3) plot(w/pi,abs(X));grid title(X(e^{j\omega})幅度谱) xlabel(\omega/\pi); ylabel(振幅); subplot(4,1,4) plot(w/pi,angle(X));grid title(相位谱arg[X(e^{j\omega})]) xlabel(\omega/\pi); ylabel(以弧度为单位的相位); 实验总结： 通过这次编程实验，感觉对MATLAB语句还不是很熟练，很多处语句都不怎么明白为什么要这样写，但还是收获蛮大的，不仅加深了对数字信号处理的离散时间信号的理解，还锻炼编程能力，感觉时间没白白浪费。 成绩 指导教师 日期 注:请用A4纸书写，不够另附纸。 第　　　页，共　　　页 GDOU-B-11-112