小波信号降噪

N=1000; t=1:1000; f=sin(0.03*t);%生成正弦信号 load noissin; e=noissin;%注入噪声 subplot(221); plot(t,f);%画出原始波形 xlabel('样本序列n'); ylabel('原始信号幅值'); grid;%画出网格 subplot(222); plot(e);%画出加噪后的波形 xlabel('样本序列n'); ylabel('含有噪声的信号幅值'); grid;%画出网格 s1=wden(e,'minimaxi','s','one',5,'db3');%选用db3小波对信号进行5层分解,并对细节系数选用minimaxi阈值模式和尺度噪声 subplot(223); plot(s1);%画出降噪后的波形 xlabel('样本序列n'); ylabel('db3小波降噪后的信号幅'); grid;%画出网格 s2=wden(e,'heursure','s','one',5,'sym8');%选用sym8小波对信号进行5层分解,并对细节系数选用sure阈值模式和尺度噪声 subplot(224); plot(s2);%画出降噪后的波形 xlabel('样本序列n'); ylabel('sym8小波降噪后的信号幅度'); grid;%画出网格