matlab的综合运用

实验四 matlab 的综合运用一、实验目的(1)加深对 MATLAB 的理解和应用。(2)掌握利用MATLAB求解信号与系统中的问题。二、实验环境安装MATLAB7.0以上版本的计算机三、实验内容与结论：1. 某离散时间 LTI 系统满足差分方程yk + 0.7 yk -1 - 0.45 yk - 2 - 0.6 yk - 3二 0.8xk - 0.44 xk -1+0.36 xk - 2 + 0.22xk 试利用impz函数求其单位脉冲响应，并画出前30点的图。k=0:30;a=1 0.7 -0.45 -0.6;b=0.8 -0.44 0.36 0.22;h=impz(b,a,k);stem(k,h);2. 某RC电路如图，试求：(1) 对不同RC值，画出该系统的幅度响应曲线(2) 信号x(t)=cos100t+cos3000t。试确定适当的RC值滤除信号中的高频分量，并画出 信号x(t)和滤波后的信号y(t)的波形。(3) 50Hz的交流信号经过全波整流后可表示为x(t)=10|sin100nt|,试取不同的RC值, 计算并画出x (t)过该系统的响应y (t),并计算两信号的直流分量。w=linspace(0,2*pi,200);b=1;a=0.3 1;H=freqs(b,a,w);plot(w,H);RC=0.3;t=linspace(-2,2,1024);w1=100;w2=3000;H1=1/(RC*j*w1+1);H2=1/(RC*j*w2+1);x=cos(100*t)+cos(3000*t);y=abs(H1)*cos(w1*t+angle(H1)+abs(H2)*cos(w2+angle(H2); subplot(2,1,1);plot(t,x);ylabel(x(t);xlabel(时间(秒);subplot(2,1,2);plot(t,y);ylabel(y(t);xlabel(时间(秒);时间（秒）RC=0.001;t=linspace(-2,2,100);H1=1/(RC*j*100*pi+1);x=10*abs(sin(100*pi*t);y=10*abs(H1)*abs(sin(100*pi*t+angle(H1); subplot(2,1,1);plot(t,x);ylabel(x(t);xlabel(时间(秒);subplot(2,1,2);plot(t,y);ylabel(y(t);xlabel( 时间(秒) );RC=0.04(3)RC=0.031050-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 时间(秒)1050-2-1.5 -1 -0.5 0 0.5 1 1.5 2时间（秒）1.510.50 0.5 1 1.5 2时间（秒）0-2-1.5 -1 -0.50.80.60.40.20-2 -1.5 -1 -0.5 0 0.5 1 1.5 2时间(秒)RC=0.0013已知H (s)二，画出该系统的零极点分布图，判断系统的稳定性，求出系统的冲激响应，阶跃响应和频率响应。零点分布图从零极点分布图可知系统的极点位于左半平面，故系统稳定。 冲激响应num=12;den=1 2 2 1;sys=tf(num,den);figure(1);pzmap(sys);t=1:0.02:10;w=0:0.02:5;h=impulse(num,den,t);figure(2);plot(t,h);xlabel(time(s);title(冲激响应);H=freqs(num,den,w);figure(3);plot(w,H);xlabel(frequencyomega);title( 频率响应);g = tf(1 2,1 2 2 1); figure(1); tep(g,10); title(阶跃响应)