序列的matlab实现样本

实验十一 离散信号时域分析的MATLAB实现 一、实验目的 1. 熟悉MATLAB 编程方法、常用语句和可视化绘图技术； 2. 掌握序列时域运算的 MATLAB 编程方法。 二、实验原理 在用 MATLAB 表示离散信号并将其可视化时，由于矩阵元素个数是有限的，因此无法表示无限长序列；另外，离散信号无法进行符号运算。在 MATLAB 中，绘制离散序列波形图的专用命令为 stem( )。其格式有： （1）stem(k,f) 在图形窗口中，绘制出样值顶部为空心圆的序列 f(k)波形图。 （2）stem(k,f,’fill’) 在图形窗口中，绘制出样值顶部为实心圆的序列 f(k)波形图。 下面介绍离散序列的 MATLAB 表示、基本运算（相加、相乘、平移、反转、尺度变换）、卷积和的实现及其图形显示方法。 1.单位序列(k ) 单位序列的定义： 0,00,1)(kkk 下面为绘制(k -k 0)波形图的子程序： fu nction impseq(k1,k2,k0) %单位序列(k-k0),k0 为时移量 k=k1:k2; %k1,k2 为序列的起止序列号 n=length(k); x=zeros(1,n); x(1,k0-k1+1)=1; %在k0 时刻信号赋值为1 stem(k,x,'fill') axis([k1,k2,0,1.1]) title('单位序列d(k-k0)') 输入如下命令，则可获得单位序列(k -3 )的波形图，如图11-1 所示。 impseq(-1,5,3) 2 .单位阶跃序列(k ) 单位序列的定义： 0,00,1)(kkk 下面为绘制(k -k 0 )波形图的MATLAB 子程序。 function stepseq(k1,k2,k0) %单位阶跃序列,k0 为时移量 k=k1:k0-1; %k1,k2 为序列的起止序列号 kk=length(k); x=zeros(1,kk); %k0 前信号赋值为0 stem(k,x,'fill') %绘出 k1~k0-1 的波形(0 值) hold on n=k0:k2; nn=length(n); 图11-1 图11-2 x=ones(1,nn); %k0 后信号赋值为1 stem(n,x,'fill') %绘出k0~k2 的波形(1 值) hold off axis([k1,k2,0,1.1]) title('单位阶跃序列') 运行如下命令，则可获得单位序列(k -3 )的波形图，如图11-2 所示。 stepseq(-1,10,3) 3 .序列的相加（减）、相乘运算 对序列向量f1(k)、f2(k)相加或相乘，可以通过补零的方式使 f1(k)、f2(k)成为具有相同维数的序列向量s 1(k)、s 2(k)，然后对s 1(k)、s 2(k)相加或相乘。因此，序列向量f1(k)、f2(k)的维数可以不同。以下函数可实现序列向量f1(k)、f2(k)的相加或相乘运算。 function [f,k]=sigadd(f1,k1,f2,k2) %实现序列f1,f2 的相加,相减,相乘,可据实际需...