一維小波變換演算法程式目錄
阿新 • • 發佈：2019-01-25
[C,L]=wavedec(s,5,'db4');
cA5=appcoef(C,L,'db4',5);
A5=wrcoef('a',C,L,'db4',3);
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);
D5=wrcoef('d',C,L,'db4',5);
subplot(3,2,1);plot(A5);title('低頻A5')
subplot(3,2,2);plot(D1);title('高頻D1')
subplot(3,2,3);plot(D2);title('高頻D2')
subplot(3,2,4);plot(D3);title('高頻D3')
subplot(3,2,5);plot(D3);title('高頻D4')
subplot(3,2,6);plot(D3);title('高頻D5')
figure(3)
A0=waverec(C,L,'db4');
subplot(3,1,1);plot(s);title('原始訊號')
subplot(3,1,2);plot(A0);title('重構訊號')
subplot(3,1,3);plot(s-A0);title('誤差訊號')
err = max(abs(s-A0))
＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃
cA5=appcoef(C,L,'db4',5);
A5=wrcoef('a',C,L,'db4',3);
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);
D5=wrcoef('d',C,L,'db4',5);
subplot(3,2,1);plot(A5);title('低頻A5')
subplot(3,2,2);plot(D1);title('高頻D1')
subplot(3,2,3);plot(D2);title('高頻D2')
subplot(3,2,4);plot(D3);title('高頻D3')
subplot(3,2,5);plot(D3);title('高頻D4')
subplot(3,2,6);plot(D3);title('高頻D5')
figure(3)
A0=waverec(C,L,'db4');
subplot(3,1,1);plot(s);title('原始訊號')
subplot(3,1,2);plot(A0);title('重構訊號')
subplot(3,1,3);plot(s-A0);title('誤差訊號')
err = max(abs(s-A0))
＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃