snr = 4;
t=0:1/10000:1-0.0001;
y=sin(2*pi*t);
>> n = randn(size(t));
>> s=y+n;
>> xd = wden(s,'heursure','s','one',3,'sym8');
>> subplot(3,1,1);
>> plot(s);
>> xlabel('樣本序號n');
>> ylabel('幅值A');
>> title('含噪信號');
>> subplot(3,1,2);
>> plot(y);
>> title('參考信號');
>> xlabel('樣本序號n');
>> ylabel('幅值');
>> subplot(3,1,3);
>> plot(xd);
>> title('消噪信號');
>> xlabel('樣本信號');
>> ylabel('幅侄A')