D08HR_2數(shù)據(jù)采集器技術(shù)說明書

1、37DHR_2電法勘探數(shù)據(jù)采集器技術(shù)說明書 DHR_2電法勘探數(shù)據(jù)采集器 技術(shù)說明書 研制者 陸煥文 2005年9月3目 (一). 電法勘探數(shù)據(jù)采集器外形照片(二). 儀器的用途 本儀器是一臺便攜式高速高靈敏度的模擬量數(shù)據(jù)采集器. 它可以單獨使用, 也可以在RS485二總線通信綱絡(luò)中作為主機和從機使用.作為單機使用時, 它可以作模擬電壓波形放大, 自動零點補償, 自動放大增益控制, 模擬量無源和有源濾波, 模/數(shù)轉(zhuǎn)換, 數(shù)據(jù)存儲, 數(shù)字信號實時處理, 作數(shù)字濾波器,DFT離散富里葉變換進行頻譜分拆得到被采集信號的振幅和相位. 儀器內(nèi)裝現(xiàn)場數(shù)據(jù)庫, 可以對現(xiàn)傷采集到的原始記錄存檔, 存檔后的數(shù)據(jù)

2、永久保存(儀器斷電也不會消失), 數(shù)據(jù)庫內(nèi)數(shù)據(jù)可以隨時查詢顯示(可以以數(shù)據(jù)表格方式顯示, 也可以以圖形曲線方式顯示). 通過機上的RS232串行通信接口把數(shù)據(jù)庫的全部記錄發(fā)送到PC系統(tǒng)機, 在PC機方用VC+偏寫的數(shù)據(jù)采集應(yīng)用工程軟件對接收到的數(shù)據(jù)采集器的數(shù)據(jù)庫全部數(shù)據(jù)進行整理存檔(存入硬磁盤內(nèi)).PC機內(nèi)的數(shù)據(jù)也可以與通用的數(shù)據(jù)庫軟件Microsoft office Access 關(guān)聯(lián). 將數(shù)據(jù)采集器的記錄存放到通用數(shù)據(jù)庫中.儀器應(yīng)用領(lǐng)域:(1). 地質(zhì)物理勘探.l 高密度直流大地電阻率法勘探: 用於尋找地下水, 工程地質(zhì)勘探中地層剖面結(jié)構(gòu)勘察(鐵路, 公路路基勘察).l 高密度時間域和頻率

3、域激發(fā)極化法電法勘探: 用於尋找有色金屬礦(銅礦, 鉛鋅礦等), 尋找地下水.l 可控源磁大地電流法CSMT: 在該方法中作野外從機使用;采集磁探頭上來的大地電流和來自不極化電極上來的電壓信號.(2). 實驗室測試儀.l 自動量程控制的數(shù)字電壓表.l 瞬態(tài)波形采集器與波形顯示器.l 頻率和周期測量儀.l 振動信號和聲音信號采集, 頻譜分拆, 幅值和相位的波德圖顯示.( 可用作機械振動測量和監(jiān)控, 音晌器件喇叭和蜂鳴器質(zhì)量檢驗)(三). 儀器的主要技術(shù)指標.(1). 儀器的量程:+ -10微伏到 + - 5000 毫伏.(2). 儀器的最高分辨率: 10 微伏.(3). 儀器的精確度: 1/10

4、00.(4). 儀器的工作壞境: 野外條件, 氣溫在零下10攝氏度到40攝氏度. 空氣濕度低於80%相對濕度. 能經(jīng)受長途運輸和野外使用中的振動.(5). 儀器用RS485串行二總線通信最大距離為1000米, 從機可尋地址的數(shù)量為128臺, 基本通信從機臺數(shù)為32臺, 超過32臺時要添加中繼器.(6). 儀器的重量和體積: 約0.5公斤,20厘米*12厘米*11厘米.(7). 儀器的備電使用指標: 儀器內(nèi)裝可充電的5號鎳氫電池14節(jié), 串聯(lián)使用, 總電壓17伏, 儀正常工作電流約120毫安, 在充足電池容量(1200MA/H) 后可速續(xù)工作8小時. 當電池用到輸出電壓降到15伏時, 要進行充電

5、, 要使用儀器配套提供的充電器連續(xù)充電8小時. 超過8小時后充電器會自動停止充電. 儀器在使用時也可補充電, 但是儀器測量指標會降低, 這是50HZ工業(yè)電干擾信號會進入電池影響高靈度放大器的工作狀態(tài)造成的, 所以提倡儀器工作時不要充電.(四). 儀器硬件電路和軟件編制說明.1 儀器硬件方框圖Vin零點偏移與輸入的直流偏移電壓(自然電位) 補償器無源50HZ陷波器緩沖放大器(增益=1, 同相放大器)位控輸出總線GND可控增益放大器增益可變范圍 1,2,4,8,16有源低通濾波器DC_10KHZ可控增益放大器增益可變范圍 1,2,4,8,16 高8位字節(jié)緩沖輸入接口(A/D的高4位數(shù)據(jù)偏移碼在低4

6、位, 高4位=0000)微控制器CPU板W77E58 CPU32KBYTERAM 62256128KBYTEFLASH MEMORYAM29F010B實時鐘DS12C88764KBYTE EPROM27C512模數(shù)轉(zhuǎn)換器MAX120分辨率12BIT采速率 500KHZ(2微秒)A/D的低8位緩沖輸入接口VCC+5VGND0VDC/DC變換器DC/DC變換器RS485通信口RS232通信口+8.5v輕觸鍵盤5*4健1,2,3,4,5,6,7,8,9,0+,-,”.”UP,ESC,RSTLEFT,RHIT,DONW,ENT點陣式液晶顯示器LG128642128*64點陣圖形顯示0v鎳氫電池14節(jié)5

7、號1200mah-8.5v2 無源50HZ陷波器電路說明.參考下圖.示波器50HZ信號發(fā)生器圖中R29(10K) 和結(jié)型場效應(yīng)晶體管J2和J3(3DJ6F) 組成了輸入信號幅度超過最大允許電壓時對本儀器的輸入放大器保護, 它把輸入信號限幅到+8.5V到-8.5V之間. 場效應(yīng)晶體管在此作反向漏電流極小的高頻二極管使用, 在輸入信號在正常范圍時J2,J3不會產(chǎn)生反向漏電流而改變輸入信號的幅度.( 在要求不高時也可用快速二極管代替).從R30_R33和C11_C24和P2_P3電位器組成了一個雙T電橋, 它可以將輸入信號中的50HZ工業(yè)電干擾澡聲波陷掉. 在調(diào)整好后可以使50HZ波衰減到原來的1/

8、100. 具體調(diào)整是調(diào)節(jié)P2和P3兩個電位器. 讓儀器輸入50HZ正弦波信號它的幅值是1000mv, 用示波器測量運算放大器的輸出, 使用示波器的高靈敏度檔測量波形, 反復(fù)調(diào)節(jié)P2和P3電位器使得示波器上觀察到的波形幅度小於10mv.3. 自動調(diào)零和輸入自然電位補償器電路工作原理與程序設(shè)計. l 參看上圖說明電路工作原理。圖IC9:X9312WP是一片數(shù)字電位器,VH到VL腳之間的電阻是10KR,VW腳是電位器的中間滑動端. 總共可滑動100步. 可向上和向下滑動, 方向由來自單片機的位輸出口OX1_6控制, 當OX1_6=”1”(5V) 可以使VW向上滑動.OX1_6=”0” 時可以向下滑動

9、. 具體滑動一步要用OX1_7輸出一個負脈沖(即”1”,”0”,”1”). 當VW中間滑動端向上或向下滑動到端點時(即VH,VL), 再給負脈沖就不再滑動. 跟據(jù)以上所述, 要想使VW滑到端點只要用OX1_7發(fā)出100個負脈沖,不管原來VW在什么位置. 到了端點后就可以用OX1_7發(fā)出N個脈沖就使VW到達要求位置N.電路IC5(OP_07) 是輸入緩沖放大器, 它的增益=1,IC015(OP_07) 是執(zhí)行零點和自電補償?shù)臏p法器. 相減點是運放的”+” 輸入腳, 減法器的輸出端是運放的6腳, 運放的增益是-1. 當要求在輸入信號中減去直流分量SPmv時, 只要在運放IC015:3腳提供SP/2

10、mv的補償電壓. 這時運放就會輸出0v.補償電壓SP/2是由數(shù)字電位器和一些電阻提供的, 通過調(diào)節(jié)P1和P2電位器, 使VSP電位跟隨數(shù)字電位器的VW滑動點向下滑動而變化. 最后達到SP在-100mv到+100mv范圍內(nèi)可補償.具體調(diào)試步驟如下:使信號發(fā)生器輸出0mv, 用儀器(system test) 系統(tǒng)自檢軟件操作使數(shù)字電位器向下滑100步到達端點VL, 這時可用數(shù)字電在表測VW腳, 應(yīng)該=0mv.用儀器(system test) 系統(tǒng)自檢軟件操作使數(shù)字電位器向上滑50步.調(diào)節(jié)P2電位器, 用數(shù)字電壓表測運放IC015的6腳, 使它 = 0mv.用儀器(system test) 系統(tǒng)自檢

11、軟件操作使數(shù)字電位器向下滑100步.調(diào)節(jié)P1電位器, 用數(shù)字電壓表測運放IC015的6腳, 使它+100mv.用儀器(system test) 系統(tǒng)自檢軟件操作使數(shù)字電位器向上滑100步.用數(shù)字電壓表測運放IC015的6腳, 它應(yīng)該 = -100mv.再重復(fù)一次以上各步驟.6儀器的電源電路原理與調(diào)試。(7)RS485串行二總線多機通訊硬件與軟件說明.本儀器是一套野外局域計算機綱絡(luò)中的一個從機.它可以成為32臺從機中的某一號從機.從機的偏號可用本機的參數(shù)設(shè)置中設(shè)本機地址操作完成.設(shè)置的地址范圍是1_32.地址一旦設(shè)定,本儀器中的串行EEPROM(24LC65)中就保存好這個地址,儀器停電也不會消

12、失.要想改變必須重新設(shè)置.本套儀器的主機是D03T電法勘探發(fā)送機.一臺主機可以用1根RS485總線與多達32臺從機連接.主機可與任意1臺從機進行半雙工通(一臺機同時只能處在發(fā)送或接收兩種狀態(tài)之一).RS485二總線是差動輸出或輸入雙向接口.輸入/輸出線有三根線:RA,RB,GND.這三根線用9蕊D型扦座在儀器面板上引出.對應(yīng)關(guān)系是:1腳是RA,4腳是RB,5腳是GND.*上圖是D05R數(shù)據(jù)采集器中隔離的RS485串行通訊接口電路.IC1是單片機W77E58,它有兩個串行接口SBUF和SBUF1,本儀器中用SBUF1作為RS485通訊接口.RXD1(P1.2)TXD1(P1.3)串行輸入和輸出線

13、經(jīng)過6N137高速光電耦合器隔離后連接到RS485接口蕊片MAX487E的1腳和4腳.RS485的方向控制是由W77E58的P1.6控制的.P1.6=0:RS485處於接收狀態(tài),P1.6=1:RS485處於發(fā)送狀態(tài).本電路的隔離電源是由DC/DC變換器電路提供,SG3525 PWM芯片把+8V_-8V電源逆變成20KHZ的方波.再用高壓隔離高頻變壓器輸出7V左右的方波,再由全波整流電路輸出7V左右的直流電壓.再經(jīng)750L05低壓降穩(wěn)壓片穩(wěn)成5V的浮空隔離電源.提供RS485接口蕊片MAX487E和光電耦合器6N137使用.*離散富里葉變換(DFT)在地質(zhì)儀器中應(yīng)用 對于大地系統(tǒng)的分析，電法勘探

14、中主要是用電流作為輸入信號，電壓作為輸出信號，大地是信號傳遞系統(tǒng)，如下圖所示： Vt發(fā)送機輸出電流時域信號方波I (t)頻域信號I (f)接收機獲得信號時域，準方波大地系統(tǒng)時域脈沖響應(yīng)函數(shù) h (t)頻域傳遞函數(shù)H (f) I If Vf fo f fo f 圖（1） 上圖中的It是時域電流的波形，If是頻域中電流的頻譜。Vt是時域的電壓波形，Vf是頻域中電壓的頻譜。 以上時域信號都是可以測量獲得的，測量的辦法時對電流和電壓的波形進行采樣和模數(shù)轉(zhuǎn)換（A/D），如下圖所示： (I) t h (K0) h (K1023) (V) t K0 K1023 t T N 圖（2） 在電法勘探中使用的波形是

15、連續(xù)的周期性的穩(wěn)態(tài)波形。這是為了用頻率域分析信號的要求，要求周期性是為了用計算機使用的離散采樣信號來分析物理的連續(xù)信號，用傅立葉級數(shù)計算來代替無窮時間的連續(xù)積分。 對于時間域的信號It輸入到大地的系統(tǒng)，經(jīng)過大地系統(tǒng)的作用輸出的時間域信號Vt。Vt與It之間的關(guān)系是： Vt = It（卷積）h(t) (1)這是一個積分方程，很難用計算機計算處理。 而在頻率域中Vf和If的關(guān)系是一個簡單的代數(shù)式： Vf = If * H（f） （2） 我們的目的是要求得大地傳遞函數(shù)H（f），即大地的電阻率頻譜函數(shù)，這里的H（f）傳遞函數(shù)是一個復(fù)數(shù)，它可以寫成用振幅 | H（f）| 和相位R 表示的e指數(shù)形式： j

16、*R H (f) = | H（f）| * e （3） 上式中j是虛部-1 利用公式（2）和（3）可得到下式： j*V Vf |Vf | * e |Vf | j*(V -I ) H (f) = = = * e (4) j*IIf |If | * e |If |上式中 |Vf | 是輸入電壓的頻譜幅值，V是輸入電壓的相位 |If |是輸入電流的頻譜幅值，I是輸入電壓的相位利用公式（3）和公式（4）可得到下式：j*R |Vf | j*(V -I ) | H（f）| * e = * e （5） |If |根據(jù)（5）式我們可得到大地系統(tǒng)的復(fù)電阻幅值： |Vf | | H（f）| = （6） |If |大

17、地系統(tǒng)的相位： R =V I （7）到此我們意義竟清楚地看出，在頻率域中求大地系統(tǒng)的復(fù)電阻的幅值和相位用到的計算方式都是簡單的代數(shù)方法，計算機很容易計算處理。但問題的關(guān)鍵是我們?nèi)绾吻蟪鰘Vf |，V，|If |，I這些量。這正是我們要在下文中討論的DFT離散傅立葉變換算法。用高速A/D對電流和電壓的時域信號進行一個完整周期的采樣，然后分別對兩個數(shù)據(jù)序列進行離散傅立葉變換，就能用加法和乘法來算出|Vf |，V，|If |，I 為了說明離散的傅立葉變換，我們先要從速讀的傅立葉變換說起，傅立葉積分定義為： -j2ftH（f）= h(t) * e dt (8) -上式中H（f）是頻率域內(nèi)的頻譜函數(shù) h

18、(t)是時間域內(nèi)的時間波形信號函數(shù)因為 -j2fte = cos(2ft)- jsin(2ft) (9)用（8）式和（9）式可得到以三角函數(shù)表達的傅立葉變換式： H（f）= h(t) * cos(2ft)- jsin(2ft) dt （10） - 對于周期性的時域信號h(t)可以用下面傅立葉級數(shù)形式來代替上式： N-1 H（f n）= h(k) * cos(2f n * k)- jsin(2f n * k) K=0 N-1 N-1= h(k) * cos(2f n * k)- j * sin(2f n * k) K=0 K=0= Af n j* B f n (11) 上式中Af n是頻譜函數(shù)H

19、（f n）的f n諧波分量的實數(shù)部分系數(shù) B f n是頻譜函數(shù)H（f n）的f n諧波分量的虛數(shù)部分系數(shù) 如果f n是基波，則f 0=1，上式變得更加簡單。我們在地質(zhì)儀器中為了提高精度只對方波的基波進行傅立葉變換，實踐證明對基波采1024個樣，A/D轉(zhuǎn)換分辨達到12bit時，可以得到幅度和相位分辨率0.1%，相位精度在2個毫弧度左右。 我們得到對方波信號的基波離散傅立葉變換的公式如下：1023 2 2H（f0）= h(k) * cos( * k)- j * h(k)sin( * k) （12） K=0 1024 1024 參看圖（2）上式中，h(k)是波形中第K個采樣點的A/D轉(zhuǎn)換得的數(shù)據(jù)，它

20、是一個二字節(jié) 2的整數(shù)，讓 =，上式變?yōu)椋?10241023 H（f0）= h(k) * cos( * k)- j * h(k)sin( * k) K=0 =A f0 j* B f0 （13）上式中的 1023 A f0 = h(k) * cos( * k) （14） K=0 1023 B f0 = h(k) * sin( * k) （15） K=0 A f0是頻譜函數(shù)的實部，B f0是頻譜函數(shù)的虛部 用e指數(shù)形式，（13）式可寫成 B f0 jarctg 2 2 A f0H（f0）= （A f0） + （B f0） * e j*f0 = | H（f0）| * e (16) 2 2 上式中 |

21、H（f0）|= （A f0） + （B f0） （17） 它是基波f0的幅值 B f0 f0 = arctg A f0 這是基波f0 的相位。注意這個相位是隨機的，但是對電流I波形和電壓V波形進行同步采樣，用傅立葉變換求出的I電流相位）和V（電壓相位），雖然它們都是隨機的，隨機的值取決與開始采樣的那時刻，但V I = 電壓對電流的相移是固定的，它是由大地系統(tǒng)的傳遞函數(shù)決定的。它就是我們要測量的物理參數(shù)，它就是復(fù)電阻率的相位。綜上所述，要用傅立葉變換對大地系統(tǒng)測量頻譜特性的方法步驟如下： 1：讓發(fā)送機想地下供電，發(fā)送機輸出的電流波形是方波，頻率是精確的（1/10萬 精確度），電流的幅度不要求很精

22、確（5%精確度），用接收機中A/D轉(zhuǎn)換器對方波采樣1024個數(shù)據(jù)，A/D的分辨率為12bit。為了使放大器部分與發(fā)送機部分的輸出端上的高壓隔離，需要一個隔離電壓在2500V以上，頻率響應(yīng)在100KHZ以上的隔離放大器。發(fā)送機的電流值用0.1的精密電阻（自制的錳銅絲電阻）轉(zhuǎn)換為電壓值送給隔離變壓器，隔離放大器輸出的電壓與發(fā)送機輸出端高壓隔離，剩下的只是I*0.1的電壓VI值。這個電壓就是發(fā)送機輸出的電流波形，如下圖所示：高壓隔離變壓器放大器A/D發(fā)送機 A微處理器 B 接收機 2:接收機對發(fā)送機輸出電流波形采樣的同時,同步地對來自大地某測點的V波形進行采樣1024個點.A/D的分辨率是12bit

23、,為了提高測量精度對來自大地M,N接收電極上的信號要進行放大,放大到A/D能接收的最高電壓的一半以上.放大器的放大倍數(shù)是可以用計算機程序控制的(1100),這樣就形成了一個浮點A/D轉(zhuǎn)換器,它的分辨率是12bit (1/4096),而電壓測量的范圍是滿度滿度值: + 50mV + 5V最高分辨率: 50 mV /2048 = 24.4 25V 對于野外信號有25V的分辨率是足夠了,再高會被電極的極差漂移,電線上的干擾,放大器的噪聲淹沒. 3:分別用離散傅立葉變換酸楚電流波形的頻譜幅值和相位得到: jII(f0)= |I(f0)| * e jVV(f0)= |V(f0)| * e 4:把電壓的幅

24、值除以電流的幅值得到地質(zhì)體的復(fù)電阻幅值. |V(f0)| R(f0)= |I(f0)| 5:把電壓的相位減去電流的相位,得到地質(zhì)體的復(fù)電阻的相位 R = V I 6:把復(fù)電阻幅值乘上野外布極的裝置系數(shù)K值,得到地質(zhì)的視復(fù)電阻率s = K * | R(f0)| 7:改變發(fā)送機的輸出電流的頻率后再重復(fù)以上6個步驟,這樣就可以得到地質(zhì)體對不同頻率的響應(yīng),得到完整的地質(zhì)體復(fù)電阻率的幅值隨頻率而變化的相位圖,如下圖所示: H, 幅值曲線 相位曲線 f 該圖可在野外接收機上顯示出來 8:以上步驟還缺少一個儀器自身放大器的頻譜特性對測量帶來的影響問題,本文用另一篇文章來說明校正方法./*/* 離散富里變換程

25、序 */*/*DFT()函數(shù)中計算sin(M*2PI/1024)使用查表法*/*存放在CORD(ROM)中的sin()表*/static float code sin_table_float260= 0.0000000, 0.0061359, 0.0122715, 0.0184067, 0.0245412, 0.0306748, 0.0368072, 0.0429383, 0.0490677, 0.0551952, 0.0613207, 0.0674439, 0.0735646, 0.0796824, 0.0857973, 0.0919090, 0.0980171, 0.1041216, 0.

26、1102222, 0.1163186, 0.1224107, 0.1284981, 0.1345807, 0.1406582, 0.1467305, 0.1527972, 0.1588582, 0.1649131, 0.1709619, 0.1770042, 0.1830399, 0.1890687, 0.1950903, 0.2011046, 0.2071114, 0.2131103, 0.2191012, 0.2250839, 0.2310581, 0.2370236, 0.2429802, 0.2489276, 0.2548656, 0.2607941, 0.2667128, 0.272

27、6214, 0.2785197, 0.2844075, 0.2902847, 0.2961509, 0.3020059, 0.3078496, 0.3136818, 0.3195020, 0.3253103, 0.3311063, 0.3368899, 0.3426607, 0.3484187, 0.3541635, 0.3598951, 0.3656130, 0.3713172, 0.3770074, 0.3826834, 0.3883450, 0.3939920, 0.3996242, 0.4052413, 0.4108432, 0.4164295, 0.4220003, 0.427555

28、1, 0.4330938, 0.4386162, 0.4441221, 0.4496113, 0.4550836, 0.4605387, 0.4659765, 0.4713967, 0.4767992, 0.4821838, 0.4875502, 0.4928982, 0.4982277, 0.5035384, 0.5088301, 0.5141028, 0.5193560, 0.5245897, 0.5298036, 0.5349976, 0.5401714, 0.5453250, 0.5504580, 0.5555702, 0.5606616, 0.5657318, 0.5707808,

29、0.5758082, 0.5808139, 0.5857978, 0.5907597, 0.5956993, 0.6006165, 0.6055111, 0.6103828, 0.6152316, 0.6200572, 0.6248595, 0.6296383, 0.6343933, 0.6391245, 0.6438316, 0.6485144, 0.6531729, 0.6578067, 0.6624158, 0.6669999, 0.6715590, 0.6760927, 0.6806010, 0.6850837, 0.6895406, 0.6939715, 0.6983762, 0.7

30、027547, 0.7071068, 0.7114322, 0.7157308, 0.7200025, 0.7242471, 0.7284644, 0.7326543, 0.7368166, 0.7409511, 0.7450578, 0.7491364, 0.7531868, 0.7572088, 0.7612024, 0.7651672, 0.7691033, 0.7730104, 0.7768885, 0.7807372, 0.7845566, 0.7883464, 0.7921066, 0.7958369, 0.7995372, 0.8032075, 0.8068476, 0.8104

31、572, 0.8140363, 0.8175848, 0.8211025, 0.8245893, 0.8280451, 0.8314696, 0.8348629, 0.8382247, 0.8415550, 0.8448536, 0.8481203, 0.8513552, 0.8545580, 0.8577286, 0.8608670, 0.8639728, 0.8670462, 0.8700870, 0.8730950, 0.8760701, 0.8790122, 0.8819213, 0.8847971, 0.8876396, 0.8904487, 0.8932243, 0.8959662

32、, 0.8986745, 0.9013488, 0.9039893, 0.9065957, 0.9091680, 0.9117060, 0.9142098, 0.9166791, 0.9191139, 0.9215140, 0.9238795, 0.9262102, 0.9285061, 0.9307669, 0.9329928, 0.9351835, 0.9373390, 0.9394592, 0.9415441, 0.9435934, 0.9456073, 0.9475856, 0.9495282, 0.9514350, 0.9533060, 0.9551412, 0.9569404, 0

33、.9587035, 0.9604305, 0.9621214, 0.9637761, 0.9653944, 0.9669765, 0.9685221, 0.9700313, 0.9715039, 0.9729400, 0.9743394, 0.9757021, 0.9770281, 0.9783174, 0.9795698, 0.9807853, 0.9819639, 0.9831055, 0.9842101, 0.9852777, 0.9863081, 0.9873014, 0.9882576, 0.9891765, 0.9900582, 0.9909027, 0.9917098, 0.99

34、24796, 0.9932119, 0.9939070, 0.9945646, 0.9951847, 0.9957674, 0.9963126, 0.9968203, 0.9972904, 0.9977230, 0.9981181, 0.9984756, 0.9987954, 0.9990777, 0.9993224, 0.9995294, 0.9996988, 0.9998306, 0.9999247, 0.9999812, 1.0000000, 0.9999812, 0.9999247, 0.9998306, ;/* cosine table */static float code cos

35、_table_float260= 1.0000000, 0.9999812, 0.9999247, 0.9998306, 0.9996988, 0.9995294, 0.9993224, 0.9990777, 0.9987954, 0.9984756, 0.9981181, 0.9977230, 0.9972904, 0.9968203, 0.9963126, 0.9957674, 0.9951847, 0.9945646, 0.9939070, 0.9932119, 0.9924796, 0.9917098, 0.9909027, 0.9900582, 0.9891765, 0.988257

36、6, 0.9873014, 0.9863081, 0.9852777, 0.9842101, 0.9831055, 0.9819639, 0.9807853, 0.9795698, 0.9783174, 0.9770281, 0.9757021, 0.9743394, 0.9729400, 0.9715039, 0.9700313, 0.9685221, 0.9669765, 0.9653944, 0.9637761, 0.9621214, 0.9604305, 0.9587035, 0.9569404, 0.9551412, 0.9533060, 0.9514350, 0.9495282,

37、0.9475856, 0.9456073, 0.9435934, 0.9415441, 0.9394592, 0.9373390, 0.9351835, 0.9329928, 0.9307669, 0.9285061, 0.9262102, 0.9238795, 0.9215140, 0.9191139, 0.9166791, 0.9142098, 0.9117060, 0.9091680, 0.9065957, 0.9039893, 0.9013488, 0.8986745, 0.8959662, 0.8932243, 0.8904487, 0.8876396, 0.8847971, 0.8

38、819213, 0.8790122, 0.8760701, 0.8730950, 0.8700870, 0.8670462, 0.8639728, 0.8608670, 0.8577286, 0.8545580, 0.8513552, 0.8481203, 0.8448536, 0.8415550, 0.8382247, 0.8348629, 0.8314696, 0.8280451, 0.8245893, 0.8211025, 0.8175848, 0.8140363, 0.8104572, 0.8068476, 0.8032075, 0.7995372, 0.7958369, 0.7921

39、066, 0.7883464, 0.7845566, 0.7807372, 0.7768885, 0.7730104, 0.7691033, 0.7651672, 0.7612024, 0.7572088, 0.7531868, 0.7491364, 0.7450578, 0.7409511, 0.7368166, 0.7326543, 0.7284644, 0.7242471, 0.7200025, 0.7157308, 0.7114322, 0.7071068, 0.7027547, 0.6983762, 0.6939715, 0.6895406, 0.6850837, 0.6806010

40、, 0.6760927, 0.6715590, 0.6669999, 0.6624158, 0.6578067, 0.6531729, 0.6485144, 0.6438316, 0.6391245, 0.6343933, 0.6296383, 0.6248595, 0.6200572, 0.6152316, 0.6103828, 0.6055111, 0.6006165, 0.5956993, 0.5907597, 0.5857978, 0.5808139, 0.5758082, 0.5707808, 0.5657318, 0.5606616, 0.5555702, 0.5504580, 0

41、.5453250, 0.5401714, 0.5349976, 0.5298036, 0.5245897, 0.5193560, 0.5141028, 0.5088301, 0.5035384, 0.4982277, 0.4928982, 0.4875502, 0.4821838, 0.4767992, 0.4713967, 0.4659765, 0.4605387, 0.4550836, 0.4496113, 0.4441221, 0.4386162, 0.4330938, 0.4275551, 0.4220003, 0.4164295, 0.4108432, 0.4052413, 0.3996242, 0.3939920, 0.3883450, 0.3826834, 0.3770074, 0.3713172, 0.3656130, 0.3598951, 0.3541635, 0.3484187, 0.3426607, 0.3368899, 0.3311063, 0.3253103, 0.3195020, 0.3136818, 0.3078496, 0.3020059, 0.2961509, 0.2902847, 0.2844075, 0.2785197, 0.2726214