關(guān)于小波分析的matlab程序-文庫吧
2025-06-01 03:46 本頁面
【正文】 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_noise,low_depose)。 %卷積l_fre_down=dyaddown(l_fre)。 %抽取,得低頻細(xì)節(jié)h_fre=conv(y_noise,high_depose)。h_fre_down=dyaddown(h_fre)。 %信號高頻細(xì)節(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。)。% 消噪處理for i_decrease=31:44。if abs(h_fre_down(1,i_decrease))=h_fre_down(1,i_decrease)=(10^7)。endendl_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。 %消噪后信號重構(gòu)%平滑處理for j=1:2for i=60:70。sig_denoise(i)=sig_denoise(i2)+sig_denoise(i+2)/2。end。end。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。)。 %原信號與消噪后信號的比較小波譜分析mallat算法經(jīng)典程序clc。clear。%% f1=50。 % 頻率1f2=100。 % 頻率2fs=2*(f1+f2)。 % 采樣頻率Ts=1/fs。 % 采樣間隔N=120。 % 采樣點(diǎn)數(shù)n=1:N。y=sin(2*pi*f1*n*Ts)+sin(2*pi*f2*n*Ts)。 % 正弦波混合figure(1)plot(y)。title(39。兩個正弦信號39。)figure(2)stem(abs(fft(y)))。title(39。兩信號頻譜39。)%% h=wfilters(39。db3039。,39。l39。)。 % 低通g=wfilters(39。db3039。,39。h39。)。 % 高通h=[h,zeros(1,Nlength(h))]。 % 補(bǔ)零（圓周卷積，且增大分辨率變于觀察）g=[g,zeros(1,Nlength(g))]。 % 補(bǔ)零（圓周卷積，且增大分辨率變于觀察）figure(3)。stem(abs(fft(h)))。title(39。低通濾波器圖39。)figure(4)。stem(abs(fft(g)))。title(39。高通濾波器圖39。)%% (圓周卷積的快速傅里葉變換實現(xiàn))sig1=ifft(fft(y).*fft(h))。 % 低通(低頻分量)sig2=ifft(fft(y).*fft(g))。 % 高通(高頻分量)figure(5)。 % 信號圖subplot(2,1,1)plot(real(sig1))。title(39。分解信號139。)subplot(2,1,2)plot(real(sig2))。title(39。分解信號239。)figure(6)。 % 頻譜圖subplot(2,1,1)stem(abs(fft(sig1)))。title(39。分解信號1頻譜39。)subplot(2,1,2)stem(abs(fft(sig2)))。title(39。分解信號2頻譜39。)%% sig1=dyaddown(sig1)。 % 2抽取sig2=dyaddown(sig2)。 % 2抽取sig1=dyadup(sig1)。 % 2插值sig2=dyadup(sig2)。 % 2插值sig1=sig1(1,[1:N])。 % 去掉最后一個零sig2=sig2(1,[1:N])。 % 去掉最后一個零hr=h(end:1:1)。 % 重構(gòu)低通gr=g(end:1:1)。 % 重構(gòu)高通hr=circshift(hr39。,1)39。 % 位置調(diào)整圓周右移一位gr=circshift(gr39。,1)39。 % 位置調(diào)整圓周右移一位sig1=ifft(fft(hr).*fft(sig1))。 % 低頻sig2=ifft(fft(gr).*fft(sig2))。 % 高頻sig=sig1+sig2。 % 源信號%% figure(7)。subplot(2,1,1)plot(real(sig1))。title(39。重構(gòu)低頻信號39。)。subplot(2,1,2)plot(real(sig2))。title(39
語文相關(guān)推薦
鄂ICP備17016276號-1