『壹』 請問如何在matlab中對信號進行去噪操作。最好用函數的形式,方便其他程序調用
去噪有很多種方法,現在小波去噪最為常用,在matlab中有自帶的函數進行小波去噪,簡單易行。最常用且簡單的是閾值去噪,用函數ddencmp()生成信號的默認閾值,然後利用函數wdencmp()進行去噪。如:
%對含噪信號s進行3層小波分解
[c,l]=wavedec(s,3,'db1'); %s為含噪信號,3為小波分解層數,db1為採用的小波基
%獲取信號默認值
[thr,sorh,keepapp]=ddencmp('den','wv',s);
%參數den代表去噪,wv代表小波,s代表含有雜訊的信號
%實現去噪過程
s2=wdencmp('gbl',c,l,'db1',3,thr,sorh,keepapp);
%參數gbl表示採用全局去噪
『貳』 編寫用均值濾波去噪的matlab程序,用兩種方法實現.(重謝)
1、雙循環語句,移動平均法。
<p>雙循環語句,移動平均法</p> <p>%均值濾波</p>
<p>clc,clear;</p>
<p>f=imread('lena.bmp');</p>
<p>subplot(121),imshow(f),title('原圖');</p>
<p>f1=imnoise(f,'gaussian',0.002,0.0008);</p>
<p>%subplot(222),imshow(f1),title('添加高斯雜訊圖');</p>
<p>k1=floor(3/2)+1;</p>
<p>k2=floor(3/2)+1;</p> <p>X=f1;</p>
<p>[M,N]=size(X);</p> <p>uint8Y=zeros(M,N);</p>
<p>funBox=zeros(3,3);</p> <p>fori=1:M-3</p>
<p>forj=1:N-3</p> <p>funBox=X(i:i+3,j:j+3);</p>
<p>s=sum(funBox(:));</p> <p>h=s/9;</p>
<p>Y(i+k1,j+k2)=h;</p> <p>end;</p>
<p>end;</p> <p>Y=Y/255;</p>
<p>subplot(122),imshow(Y),title('均值濾波');</p>
<p>實現圖:</p>
2、filter2。
<p>filter2</p>
<p>clearall;</p>
<p>I=imread('lena.bmp');</p>
<p>%讀入預處理圖像</p>
<p>imshow(I)</p> <p>%顯示預處理圖像</p>
<p>K1=filter2(fspecial('average',3),I)/255;</p> <p>%進行3*3均值濾波</p>
<p>K2=filter2(fspecial('average',5),I)/255;</p> <p>%進行5*5均值濾波</p>
<p>K3=filter2(fspecial('average',7),I)/255;</p> <p>%進行7*7均值濾波</p>
<p>figure,imshow(K1)</p> <p>figure,imshow(K2)</p>
<p>figure,imshow(K3)</p>