数字信号处理及应用卢光跃.docxVIP

??求序列y(n)=x(n)+h(n)。?解：MATLAB程序如下：?? ???n=-10:10;? ????x=[zeros(1,10),1,2,3,4,5,zeros(1,6)];?????h=[zeros(1,10),1,2,1,2,zeros(1,7)];?????y=x+h;? ????stem(n,y,*k);??? ? ? ? grid?on??? ? ? ? ? ? ? ? ? ? ? ??? 图13?两序列的和? ??4.两序列的积? 例：已知x(n)=?δ(n)+2δ(n-1)+3δ(n-2)+4δ(n-3)+5δ(n-4)，h(n)=?δ(n)+2δ(n-1)+δ(n-2)+2δ(n-3)?????? ? ?求序列y(n)=x(n).*h(n)。?解：MATLAB程序如下：? ????n=-10:10;? ????x=[zeros(1,10),1,2,3,4,5,zeros(1,6)];?????h=[zeros(1,10),1,2,1,2,zeros(1,7)];?????y=x.*h;? ????stem(n,y,*k);??????grid?on???? ? ? ? 运行结果：? ????? ? ? ? ? ? ? ??图14??两序列的积? ??5.两序列的标乘（MATLAB指令为dot）? ? 例：已知x(n)=?δ(n)+2δ(n-1)+3δ(n-2)+4δ(n-3)+5δ(n-4)，h(n)=?δ(n)+2δ(n-1)+δ(n-2)+2δ(n-3)?????? ? ? 求序列x(n)、?h(n)的标乘??解：MATLAB程序如下：? ????n=-10:10;? ????x=[zeros(1,10),1,2,3,4,5,zeros(1,6)];?????h=[zeros(1,10),1,2,1,2,zeros(1,7)];?????y=dot(x,h)?? ? ? ? 运行结果：? ????y?=??????? ??16? ??6.序列的累加? ? 例：（1）求x(n)=?δ(n)+2δ(n-1)+3δ(n-2)+4δ(n-3)+5δ(n-4)的累加? ?解：MATLAB程序如下：?? ???n=-10:20;? ????x=[zeros(1,10),1,2,3,4,5,zeros(1,16)];?????y=cumsum(x);?????stem(n,y,*k);?????grid?on???? ? ? ? 运行结果：? ?????? ?? ???? ? ? ??图15??序列（1）的累加? ???例：（2）求h(n)=?δ(n)+2δ(n-1)+δ(n-2)+2δ(n-3)的累加??解：MATLAB程序如下：?????n=-10:20;? ????h=[zeros(1,10),1,2,1,2,zeros(1,17)];?????y=cumsum(h);?????stem(n,y,*k);?????grid?on???? ? ? ? 运行结果：? ????? ???? ? ? ?图16??序列（2）的累加? ? ??2.3序列的卷积运算? ??例：?x(n)=?δ(n)+2δ(n-1)+3δ(n-2)+4δ(n-3)+5δ(n-4)? h(n)=?δ(n)+2δ(n-1)+δ(n-2)+2δ(n-3)? 解：（1）序列卷积计算公式：y(n)=x(n)*h(n)=?x(m)h(n-m)? （2）根据定义计算y(n)?MATLAB程序如下：? N=10;?M=10;?L=N+M-1;? x=[1,2,3,4,5,zeros(1,5)];??h=[1,2,1,2,zeros(1,6)];??y=conv(x,h);?n=0:L-1;?stem(n,y,*k);??grid?on??运行结果：? ? ????????????图17?卷积和结果序列? 通用程序：? ???function[y,n]=dconv(x,nx,h,nh);????x=input(请输入序列（1）：x（n）=);????nx=input(请输入序列（1）长度（如0：5）=);????h=input(请输入序列（2）：h（n）=);????nh=input(请输入序列（2）长度（如0：5）=);????y=conv(x,h)????k1=nx(1)+nh(1);? ???k2=length(x)+length(h)-2;????n=k1:k1+k2? ???subplot(131);?stem(nx,x,*r);?xlabel(nx);????ylabel(x(n));?grid?on?? ???subplot(132);?stem(nh,h,*b);?xlabel(nh);????ylabel(h(n));?grid?on