examp13_6.m

来自「matlab7小波分析书中的例子」· M 代码 · 共 19 行

clc;
clear;
%设置信噪比和随机种子值
snr=4;
init=2055615866;
%产生原始信号sref和被高斯白噪声污染的信号s
[sref,s]=wnoise(1,11,sqrt(snr),init);
%用sym8小波对信号s进行3层分解并对细节系数
%选用SURE阈值模式和尺度噪声
xd1=wden(s,'heursure','s','one',3,'sym8');
xd2=wden(s,'sqtwolog','s','one',3,'sym8');
xd3=wden(s,'minimaxi','s','one',3,'sym8');
%对上述信号进行图示
subplot(5,1,1);plot(sref);title('参考信号');
subplot(5,1,2);plot(s);title('染噪信号');
subplot(5,1,3);plot(xd1);title('heusure消噪信号');
subplot(5,1,4);plot(xd2);title('sqtwolog消噪信号');
subplot(5,1,5);plot(xd3);title('minimaxi消噪信号');