数字信号产生和运算.docx

电气系电子信息工程  PAGE \* MERGEFORMAT 42 实验一基本信号的产生和运算 1.1实验目的 利用MATLAB产生基本信号、绘制信号波形； 实现信号的基本运算，重点掌握卷积和运算。 1.2实验原理 MATLAB的基本应用方法 命令窗（Command window）的使用： 输入各类变量或函数名称，按回车即得到当前变量或函数值； 输入各类命令，按回车即得到该命令执行结果； 若需要输入多行命令或程序，各行间用“；”间隔； 2、M文件的编制与调制执行： 一行中“%”以后内容为注释部分，不影响程序执行； 程序编制完毕后，如果出现错误，可在命令窗口看到错误类型及位置，根据错误检测信息对程序进行调制； MATLAB命令及函数 在MATLAB中，任何变量或函数均表现为向量，任何向量的元素编号均从1开始；设定坐标向量n和信号向量x；x和n为长度相同的向量，向量的编号从1开始：；坐标向量可以直接逐点写出：；也可以采用起点，终点和步长的形式写出：；信号向量可以直接逐点写出：；也可以采用与n有关的函数运算形式写出：。 作图时主要通过合理设置n的范围及步长来保证变量坐标的正确性；可以利用title,axis等函数为图形设置说明和坐标范围。特别注意，作图时必须保证坐标向量与信号向量长度完全一致； 0101：离散序列的作图 0102：将图形表现为连续曲线 直接表现离散序列 n=[2 3 4 5 6 7]； n=[2 3 4 5 6 7]； x=[1 2 3 4 3 2]； x=[1 2 3 4 3 2]； plot(n,x)； stem（n,x）； 0203：信号表现为坐标向量的函数 n=[2 3 4 5 6 7]； x=exp(j*(pi/8)*n)； plot(n,x)； 0204：图形说明和坐标范围的设置 n=[-20:0.5:20]； x=exp(j*(pi/8)*n)； plot(n,x),title(n=[-20:0.5:20]；x=exp(j*(pi/8)*n)；plot(n,x)）； axis([-20,20,-2,2])； 基本信号的产生 单位抽样序列： 单位阶跃序列： 例一： 代码： clear all N=32; x=zeros(1,N); x(1)=1; xn=0:N-1; stem(xn,x); axis([-1 33 0 1.1]) 结果：见例一图 例一图 例二图 例三图 例二： 代码： clear all N=32; K=20; x=zeros(1,N); x(K)=1; xn=0:N-1; stem(xn,x); axis([-1 33 0 1.1]) 结果：见例二图 例三： 代码： clear all N=32; x=ones(1,N); x(1)=1; xn=0:N-1; stem(xn,x); axis([-1 32 0 1.1] 结果:见例三图 例四： 代码： n=[0:10],x(n)=(0.9)^n n=[0:10];x=(0.9).^n; stem(n,x) 结果：见例四图 例四图 例五图 例五： 代码： clear; N=5; M=6; L=N+M-1; x=[1,2,3,4,5]; h=[6,2,3,6,4,2]; y=conv(x,h); nx=0:N-1; nh=0:M-1; ny=0:L-1; subplot(231); stem(nx,x,.k);xlabel(n);ylabel(x(n));grid on; subplot(232); stem(nh,h,.k);xlabel(n);ylabel(h(n));grid on; subplot(233); stem(ny,y,.k);xlabel(n);ylabel(y(n));grid on; 结果：见例五图 例六： 代码 ： clear; x=ones(100); t=1:100; b=[.001836,.007344,.011016,.007374,.001836]; a=[1,-3.0544,3.8291,-2.2925,.55075]; y=filter(b,a,x); plot(t,x,r.,t,y,k-);grid on; ylabel(x(n) and y(n)) xla