數(shù)學(xué)之美:嵌入式編程凹凸性之妙用(附C代碼)
掃描二維碼
隨時隨地手機看文章
[導(dǎo)讀] 咦,你已被成功吸引進來了,不是你想的那樣哈~~~
皮一下哈,言歸正傳,今天遇到一個網(wǎng)友問一個問題,他有一個傳感器測量一個物理量,需要判斷其變化趨勢,我給了一些建議,這里將這個建議展開做些深入分析,并分享給大家。
本文想借此表達(dá)一下個人的一個觀點,做開發(fā)如果遇到無法解決的難題,可以試著從數(shù)序的角度出發(fā),看能否找到答案。
注:文中配圖只為閱讀輕松一點,本人數(shù)學(xué)也是半吊子,有錯誤幫忙指正。
一個項目中用到一個傳感器測量一物理量,這里假定測量溫度吧。需要判斷其變化趨勢,利用這個變化趨勢去做一些應(yīng)用。
那么要怎么判斷一個物理量的變化趨勢呢?我們能自然能想到去求取該隨機序列的變化率。這里涉及到一些數(shù)序定義。
這樣將S(t)信號轉(zhuǎn)換為離散信號序列S(n),那么對于當(dāng)前時刻其斜率怎么求取呢?(這里忽略中間的過度態(tài),僅將其看為線段相連,當(dāng)然現(xiàn)實應(yīng)用中如果有更高要求,可以做曲線擬合)
但是如果只判斷,斜率極容易誤判,比如下面這樣的情況:
則稱函數(shù)f為l上凹函數(shù),有的書上也稱為下凸函數(shù)。
如果把上述條件中的“≥”改成“>”,則叫做嚴(yán)格上凹函數(shù),或叫做嚴(yán)格下凸函數(shù)。
上面是一維函數(shù)情況,這里來個2維函數(shù)的圖,剛方便理解
設(shè)f為定義在區(qū)間I上的函數(shù),若對I上的任意兩點x1
可見,凹凸是相對的,如f(x)在某區(qū)間為凹,則-f(x)則在該區(qū)間為凸。
證明,這里就不推導(dǎo)了,可以利用拉格朗日中值定理可以推導(dǎo)出上面這個性質(zhì)。
來看一下會動的圖,加深一下理解:
函數(shù) 從 到 切線為藍(lán)色,曲線向上凹,綠色表示曲線是向下凹的,紅色表示曲線的拐點。
是否可以利用離散序列的求導(dǎo)數(shù)來判斷傳感器的變化趨勢?
S[n]表示當(dāng)前測量點,S[n-1]表示前一個測量點,S[n-2]表示前第2個測量點。
利用excel生成曲線:
再進一步:
一階導(dǎo)數(shù)與二階導(dǎo)數(shù)結(jié)合起來看,就可以看出測量值變化趨勢的趨勢,比如在前1/4周期,此區(qū)間變換趨勢為增,也即一階導(dǎo)數(shù)為正,而其二階導(dǎo)數(shù)為負(fù),也可以看出遞增的趨勢是逐漸減小到0的。
如果只是做定性判斷,上述函數(shù),完全沒必要與采樣周期做除法,只需要考察其增量即可,代碼可優(yōu)化如下:
這里意外引入一個可能很多人沒注意的知識點NaN,在計算中,NaN代表非數(shù)字,是數(shù)字?jǐn)?shù)據(jù)類型的成員,可以將其解釋為不確定的或無法表示的值,尤其是在浮點運算中。1985年,IEEE 754浮點標(biāo)準(zhǔn)引入了NaN的系統(tǒng)使用,并表示了其他無限量(如無窮大)。
前述函數(shù)返回0x7FBFFFFF,也就是表示無窮大。
不同的操作系統(tǒng)和編程語言可能具有NaN的不同字符串表示形式:
實際上,由于編碼的NaN具有符號,因此通常也可以在NaN的字符串表示中找到它們,例如:
這里給出我的建議方案:
將傳感器信號經(jīng)由電路處理,模數(shù)采樣,在進入前級數(shù)字濾波器,濾除不必要的噪聲,在進行一階/二階求導(dǎo)。對于一階和二階求導(dǎo)再做一級移動平均濾波,最后在按照上面描述進行判別變化趨勢,則個人認(rèn)為基本就比較健壯了。實際移動均值濾波長度不宜選擇過長,否則響應(yīng)就比較滯后了。不能對傳感器的變化趨勢做出實時的判別。加了后級均值濾波器,則會消除由于波形忽上忽下的隨機噪聲干擾影響,使得系統(tǒng)判別更為健壯,實際濾波器長度需根據(jù)不同的場合進行調(diào)試優(yōu)化?;蛘咭部梢赃x擇別的IIR/FIR濾波器形式實現(xiàn)。
具體實現(xiàn)可參考(點擊可閱讀):
手把手教系列之IIR數(shù)字濾波器設(shè)計實現(xiàn)
做為嵌入式er編程,有時候有必要去看看數(shù)學(xué)書,了解一下數(shù)學(xué)原理的背后故事,可能會給你帶來意想不到的作用哦。
本文辛苦原創(chuàng),也是為了幫助網(wǎng)友,如有錯誤之處,也請聯(lián)系我指出。希望大家不吝點個贊,隨心打個賞也是可以的,哈哈。也請幫忙轉(zhuǎn)發(fā)分享支持一下~~~
—END—
是個啥坑?
函數(shù)的凹凸性
凹函數(shù)
凹函數(shù)是一個定義在某個向量空間的凸集C(區(qū)間)上的實值函數(shù)f。設(shè)f為定義在區(qū)間I上的函數(shù),若對I上的任意兩點x1
凸函數(shù)
性質(zhì)
回到坑里
上代碼
#include
代碼優(yōu)化
typedef struct _T_2ND_DRV { float xn1; float xn2;
}t_2ND_DRV; typedef struct _T_1ST_DRV { float xn1;
}t_1ST_DRV; void init_second_derivative(t_2ND_DRV *pSndDrv) {
pSndDrv->xn1 = 0;
pSndDrv->xn2 = 0;
} float second_derivative(t_2ND_DRV *pSndDrv, float xn) { float result=0.0f;
result = xn-2*pSndDrv->xn1-pSndDrv->xn2;
pSndDrv->xn2 = pSndDrv->xn1;
pSndDrv->xn1 = xn; return result;
} void init_fisrt_derivative(t_1ST_DRV *p1stDrv) {
p1stDrv->xn1 = 0;
} float fisrt_derivative(t_1ST_DRV *p1stDrv, float xn) { float result=0.0f;
result = xn-p1stDrv->xn1;
p1stDrv->xn1 = xn; return result;
}
意外收獲
nan
NaN
NaN%
NAN
NaNQ
NaNS
qNaN
sNaN 1.#SNAN 1.#QNAN -1.#IND
-NaN
NaN12345
-sNaN12300
-NaN(s1234)
工程應(yīng)用
總結(jié)一下
免責(zé)聲明:本文內(nèi)容由21ic獲得授權(quán)后發(fā)布,版權(quán)歸原作者所有,本平臺僅提供信息存儲服務(wù)。文章僅代表作者個人觀點,不代表本平臺立場,如有問題,請聯(lián)系我們,謝謝!