数字信号的数字相关和数字卷积运算《生物医学信号处理》实习报告2.docVIP

《生物医学信号处理》实习报告 学生姓名： 学号： 实验室名称： 项目名称：数字信号的数字相关和数字卷积运算 项目内容： 编写函数实现两个随机序列的线性、循环相关和线性、循环卷积。 原理（写出具体的计算公式） 一．数字相关和数字卷积的定义 1.线性相关：是相关的一种运算，讨论两信号之间的同步性、相似性或相同性，或两信号的变化规律是否具有线性关系或接近线性关系的程度。 2.循环相关：是针对序列的循环移位的一种相关运算。 3.线性卷积：是线性时不变系统分析的重要工具，很多滤波器的设计中都要用到卷积运算。 4.循环卷积：是针对序列的循环移位的一种运算。 二．数字相关和数字卷积的计算公式 1.线性相关函数： 2.循环相关函数： 3.线性卷积函数： 4.循环卷积函数： 编写的源程序： function v1=linel(x,y)%线性相关函数 N=length(y); M=2*N-1; v=zeros(N,M);% 先构建N*M的矩阵，实现线性左移 for i=1:N for j=N:M v(i,j-i+1)=y(j-N+1); end end v1=x*v; end function v2=circlel(x,y)%循环相关函数 N=length(y); v=zeros(N,N);%先构建N*N的矩阵，实现循环左移 for i=1:N for j=1:N v(i,j)=y(j); end L=y(1); for k=1:N-1 y(k)=y(k+1); end y(N)=L; end v2=x*v; end function v3=liner(x,y)%线性卷积函数 N=length(y); M=2*N-1; v=zeros(N,M);%先构建N*M的矩阵，实现线性右移 for i=1:N for j=1:N v(i,j+i-1)=y(j); end end v3=x*v; end function v4=circler(x,y)%循环卷积函数 N=length(y); v=zeros(N,N);% 先构建N*N的矩阵，实现循环右移 for i=1:N for j=1:N v(i,j)=y(j); end L=y(N); for k=1:N-1 y(N-k+1)=y(N-k); end y(1)=L; end v4=x*v; end 编写的源程序： A = textread(ECG1.txt); % 导入数据 subplot(8,1,1); %页面分八行一列，画在第一个区域 plot(A); %画ECG1 title(ECG1); grid on; %显示网格 B = textread(ECG2.txt); %导入数据 subplot(8,1,2); %页面分八行一列，画在第二个区域 plot(B); %画ECG2 title(ECG2); grid on; %显示网格 subplot(8,1,3); %页面分八行一列，画在第三个区域 plot(linel(A,B)); %计算线性相关，并画曲线 title(线性相关); grid on; %显示网格 subplot(8,1,4); %页面分八行一列，画在第四个区域 plot(xcorr(A,B)); %计算线性相关，并画曲线 title(MATLAB自带线性相关); grid on; %显示网格 subplot(8,1,5); %页面分八行一列，画在第五个区域 plot(circlel(A,B)); %计算循环相关，并画曲线 title(循环相关); grid on; %显示网格 subplot(8,1,6); %页面分八行一列，画在第六个区域 plot(liner(A,B)); %计算线性卷积，并画曲线 title(线性卷积); grid on; %显示网格 subplot(8,1,7); %页面分八行一列，画在第七个区域 plot(conv(A,B)); %计算线性卷积，并画曲线 title(MATLAB自带线性卷积); grid on; %显示网格 subplot(8,1,8); %页面分八行一列，画在第八个区域 plot(circler(A,B)); %计算循环卷积，并画曲线 title(循环卷积); grid on; %显示网格 结论（画出要求的图形） 图