实验一 离散时间信号与系统的时域分析.pptVIP

实验一 离散时间信号与系统的 时域分析 一、实验目的 1 、掌握离散时间信号的 MATLAB 表示； 2 、信号运算； 3 、离散时间信号的卷积运算。 二、实验原理 离散时间信号 ? 离散时间信号只在某些离散的瞬时给出函数值，而在 其他时刻无定义。它是时间上不连续按一定先后次序 排列的一组数的集合，称为时间序列，用 x(n) 表示， n 取整数代表时间的离散时刻。 ? 在 matlab 中用向量来表示一个有限长度的序列。 序列的类型 ? 为了分析的方便，在数字信号处理中规定了一些基本 的序列。 单位采样序列 ? function [x,n]=impseq(n 1 ,n 2 ,n 0 ) ? n=[n1:n2]; ? x=[(n-n 0 )==0]; 0 0 0 1 0 ( ) 0 0 1 ( ) 0 n n n n n n n n n ? ? ???? ? ? ? ? ???? ? ? ???? ? ? ? ? ???? ? ? ? 调用该函数 ? [x,n]=impseq(-2,8,2); ? stem(n,x) 单位采样序列的另一种生成方法 ? n0=-2; ? n=[-10:10]; ? nc=length(n); ? x=zeros(1,nc); ? for i=1:nc ? if n(i)==n0 ? x(i)=1 ? end ? end ? stem(n,x) 单位阶跃序列 ? function [x,n]=stepseq(n 1 ,n 2 ,n 0 ) ? n=[n1:n2]; ? x=[(n-n 0 )=0]; 0 0 0 1 0 ( ) 0 0 1 ( ) 0 n n n n n n n n n ? ? ???? ?? ? ? ? ???? ? ? ???? ?? ? ? ? ???? ? ? ? 调用该函数 ? [x,n]=stepseq(-2,8,2); ? stem(n,x) 实数指数序列 x(n)=a n ( 运算符“ .^”) ? n=[0:10]; ? x=0.9.^n; ? stem(n,x) 复数指数序列 ? n=[-10:10]; ? alpha=-0.1+0.3*j; ? x=exp(alpha*n); ? real_x=real(x); image_x=imag(x); ? mag_x=abs(x); phase_x=angle(x); ? subplot(2,2,1); stem(n,real_x) ? subplot(2,2,2); stem(n,image_x) ? subplot(2,2,3); stem(n,mag_x) ? subplot(2,2,4); stem(n,phase_x) ( ) ( ) j n x n e ? ? ? ? ( 0.1 j0.3)n x(n) e ( 10 n 10) ? ? ? ???? ? ? ? 正、余弦序列 0 ( ) sin( ) x n n ? ? ? ? 0 ( ) cos( ) x n n ? ? ? ? ? n=[0:10]; ? x=3*cos(0.1*pi*n+pi/3); ? stem(n,x) 随机序列 ? rand(1,N) 产生其元素在 [0 ， 1] 之间均匀分布长度为 N 的随机序列 ? randn(1,N) 产生均值为 0 ，方差为 1 ，长度为 N 的高斯随 机序列 周期序列 ? 如果序列 x(n)=x(n+N) ? 如何生成周期序列 1 、 将一个周期复制 p 次； 2 、借助矩阵运算、 matlab 下标能力。先生成一个包含 p 列 x(n) 值的矩阵，然后用结构（ : ）来把 p 列串接成一个 长周期序列。因为这个结构只能用于列向，最后还需要 做矩阵转置获得所需序列。 ? 例：生成 3 个周期的周期序列 ? x=[1,2,3]; 一个 x(n) ? xn=x ‘ *ones(1,3) 生成 p 列 x(n) ? xn=xn(:) ‘ 将 p 列串接成长列序列并转置 信号运算 ? 信号时移 ? 信号倒置 ? 信号加 ? 信号乘 ? 信号微分 ? 信号积分等 序列运算 ? 在 matlab 中进行序列运算要求参与运算的序列的长度要一样，如 果出现长度不一样或者长度相同但采样位置不同时，不能直接利 用加减运算符，利用 matlab 中的下标运算使其具有相同的长度。 ? function [y,n]=sigadd(x 1 ,n 1 ,x 2 ,n 2 ) ? n=min(min(n 1 ),min(n 2 )) : max(max(n 1 ),max(n 2 )); ? y1=zeros(1,length(n)); ? y2=y1; ? y1(find((n=min(n1))(n=ma