【语音处理】基于加权压力匹配方法（WPMM）的私人声音系统研究（Matlab代码实现）

%% Algorithm 2
function [psi_est, q_temp] = Alg2_RecursiveBisectionSearch(zBT, zBT_zBconj, betaI, ZDH_ZD, zBconj, pBmin, Psi_Max, Psi_Min, thresh)
% Implementation of Alg 4
% Author: F Olivieri
psi_est = (Psi_Max - Psi_Min)/2 + Psi_Min;
q_temp = (zBT_zBconj + psi_est*ZDH_ZD + betaI)\zBconj;
pB = real(zBT*q_temp);
condition1 = abs(pB - pBmin);
if condition1 <= thresh ||  Psi_Max <= thresh % condition is met: Stop!
    return;
else % search for the parameter
    if pB < pBmin % yes: need to decrease psi_est
        Psi_Max =  psi_est;
    else % no: need to increase psi_est
        Psi_Min = psi_est;
    end
    psi_est = Alg2_RecursiveBisectionSearch(zBT, zBT_zBconj, betaI, ZDH_ZD, ...
        zBconj, pBmin, Psi_Max, Psi_Min, thresh);
end
end
%% Iterative search
function [q_fsp, psi_D_fsp] = fsp(q_fsp, psi_D_fsp, MatTemp, ZDH_ZD, zBconj, zBT, pB_min, delta_psi_Alg_1, lengthfreqbins, L)
% Implementation of Full-search of psi_D
% Author: F Olivieri
for freq_idx = 2:lengthfreqbins
    psi_D_temp = 1;
    q_temp = (MatTemp(:, :, freq_idx) + psi_D_temp*ZDH_ZD(:, :, freq_idx))\ zBconj(:, freq_idx);
    while zBT(freq_idx, :)*q_temp < pB_min && psi_D_temp > 0
        psi_D_temp = psi_D_temp - delta_psi_Alg_1;
        q_temp = (MatTemp(:, :, freq_idx) + psi_D_temp*ZDH_ZD(:, :, freq_idx))\ zBconj(:, freq_idx);
    end
    q_fsp(:, freq_idx) = q_temp;
    psi_D_fsp(freq_idx) = psi_D_temp;
end
end