基于MATLAB线性与圆周卷积计算及分析.pdfVIP

己知两个有限长序列：

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）实验前，预先笔算好这两个序列的线性卷积及下列几种情况的圆周卷积

x(n)⑤h(n)②x(n)⑥h(n)②x(n)⑨h(n)③x(n)⑩h(n)

①

（2）编制一个计算两个序列线性卷积的通用程序，计算x(n)*h(n)。

（3）编制一个计算圆周卷积的通用程序，计算上述4种情况下两个序列x(n)与h(n)的

圆周卷积。

首先，要建立两个.m文件的函数，以供计算不同x(n)*h(n)时调用！

函数一：序列x的圆周移位函数cirshifted

functiony=cirshifted(x,m,N)

iflength(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);

函数二：圆周卷积函数yzjj

functionyc=yzjj(x1,x2,N)

iflength(x1)N

error;

end

iflength(x2)N

error;

end

x1=[x1,zeros(1,N-length(x1))];%填充序列x1(n)使其长度为N，序列h(n)的长度为N1,序列x(n)

的长度为N2

x2=[x2,zeros(1,N-length(x2))];

n=[0:1:N-1];

x2=x2(mod(-n,N)+1);

H=zeros(N,N);forn=1:1:N

H(n,:)=cirshifted(x2,n-1,N);end

yc=x1*H;%计算圆周卷积

己知两个有限长序列:)4(

=δ+2δ−1+3δ−2+4δ−3+5δ−)()()()

(xnnnnnn)3(=δ−2+2δ−)

()+2δ−1+δ

()()()(hnnnnnn

(1)实验前，预先笔算好这两个序列的线性方程及下列几种情况的非线性方程

1)(NX5)(Hn2)(NX6)(HN2)(NX9)(HN3)(NX10)(HN

(2)编制一个计算两个序列线性方程的通用程序，计算)(*)(xNHN。

(3)编制一个计算波形波形的通用程序，计算上述4种情况下两个序列)(NX与)(Hn的波

形波形。

首先，要建立两个.m文件的函数，以供计算不同)(*)(xNHN时调用!

函数一:序列x的圆周升降函数cirshifted

=循环移位(x,m,n)的函数if

length(x)nerror(x的长度

必须小于n);结束x=[x，零

（1，N‑长度（x））]；n=[0：

1：n‑1]；y=x(mod(n‑m,

n)+1);

函数二:波形函数YZJJ

函数ycyzjj(x1,x2,n)如果长度(x1)n错误；如果Length(x2)n错误，则结束；END

1=[x1,Zeros(1,N‑Length(x1))]