我现在正在做mtech项目。实际上我正在研究离散小波变换。我已经完成了小波分解。我一直坚持将信号重建回原来的状态。请帮帮我。

代码 -

f=10;
fs=200;
amp_x = 14;
% amplitude for sinusodial 2
amp_y = 9;
% Time vector 
t=(0:1:5000);
% Create a sine wave 
x = amp_x * (sin(2*pi*f/fs*t));
%filter coefficients
LO_D =[ 0.7071    0.7071];
HI_D =[-0.7071    0.7071];

C = conv(LO_D,x);
k = 2;
CA= C(2:k:length(C));

D = conv(HI_D,x);
g=2;
CD= D(2:g:length(D));

CA 和 CD 是我的近似和详细系数。现在我对 CD 系数进行平方以去除伪影。

for i=1:length(CD)
    cd1(i)=CD(i).^2;
end

平方后，我设置了一个阈值并删除了大于阈值的值。所以我得到一个去噪信号。

现在我需要将去噪信号重建为原始信号。如何在不使用内置命令的情况下做到这一点。请回答我。