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

用MATLAB实现常用离散时间信号及其时域运算 用MATLAB实现常用离散时间信号及其时域运算 摘要:本实验实现了运用MATLAB的绘图指令绘制离散时间信号，用MATLAB语言实现离散时间信号的时域运算。达到了熟悉MATLAB基本使用方法以及了解典型信号的MATLAB表示方法，掌握MATLAB函数的应用的目的。 关键词: 离散时间信号 时域运算 MATLAB 正文: 离散时间信号在时间上是离散的，只在某些不连续的规定瞬间给出函数值，在其他时间没有定义 典型离散时间信号: 1n,0单位样值信号: , ,(n),, 0n,0, 单位阶跃序列: ,1nn,0 ,,u(nn),0 ,0nn0, 矩形序列: 1(0,n,N,1), R(n),,N 0(n,0orn,N), 离散时间信号的时域运算 信号的时域运算包括信号的平移、反褶、尺度及信号的相加、相乘.用MATLAB 的符号运算命令来表示离散时间信号的时域运算，然后用ezplot命令绘制出其结 果波形图。 相加和相乘: Function[f,k]=lsxj(f1,f2,k1,k2) K=min(min(k1),min(k2)):max(max(k1),max(k2)); S1=zeros(1,,length(k));s2=s1; S1(find((k=min(k1))(k=max(k1))==1))=f1; S2(find((k=min(k2))(k=max(k2))==1))=f2; F=s1+s2; Stem(k,f,’filled’) Axis([(min(min(k1),min(k2))+1)，min(f)-0.5),(max(f)+0.5)]) 两序列相乘可参照此例 平移 Function[f,k]=lspy(ff,kk,k0) K=kk+k0 F=ff 反褶 Function[f,k]=lsfz(f1,k1) F=fliplr(f1); K=-fliplr(k1) 任务 1. 编制用于产生下列信号的通用程序，要求对于任意给定的参数都能实现所要求信号，调试并运行这些通用程序，具体产生有指导教师指定的信号并绘制信号的波形。 (1)x(n)=Aδ(n-n0) (2) x(n)=A[u(n-m1)-u(n-m2)] 2已知离散序列波形如图所示，实用MATLAB绘出满足下列要求的序列波形 (1)f(k-2)u(k) (2)f(-k+2) (3)f(k-2)u(k-2) 实验结论 MATLAB能对信号进行时域分析，将信号随时间的变化规律用二维曲线表示出来就可 以描述，而且MATLAB具有强大的图形处理功能及符号运算功能，为我们直观的看 到信号的模型提供了方便。 附录: 任务一 (1)x(n)=Aδ(n-n0) clear t1=input(t1=); t0=input(n0=); t2=input(t2=); A=input(A=); t=t1:t2; n=length(t); tt=t1:t0; n1=length(tt); f=zeros(1,n); f(t0-1)=A; stem(t,f,filled),grid on title(????DòáD X(n)=A(n-n0)) (2) x(n)=A[u(n-m1)-u(n-m2)] clear m1=input(μúò???ì???μ?m1=); m2=input(μú?t??ì???μ?m2=); A=input(ê?è??ù?è?μA=); if m1m2 t=m1:1:m2-1; x=ones(1,m2-m1); else t=m2:1:m1-1; x=-ones(1,m1-m2); end stem(t,A.*x,filled,r),grid on title( X(n)=A[u(n-m1)-u(n-m2)]) 任务二 clear,i=1; for k=-3:3 if k0 f(i)=k+3; else f(i)=3; end i=i+1; f(i)=0; end k=-3:4; subplot(221),stem(k,f,filed),grid on,title(f(k)) axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]) t=-3:1:6; u=[zeros(1,3),ones(1,7)]; f1=f; k1=k+2; h=min(min(k1),min(t)):max(max(k1),max(t)); s1=zeros(1,length(h));s2=s1; s1(find((h=min(k1))(h=max(k1))))=f1