在学习数字信号处理时遇到需要进行高点数的循环卷积处理时,不可能进行人工计算。这里借用matlab进行矩阵循环卷积。这里放出我自用的matlab函数。
下面展示 matlab源代码
。
function yn=circonv(hn,xn,N)
%先要根据N点数的长度对矩阵进行补零
%然后对hn进行处理获得其循环矩阵hn_f
%最后将hn_f与xn相乘获得yn
%YSC 2021.1.7 数字信号处理
x1 = length(hn); %先获得hn与xn的长度值
x2 = length(xn);
hn(x1+1:N)=0; %对其补零获得hn与xn
xn(x2+1:N)=0;
hn_f=ones(N); %生成一个N*N矩阵
%要先获得hn的循环矩阵
%总共赋值N次
for i=1:N
hn_f(:,i)=hn;
hn=circshift(hn,1); %将hn矩阵向右移位一位;
end %经过for循环后获得hn_f矩阵
%将hn_f与xn相乘便可
xn =xn';%要相乘必须先行列互换。
yn =hn_f*xn;%计算出结构
yn = yn'; %再按习惯改成行
end