陈后金信号与系第二版 第二章matlab.docxVIP

52页例2-13 %program2-8 t=-3:0.001:3; ft=tripuls(2*t,4,0.5); ft1=tripuls(2*t,4,0.5); subplot(2,1,1) plot(t,ft1) title(f(2t)) ft2=tripuls((2-2*t),4,0.5); subplot(2,1,2) plot(t,ft2) title(f(2-2*t)) 53页例2-14 %program2_9 the energy of exponential sequence k=0:10; A=1;a=-0.6; fk=A*a.^k; W=sum(abs(fk).^2) W = 1.5625 54页例2-15 function yt = f2_2(t) yt=tripuls(t,4,0.5); %program2_10 differentiation h=0.001;t=-3:h:3; y1=diff(f2_2(t))*1/h; plot(t(1:length(t)-1),y1) title(df(t)/dt) %program2_11 integration t=-3:0.1:3; for x =1:length(t) y2(x)=quad(f2_2,-3,t(x)); end plot(t,y2) title(integral of f(t)) 48页图2-2 %program2_1 decaying exponential signal A=1;a=-0.4; t=0:0.01:10; ft=A*exp(a*t); plot(t,ft) 48页图2-3 %program2_2 sinusoidal signal A=1; w0=2*pi; phi=pi/6; t=0:0.001:8; ft=A*sin(w0*t+phi); plot(t,ft) %program2_3sample function t=-3*pi:pi/100:3*pi; ft=sinc(t/pi); plot(t,ft)