FxNewton.m

function [e W y]=FxNewton（JKLfsf0xdSS_miu）
    W=zeros（LJ）;
    e=zeros（size（d））;
    v=zeros（size（d））;
    y=zeros（Jlength（e（1:）））;
    X=zeros（L1）;
    P=mean（x.^2）;
    XB=X;
    XQB=XB;
    M=length（S（:1））;
    Y=zeros（MJ）;
    Sf_=zeros（KJ）;
    for j=1:J
        for k=1:K
            B_=freqz（S_（:（j-1）*K+k）1[f0 2*f0]fs）;
            Sf_（kj）=B_（1）;
        end
    end
    Sf_inv=pinv（Sf_）;
    p=real（Sf_inv）;
    q=imag（Sf_inv）;
    x_=hilbert（x）;
    xB=real（x_）;
    xQB=imag（x_）;
    for