hunxiang

matlab上实现的音频处理数字混响算法

function [y] = schroeder（xdelaysrt60fs）
% function [y] = schroeder（xdelaysrt60fs）
% x = input signal （monophonic sound or impulse）
% delays = delay line lengths （of comb filters） in samples
% rt60 = reverberation time at zero-frequency （e.g. 1.5）
% fs = sampling rate （e.g. 44100）
%
% This code was created by
% Riitta V滗n鋘en
% Helsinki UNiversity of Technology
% Laboratory of Acoustics and Audio Signal Processing
% email: riitta.vaananen@hut.fi
% 
% Copyright 10.5.2000 Riitta V滗n鋘en
%
%

% example delays:
% delays = [2191 2971 3253  3307];

% allpass filter delays:
 apdelays = [441 713];
 apgains = [0.5 0.5];
%
% fs = 44100
% lowpass filter coefficients
%
% Feedback gain:
%Kp = zeros（1length（delays））;
Kp = 10.^（（-3.0.*delays）./（rt60*fs））

% Lowpass fi

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3096  2009-09-25 21:23  混响\schroeder.m

     文件        309  2009-10-10 13:10  混响\use_schroeder.m

     文件       2009  2009-10-15 12:38  混响\混响C1.0.txt

     文件       2342  2009-12-02 00:07  混响\混响C1.1.txt

     文件       2530  2009-12-31 21:17  混响\混响1.2.txt

     目录          0  2009-11-22 22:39  混响

----------- ---------  ---------- -----  ----

                10286                    6