信号与系统实验MATLAB实验代码.docx

1-1clear all;close all;t=0:0.01:10;A= 1;a=-0.4;ft=A*exp(a*t);plot(t,ft);1-2clear all;close all;t=-3:0.001:3;ft=tripuls(t,4,0.5);plot(t,ft);ft1=tripuls(t,4,1);figure,plot(t,ft1);1-3clear all;close all;t=0:0.001:4;T= 1;Ft=rectpuls(t-2*T,T);plot(t,Ft);1-4clear all;close all;k=-50:50;delta=zeros(1,50),1,zeros(1,50);stem(k,delta);1-5clear all;close all;k=-50:50;uk=zeros(1,50),ones(1,51);stem(k,uk);2-1clear all;close all;t=-3:0.001:3;ft=tripuls(t,4,0.5);subplot(3,1,1);plot(t,ft);title(f(t);ft1=tripuls(2*t),4,0.5);subplot(3,1,2);plot(t,ft1);title(f(2t);ft2=tripuls(2-2*t),4,0.5);subplot(3,1,3);plot(t,ft2);title(f(2-2t);2-2clear all;close all;h=0.001;t=-3:h:3;ft=tripuls(t,4,0.5);subplot(3,1,1);plot(t,ft);title(f(t);y1=diff(ft)*1/h;subplot(3,1,2);plot(t(1:length(t)-1),y1);title(the differentiation of f(t);3-1clear all;close all;t=-10:0.01:10;A= 1;a=-0.4;ft=A*exp(a*abs(t);plot(t,ft);3-2clear all;close all;r=0.02;t=-5:r:5;f1=1/2*exp(-2*t).*stepfun(t,0);f2=1/2*exp(-2*(t-1).*stepfun(t-1,0);ft1=diff(f1)*1/r;ft2=diff(f2)*1/r;subplot(3,1,1);plot(t,f1);grid on;xlabel(t);ylabel(f(t);title(f(t);subplot(3,1,2);plot(t(1:length(t)-1),ft1);title(diff1);subplot(3,1,3);plot(t(1:length(t)-1),ft2);title(diff2);4-1clear all;close all;syms t f;f= fourier(exp(-2)*abs(t);ezplot(f);4-2clear all;close all;syms t w;f=ifourier(1/(1+(w2),t);ezplot(f);a ns1/2*exp(-t)*heaviside(t)+1/2*exp(t)*heaviside(-t)4-3clear all;close all;r=0.02;t=-5:r:5;N=200;w=2*pi;k=-N:N;w=k*w/N;f1=1/2*exp(-2*t).*stepfun(t,0);F=r*f1*exp(-j*t*w);F1=abs(F);P1=angle(F);subplot(3,1,1);plot(t,f1);grid on;xlabel(t);ylabel(f(t);title(f(t);subplot(3,1,2);plot(w,F1);xlabel(w);grid on;ylabel(F(jw);subplot(3,1,3);plot(w,P1*180/pi);grid;xlabel(w);ylabel(W);4-4clear all;close all;R=0.02;t=-2:R:2;f=stepfun(t,-1)-stepfun(t,1);f1=f.*exp(-j*5*t);f2=f.*exp(j*5*t);W1=pi*5;N=500;k=-N:N;W=k*W1/N;F1=f1*exp(-j*t*W)*R;F2=f2*exp(-j*t*W)*R;F1=real(F1);F2=real(F2);subplot(2,1,1);plot(W,F1);xla