%書本實(shí)例《實(shí)用小波分析入門》P114
%
clcclear;
%第三種信號(hào)源
f1?=?100;?
f2?=?500;?
f3?=?400;
Fs?=?2000;
N?=?1000;
t?=?（0:N-1）/Fs;
s?=?0.3*sin（2*pi*f1*t）?+?0.5*sin（2*pi*f2*t）?+?0.2*sin（2*pi*f3*t）?+?0.1*randn（1N）;

subplot（221）plot（s）title（‘原始信號(hào)‘）;
%?2.?利用coif3小波對(duì)信號(hào)進(jìn)行多尺度分解
%定義小波的名稱和分解尺度
w?=?‘db8‘;maxlev?=?4;
[c?l]?=?wavedec（s?maxlev?w）;
newc?=?c;
subplot（222）plot（c）;%axis（[0?1000?-5?5]）title（‘原始信號(hào)在4個(gè)尺度上細(xì)節(jié)分量‘）;

%?3.?將個(gè)尺度上的細(xì)節(jié)分量進(jìn)行修改
%（1）將尺度3和尺度4上的分量變?yōu)榱悖瑢?duì)小波分解后的數(shù)據(jù)進(jìn)行后續(xù)處理
newc?=?wthcoef（‘d‘?c?l?[34]）;
%（2）將尺度1上的部分細(xì)節(jié)分量變?yōu)榱恪４_定尺度1上的細(xì)節(jié)分量的位置。
k?=?maxlev?+?1;
first?=?sum（l（1:k-1））?+1;?last?=?first?+?l（k）?-1;
indd1?=?first:last;

%將尺度1上的其他分量變?yōu)樵瓉淼?/3
newc（indd1）?=?c（indd1）/3;

%確定尺度1上細(xì)節(jié)分量在時(shí)間段[400:600]的細(xì)節(jié)分量
%必須注意的是，尺度k上的細(xì)節(jié)分量的時(shí)間t的氣勢(shì)點(diǎn)
%與原始信號(hào)的時(shí)間t的對(duì)應(yīng)關(guān)系為t/2^k這里k?=?1;
indd1?=?first?+?400/2:first?+?600/2;

%（3）將尺度1在[400:600]時(shí)間段上的細(xì)節(jié)分量變?yōu)榱恪＠胢atlab上的zeros命令將小波分解后的尺度為1上的數(shù)據(jù)中時(shí)間段為[400:600]的細(xì)節(jié)分量變?yōu)榱恪?br/>newc（indd1）?=?zeros（size（indd1））;
%（4）將尺度2在原始信號(hào)上的時(shí)間t?=?500處的上的細(xì)節(jié)分量變?yōu)?.利用比例關(guān)系找到時(shí)間?t?=?500的數(shù)據(jù)在頻率域中的相應(yīng)位置。
k?=?maxlev;first?=?sum（l（1:k-1））+1;
newc（first?+?500/2^2）?=?4;
subplot（224）plot（newc（1:1000））axis（[0?1000?-5?5]）title（‘修改4個(gè)尺度上的細(xì)節(jié)分量‘）;

%?4.?重構(gòu)修改后的信號(hào)
synth?=?waverec（newc?l?w）;
subplot（223）;plot（synth（1:1000））title（‘更新后的信號(hào)‘）;

%我的測(cè)試分析
figure（2）
subplot（211）plot（c）axis（[0?1000?-5?5]）title（‘原始信號(hào)在4個(gè)尺度上細(xì)節(jié)分量‘）;
subplot（212）plot（l）title（‘原始信號(hào)在4個(gè)尺度上近似分量‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1748??2013-04-10?22:07??小波去噪.m

-----------?---------??----------?-----??----

?????????????????1748????????????????????1