hc371.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第三章例3.7.1程序hc371
% 循环卷积与线性卷积
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x1=[1,2,2,1];x2=[1,-1,-1,1];    	% 给出x1,x2序列	
y=conv(x1,x2)                  % 调用线性卷积函数求线性卷积
N=7;
x3 = circonv(x1,x2,N)       % 调用循环卷积函数求循环卷积
for N=[6,5,4]               % 对N取三种值
    x3 = circonv(x1,x2,N)   % 调用循环卷积函数求循环卷积 
    e=x3-y(1:N)             % 求线性卷积与循环卷积之差
end