时域与频域的matlab程序

一．典型连续信号和离散信号的时域波形。

1．单边指数信号)()(t u Ae t y t α=；

2．单位冲激信号)()(0t t t y +=δ；

3．单位阶跃信号)()(0t t u t y +=；

4．矩形脉冲信号)]()([)(21t t u t t u A t y +-+?=；

5．正弦信号)()sin()(t u t A t y ω?=；

6．单位序列)()(0n n n y +=δ；

7．单位阶跃序列)()(0n n u n y +=；

8．单位矩形序列)()()(21n n u n n u n y +-+=；

9．指数序列)()(n u a A n y n ?=；

10．正弦序列)()sin()(n u n A n y ω?=。

单边指数信号

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('单位冲激信号')

单位阶跃信号

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;

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)')

单位序列

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);

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('单位矩形序列')

指数序列

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)')

正弦序列

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)')