资源描述:
《matlab的综合运用》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实验四matlab的综合运用一、实验目的(1)加深对MATLAB的理解和应用。(2)掌握利用MATLAB求解信号与系统中的问题。二、实验环境安装MATLAB7.0以上版本的计算机三、实验内容与结论:1.某离散时间LT1系统满足差分方程y[幻+O.ly[k-1J-0A5y[k-2J-0.6)仏—3]=O.Sx[k]-O.44x[k-1]+036x[k-2]+Q.22x[k]试利用impz函数求其单位脉冲响应,并画出前30点的图。><>-0.5-11O15202530k=0:30;a=[l0.7-0.45-0.6];b=[0.8-0.440.3
2、60.22];h=impz(bza,k);stem(k,h);2.某RC电路如图,试求:(1)对不同RC值,画出该系统的幅度响应曲线(2)信号x(t)=cosl00t+cos3000to试确定适当的RC值滤除信号中的高频分量,并画出信号x(t)和滤波后的信号y(t)的波形。(3)50Hz的交流信号经过全波整流后可表示为x(t)=10
3、sinl00nt
4、,试取不同的RC值,计算并画出x(t)过该系统的响应y(t),并计算两信号的直流分量。w=linspace(0,2*pi,200);b=[l];a=[0.31];H=freqs(bza,w);pl
5、ot(w*H);RC=0・3;t=linspace(-2,2,1024);wl=100;w2=3000;Hl=l/(RC*j*wl+l);H2=l/(RC*j*w2+l);x=cos(100*t)+cos(3000壮);y=abs(Hl)*cos(wl*t+angle(Hl))+abs(H2)*cos(w2+angle(H2));subplot(2,1,1);plot(tzx);ylabel('x(t)');xlabel('时间(秒)*);subplot(2,1,2);plot(t,y);ylabel('y(t)');xlabel(,时间(秒)
6、T);ROO.001;t=linspace(-2,2,100);Hl=l/(RC*j*100*pi+l);x=10*abs(sin(100*pi*t));y=10*abs(Hl)*abs(sin(100*pi*t+angle(Hl)));subplot(2,1,1);plot(tzx);ylabel(*x(t)*);xlabel(*时间(秒)f;subplot(2,1,2);plot(t,y);ylabel(*y(t)');xlabel(,时间(秒)1);RC=0.04⑶RC=0.03RC=O.OO11O-0.500.5时血J(秒》-2-1.5
7、-0.500.5时问(秒》1.53.s+2s'+2厂+2s+1画出该系统的零极点分如图,判断系统的稳定性,求岀系统的冲激响应,阶跃响应和频率响应。零点分布图FAz>lo-J£oroO.«0.00.40.20.2C.4o.e-o.e--i.e-o.e-o.*i-0.2从零极点分布图可知系统的极点位于左半平面,故系统稳定。冲激响应2.5卩介跃
8、轲应num=[12];den=[l221];sys=tf(num,den);figure(1);pzmap(sys);t=l:0•02:10;w=0:0.02:5;h=impulse(num,den,t);f
9、igure(2);plot(t,h);xlabel('time(s)f;titled冲激响应H=freqs(num,den,w);figure(3);plot(wzH);xlabel(1frequency'omega1);titled频率响应J;g=tf([12],[1221]);figure(1);tep(g,10);title(1阶跃响应,