信号与系统课程设计09级.pptVIP

二、信号基本运算的MATLAB实现*t=-3:0.001:3;ft1=tripuls(2*t,4,0.5);subplot(2,1,1)plot(t,ft1)title(‘x(2t))ft2=tripuls((2-2*t),4,0.5);subplot(2,1,2)plot(t,ft2)title(‘x(2-2t))1.信号的尺度变换、翻转、时移（平移）[例1]已知三角波x(t)，用MATLAB画出的x(2t)和x(2-2t)波形。*解：二、信号基本运算的MATLAB实现2.信号的相加与相乘相加用算术运算符“+”实现相乘用数组运算符“.*”实现[例2]画信号Aeatcos(w0t+f)的波形。解：t=0:0.001:8;A=1;a=-0.4;w0=2*pi;phi=0;ft1=A*exp(a*t).*sin(w0*t+phi);plot(t,ft1)*二、信号基本运算的MATLAB实现*3.离散序列的差分与求和连续信号的微分与积分差分 y=diff(f);求和 y=sum(f(k1:k2));微分y=diff(f)/h;h为数值计算所取时间间隔定积分quad(‘function_name’,a,b);function_name为被积函数名，a和b指定积分区间。二、信号基本运算的MATLAB实现*3.离散序列的差分与求和连续信号的微分与积分[例3]已知三角波x(t)，画出其微分与积分的波形。解：%differentiationh=0.001;t=-3:h:3;y1=diff(f2_2(t))*1/h;plot(t(1:length(t)-1),y1)%integrationt=-3:0.1:3;f2=@(t)tripuls((t),4,0.5);forx=1:length(t)y2(x)=quad(f2,-3,t(x));endplot(t,y2)三角波x(t)微分与积分的波形*利用MATLAB进行系统的时域分析*01连续时间系统零状态响应的求解02连续时间系统冲激响应和阶跃响应的求解03离散时间系统零状态响应的求解04离散时间系统单位脉冲响应的求解05离散卷积的计算y=lsim(sys,x,t)一、连续时间系统零状态响应的求解*t表示计算系统响应的抽样点向量；a=[a3,a2,a1,a0];b=[b3,b2,b1,b0];sys=tf(b,a)123654x是系统输入信号向量；sys是LTI系统模型，借助tf函数获得b和a分别为微分方程右端和左端各项的系数向量。sys=tf(b,a)二、连续系统冲激响应和阶跃响应求解*连续时间系统阶跃响应可用step函数直接求出，其调用形式为连续时间系统冲激响应可用impulse函数直接求出，其调用形式为y=impulse(sys,t)t表示计算系统响应的抽样点向量；123sys是LTI系统模型。y=step(sys,t)4可用MATLAB表示为三、离散时间系统零状态响应的求解*b,a分别是差分方程左、右端的系数向量；b=[b0,b1,b2,??,bM];a=[a0,a1,a2,??,aN];x表示输入序列；y表示输出序列。y=filter(b,a,x)0102030405四、离散时间系统单位脉冲响应的求解*b,a分别是差分方程左、右端的系数向量；k表示输出序列的取值范围；h就是单位脉冲响应。h=impz(b,a,k)五、离散卷积的计算*conv函数也可用于计算两个多项式的积a=[1,0,2,3];b=[1,3,2];c=conv(a,b)例如(s3+2s+3)(s2+3s+2)可用下面MATLAB语句求出。c=conv(a,b)式中a,b为待卷积两序列的向量表示，c是卷积结果。[例1]求系统y(t)+2y(t)+100y(t)=10x(t)的零状态响应，已知x(t)=sin(2pt)u(t)。解：%program3_1微分方程求解ts=0;te=5;dt=0.01;sys=tf([10],[12100]);t=ts:dt:te;x=sin(2*pi*t);y=lsim(sys,x,t);plot(t,y);xlabel(Time(sec))ylabel(y(t))*[例2]求系统y(t