用MATLAB实现常用的离散时间信号及其时域运算

1、用MATLAB实现常用的离散时间信号及其时域运算摘要:1. 运用MATLAB勺绘图指令绘制离散时间信号；2. 用MATLA实现离散时间信号的时域运算。关键词：MATLAB 离散时间信号时域运算实验目的：1. 复习MATLA语言以及常用指令2. 学习和掌握用MATLA语言产生离散时间信号的编程方法3. 通过编程绘制出离散时间信号波形，加深理解信号的时域运算实验原理：1.1单位阶跃序列和单位样值序列离散时间信号只在某些离散瞬间时给出信号的值，因此，它是时间上不连续的序列。单位阶跃序列和单位样值序列在离散信号与系统的分析中是两个非常典型的序列，分别记为u(n)和3(n),他们的定义分别如下：0 (n

2、 0)u(n)=U (n 工 0)(n)二1 (n = 0)0 (n = 0)若单位节约序列的起始点为 no,单位样值序列出现在no时刻，则表达式为:u(n - no)0 (n no)1 (n 一 no)1 (n = no)(n - no)二io (n 式 no)2离散时间信号的时域运算与连续时间信号的研究类似，在离散系统分析中，经常遇到离散时间的运算，包括 两信号的相加、相乘以及序列本身的移位、反褶、尺度等等，也需要了解在运算过程 中序列的以及相应的波形的变化。序列x(n)的反褶:x( _ n)序列x(n)的移位：x( n - no)序列x(n)的尺度：x ( a n )两序列x(n)和xn

3、)相乘:x(n) x£n)两序列x(n)和xn)相乘:x(n) x£n)两序列x(n)和x(n)相加减:x(n) _ X2(n)两序列x(n)和xn)相乘:x(n) x£n)试验任务:1绘制用于产生以下信号的通用程序1 x( n)二 A ( nn o)2 x( n) = Au( n - mi) - u( n - m2)1.1程序和图形见附录11.2程序和图形见附录22已知波形，求表达式并绘出各信号波形f ( k - 2)u( k)f(k2)f ( k - 2)u( k - 2)2.程序和图形见附录3附录1:% n1:起始时刻；n0:跳变时刻；n2 :终止时刻;t1

4、=i nput('请输入起始时刻n1=');t0=i nput('请输入跳变时刻n0=');t2=i nput('请输入终止时刻n2=');A=i nput('请输入幅度A=');t=t1:t2; n=len gth(t); tt=t1:t0;n 1=le ngth(tt); x=zeros(1, n); x(n 1)=A*1; stem(t,x),gridoncleartitle('x(n)')axis(t1 t2 -0.2 A+1)n1=-2, n0=1, n2=3, A=2附录2:Clear%n1:起始时刻；

5、m1第一次跳变时刻；m2:第二次跳变时刻；n2 :终止时刻;t1=input('请输入起始时刻n1=');t0=in put('请输入第一次跳变时刻m仁');t3=input('请输入第二次跳变时刻m2=');t2=input('请输入终止时刻n2=');A=input('请输入幅度A=');t=t1:t2;n=len gth(t);tt=t0:t2;n 1=le ngth(tt);x=zeros(1, n-n 1),o nes(1, n1);ttt=t3:t2;n 2=le ngth(ttt);y=zeros(1

6、, n-n 2),o nes(1, n2);z=A*(x-y);stem(t,z),gridontitle( 'x ( n)')axis(t1, t2 ,min( z)-1, max(z)+1)n1=-3， n2=3， m仁-1，m2=1，A=2附录3:1. 离散序列反褶函数(Isfz.m)fun ctionf,k=lsfz(f1,k1)f=fliplr(f1);k=-fliplr(k1);2. 离散序列移位函数(lspy.m)fun cti onf,k=lspy(ff,kk,kO)k=kk+kO;f=ff;3. 离散序列相乘函数(lsxc.m)fun ctionf,k=lsx

7、c(f1,f2,k1,k2)k=mi n(min (k1),mi n(k2):max(max(k1),max(k2); s仁 zeros(1,le ngth(k);s2=s1;s1(fi nd(k>=mi n(k1)&(k<=max(k1)=1)=f1;s2(fi nd(k>=mi n( k2)&(k<=max(k2)=1)=f2; f=s1.*s2;Retur n4.主程序:'f(k)')'f(k-2)u( n)')'f(-k+2)')cleari=1;for m=-3:4if m<0f(i)=m+

8、3;elseif (m=4)f(i)=0；elsef(i)=3；endi=i+1;end%f(k)m=-3:4subplot(221),stem(m,f,'filled' ),grid on ,title(axis(mi n( m)-1,max(m)+1,mi n(f)-0.5,max(f)+1)i=1;for n=-4:6if n<0u(i)=0;elseu(i)=1;endi=i+1;endn=-4:6m=-3:4%f(k-2)y1,k1=lspy(f,m,2);%f(-k)y2,k2=lsfz(f,m);%u(k-2)y3,k3=lspy(u, n,2);%f(k-2)u(k)z1,k=lsxc(y1,u,k1, n)subplot(222),ste m( k,z1,'filled' ),grid on ,title(axis(mi n(k)-1,max(k)+1,mi n(z1)-0.5,max(z1)+1)%f(-k+2)z2,w2=lspy(y2,k2,2)subplot(223),stem(w2,z2,'filled' ),grid on ,title(axis(mi n(w2)-1,max(w2)+1,mi n(z2)-1,max(z2)+1)%f(k-2)u(k-2)z3,w3=lsxc(y1,y