滑動(dòng)平均系統(tǒng)的去噪處理

1、信號(hào)與系統(tǒng)實(shí)驗(yàn)大作業(yè)滑動(dòng)平均系統(tǒng)的去噪處理題目：(1)畫出(k=1:50)原始彳S號(hào)fk波形圖，噪聲干擾信號(hào)dk及加噪后的信號(hào)時(shí)域圖。(其中，dk可以下列語(yǔ)句產(chǎn)生：d=2*(rand(1,R)(2)畫出去噪后yk(k=1:50)的信號(hào)時(shí)域圖(3)討論M=2時(shí)，去噪效果如何內(nèi)容：設(shè)輸入信號(hào)為受噪聲干擾的信號(hào)為fk=sk+dk,其中sk2*k*sin(k/400)是原始信號(hào)，dk是噪聲。已知M點(diǎn)滑動(dòng)平均系統(tǒng)的輸入輸出關(guān)系為:1M1ykMn0fkn試?yán)肕ATLAB編程實(shí)現(xiàn)(k=1:50)用M=5點(diǎn)滑動(dòng)平均系統(tǒng)對(duì)受噪聲干擾的信號(hào)去噪。程序：%program8_3Signalsmoothingbymo

2、vingaveragefilterR=51;%lengthofinputsignal%generate,uniformlydistributedrandomnumberd=rand(1,R);k=0:R-1;s=2*k.*sin(k.*pi./400);f=s+d;figure(1);plot(k,d,'r-.',k,s,'b:',k,f,'k-');xlabel('Timeindexk');legend('dk','sk','fk');M=5;b=ones(M,1)/M;a=1;y

3、=filter(b,a,f);figure(2);plot(k,s,'b:',k,y,'r-');xlabel('Timeindexk');legend('sk','yk');M=2;b=ones(M,1)/M;a=1;y=filter(b,a,f);figure(3);plot(k,s,'b:',k,y,'r-');xlabel('Timeindexk');legend('sk','yk');圖形：題目1圖（a）題目2-&k4035-30-252015-1Q-5-Q-5-a5101520253035404550Timeindexk圖（b）圖（c）題目3分析：1 .圖（a）中三條曲線分別為噪聲信號(hào)dk、原始信號(hào)sk和受噪聲干擾的輸入信號(hào)fk。2 .圖（b）中sk為有用信號(hào)，yk是經(jīng)過(guò)5點(diǎn)滑動(dòng)平均系統(tǒng)去噪的結(jié)果。比較這兩條曲線可以看出，yk基本將所有噪聲去除，但有一部分有用信號(hào)也被去除，曲線雖然很光滑，但與標(biāo)準(zhǔn)曲線偏離較大。3 圖（C）是經(jīng)過(guò)2點(diǎn)滑動(dòng)平均系統(tǒng)去噪的結(jié)果，比較兩條曲線可以看出，yk不但去除了大部分的噪聲信號(hào)，而且很好地保留了有用信號(hào)，曲線與標(biāo)準(zhǔn)曲線很接近。4 所以綜上可以得出，當(dāng)M取值越