资源描述：
《哈工大数字信号处理大作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数字信号处理上机实验报告学号：姓名:实验题目一1.实验要求：序列卷积计算（1）编写序列基本运算函数，序列相加、相乘、翻转、求和；（2）使用自定义函数计算序列线性卷积，并与直接计算结果相比较。两个序列分别为：，2.实验过程和步骤：包含题目分析，实验程序和流程图（程序要有必要的注释）3.实验结果和分析：包含程序运行结果图，结果分析和讨论(一)基本运算函数1.原序列2.序列相加序列相加程序function[y,n]=sigadd(x1,n1,x2,n2)%implementsy(n)=x1(n)+x2(n)%-------------------------
2、--------------------%[y,n]=sigadd(x1,n1,x2,n2)%y=sumsequenceovern,whichincludesn1andn2%x1=firstsequenceovern1%x2=secondsequenceovern2(n2canbedifferentfromn1)%n=min(min(n1),min(n2)):max(max(n1),max(n2));%durationofy(n)y1=zeros(1,length(n));y2=y1;y1(find((n>=min(n1))&(n<=max(n1))=
3、=1))=x1;%x1withdurationofyy2(find((n>=min(n2))&(n<=max(n2))==1))=x2;%x2withdurationofyy=y1+y2;%sequenceaddition1.序列相乘序列相乘程序function[y,n]=sigmult(x1,n1,x2,n2)%implementsy(n)=x1(n)*x2(n)%---------------------------------------------%[y,n]=sigmult(x1,n1,x2,n2)%y=productsequenceover
4、n,whichincludesn1andn2%x1=firstsequenceovern1%x2=secondsequenceovern2(n2canbedifferentfromn1)%n=min(min(n1),min(n2)):m(min(n1),min(n2))%durationofy(n)y1=zeros(1,length(n));y2=y1;y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;%x1withdurationofyy2(find((n>=min(n2))&(n<=max(n2))==1))=x2
5、;%x2withdurationofyy=y1.*y2;%sequencemultiplication1.序列翻转序列翻转程序function[y,n]=sigfold(x,n)%implementsy(n)=x(-n)%---------------------------------------------%[y,n]=sigfold(x,n)%y=fliplr(x);n=-fliplr(n);1.序列移位序列移位程序function[y,n]=sigshift(x,m,n0)%implementsy(n)=x(n-n0)%------------
6、---------------------------------%[y,n]=sigshift(x,m,n0)%n=m+n0;y=x;主程序x1=[0:5];x2=[0,1,2,3];n1=0:5;n2=0:3;%N=n1+n2-1;figure(1)subplot(211)stem(x1)xlabel('x1')subplot(212)stem(x2)xlabel('x2')title('原序列')x=sigadd(x1,n1,x2,n2);figure(2)stem(x)xlabel('x1+x2')title('序列相加')figure(3)[
7、x,n]=sigfold(x1,n1);stem(n,x)xlabel('x1(-n)')title('序列翻转')[x,n]=sigshift(x,n,2);figure(4)stem(n,x)xlabel('x1(-n+2)')title('序列移位')x=sigmult(x1,n1,x2,n2);figure(5)stem(x)title('序列相乘')xlabel('x1*x2')(一)自定义函数计算线性卷积1.题目分析使用上一题中的序列相乘、翻转和求和子函数计算线性卷积，并与这直接用conv函数计算的线性卷积结果相比较。1.流程图1.实验程序
8、clearx1=[0:5];x2=[0:3];n1=[0:5];n2=[0:3];leng1=
