用matlab實現4種濾波降噪處理演算法,訊號為正弦訊號疊加上白雜訊,引數可自主定製。
%matlab
clear all;clc;
fs=300;
time=128;
t=0:1/fs:time;
f1=1; t1=1/f1;
y1=sin(2*pi*f1*t);
figure;
plot(t,y1);axis([0
4 -3
3]);
title('原始正弦訊號');
y2=y1+randn(1,length(y1));
figure;
plot(t,y2);axis([0
4 -3
3]);
title('加噪訊號');
n=time/t1-3;%累加中引數最大值不能超過訊號序列否則沒有數值,每個tt都要往後推n-1個週期,準備展示4個週期,故n最大為總週期數-3
t2=t1*fs;
%線性累加平均演算法
a1=linspace(0,4*t1*fs,4*t1*fs-1);%建立序列用來儲存結果,由於是序列序號,故還需乘上fs
for tt=0:4*t1*fs
x1=0;
fori=0:n-1
%進行n次累加用以平均
x1=x1+y2(1+tt+i*t2);
enda1(1+tt)=x1/n;
endfigure;
plot(a1);axis([0
1200 -3
3]);grid;
title('線性累加演算法');
%遞推式平均演算法
a2=linspace(0,length(y2),length(y2)-1);
ap1=linspace(0,n,n-1);
for tt=0:4*t1*fs
ap1(1)=y2(1+tt);
for n=2:n
ap1(n)=(n-1)/n*ap1(n-1)+y2(1+tt+(n-1)*t2)/n;
enda2(1+tt)=ap1(n);
endfigure;
plot(a2);axis([0
1200 -3
3]);grid;
title('遞推式平均演算法');
%指數加權平均演算法
a3=linspace(0,length(y2),length(y2)-1);
ap2=linspace(0,n,n-1);
alpha=30;
beta=(alpha-1)/alpha;
for tt=0:4*t1*fs
ap2(1)=y2(1+tt)*(1-beta);
for n=2:n
ap2(n)=beta*ap2(n-1)+y2(1+tt+(n-1)*t2)*(1-beta);
enda3(1+tt)=ap2(n);
endfigure;
plot(a3);axis([0
1200 -3
3]);grid;
title('指數加權平均演算法');
%五點移動平均演算法
a4=linspace(0,length(y2),length(y2)-1);
for tt=0:4*t1*fs
if tt<=t2
a4(1+tt)=(2*y2(1+tt)+y2(1+tt+t2)+y2(1+tt+2*t2))/4;
elseif t22*t2
a4(1+tt)=(y2(1+tt-t2)+2*y2(1+tt)+y2(1+tt+t2)+y2(1+tt+2*t2))/5;
else
a4(1+tt)=(y2(1+tt-2*t2)+y2(1+tt-t2)+2*y2(1+tt)+y2(1+tt+t2)+y2(1+tt+2*t2))/6;
endendfigure;
plot(a4);axis([0
1200 -3
3]);grid;
title('五點移動平均演算法');
執行後即可得到以下影象:
對比處理前後的影象可見不同演算法對應的去噪效果(不同的取樣點演算法表現是不同的)。
注:演算法公式及降噪原理詳見 高晉佔《微弱訊號檢測》,清華大學出版社。
運動目標檢測的四種方法
2.背景差分法 背景差分法通過輸入影象與背景影象進行比較從而分割出運動目標,在運用背景差分法時需要有一定的限制 要求前景 運動物體 象素的灰度值和背景象素的灰度值存在一定的差別,同時要求攝像機是靜止的改善了連續幀間查分法所存在的問題。背景差分法 background subtraction 是目前運...
js中this的四種使用方法
0x00 js中this的四種呼叫模式 1,方法呼叫模式 2,函式呼叫模式 3,構造器呼叫模式 0x01 第一種 方法呼叫模式 也就是用.呼叫的 this就是這個物件本身 當乙個函式被儲存為物件的乙個屬性時,我們稱它為乙個方法,當乙個方法被呼叫時,this被繫結到該物件,如果乙個呼叫表示式包含乙個屬...
筆記 圓角四種方法的對比以及效能檢測
這篇文章是繼筆記 ios設定圓角方法以及指定位置設圓角文章而寫的,因為上篇文章發出來後,沒有驗證,也有同行的朋友讓我給出一些測試資料來證實一下,所以這裡就給出一下我個人的一些測試資料,正確是否,還請大家作為參考。另外,我寫這個僅僅只是自己作為筆記使用,原來都是放在草稿裡的,但是手機版的無法檢視草稿,...