單片機(jī)浮點(diǎn)數(shù)的實(shí)用快速降法
在較為復(fù)雜的單片機(jī)系統(tǒng)中,為擴(kuò)大取值范圍,實(shí)現(xiàn)復(fù)雜的計(jì)算和控制,一般都要涉及浮點(diǎn)數(shù)的運(yùn)算。而一般單片機(jī)是沒有浮點(diǎn)數(shù)運(yùn)算指令的,必須自行編制相應(yīng)軟件。在進(jìn)行除法計(jì)算時(shí),通常使用的方法是比較除法[1],即利用循環(huán)移位和減法操作來得到24~32 位商,效率很低。有些文獻(xiàn)給出了一些改進(jìn)方法[2],但思想不清晰,很難推廣使用。這里給出一種浮點(diǎn)數(shù)除法運(yùn)算的實(shí)用快速算法。該方法以數(shù)值計(jì)算中的預(yù)估 -修正方法為指導(dǎo),充分利用了16位單片機(jī)的乘除法功能,很輕易地實(shí)現(xiàn)了浮點(diǎn)數(shù)的除法。
1 浮點(diǎn)數(shù)格式
IEEE的浮點(diǎn)數(shù)標(biāo)準(zhǔn)規(guī)定了單精度(4字節(jié))、雙精度(8字節(jié))和擴(kuò)展精度(10字節(jié))三種浮點(diǎn)數(shù)的格式。最常用的是單精度浮點(diǎn)數(shù),格式如圖1所示。但是這種格式的階碼不在同一個(gè)字節(jié)單元內(nèi),不易尋址,從而會(huì)影響運(yùn)算速度。
通常在單片機(jī)上采用的是一種變形格式的浮點(diǎn)數(shù),如圖2所示。其中的23位尾數(shù)加上隱含的最高位1,構(gòu)成一個(gè)定點(diǎn)原碼小數(shù),即尾數(shù)為小于1大于等于0.5的小數(shù)。有關(guān)浮點(diǎn)數(shù)格式的詳細(xì)內(nèi)容請參考有關(guān)文獻(xiàn)[1][2]。
2 快速除法的算法原理
在16 位單片機(jī)中只有16位的乘除法,而浮點(diǎn)數(shù)的精度(即尾數(shù)的有效位數(shù))達(dá)24位,因此無法直接相除,但依然可以利用16位的乘除法指令來實(shí)現(xiàn)24位除法。不過,如果只進(jìn)行一次16位的除法必定會(huì)帶來很大誤差,因此問題的關(guān)鍵在于如何消除這個(gè)誤差,從而達(dá)到要求的精度。這其實(shí)就是通常數(shù)值計(jì)算中所采用的預(yù)估- 修正方法。
假設(shè)兩個(gè)浮點(diǎn)數(shù)經(jīng)過預(yù)處理后,被除數(shù)和除數(shù)尾數(shù)擴(kuò)展為32位(末8位為0)分別放入X和Y中。鄰YL為Y的低16位,并記YH=Y-YL。顯然YH≈Y,X/Y與Y/YH相差不多:
(X/Y)/(X/YH)=(YH/Y)
=YH/(YH+YL)
=1/(1+YL/YH)
≈1-YL/YH
=(YH-YL)/YL
評論