3、代碼

%含噪聲的三角波與正弦波的組合

%利用db5小波對信號進行7層分解

%生產(chǎn)正弦信號

clc;close all;clear all;

N=1000;

t=1:N;

sig1=sin(0.3*t);

%生成三角形波形

sig2(1:500)=((1:500)-1)/500;

sig2(501:N)=(1000-(501:1000))/500;

figure(1);

subplot(211);

plot(t，sig1，‘linewidth’，2);

xlabel(‘樣本序號 N’);

ylabel(‘幅值A(chǔ)’);

subplot(212);

plot(t，sig2，‘linewidth’，2);

xlabel(‘樣本序號 N’);

ylabel(‘幅值A(chǔ)’);

%疊加信號

x=sig1+sig2+randn(1，N);

figure(2);

plot(t，x，‘linewidth’，2);

xlabel(‘樣本序號 N’);

ylabel(‘幅值A(chǔ)’);%一維小波分解

［c，l］=wavedec(x，7，‘db5’);%重構(gòu)第1-7層逼近系數(shù)

a7=wrcoef(‘a(chǎn)’，c，l，‘db5’，7);

a6=wrcoef(‘a(chǎn)’，c，l，‘db5’，6);

a5=wrcoef(‘a(chǎn)’，c，l，‘db5’，5);

a4=wrcoef(‘a(chǎn)’，c，l，‘db5’，4);

a3=wrcoef(‘a(chǎn)’，c，l，‘db5’，3);

a2=wrcoef(‘a(chǎn)’，c，l，‘db5’，2);

a1=wrcoef(‘a(chǎn)’，c，l，‘db5’，1);%顯示逼近系數(shù)

figure(3)

subplot(711)

plot(a7，‘linewidth’，2);

ylabel(‘a(chǎn)7’);

subplot(712)

plot(a6，‘linewidth’，2);

ylabel(‘a(chǎn)6’);

subplot(713)

plot(a5，‘linewidth’，2);

ylabel(‘a(chǎn)5’);

subplot(714)

plot(a4，‘linewidth’，2);

ylabel(‘a(chǎn)4’);

subplot(715)

plot(a3，‘linewidth’，2);

ylabel(‘a(chǎn)3’);

subplot(716)

plot(a2，‘linewidth’，2);

ylabel(‘a(chǎn)2’);

subplot(717)

plot(a1，‘linewidth’，2);

ylabel(‘a(chǎn)1’);

xlabel(‘樣本序號 N’);%重構(gòu)第1-7層細節(jié)系數(shù)

d7=wrcoef(‘d’，c，l，‘db5’，7);

d6=wrcoef(‘d’，c，l，‘db5’，6);

d5=wrcoef(‘d’，c，l，‘db5’，5);

d4=wrcoef(‘d’，c，l，‘db5’，4);

d3=wrcoef(‘d’，c，l，‘db5’，3);

d2=wrcoef(‘d’，c，l，‘db5’，2);

d1=wrcoef(‘d’，c，l，‘db5’，1);

%顯示細節(jié)系數(shù)

figure(4)

subplot(711)

plot(d7，‘linewidth’，2);

ylabel(‘d7’);

subplot(712)

plot(d6，‘linewidth’，2);

ylabel(‘d6’);

subplot(713)

plot(d5，‘linewidth’，2);

ylabel(‘d5’);

subplot(714)

plot(d4，‘linewidth’，2);

ylabel(‘d4’);

subplot(715)

plot(d3，‘linewidth’，2);

ylabel(‘d3’);

subplot(716)

plot(d2，‘linewidth’，2);

ylabel(‘d2’);

subplot(717)

plot(d1，‘linewidth’，2);

ylabel(‘d1’);

xlabel(‘樣本序號 N’);

總結(jié)

以上是生活随笔為你收集整理的matlab db5是什么小波,3、代码 - matlab小波分析步骤是什么的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。