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