广石化信号与系统实验2

《广石化信号与系统实验2》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验代码与图片(一)代码：t0=-5;t仁5;dt=0.01;t=tO:dt:t1;f1=exp(-2*t).*heaviside(t);f2=heaviside(t)-heaviside(t-4);y=dt*conv(f1,f2);%卷积subplot(2,2,1);plot(t,f1,'k1,'linewidth',2);title(,eA(-2t)[u(t)],);xlabel(,x轴');ylabel('y轴TRotatiorTQ);set(gca,'xtick：[・5:1:5]);gridon;su

2、bplot(2,2,2);plot(tJ2,,k7linewidth,,2);title(l(t)・u(t・4));xlabel「x轴^ylabelfy轴TRotatiorTQ);set(gca,xtick：卜5:1:5]);gridon;t=2*t0:dt:2*t1;%重新定义卷积的范围，且要等于两个被卷积函数的范围的两倍。subplot(212);plot(t,y,k」inewidth：2);titleCy=f1*f2');xlabelCx轴^ylabelfy^•/Rotation*,0);gridon;

3、eL2t)[u(t)]u(t)-u(t-4)5-4-3-2-1012345-5-4-3-2-1012345x轴x轴%使用Isim函数来获得零状态响应。t=(0:0.01:10);a=[134];%微分方程的左端系数向量。b=[013];%微分方程的右端系数向量。sys=tf(b,a);%LTI系统模型f=exp(-2*t).*heaviside(t);y=lsim(sys,f,t);%f是紊统输入函数，t是抽样点向量plotfty/k','linewidth',2);gridon;(三)代码：%使用impul

4、se()还有step()函数来计算冲激函数和阶跃函数t=0:0.01:10;a=[218];b=[001];y=impulse(b,a,t);subplot(2X);plot(t,y,kTlinewidth：2);title(‘冲激响应')；y1=step(b,a,t);subplot(212);plot(t,y/klinewidth,,2);titlef阶跃响应')；冲激响应0.311111-0.2012345678910(四)代码：%验证x(t)*delta(t)=x(t)t=-5:0.01:5;x=s

5、in(2*t);subplot(221);plot(t,x,'k','linewidth1,2);title(,x(t),)；set(gca/xtick,,[-10:1:10]);delta=(heaviside(t)-heaviside(t-0.01))/0.02;subplot(222);plot(t,delta,'linewidth',2);axis([-5501]);title('delta')；%set(gca；xtick或者ytick：[差坐标值：刻度：右坐标值]),这个函数用于设置坐标轴的刻度。

6、y=0.01*conv(x,delta);t=/0:0・01:10;subplot(212);plot(t,y,kTlinewidth：2);set(gca,'xtickH0:1:10]);axis([-55-11]);title(y);效果图：5-4-3-2-101234510.80.60.4050.2°5可见验证了x(t)*6(t)=x(t)o(五)length()函数即数组长度，b=length(A)作用是：如果A为非空数组，返回行数和列数两者之间数值较大的那一个值，即相当于执行了max(size(A)

7、);如果A为空数组，则返回0;如果A是一个向量则返回A的长度。调用格式是n=length(M)。ones()函数的作用是产生mxn的全1矩阵。如果调用方式为ones(n),贝

8、」表示产生nm的方阵。调用格式是ones(m,n)o(六)代码：%首先构建一个LC低通滤波的传递函数%H(s)=(2500)/((sA2)+70.7s+2500)%%对应的微分方程是：yH(t)+70.7y,(t)+2500y(t)=2500f(t)-—%t=(0:0.01:5);a=[170.72500];b=[002500];sys

9、=tf(b,a);%LTI系统模型f=sin(2*pi*t)+0.1*sin(50*pi*t);%输入函数y=lsim(sys,f,t);%零状态响应函数,即输出subplot(211);plot(t,f,'k','linewidth',2);gridon;title('f=sin(2*pi*t)+0.1*sin(50*pi*t)');subplot(212);plot(t,y,'k','linew