声音反卷积matlab,用MATLAB做反卷积

声⾳反卷积matlab,⽤MATLAB做反卷积

关键词：反卷积 MATLAB fft 频移分母中频谱零点卷积核

%==============================================================

代码如下：

clear all;clc;

h = [1 1 1 1] % 要求

f = [1 -2 3 -2] % 已知

g = conv(h,f) % 已知 g = h*f 这⾥卷积结果g知道，f知道，f视作卷积核，反卷积求h

Nf = length(f);

Nh = length(h);

N = 16;

% N = input('请输⼊⽐Nf*Nh⼤的数：');

Gk = fftshift(fft(g,N));

Fk = fftshift(fft(f,N));

container = repmat(zeros(size(Fk)),Nh-1,1);

for i = 1:(Nh-1)

棱镜常数

k = i*Nf;

container(i,:) = Fk([end-k+1:end 1:end-k]); %此处是做频移，算的是上⾯的Fpk那个连乘 Nh-1 次的公式；container的每⼀⾏保存的都是将要连乘的频移后结果

end

Fpk = prod(container) % 此处连乘

Fek = Fk.*Fpk

hpn = ifft(ifftshift(Gk.*Fpk))

fen = ifft(ifftshift(Fek))

凝聚力工程

n = (0:15)';

割草机设计

fen = fen';

艾滋病实名制hpn = hpn';

table = real([n fen hpn])

利比亚战争运⾏结果如下：

和论⽂中的结果是⼀致的

本文发布于:2024-09-21 10:59:42，感谢您对本站的认可！

标签：卷积知道频谱禁播

留言与评论（共有 0 条评论）