时域与频域的matlab程序

1 一．典型连续信号和离散信号的时域波形。 1．单边指数信号)()(tuAetyt； 2．单位冲激信号)()(0ttty ； 3．单位阶跃信号)()(0ttuty； 4．矩形脉冲信号)]()([)(21ttuttuAty； 5．正弦信号)()sin()(tutAty； 6．单位序列)()(0nnny ； 7．单位阶跃序列)()(0nnuny； 8．单位矩形序列)()()(21nnunnuny； 9．指数序列)()(nuaAnyn； 10．正弦序列)()sin()(nunAny。 2 单边指数信号 function zhishu(A,a,t1,t2,dt) t1=0 t2=10 A=1 A=-0.4 dt=0.01 t=t1:dt:t2; y=A*exp(a*t); plot(t,y) axis([t1,t2,0,1.2]) xlabel('t') ylabel('y(t)') title(' 单边指数信号') 单位冲激信号 function chongji(t1,t2,t0) dt=0.01; t1=10; t2=-5; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(-t0-t1)/dt+1)=1/dt; stairs(t,x); axis([t1,t2,0,1.2/dt]) xlabel('t') ylabel('y(t)') title('单位冲激信号') 3 单位阶跃信号 function jieyao(t1,t2,t0) t1=0;t2=10;t0=-4 t=t1:0.01:-t0; tt=-t0:0.01:t2; n=length(t); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); plot(tt,uu) hold on plot(t,u) plot([-t0,-t0],[0,1]) hold off title('单位阶跃信号 y(t)') axis([t1,t2,-0.2,1.5]) 矩形脉冲信号 function jxmcxh(A,width,T1,T2,dt,T0) A=3;width=2; 4 T1=-3;T2=3; T0=0;dt=0.01 t=T1:dt:T2; ft=A*rectpuls(t-T0,width); plot(t,ft); xlabel('t') ylabel('y(t)') title('矩形脉冲信号') axis([t1,t2,0,4]); 正弦信号 function zhengxian(A,w,t1,t2,dt) A=5;w=0.5*pi;t1=0;t2=15;dt=0.01 t=t1:dt:t2; f=A*sin(w*t); plot(t,f) title('正弦信号') xlabel('t') ylabel('y(t)') 5 单位序列 function dwxulie(k1,k2,k0) k1=-8;k2=12;k0=-2; k=k1:k2; n=length(k); f=zeros(1,n); f(1,-k0-k1+1)=1; stem(k,f,'filled') axis([k1,k2,0,1.5]) title('单位冲序列') 单位阶跃序列 function jyxulie(k1,k2,k0) k1=-10;k2=10;k0=4; k=k1:-k0-1; kk=-k0:k2; n=length(k); nn=length(kk); u=zeros(1,n); 6 uu=ones(1,nn); stem(kk,uu,'filled') hold on stem(k,u,'filled') hold off title('单位阶跃序列') axis([k1,k2,0,1.5]) 单位矩形序列 function jyxulie(k1,k2,k0) k1=-8;k2=12;k0=1; axis([k1,k2,0,1.5]); k=k1:-k0-1; kk=-k0:6; kkk=7:k2 n=length(k); nn=length(kk); nnn=length(kkk); u=zeros(1,n); uu=ones(1,nn); uuu=zeros(1,nnn); stem(kk,uu,'filled') hold on stem(k,u,'filled') stem(kkk,uuu,'filled') hold off title('单位矩形序列') 7 指数序列 function dszsu(c,a,k1,k2) %c: 指数序列的幅度 %a: 指数序列的底数 %k1: 绘制序列的起始序号 %k2: 绘制序列的终止序号 c=1;a=2;k1=-2;k2=10; k=k1:k2; x=c*(a.^k); stem(k,x,'filled') hold on plot([k1,k2],[0,0]) hold off title('指数序列') xlabel('n') ylabel('f(n)') 8 正弦序列 function zxxulie(A,w,k1,k2) k1=-30;k2=30;a=2;w=0.25 k=k1:k2; stem(k,A*sin(k*w),'filled') title('离散时间正弦序列f(n)=Asin(wn)') xlabel('n') ylabel('f(n)')