關于小波分析的matlab程序(已修改)
2025-06-28 03:46 本頁面
【正文】 （個人收集關于小波分析的matlab程序）小波濾波器構造和消噪程序 1小波譜分析mallat算法經(jīng)典程序 7小波包變換分析信號的MATLAB程序 9利用小波變換實現(xiàn)對電能質(zhì)量檢測的算法實現(xiàn) 15基于小波變換的圖象去噪 Normalshrink算法 17小波濾波器構造和消噪程序1． 重構% % 此函數(shù)用于研究Mallet算法及濾波器設計% 此函數(shù)僅用于消噪a=pi/8。 %角度賦初值b=pi/8。%低通重構FIR濾波器h0(n)沖激響應賦值h0=cos(a)*cos(b)。h1=sin(a)*cos(b)。h2=sin(a)*sin(b)。h3=cos(a)*sin(b)。low_construct=[h0,h1,h2,h3]。L_fre=4。 %濾波器長度low_depose=low_construct(end:1:1)。 %確定h0(n)，低通分解濾波器for i_high=1:L_fre。 %確定h1(n)=(1)^n,高通重建濾波器if(mod(i_high,2)==0)。coefficient=1。elsecoefficient=1。endhigh_construct(1,i_high)=low_depose(1,i_high)*coefficient。endhigh_depose=high_construct(end:1:1)。 %高通分解濾波器h1(n)L_signal=100。 %信號長度n=1:L_signal。 %信號賦值f=10。t=。y=10*cos(2*pi*50*n*t).*exp(20*n*t)。figure(1)。plot(y)。title(39。原信號39。)。check1=sum(high_depose)。 %h0(n)性質(zhì)校驗check2=sum(low_depose)。check3=norm(high_depose)。check4=norm(low_depose)。l_fre=conv(y,low_depose)。 %卷積l_fre_down=dyaddown(l_fre)。 %抽取,得低頻細節(jié)h_fre=conv(y,high_depose)。h_fre_down=dyaddown(h_fre)。 %信號高頻細節(jié)figure(2)。subplot(2,1,1)plot(l_fre_down)。title(39。小波分解的低頻系數(shù)39。)。subplot(2,1,2)。plot(h_fre_down)。title(39。小波分解的高頻系數(shù)39。)。l_fre_pull=dyadup(l_fre_down)。 %0差值h_fre_pull=dyadup(h_fre_down)。l_fre_denoise=conv(low_construct,l_fre_pull)。h_fre_denoise=conv(high_construct,h_fre_pull)。l_fre_keep=wkeep(l_fre_denoise,L_signal)。 %取結(jié)果的中心部分,消除卷積影響h_fre_keep=wkeep(h_fre_denoise,L_signal)。sig_denoise=l_fre_keep+h_fre_keep。 %信號重構pare=sig_denoisey。 %與原信號比較figure(3)。subplot(3,1,1)plot(y)。ylabel(39。y39。)。 %原信號subplot(3,1,2)。plot(sig_denoise)。ylabel(39。sig\_denoise39。)。 %重構信號subplot(3,1,3)。plot(pare)。ylabel(39。pare39。)。 %原信號與消噪后信號的比較2． 消噪、% 此函數(shù)用于研究Mallet算法及濾波器設計% 此函數(shù)用于消噪處理%角度賦值%此處賦值使濾波器系數(shù)恰為db9%分解的高頻系數(shù)采用db9較好,即它的消失矩較大%分解的有用信號小波高頻系數(shù)基本趨于零%對于噪聲信號高頻分解系數(shù)很大，便于閾值消噪處理[l,h]=wfilters(39。db1039。,39。d39。)。low_construct=l。L_fre=20。 %濾波器長度low_depose=low_construct(end:1:1)。 %確定h0(n)，低通分解濾波器for i_high=1:L_fre。 %確定h1(n)=(1)^n,高通重建濾波器if(mod(i_high,2)==0)。coefficient=1。elsecoefficient=1。endhigh_construct(1,i_high)=low_depose(1,i_high)*coefficient。endhigh_depose=high_construct(end:1:1)。 %高通分解濾波器h1(n)L_signal=100。 %信號長度n=1:L_signal。 %原始信號賦值f=10。t=。y=10*cos(2*pi*50*n*t).*exp(30*n*t)。zero1=zeros(1,60)。 %信號加噪聲信號產(chǎn)生zero2=zeros(1,30)。noise=[zero1,3*(randn(1,10)),zero2]。y_noise=y+noise。figure(1)。subplot(2,1,1)。plot(y)。title(39。原信號39。)。subplot(2,1,2)。plot(y_noise)。
語文相關推薦
鄂ICP備17016276號-1