实验七 MATLAB在信号处理中的应用(一).docVIP

Matlab在电子信息类中的应用实验 班级： 成绩： 姓名： 学号： 实验七 MATLAB在信号处理中的应用(一) 实验目的: 1. 掌握信号及其表示方法、信号的基本运算； 2. 掌握线性时不变系统的响应 实验内容: 用MATLAB命令绘出连续时间信号xt=e-0.707tsin(t)*cos(t) 关于t的曲线，t的范围为0~30s,并以0.1s递增。请给图形加上标题， 编写脚本产生三角波函数，方波函数，sinc函数，以及非周期方波和非周期三角波函数，并作图。 P182第二、五、六题； 有二阶系统： 当输入周期为 4s的方波时的输出响应。 有二阶系统： 求系统对100点随机噪声的输出响应曲线。 实验结果： 1. t=0:0.1:30; x=exp(-0.707*t).*sin(t).*cos(t); plot(t,x); grid; xlabel(time(sec)); ylabel(x(t)) 2. Fs=10000;t=0:1/Fs:1; x1=sawtooth(2*pi*50*t,0); x2=square(2*pi*50*t,80); subplot(2,1,1),plot(t,x1),axis([0,0.2,-1,1]); subplot(2,1,2),plot(t,x2),axis([0,0.2,-1.5,1.5]); x=linspace(-4,4); y=sinc(x); plot(x,y) t=0:0.01:1; y1=rectpuls(t);y2=rectpuls(t,0.7); y3=tripuls(t);y4=tripuls(t,0.7); subplot(4,1,1);plot(t,y1),grid; subplot(4,1,2);plot(t,y2),grid; subplot(4,1,3);plot(t,y3),grid; subplot(4,1,4);plot(t,y4),grid; P182. 2.. clc; clear all; clf; m = 10; e = 0:19; c = 0.6; k = 5; a = k*c.^e; a =a; b = circshift(a,m); L = length(a)-1; n = 0:L; subplot(2,1,1); stem(n,a); axis([0,L,min(a),max(a)]); title(?-DòáD); subplot(2,1,2); stem(n,b); axis([0,L,min(a),max(a)]); title([?2?üò???,num2str(m),???ù±?DòáD]); 5.先把分子分母同乘以z^2等价 num=[2,3,0]; den=[1,0.4,1]; disp(?μí3′?μYoˉêyH(z)); printsys(num,den,z); disp(×a?aá???μ???ò??￡Dí); [z1,p1,k1] = tf2zp(num,den) 系统传递函数H(z) num/den = 2 z^2 + 3 z --------------- z^2 + 0.4 z + 1 转为零极点增益模型 z1 = 0 -1.5000 p1 = -0.2000 + 0.9798i -0.2000 - 0.9798i k1 = 2 本题直接调用了一个化因式的指令z’sys，是我在网上查的，还有这道题用了上道题的零极点来画图。以及pzmap指令。 程序：num=[4,15.6,6,2.4,-6.4]; den=[3,3.4,6.3,-11.4,6]; disp(?μí3′?μYoˉêyH(z)); g=tf(num,den); printsys(num,den,z); disp(×a?aòòê?D?ê?); zsys=zpk(num./den) [z1,p1,k1] = tf2zp(num,den) g1=zpk(g); [z1,p1,k1] = tf2zp(num,den); pzmap(g1) grid on 系统传递函数H(z) num/den = 4 z^4 + 15.6 z^3 + 6 z^2 + 2.4 z - 6.4 -------------------------------------- 3 z^4 + 3.4 z^3 + 6.3 z^2 - 11.4 z + 6 转为因式形