關(guān)于小波分析的matlab程序-wenkub
2023-07-01 03:46:33 本頁面
【正文】 re_denoise=conv(low_construct,l_fre_pull)。小波分解的高頻系數(shù)39。)。 %信號高頻細(xì)節(jié)figure(2)。l_fre=conv(y,low_depose)。check1=sum(high_depose)。plot(y)。 %信號賦值f=10。endhigh_construct(1,i_high)=low_depose(1,i_high)*coefficient。 %確定h0(n)，低通分解濾波器for i_high=1:L_fre。h3=cos(a)*sin(b)。 %角度賦初值b=pi/8。%低通重構(gòu)FIR濾波器h0(n)沖激響應(yīng)賦值h0=cos(a)*cos(b)。low_construct=[h0,h1,h2,h3]。 %確定h1(n)=(1)^n,高通重建濾波器if(mod(i_high,2)==0)。endhigh_depose=high_construct(end:1:1)。t=。title(39。 %h0(n)性質(zhì)校驗check2=sum(low_depose)。 %卷積l_fre_down=dyaddown(l_fre)。subplot(2,1,1)plot(l_fre_down)。subplot(2,1,2)。)。h_fre_denoise=conv(high_construct,h_fre_pull)。 %信號重構(gòu)pare=sig_denoisey。y39。ylabel(39。plot(pare)。 %原信號與消噪后信號的比較2． 消噪、% 此函數(shù)用于研究Mallet算法及濾波器設(shè)計% 此函數(shù)用于消噪處理%角度賦值%此處賦值使濾波器系數(shù)恰為db9%分解的高頻系數(shù)采用db9較好,即它的消失矩較大%分解的有用信號小波高頻系數(shù)基本趨于零%對于噪聲信號高頻分解系數(shù)很大，便于閾值消噪處理[l,h]=wfilters(39。)。 %確定h0(n)，低通分解濾波器for i_high=1:L_fre。endhigh_construct(1,i_high)=low_depose(1,i_high)*coefficient。 %原始信號賦值f=10。 %信號加噪聲信號產(chǎn)生zero2=zeros(1,30)。subplot(2,1,1)。)。受噪聲污染的信號39。check3=norm(high_depose)。 %抽取,得低頻細(xì)節(jié)h_fre=conv(y_noise,high_depose)。title(39。plot(h_fre_down)。% 消噪處理for i_decrease=31:44。l_fre_denoise=conv(low_construct,l_fre_pull)。sig_denoise=l_fre_keep+h_fre_keep。end。ylabel(39。plot(sig_denoise)。 %消噪后信號subplot(3,1,3)。)。 % 頻率1f2=100。 % 采樣點數(shù)n=1:N。兩個正弦信號39。)%% h=wfilters(39。)。h39。 % 補零（圓周卷積，且增大分辨率變于觀察）figure(3)。)figure(4)。)%% (圓周卷積的快速傅里葉變換實現(xiàn))sig1=ifft(fft(y).*fft(h))。title(39。分解信號239。分解信號1頻譜39。)%% sig1=dyaddown(sig1)。 % 2插值sig1=sig1(1,[1:N])。 % 重構(gòu)高通hr=circshift(hr39。 % 位置調(diào)整圓周右移一位sig1=ifft(fft(hr).*fft(sig1))。subplot(2,1,1)plot(real(sig1))。subplot(2,1,2)plot(real(sig2))。figure(8)。)。)。linewidth39。legend(39。)title(39。s1=sin(2*pi*50*t*)+sin(2*pi*120*t*)+rand(1,length(t))。s2(t)=sin(2*pi*200*t*)+sin(2*pi*120*t*)+rand(1,length(t))。S139。S239。shannon39。s131=wprcoef(wpt,[3,1])。s135=wprcoef(wpt,[3,5])。s11=norm(s131)。s15=norm(s135)。st11=std(s131)。st15=std(s135)。正常信號的特征向量39。ylabel(39。plot(s131)。subplot(9,2,7)。)。S13339。ylabel(39。plot(s135)。subplot(9,2,15)。)。S13739。,39。s230=wprcoef(wpt,[3,0])。s234=wprcoef(wpt,[3,4])。s20=norm(s230)。s24=norm(s234)。st20=std(s230)。st24=std(s234)。disp(39。plot(s230)。subplot(9,2,6)。)。S23239。ylabel(39。plot(s234)。subplot(9,2,14)。)。S23639。ylabel(39。py1=y1.*conj(y1)/1024。py130=y130.*conj(y130)/1024。py132=y132.*conj(y132)/1024。py134=y134.*conj(y134)/1024。py136=y136.*conj(y136)/1024。py230=y230.*conj(y230)/1024。py232=y232.*conj(y232)/1024。py23
語文相關(guān)推薦
鄂ICP備17016276號-1