hthr = 2logn^05; eta = sumabsx^2-n1n; crit = lognlog2^15n^05; thr = thselectxrigrsure; thrthr hthr = hthr; threta crit = hthr;对上述代码进行改进写一段改进的阈值函数的代码

function [thr] = improvedThreshold(x, n) hthr = (2*log(n)).^0.5; eta = sum(abs(x).^2-n,1)./n; crit = (log(n)/log(2))^(1.5)/(n.^0.5); thr = thselect(x, 'rigrsure'); thr(thr > hthr) = hthr; thr(eta < crit) = hthr; thr(eta < hthr) = (hthr + thr(eta < hthr)) / 2; % added line end