信号与系统实验2-9.docVIP

信号与系统 ***************************注意*************************** （1）这些程序是基于MATLAB7.0编写，在2013.a中会有不兼容的现象程序之前建议清空（）命令窗口 ********************************************************* 2.6-1 %绘制单位阶跃信号 t=-5:0.01:5; y1=(t0); subplot(2,2,1); plot(t,y1); axis([-5,5,-0.5,1.5]); xlabel(t); ylabel(u(t)); title(单位阶跃信号) %绘制指数信号 t=0:0.01:10; A=1; a=-0.4; y2=A*exp(a*t); subplot(2,2,2); plot(t,y2); axis([0,10,-0.5,1.5]); xlabel(t); ylabel(exp(-0.4t)); title(指数信号) %绘制正弦信号 A=5; w=0.5*pi; t=0:0.01:16; y3=A*exp(a*t); subplot(2,2,3); plot(t,y3); title(正弦信号) xlabel(t); ylabel(5sin(Omegat)); axis([0,16,-5,5]); line([0,16],[0,0]); %绘制抽样信号 t=-15:0.01:15; t1=t/pi; y4=sinc(t1); subplot(2,2,4); plot(t,y4); title(抽样信号) xlabel(t); ylabel(sinc(t)); axis([-15,15,-0.3,1.1]); line([-15,15],[0,0]); 2.6-2 t=-6:0.01:6; u1=0.5+0.5.*sign(t); u2=0.5+0.5.*sign(t-4); f1=0.5.*t.*(u1-u2); subplot(2,2,1); %figure(1); plot(t,f1); title(f1(t)=0.5t[u(t)-u(t-4)]); xlabel(t); axis([-1,6,-0.2,2.2]); f2=sin(4*pi*t) subplot(2,2,2); %figure(2); plot(t,f2); title(f2(t)=sin(4*pi*t)); xlabel(t); axis([-6,6,-1.5,1.5]) f3=f1+f2 subplot(2,2,3); %figure(3); plot(t,f3); title(f1(t)+f2(t)); xlabel(t); axis([-6,6,-1.5,3]) f4=f1.*f2 subplot(2,2,4); %figure(4); plot(t,f4); title(f1(t)*f2(t)); xlabel(t); axis([-6,6,-2,2]); 2.6-3 先新建一个M文件，写一个Heaviside(t)子函数 function f=Heaviside(t) f=(t0); 再在命令窗口输入指令syms t; f=sym(Heaviside(t+1)-t*Heaviside(t)+(t-1)*Heaviside(t-1)); f1=2*subs(f,t,0.5*t); subplot(2,2,1); ezplot(f1,[-3,3]); title(f1(t)=2f(0.5t)); f2=subs(f,t,-2*t+3); subplot(2,2,2); ezplot(f2,[0,3]); title(f2(t)=f(-2t+3)); f3=sym(-Heaviside(t)+Heaviside(t-1)); subplot(2,2,3); ezplot(f3,[-2,2]); line([-1,-1],[0,1]); axis([-2,2,-1.5,1.5]); title(f3(t)=df(t)/dt); f4=int(f); subplot(2,2,4); ezplot(f4,[-2,3]); title(f4(t)=∫f(τ)dτ); 3.10-1 n=-10:10; w1=0.4*pi; n1=-10:-1; ft1=sin(0.2*pi*n1)./(pi*n1); n2=1:10;ft2=sin(0.2*pi*n2)./(pi*n2); ft=[ft1,0.2,ft2]; fn=abs(ft); phase=angle(ft); subplot(2,2,1); stem(n,fn); title(幅度谱); subplot(2,2,2); stem(n,