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

1、用MATLAB实现常用离散时间信号及其时域运算摘要：本实验实现了运用MATLAB的绘图指令绘制离散时间信号，用MATLAB语言实现离散时间信号的时域运算。达到了熟悉MATLAB基本使用方法以及了解典型信号的MATLAB表示方法，掌握MATLAB函数的应用的目的。关键词:离散时间信号 时域运算 MATLAB 正文:离散时间信号在时间上是离散的，只在某些不连续的规定瞬间给出函数值，在其他时间没有定义典型离散时间信号：单位样值信号： 单位阶跃序列：矩形序列：离散时间信号的时域运算信号的时域运算包括信号的平移、反褶、尺度及信号的相加、相乘.用MATLAB的符号运算命令来表示离散时间信号的时域运算，然后

2、用ezplot命令绘制出其结果波形图。相加和相乘：Functionf,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)两序列相乘可参

3、照此例平移Functionf,k=lspy(ff,kk,k0)K=kk+k0F=ff反褶Functionf,k=lsfz(f1,k1)F=fliplr(f1);K=-fliplr(k1)任务1. 编制用于产生下列信号的通用程序，要求对于任意给定的参数都能实现所要求信号，调试并运行这些通用程序，具体产生有指导教师指定的信号并绘制信号的波形。（1）x(n)=A(n-n0)(2) x(n)=Au(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能对信号进

4、行时域分析，将信号随时间的变化规律用二维曲线表示出来就可以描述，而且MATLAB具有强大的图形处理功能及符号运算功能，为我们直观的看到信号的模型提供了方便。附录: 任务一(1)x(n)=A(n-n0)cleart1=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 ont

5、itle('³å¼¤ÐòÁÐ X(n)=A&(n-n0)')(2) x(n)=Au(n-m1)-u(n-m2)clearm1=input('µÚÒ»¸öÌø±äµãm1=');m2=input('µÚ¶þ¸öÌø±äµãm2=

6、9;);A=input('ÊäÈë·ù¶ÈÖµA=');if m1<m2 t=m1:1:m2-1; x=ones(1,m2-m1);else t=m2:1:m1-1; x=-ones(1,m1-m2);endstem(t,A.*x,'filled','r'),grid ontitle(' X(n)=Au(n-m1)-u(n-m2)')任务二clear,i=1;for k=-3:3 if k<0 f(i)=k+3; els

7、e f(i)=3; end i=i+1; f(i)=0; endk=-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<=m

8、ax(k1)=f1;s2(find(h>=min(t)&(h<=max(t)=u;f21=s1.*s2;subplot(222),stem(h,f21,'filed'),grid on,title('f(k-2)u(k)')axis(min(k)-1,max(k)+1,min(f21)-0.5,max(f21)+0.5) f22=fliplr(f);k22=-fliplr(k-2);subplot(223),stem(k22,f22,'filed'),grid on,title('f(-k+2)') axis(min(k)-1,max(k)+1,min(f22)-0.5,max(f22)+0.5)f1=f;k1=k+2;t=t+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;s2(find(h>=min(t)&(h<=max(t)=u;