文档介绍：小波分析实验报告
姓名：
班级:
学号:二
成绩：
教师签名:
头验一 小波函数的Fourier变换和Fourier逆变换
实验目的
用Matlab实现函数的Fourier变换和Fourier逆变换
实验内容
用Ma;
>> f=sin(2*pi*t)+sin(4*pi*t)+sin(10*pi*t);
>>plot(t,f);
Th— I
r, t- 一
JJI.
f=si n(2nt)+ si n(4nt)+s in (10nt)
flHn
RnHn
An
■5
JI
-111.
Muv
Ann
AMnH
JI
I — I— a
Ann
AnuHMP
」一
Hh-
ahh
■5
0
■5
4
nu
mi J n” —
UHy-
■5
1
Hu
JI
■5
1
实验三：一维离散小波变换
实验目的
实验内容
尺度小波分解重构层数为4并显示各层低频高频图形，加以比较。具体格式如上 例
单尺度：
>> load noissin;
>> s=noissin(1:1000);
>> [cA1,cD1]=dwt(s,'db4');
>> A1=upcoef('a',cA1,'db4',1);
>> D1=upcoef('d',cD1,'db4',1);
>> subplot(4,1,1);plot(s);title('原始信号')
>> subplot(4,1,2);plot(A1);title('低频')
>> subplot(4,1,3);plot(D1);title('高频')
>> sO=idwt(cA1,cD1,'db4');
>> subplot(4,1,4);plot(s0);title('重构信号')
1
0
0
-1 0
2
-2
100 200 300 400 500 600 700 800 900 1000
200
400 重构60信号 800
1000
1200
多尺度:
>> sO=idwt(cAl,cDl,'db4');
>> [C,L]=wavedec(s,4,'db4');
>> cA5=appcoef(C,L,'db4',4);
>> A4=wrcoef('a',C,L,'db4',4);
>> D1=wrcoef('d',C,L,'db4',1);
>> D2=wrcoef('d',C,L,'db4',2);
>> D3=wrcoef('d',C,L,'db4',3);
>>D4=wrcoef('d',C,L,'db4',4)
>> figure(2);
>> subplot(5,1,1);plot(A4);title('第四层低频') >>subplot(5,1,2);plot(D4);title('第四层高频')
>> subplot(5,1,3);plot(D3);title('第三层高频')
>> subplot(5,1,4);plot(D2);title('第二层高频')
>> subplot(5,1,5);plot(D1);title('第一层高频')
>> figure(3);
>> s1=waverec(C,L,'db4');
>> subplot(3,1,1);plot(s);title('原始信号')
>> subplot(3,1,2);plot(s1);title('重构信号')
>> subplot(3,1,3);plot(s-s1);title('误差信号')
第四层低频
100 200 300 40第四层高频600 700
-2」
-」
°
-■- v^- ■•" x'V^'V'v 's/j V ■ ■. iJW' ■■/'1 '\/'ilJ-j
40第三层高频C L
‘1制卜l 7■叫』［怖!伽,
100 200 300 40第 二 层高 频600 700 800 900 10,00
100 200 300 4Q
超‘
300__700__800__900__1000
-"
-」
10
100 200 300 40第一层高
1 J I I I I I I I I I I
"0 100 200