循环卷积函数(Matlab) - 梦想时空 - dreamland2006

载入中

自定义HTML载入中...

快速链接

循环卷积函数(Matlab) [转贴 2008-01-25 16:22:48]

文章来源: 小刚的个人空间

function fn=circonvt(x1,x2,N)
%circonvt函数实现输入序列x1和x2的循环卷积,fn为输出序列
%N 为循环卷积长度
%实现方法:fn=sum(x1(m)*x2((n-m) mod N))
if (length(x1)>N|length(x2)>N)%判断输入信号的长度
error('N的长度必须大于输入数据的长度');
end
x1=[x1,zeros(1,N-length(x1))];
x2=[x2,zeros(1,N-length(x2))];
m=0:N-1;
x=zeros(N,N);
for n=0:N-1
x(:,n+1)=x2(mod((n-m),N)+1)';
end;
fn=x1*x; %循环计算卷积