幅頻均衡(F題)

1、數(shù)字幅頻均衡放大器摘要：本系統(tǒng)采用LM3S811作為主控制器，通過前置放大、濾波、阻帶網(wǎng)絡(luò)，再經(jīng)過幅頻均衡，對信號進(jìn)行處理，使得數(shù)字濾波處理后波形與阻帶網(wǎng)絡(luò)前波形相同，再對信號進(jìn)行功率放大，盡量使得所得功率最大。關(guān)鍵詞：LM3S811、前置放大、帶通濾波、數(shù)字均衡、低頻功放一、系統(tǒng)方案分析隨著數(shù)字信號處理技術(shù)的發(fā)展，單片機(jī)技術(shù)已廣泛應(yīng)用于各個(gè)領(lǐng)域。借助于現(xiàn)代數(shù)字電子及數(shù)字信號處理技術(shù)，需對不規(guī)則的信號進(jìn)行還原。圍繞這一課題我們進(jìn)行方案選擇與論證、系統(tǒng)的軟硬件設(shè)計(jì)與調(diào)試，基本實(shí)現(xiàn)了課目的某些項(xiàng)指標(biāo)和要求。并在此基礎(chǔ)上，撰寫了本報(bào)告的。整個(gè)系統(tǒng)分為前置放大、帶通濾波、阻帶網(wǎng)絡(luò)、數(shù)字均衡及功率放大幾

2、個(gè)部分，系統(tǒng)總方案和實(shí)現(xiàn)方框圖如下： 由于信號比較微弱，放大倍數(shù)比較大，因此可采用低噪聲窄運(yùn)放對信號進(jìn)行放大，只要在20HZ20KHZ頻帶內(nèi)保持平坦，紋波小于1dB就可以滿足要求，此外，要求輸出阻抗為600，可以在運(yùn)放的輸出端串聯(lián)600的電阻直接接入下一級;在電路設(shè)計(jì)的每一個(gè)模塊都對前后的電路都有影響，故用電壓跟隨器分隔開?？紤]到信號經(jīng)過帶阻網(wǎng)絡(luò)后信號沒有直流偏置，而LM3S811數(shù)字濾波采樣只能對正電壓進(jìn)行濾波處理，故在數(shù)字濾波前時(shí)采用鉗位電路對信號進(jìn)行抬高。功率放大模塊采用甲乙類功放，因?yàn)樗睦碚撔士蛇_(dá)78.5%，并且發(fā)熱相對較小，失真率比小,符合題目的要求。二、前置放大器的設(shè)計(jì)2.1

3、前置放大的硬件設(shè)計(jì) 可控增益寬帶放大器由芯片OP37構(gòu)成?？煽卦鲆娣糯笃饔蒾p37運(yùn)放構(gòu)成，op37是低漂移，高速，增益帶寬高，允許精確的高增益。本題中采用一片此芯片構(gòu)成放大級，電路圖如圖1所示。圖1前級放大電路 由圖可知，最大放大倍數(shù)可達(dá)650倍，可根據(jù)電位器來調(diào)節(jié)，符合題目要求。實(shí)驗(yàn)中輸入10mv(Vpp）電壓，輸出5Vpp，放大倍數(shù)為500倍。2.2帶通濾波器的設(shè)計(jì)方案一：采用無源濾波 由于無源濾波通過電感和電容的匹配對某次諧波并聯(lián)低阻（調(diào)諧濾波）狀態(tài),給某次諧波電流構(gòu)成一個(gè)低阻態(tài)通路。這樣諧波電流就不會(huì)流入系統(tǒng)，缺點(diǎn)為諧波濾除率一般只有80，對基波的無功補(bǔ)償也是一定的 。方案二：有源濾

4、波 有源濾波自身就是諧波源。在檢測到系統(tǒng)諧波的同時(shí)產(chǎn)生一組和系統(tǒng)幅值相等，相位相反的諧波向量，這樣可以抵消掉系統(tǒng)諧波，使其成為正弦波形。有源濾波除了濾除諧波外，同時(shí)還可以動(dòng)態(tài)補(bǔ)償無功功率。其優(yōu)點(diǎn)是反映動(dòng)作迅速，濾除諧波可達(dá)到95以上，補(bǔ)償無功細(xì)致。故采用此方案。2.3帶阻網(wǎng)絡(luò)設(shè)計(jì) 根據(jù)本題中要求，所制作的帶阻網(wǎng)絡(luò)對前置放大電路所輸出的信號進(jìn)行帶阻濾波，根據(jù)題目要求，制作帶阻網(wǎng)絡(luò)對前置放大電路輸出信號進(jìn)行濾波，以10kHz時(shí)輸出信號電壓幅度為基準(zhǔn)，要求最大衰減10dB。2.14鉗位電路設(shè)計(jì)由于ARM系統(tǒng)對信號的處理只能是正的，不能處理負(fù)的數(shù)值，故需要對輸入波形進(jìn)行直流偏置，使得電壓抬高都為正值，

5、故采用鉗位電路設(shè)計(jì)方案。鉗位電路的作用是將周期性變化的波形的頂部或底部保持在某一確定的直流電平上。圖上部為常見的二極管鉗位電路。二極管的鉗位作用是指利用二極管正向?qū)▔航迪鄬Ψ€(wěn)定，且數(shù)值較小(有時(shí)可近似為零)的特點(diǎn)，來限制電路中某點(diǎn)的電位。鉗位電路設(shè)計(jì)如圖5所示。圖2鉗位電路在上面電路中，Vcc為輸入直流電壓，已達(dá)到電路的完全偏移。電路中電阻設(shè)的要很大，以便延遲充放電時(shí)間。實(shí)際測量中，Vcc為0.5V，R=200K，前級噪聲24.4mv(Vpp)。三、數(shù)字均衡方法比較與選擇3.1系統(tǒng)方案論述3.1.1數(shù)字信號處理器選擇方案一：高速FPGA可用VerilogHDL語言設(shè)計(jì)專用的數(shù)字信號處理電路，

6、由硬件實(shí)現(xiàn)信號處理，速度得到了保證，也就可以得到很好的帶寬平坦度，但FPGA價(jià)格昂貴，開發(fā)周期長，且本課題僅要求對20Hz-20kHz的音頻信號進(jìn)行處理，采用高速FPGA過于浪費(fèi)資源。方案二：DSP處理器具有合適于數(shù)據(jù)進(jìn)行快速運(yùn)算的硬件結(jié)構(gòu)，工作頻率高，算法成熟，是實(shí)現(xiàn)數(shù)字信號運(yùn)算與處理的理想方案，完全可勝任系統(tǒng)要求，但DSP處理器同樣存在成本高的問題。方案三：ARM微控制器通過鎖相環(huán)倍頻可使主頻達(dá)到相當(dāng)高的工作頻率，其成本相對FPGA、DSP更低，只要數(shù)字信號處理算法設(shè)計(jì)合理，用于處理音頻段范圍內(nèi)的信號問題不大。綜上所述，采用ARM系列的LM3S811微控制器作為數(shù)字信號處理器較為合適。3.

7、1.2數(shù)字幅頻均衡方法選擇方案一：利用FFT快速運(yùn)算的辦法對輸入信號進(jìn)行離散傅立葉變換，分析其頻譜，然后根據(jù)所希望的頻率特性進(jìn)行濾波，再利用傅立葉反變換將頻域信號恢復(fù)至?xí)r域信號，但FFT是塊計(jì)算，對速度較高的FPGA和DSP處理器可行，若采用單片機(jī)或ARM微控制器作為數(shù)字信號處理器，則難以滿足設(shè)計(jì)要求。方案二：采用無限沖激響應(yīng)濾波器IIR，IIR濾波器設(shè)計(jì)簡單，實(shí)現(xiàn)的階數(shù)較低，但它具有非線性相位，且由于其為反饋型結(jié)構(gòu)，對濾波器的參數(shù)精度要求較高，否則可能引起震蕩或發(fā)散。方案三：采用有限沖激響應(yīng)FIR濾波器，F(xiàn)IR濾波器采用非遞歸結(jié)構(gòu)，可以得到嚴(yán)格的線性相位，運(yùn)算誤差較小，且傳遞函數(shù)不存在極點(diǎn)，

8、穩(wěn)定性好，但與IIR相比，相同條件下需要的階數(shù)更高，導(dǎo)致延時(shí)時(shí)間較長。由于本次設(shè)計(jì)采用ARM控制器，運(yùn)算速度相對DSP和FPGA較低，所以濾波器階數(shù)不能太高，故采用了IIR數(shù)字濾波器。四、低頻功放的硬件設(shè)計(jì)按功放中功放管的導(dǎo)電方式不同功率放大器分有甲類（又稱A類）、乙類（又稱B類）、甲乙類（又稱AB類）、和丁類（又稱D類）。方案一：甲類功放。甲類功放是指在信號的整個(gè)周期內(nèi)（正弦波的正負(fù)兩個(gè)半周），放大器的任何功率輸出元件都不會(huì)出現(xiàn)電流截至（即停止輸出）的一類放大器。甲類放大器工作時(shí)會(huì)產(chǎn)生高熱，效率很低，但固有的有點(diǎn)是不存在就交越失真。單端放大器都是甲類工作方式，推免放大器可以是甲類，也可以是乙

9、類或甲乙類。因此甲類在輸入信號的一個(gè)周期內(nèi)都導(dǎo)通，晶體管導(dǎo)通角度為360°，特點(diǎn)是線性好、失真小且失真成分以偶次諧波為主。效率小于50%，所以不符合題目中60%的要求。方案二：乙類功放。乙類功放是指正弦信號的正負(fù)兩個(gè)半周分別由推免輸出級的兩“臂”輪流放大輸出的乙類放大器，每一“臂”的到點(diǎn)時(shí)間為信號的半個(gè)周期，乙類放大器的優(yōu)點(diǎn)是效率高，缺點(diǎn)是會(huì)產(chǎn)生交越失真。乙類功放晶體管只在輸入信號正半周工作在放大區(qū)（導(dǎo)通）晶體管導(dǎo)通角度為180°，乙類功放的效率平均約為75%，產(chǎn)生的熱量較甲類機(jī)低，容許使用較小的散熱器。方案三：甲乙類功放。甲乙類功放界于甲類和乙類之間，推免放大的每一個(gè)“臂

10、”導(dǎo)通時(shí)間大于信號的半個(gè)周期而小于一個(gè)周期，甲乙類放大有效解決了乙類放大器的交越失真問題，效率又比甲類放大器高，因此獲得了極為廣泛的應(yīng)用。甲乙類功放在信號的大半周期間導(dǎo)通，導(dǎo)通角度180°<<360°。方案四：丁類功放。D類功率放大器工作于開關(guān)狀態(tài)，理論效率可達(dá)100%，實(shí)際的運(yùn)用也可達(dá)80%以上。功率器件的耗散功率小，產(chǎn)生熱量少，可以大大減小散熱器的尺寸，連續(xù)輸出功率很容易達(dá)到數(shù)百瓦。功率MOS有自保護(hù)電路，可以大大簡化保護(hù)電路，而且不會(huì)引入非線性失真。但是相對設(shè)計(jì)電路比較復(fù)雜。 根據(jù)題目的要求，功率放大器接8電阻負(fù)載（一端接地）時(shí)，要求輸出功率10W，輸出電

11、壓波形無明顯失真，功率放大電路的-3dB通頻帶為20Hz20kHz，功率放大電路的效率60。綜上四個(gè)方案的比較，我們決定采用甲乙類功放，效率高于60，且所需的電路并不難。功放的電路流程圖如下：圖3 功放硬件電路設(shè)計(jì)流程圖 音頻放大器輸出失真度小于某一數(shù)值(r<1%)的最大功率稱為額定功率，表達(dá)式 Po= Uo2/RL Pi=Ui2*I Uo為負(fù)載兩端的最大不失真電壓，RL為額定負(fù)載阻抗。五、軟件設(shè)計(jì)5.1軟件設(shè)計(jì)本次設(shè)計(jì)，利用LM3S811的片內(nèi)的十位AD，很方便地實(shí)現(xiàn)AD轉(zhuǎn)換，經(jīng)過IIR濾波器后，DA輸出。由于沒有片內(nèi)DA，所以DA需要外擴(kuò)，采用了12位高速DA轉(zhuǎn)換芯片MAX7541，

12、只要不斷地往其12個(gè)數(shù)據(jù)口發(fā)送數(shù)據(jù)，就能實(shí)現(xiàn)DA轉(zhuǎn)換。5.1.1基于matlab的IIR濾波器設(shè)計(jì) 設(shè)IIR數(shù)字濾波器的輸入為x(k)，輸出為y(k),則該濾波器的可用常系數(shù)差分方程表示為其中ai、bi為常系數(shù)，N、M為整數(shù)，i=0,1,2,，k=0,1,2,，其傳遞函數(shù)為IIR濾波器的設(shè)計(jì)過程主要是尋找濾波器的各個(gè)系數(shù)，使其逼近所要求的特性指標(biāo)。其設(shè)計(jì)方法有兩種：一是直接法，也稱最優(yōu)化設(shè)計(jì)法，先確定一種最優(yōu)化準(zhǔn)則（均方誤差最小準(zhǔn)則或最大誤差最小化準(zhǔn)則），再求此最佳準(zhǔn)則下的濾波器系數(shù)；二是間接法，先根據(jù)已有模擬濾波器的設(shè)計(jì)理論設(shè)計(jì)一個(gè)合適的濾波器，然后變換成滿足技術(shù)指標(biāo)的數(shù)字濾波器。間接設(shè)計(jì)法

13、也稱經(jīng)典設(shè)計(jì)法，如果設(shè)計(jì)的IIR數(shù)字濾波器頻率特性有特殊要求，經(jīng)典設(shè)計(jì)往往顯得無能為力，而最優(yōu)設(shè)計(jì)法在這些頻率特性有要求的場合很有用場。MATLAB提供了yulewalk和maxflat兩個(gè)最優(yōu)設(shè)計(jì)函數(shù)Error! Reference source not found.。5.1.2IIR具體設(shè)計(jì)本系統(tǒng)處理信號的頻率特性比較特殊，故利用yulewalk最優(yōu)設(shè)計(jì)函數(shù)設(shè)計(jì)IIR濾波器，首先測出帶阻網(wǎng)絡(luò)的幅頻特性，fs為采樣頻率，f為0到fs/2的頻率點(diǎn)，M為各頻率點(diǎn)對應(yīng)的放大倍數(shù)，然后對M取倒數(shù)，得到m。本次設(shè)計(jì)濾波器的階數(shù)為5，則調(diào)用函數(shù)b,a=yulewalk(4,f,m)得到濾波器的參數(shù)如下b

14、= 1.2083,-0.9642,-0.3355,0.3357 ,0.0567a= 1.0000,-1.2246,0.0116,0.5014 ,-0.1455由系數(shù)b和a可得常系數(shù)差分方程其中b= 1.2083,-0.9642,-0.3355,0.3357 ,0.0567, a= 1.0000,-1.2246,0.0116,0.5014 ,-0.1455, k=0,1,2,3,4。根據(jù)上述差分方程進(jìn)行軟件編程并固化程序，即可使微控制器LM3S811實(shí)現(xiàn)IIR濾波的功能。圖中紅線為理想濾波器的幅頻特性，而藍(lán)線為得到的濾波器的幅頻特性。六、附錄（一）硬件電路附錄6.1帶通濾波電路圖與測試結(jié)果：圖2

15、 -1dB（20HZ-20KHZ)帶通濾波電路 在電路中，虛線左端是20Hz高通濾波器，虛線右端是20KHz低通濾波器。此電路是帶通可調(diào)的濾波器，可以調(diào)節(jié)電位器的值來調(diào)節(jié)帶通范圍，具有更實(shí)用性。頻率(KHz)0.020.030.040.080.20.61.05.010.014.018.019.020.0輸出Vpp(v)4.484.804.965.125.125.125.125.045.004.964.804.644.45-1dB時(shí)，即衰減到0.89時(shí)，如輸出通頻帶為5V，則截止時(shí)為4.45。根據(jù)表格可知，符合帶通濾波設(shè)計(jì)要求。6.2 前級輸出電阻測量空載輸出電壓加入600歐電阻負(fù)載電壓2.74

16、V1.37V可知符合600歐輸出電阻的要求。6.3本次制作的帶阻網(wǎng)絡(luò)電路圖帶阻網(wǎng)絡(luò)測試測試輸入信號 VPP=5V頻率/HZ20401003005006008009001K輸出1.211.050.660.450.410.410.390.380.4頻率/HZ4K6K8K10K12K14K16k18k20k輸出0.520.620.720.820.920.961.021.041.06圖3 帶阻網(wǎng)絡(luò)電路圖根據(jù)題目所給的帶阻網(wǎng)絡(luò)結(jié)構(gòu)，我們采用Multism進(jìn)行了輔助分析與設(shè)計(jì)，其幅頻性的分析結(jié)果如下圖4所示。 圖4 仿真圖由圖可知，在10KHZ時(shí)，輸出信號大約在-10dB,信號的最大衰減為-25dB，故滿

17、足題目要求。但特別注意，仿真都是理想狀態(tài)，要想實(shí)際達(dá)到要求，必須提高電路中電容、電感、電阻的精度。6.4功率放大電路圖及測試圖7 功放前置20Hz到20KHz的帶通濾波圖8 前級放大和功率放大電路功放測試：測試條件：帶8歐姆負(fù)載，由信號發(fā)生器輸入5V VPP信號頻率205010020050010002000500010k15k20k輸入電壓/V10.610.610.610.610.610.610.610.610.610.610.6波形失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真無明顯失真功率輸出/W10.210.312.112.110

18、.211.011.111.711.111.810.8效率6.5外擴(kuò)D/A電路圖4 MAX7541引腳圖則DA外擴(kuò)電路如下：圖5 MAX7541外接電路圖（二）程序設(shè)計(jì)附錄程序流程圖如下主函數(shù)流程圖 中斷服務(wù)程序流程圖源程序代碼/*PD0PD7接B0B7，PA0PA3接B8B11*/#include <hw_types.h>#include <hw_memmap.h>#include <hw_ints.h>#include <interrupt.h>#include <sysctl.h>#include <gpio.h>#i

19、nclude <timer.h>#include <adc.h>/ 將較長的標(biāo)識符定義成較短的形式#define SysCtlPeriEnable SysCtlPeripheralEnable#define SysCtlPeriDisable SysCtlPeripheralDisable#define GPIOPinTypeIn GPIOPinTypeGPIOInput#define GPIOPinTypeOut GPIOPinTypeGPIOOutput#define GPIOPinTypeOD GPIOPinTypeGPIOOutputOD#define ADCS

20、equEnable ADCSequenceEnable#define ADCSequDisable ADCSequenceDisable#define ADCSequConfig ADCSequenceConfigure#define ADCSequStepConfig ADCSequenceStepConfigure#define ADCSequDataGet ADCSequenceDataGet#define KEY_PERIPH SYSCTL_PERIPH_GPIOB#define KEY_PORT GPIO_PORTB_BASE#define KEY_PIN GPIO_PIN_5#de

21、fine DATA_PERIPH_8 SYSCTL_PERIPH_GPIOD#define DATA_PORT_8 GPIO_PORTD_BASE#define DA_PINS_8GPIO_PIN_0|GPIO_PIN_1|GPIO_PIN_2|GPIO_PIN_3|GPIO_PIN_4|GPIO_PIN_5|GPIO_PIN_6|GPIO_PIN_7#define DATA_PERIPH_4 SYSCTL_PERIPH_GPIOA#define DATA_PORT_4 GPIO_PORTA_BASE#define DATA_PINS_4 GPIO_PIN_0|GPIO_PIN_1|GPIO_

22、PIN_2|GPIO_PIN_3#define N 7/ IIR濾波器階數(shù)int bN+1 = 6656,-9185,1559,5564,-4080,168,582,-92;/ IIR濾波器系數(shù)int aN+1 = 10000,-14301,886,9744,-5991,-639,1146,-25;unsigned long xN+1=0;/采樣序列unsigned long yN+1=0;/輸出序列unsigned long TheSysClock;/系統(tǒng)時(shí)鐘unsigned long ADC_DATA;int ADC_EndFlag = 0;/AD采樣階數(shù)標(biāo)志/*防止jtag鎖死函數(shù)*/v

23、oid jtagWait(void) SysCtlPeriEnable(KEY_PERIPH); / 使能KEY所在的GPIO端口GPIOPinTypeIn(KEY_PORT, KEY_PIN); / 設(shè)置KEY所在管腳為輸入if (GPIOPinRead(KEY_PORT, KEY_PIN) = 0x00) / 若復(fù)位時(shí)按下KEY，則進(jìn)入for(;); / 死循環(huán)，以等待JTAG連接SysCtlPeriDisable(KEY_PERIPH); / 禁止KEY所在的GPIO端口/*系統(tǒng)時(shí)鐘初始化函數(shù)*/void clockInit(void) /*SysCtlLDOSet(SYSCTL_LDO

24、_2_50V); / 設(shè)置LDO輸出電壓SysCtlClockSet(SYSCTL_USE_OSC | / 系統(tǒng)時(shí)鐘設(shè)置SYSCTL_OSC_MAIN| / 采用主振蕩器SYSCTL_XTAL_6MHZ| / 外接6MHz晶振SYSCTL_SYSDIV_1); / 不分頻*/SysCtlLDOSet(SYSCTL_LDO_2_75V); / 配置PLL前須將LDO設(shè)為2.75VSysCtlClockSet(SYSCTL_USE_PLL| / 系統(tǒng)時(shí)鐘設(shè)置，采用PLLSYSCTL_OSC_MAIN| / 主振蕩器SYSCTL_XTAL_4MHZ| / 外接4MHz晶振SYSCTL_SYSDIV_

25、4); /4分頻TheSysClock = SysCtlClockGet( ); / 獲取當(dāng)前的系統(tǒng)時(shí)鐘頻率/*Timer初始化*/void timerInit(void)SysCtlPeriEnable(SYSCTL_PERIPH_TIMER2); / 使能Timer模塊TimerConfigure(TIMER2_BASE, TIMER_CFG_32_BIT_PER); / 配置Timer為32位周期定時(shí)TimerControlTrigger(TIMER2_BASE, TIMER_A, true); / 使能內(nèi)部觸發(fā)脈沖的產(chǎn)生TimerControlStall(TIMER2_BASE, TI

26、MER_A, true); / 調(diào)試時(shí)暫停計(jì)數(shù)（必要！）TimerLoadSet(TIMER2_BASE, TIMER_A, 600); / 設(shè)置Timer初值,SysCtlClockGet()是定時(shí)1s,SysCtlClockGet()/2是定時(shí)0.5sTimerEnable(TIMER2_BASE, TIMER_A); / 使能Timer計(jì)數(shù)/*ADC初始化*/void adcInit(void)SysCtlPeriEnable(SYSCTL_PERIPH_ADC); / 使能ADC模塊SysCtlADCSpeedSet(SYSCTL_ADCSPEED_250KSPS); / 設(shè)置ADC采

27、樣率ADCSequDisable(ADC_BASE, 1); / 禁止采樣序列ADCSequConfig(ADC_BASE,1,ADC_TRIGGER_TIMER,0);/ 采樣序列配置：ADC基址，采樣序列0，定時(shí)器觸發(fā)，采樣優(yōu)先級0ADCSequStepConfig(ADC_BASE,1,0,ADC_CTL_CH1|ADC_CTL_END|ADC_CTL_IE);/ 采樣步進(jìn)設(shè)置：ADC基址，采樣序列0，步值0，采樣ADC0后停止并申請中斷ADCIntEnable(ADC_BASE,1); / 使能ADC中斷IntEnable(INT_ADC1); / 使能ADC采樣序列中斷IntMast

28、erEnable(); / 使能處理器中斷ADCSequEnable(ADC_BASE, 1); / 使能采樣序列/* ADC采樣*/unsigned long adcSample(void)unsigned long ulValue;while (!ADC_EndFlag); / 等待采樣結(jié)束ADC_EndFlag = 0; / 清除ADC采樣結(jié)束標(biāo)志ADCSequDataGet(ADC_BASE, 1, &ulValue); / 讀取ADC轉(zhuǎn)換結(jié)果return(ulValue);/* *DAC輸出函數(shù)*/void DAC(unsigned long data) GPIOPinWri

29、te(DATA_PORT_4,DATA_PINS_4,(data&0xF00)>>8);/送高四位數(shù)據(jù) GPIOPinWrite(DATA_PORT_8,DATA_PINS_8,(data&0x0FF);/送低八位數(shù)據(jù)/*將上次存入數(shù)組元素整體左移一次,丟棄第一個(gè)元素*/void MoveNext(void) int i=0; for(i=0; i<N; i+) xi=xi+1; yi=yi+1; /*IIR濾波*/void IIR() int i; unsigned long sum=0; for(i=0;i<N;i+) sum+=bi*xN-i-yi*aN-i; yN=sum+bN*x0;main() jtagWait();/系統(tǒng)初始化 clockIni