[哈工大数字信号处理大作业.docVIP

数字信号处理上机实验报告 学号： 姓名: 实验题目一 1. 实验要求：序列卷积计算 （1）编写序列基本运算函数，序列相加、相乘、翻转、求和； （2）使用自定义函数计算序列线性卷积，并与直接计算结果相比较。两个序列分别为： ， 2. 实验过程和步骤：包含题目分析，实验程序和流程图（程序要有必要的注释） 3. 实验结果和分析：包含程序运行结果图，结果分析和讨论 基本运算函数 原序列 序列相加 序列相加程序 function [y,n]=sigadd(x1,n1,x2,n2) %implements y(n)=x1(n)+x2(n) %--------------------------------------------- % [y,n] = sigadd(x1,n1,x2,n2) % y = sum sequence over n, which includes n1 and n2 % x1 = first sequence over n1 % x2 = second sequence over n2 (n2 can be different from n1) % n=min(min(n1),min(n2)):max(max(n1),max(n2)); %duration of y(n) y1=zeros(1,length(n)); y2=y1; y1(find((n=min(n1))(n=max(n1))==1))=x1; %x1 with duration of y y2(find((n=min(n2))(n=max(n2))==1))=x2; %x2 with duration of y y=y1+y2; %sequence addition 序列相乘 序列相乘程序 function [y,n]=sigmult(x1,n1,x2,n2) %implements y(n)=x1(n)*x2(n) %--------------------------------------------- % [y,n] = sigmult(x1,n1,x2,n2) % y = product sequence over n, which includes n1 and n2 % x1 = first sequence over n1 % x2 = second sequence over n2 (n2 can be different from n1) % n=min(min(n1),min(n2)):m(min(n1),min(n2)) %duration of y(n) y1=zeros(1,length(n)); y2=y1; y1(find((n=min(n1))(n=max(n1))==1))=x1; %x1 with duration of y y2(find((n=min(n2))(n=max(n2))==1))=x2; %x2 with duration of y y=y1.*y2; %sequence multiplication 序列翻转 序列翻转程序 function [y,n]=sigfold(x, n) %implements y(n)=x(-n) %--------------------------------------------- % [y,n] = sigfold(x,n) % y=fliplr(x); n=-fliplr(n); 序列移位 序列移位程序 function [y,n]=sigshift(x,m,n0) %implements y(n)=x(n-n0) %--------------------------------------------- % [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) [x,n] = sigfold(x1,n1); stem(n,x) xlabel(x1(-n)) title(序列翻转) [x,n] = sigshif