实验二程序示例.doc

《实验二程序示例.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、四、实验内容1、编写程序Q2_1，完成与两函数的卷积运算。其中：在一个图形窗口中，画出、以及卷积结果。要求每个坐标系有标题、坐标轴名称。clear;closeall;t0=-2;t1=4;dt=0.01;t=t0:dt:t1;f1=exp(-2*t).*ut(t);f2=ut(t)-ut(t-4);f=dt*conv(f1,f2);subplot(2,2,1)plot(t,f1),gridon,title('Signalf1(t)'),xlabel('Timetsec'),axis([t0,t1,-0.2,1.2])subplot(2,2,2

2、)plot(t,f2),gridon,title('Signalf2(t)'),axis([t0,t1,-0.2,1.2])subplot(2,1,2)t=2*t0:dt:2*t1;plot(t,f),gridon,title('Theconvolutionoff1(t)andf2(t)'),axis([2*t0,2*t1,-0.1,0.6]),xlabel('Timetsec')2、编写程序Q2_2。若系统模型为：其中求零状态响应，画出波形。并与手工计算的零状态响应比较。ts=0;te=5;dt=0.01; b=[1,3];a=[1,3,4

3、]; sys=tf(b,a);  t=ts:dt:te;       f=exp(-t).*ut(t);    y=lsim(sys,f,t);      plot(t,y);       xlabel('Time(sec)');     ylabel('y(t)'); 3、编写程序Q2_3。已知某连续系统的微分方程为：，绘出该系统的冲激响应和阶跃响应的波形。a=[2,1,8];  b=[1];t=0:0.1:10; y=impulse(b, a, t);   subplot(2,1,1) plot(t,y)         title('单

4、位冲激响应') y=step(b, a, t);   subplot(2,1,2) plot(t,y)     title('单位阶跃响应 ') xlabel('t') 4、编写程序Q2_4。验证卷积的相关性质。验证性质：选择信号x(t)的数学表达式为：x(t)、δ(t)和x(t)*δ(t)的波形t=[0:0.01:3]; x=exp(-2*t).*ut(t); n=length(t); chongji=zeros(1,n); chongji(1,1/0.01+1)=1/0.01; subplot(2,2,1); y=0.01*conv(x,

5、chongji); plot(t,x); xlabel('t'); ylabel('x(t)'); subplot(2,2,2); plot(t,chongji); xlabel('t'); ylabel('chongji()'); t=0:0.01:6; subplot(2,1,2);plot(t,y); xlabel('t'); ylabel('x(t)*chongji()'); 5、length()、ones()两条指令的作用是什么？说明其格式。Length():是求出向量数组包合含有多少个元素，例如：t=0:0.1:5;n=lengt

6、h(t);n;就输出n=51。    Ones():作用是将向量数组元素初始化为1，例如：ones(1,n)就是就会输出一个一行n列的数组，而里面的元素全部是1。