線性卷积与圆周卷积的计算.docxVIP

线性卷积与圆周卷积的计算 一、基本原理 1. 线性卷积 当系统输入序列为x(n),系统的单位冲击响应为h(n)，输出序列为y（n），则线性时不变系统输入、输出间的关系为： Y（n）=h（n）*x（n） 2. 圆周卷积 设两个有限长序列和，均为点，其点的DFT分别为和，如果，则 eq \o\ac(○,N) 实验内容与要求 已知两个有限长序列： 实验前，预先笔算好这两个序列的线性卷积及下列几种情况的圆周卷积 = 1 \* GB3 \* MERGEFORMAT ①⑤ = 2 \* GB3 \* MERGEFORMAT ②⑥ = 3 \* GB3 \* MERGEFORMAT ③⑨ = 4 \* GB3 \* MERGEFORMAT ④⑩ （2）编制一个计算两个序列线性卷积的通用程序，计算。 （3）编制一个计算圆周卷积的通用程序，计算上述4种情况下的两个序列的圆周卷积。 （4）上机调试并记录实验结果 （5）将实验结果和预先笔算的结果比较，验证其正确性。 三、实验过程 function yc=circonv(x1, x2, N) if length(x1)N error(N必须大于等于x1的长度); end if length(x2)N error(N必须大于等于x2的长度); end x1=[x1, zeros(1,N-length(x1))]; x2=[x2, zeros(1,N-length(x2))]; n=[0:1:N-1]; x2=x2(mod(-n, N)+1); H=zeros(N, N); for n=1:1:N H(n,:)=cirshiftd(x2, n-1, N); end yc=x1*H; %计算圆周卷积 function y=cirshiftd(x, m, N) if length(x)N error(x的长度必须小于N); end x=[x, zeros(1, N-length(x))]; n=[0:1:N-1]; y=x(mod(n-m,N)+1); = 1 \* GB3 \* MERGEFORMAT ①x(n)⑤y(n) clear all xn=[1 2 3 4 5]; hn=[1 2 1 2]; N1=length(xn); N2=length(hn); y1n=conv(xn, hn); ycn=circonv(xn, hn, 5); ny1=[0:1:length(y1n)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1, y1n); subplot(2,1,2); stem(ny2, ycn); = 2 \* GB3 \* MERGEFORMAT ②x(n)⑥y(n) clear all; N1=5; N2=4; xn=[1 2 3 4 5]; hn=[1 2 1 2]; yln=conv(xn,hn); ycn=circonv(xn,hn,6); ny1=[0:1:length(yln)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,yln); ylabel(线性卷积); subplot(2,1,2); stem(ny2,ycn); ylabel(圆周卷积); = 3 \* GB3 \* MERGEFORMAT ③x(n)⑨y(n) clear all; N1=5; N2=4; xn=[1 2 3 4 5]; hn=[1 2 1 2]; yln=conv(xn,hn); ycn=circonv(xn,hn,9); ny1=[0:1:length(yln)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,yln); ylabel(线性卷积); subplot(2,1,2); stem(ny2,ycn); ylabel(圆周卷积); = 4 \* GB3 \* MERGEFORMAT ④x(n)⑩y(n) clear all; N1=5; N2=4; xn=[1 2 3 4 5]; hn=[1 2 1 2]; yln=conv(xn,hn); ycn=circonv(xn,hn,10); ny1=[0:1:length(yln)-1]; ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,yln); ylabel(线性卷积); subplot(2,1,2); stem(ny2,ycn); ylabel(圆周卷积); 思考题解答 （4