實(shí)戰(zhàn)：存儲(chǔ)架構(gòu)、DMA與實(shí)時(shí)性優(yōu)化入門)
1. 先別急著寫代碼搞清楚DSP到底指什么剛接觸這塊的時(shí)候我被問得最多的一句話是「學(xué)DSP是不是得先把信號(hào)與系統(tǒng)啃一遍」。問這話的人里有一半連示波器都沒怎么摸過另一半則是從單片機(jī)轉(zhuǎn)過來(lái)、看到「哈佛架構(gòu)」「流水線」「EDMA」這幾個(gè)詞就頭皮發(fā)麻。我特別理解這種感受因?yàn)槲易约寒?dāng)年也是從一份全英文的數(shù)據(jù)手冊(cè)開始翻了三頁(yè)就合上了。先把最容易混淆的概念拆開。DSP這三個(gè)字母在不同語(yǔ)境下指三樣完全不同的東西第一種是芯片一類專門為乘加運(yùn)算和實(shí)時(shí)數(shù)據(jù)流設(shè)計(jì)的處理器第二種是學(xué)科數(shù)字信號(hào)處理這套數(shù)學(xué)理論第三種是開發(fā)范式一種以確定性時(shí)序、數(shù)據(jù)流驅(qū)動(dòng)、資源極度受限為前提的工程方法。小白最容易踩的坑就是把三者當(dāng)成一件事結(jié)果去啃了三個(gè)月的濾波器推導(dǎo)回頭連開發(fā)環(huán)境都裝不利索。這篇文章講的是第一種和第三種——怎么在一塊DSP芯片上把程序跑起來(lái)、跑對(duì)、跑快第二種只講夠用的部分。適合誰(shuí)看如果你有單片機(jī)基礎(chǔ)能看懂C語(yǔ)言和基本的電路圖想在音頻、電機(jī)控制、雷達(dá)前端、工業(yè)檢測(cè)這類場(chǎng)景里做實(shí)時(shí)信號(hào)處理那這篇內(nèi)容就是給你準(zhǔn)備的。零基礎(chǔ)也能跟我會(huì)把每個(gè)術(shù)語(yǔ)都落到具體操作上。2. 選型和硬件準(zhǔn)備別在第一步就把錢浪費(fèi)掉2.1 為什么不用普通單片機(jī)硬扛很多人第一反應(yīng)是「我用一顆主頻兩三百兆的單片機(jī)不行嗎」。短答案看你的算法預(yù)算。長(zhǎng)答案得從架構(gòu)說(shuō)起。通用單片機(jī)的核心思路是「什么活都能干」指令集通用、分支預(yù)測(cè)復(fù)雜、緩存層級(jí)深。而DSP的核心思路是「把乘加做到極致」單周期完成一次乘加MAC、硬件零開銷循環(huán)、位反轉(zhuǎn)尋址專門服務(wù)FFT、飽和運(yùn)算防止溢出、多總線哈佛結(jié)構(gòu)讓取指和取數(shù)互不打架。舉個(gè)實(shí)測(cè)量的對(duì)比同樣做1024點(diǎn)復(fù)數(shù)FFT一顆200MHz左右的通用MCU用CMSIS-DSP庫(kù)跑下來(lái)大概在20萬(wàn)到30萬(wàn)個(gè)時(shí)鐘周期量級(jí)而一顆600MHz的浮點(diǎn)DSP能壓到1萬(wàn)個(gè)周期出頭。差了將近兩個(gè)數(shù)量級(jí)這不只是主頻的差距是架構(gòu)層面的差距。更關(guān)鍵的是確定性。實(shí)時(shí)系統(tǒng)里最怕的不是慢是「有時(shí)候快有時(shí)候慢」。通用MCU的緩存命中與否會(huì)讓同一段代碼的執(zhí)行時(shí)間浮動(dòng)好幾倍而DSP可以把L1配置成SRAM而不是cache讓每次訪問的周期數(shù)固定下來(lái)。這個(gè)特性在電機(jī)環(huán)路、音頻采樣這種硬實(shí)時(shí)場(chǎng)景里是剛需。注意不要因?yàn)椤窪SP聽起來(lái)高級(jí)」就無(wú)腦上。如果你的算法是簡(jiǎn)單的PID加幾個(gè)濾波幾百兆的MCU完全夠用開發(fā)成本還低得多。DSP的價(jià)值在算法復(fù)雜度真正壓到通用核喘不過氣的時(shí)候才體現(xiàn)。2.2 開發(fā)板怎么挑看三個(gè)硬指標(biāo)第一塊板子選錯(cuò)后面全是坑。我總結(jié)下來(lái)就看三樣考察項(xiàng)具體看什么為什么重要外設(shè)閉環(huán)是否同時(shí)有ADC/DAC或音頻編解碼器、至少一個(gè)外部存儲(chǔ)接口、串口、若干GPIO沒有輸入輸出你只能跑純計(jì)算驗(yàn)證不了實(shí)時(shí)性資料完整度是否提供原理圖不是框圖是完整原理圖、BOM、成套例程出問題時(shí)能對(duì)著電路查信號(hào)這是救命的東西調(diào)試鏈路仿真器是板載還是外接、接口定義是否標(biāo)準(zhǔn)便宜的板子經(jīng)常省掉這部分自己配很折騰我并不推薦一上來(lái)就自己畫板。DSP的電源時(shí)序、內(nèi)核電壓、PLL環(huán)路濾波、去耦電容布局任何一個(gè)細(xì)節(jié)出錯(cuò)都可能導(dǎo)致「程序能下載但一跑就飛」。先在成熟評(píng)估板上把軟件跑通等你能獨(dú)立定位「是硬件問題還是軟件問題」了再考慮自制最小系統(tǒng)。2.3 自制最小系統(tǒng)時(shí)的三個(gè)致命細(xì)節(jié)如果你確實(shí)要自己畫這三處必須盯死。上電時(shí)序。多數(shù)DSP要求內(nèi)核電壓先于或同時(shí)于IO電壓建立反過來(lái)的話內(nèi)部ESD二極管可能倒灌長(zhǎng)期會(huì)損傷芯片。穩(wěn)妥做法是用帶時(shí)序控制的電源管理芯片或者至少讓兩路電源的使能信號(hào)有明確的先后。我見過一塊板子因?yàn)閮陕冯娫赐瑫r(shí)上電、上升沿斜率不同導(dǎo)致大約每十次上電有一次啟動(dòng)失敗查了兩周才定位。復(fù)位保持時(shí)間。復(fù)位信號(hào)必須在電源穩(wěn)定之后還要保持一段時(shí)間具體數(shù)值查手冊(cè)的電氣特性表。太短會(huì)導(dǎo)致內(nèi)部PLL沒鎖住就釋放復(fù)位表現(xiàn)是「偶爾能啟動(dòng)大部分時(shí)間卡死」。時(shí)鐘與PLL。外部晶振的負(fù)載電容要按晶振廠家的規(guī)格配不是隨便放兩個(gè)22pF。PLL的環(huán)路濾波電阻電容更要嚴(yán)格按照手冊(cè)給的參數(shù)這顆電容選錯(cuò)輸出時(shí)鐘抖動(dòng)會(huì)大到讓高速外設(shè)工作異常。心得新板子第一次上電先不裝芯片或者先只焊電源部分用示波器把每一路電壓的上升時(shí)間和紋波看一遍。紋波超過手冊(cè)要求的一半就要改別等調(diào)試階段再回頭懷疑電源。3. 開發(fā)環(huán)境落地從裝軟件到點(diǎn)亮第一個(gè)中斷3.1 工具鏈的四個(gè)組成部分一套完整的DSP開發(fā)環(huán)境包含四塊集成開發(fā)環(huán)境編輯器項(xiàng)目管理調(diào)試前端、C編譯器含優(yōu)化器和匯編器、鏈接器配合鏈接命令文件工作、以及芯片廠商提供的外設(shè)驅(qū)動(dòng)庫(kù)和算法庫(kù)。這四塊缺一不可尤其后面兩個(gè)是新手最容易忽略的。我的建議是第一次安裝全部用默認(rèn)路徑、默認(rèn)版本。工具鏈版本不匹配導(dǎo)致的詭異問題比如某個(gè)庫(kù)編譯通過但鏈接報(bào)符號(hào)沖突排查成本極高新手根本扛不住。工程目錄我習(xí)慣這樣組織project/ ├── src/ 源文件 ├── include/ 頭文件 ├── cmd/ 鏈接命令文件 ├── lib/ 第三方庫(kù)和廠商庫(kù) ├── obj/ 編譯中間產(chǎn)物不進(jìn)版本庫(kù) └── tools/ 燒寫腳本、輔助工具3.2 鏈接命令文件DSP開發(fā)的第一道門檻這是從單片機(jī)轉(zhuǎn)過來(lái)的人最先懵的地方。通用MCU通常由IDE自動(dòng)生成分散加載文件你基本不用管而DSP的鏈接命令文件.cmd必須自己寫因?yàn)樗苯記Q定了代碼和數(shù)據(jù)放在哪塊存儲(chǔ)上而這個(gè)選擇對(duì)性能的影響可能有十倍。MEMORY { L2SRAM : o 0x00800000, l 0x00040000 /* 片上 256KB */ DDR : o 0x80000000, l 0x04000000 /* 外部 64MB */ } SECTIONS { .text L2SRAM /* 代碼放片上取指快且確定 */ .cinit L2SRAM .const L2SRAM .switch L2SRAM .bss L2SRAM .stack L2SRAM .data DDR /* 大塊數(shù)據(jù)放外部 */ .far DDR }為什么這么分.text代碼段放在片上SRAM取指延遲固定中斷響應(yīng)時(shí)間才可預(yù)測(cè).data、.far這些放大數(shù)組的段丟到外部DDR容量大但訪問慢正好用DMA去搬運(yùn)。你要是把代碼放DDR里跑每次取指都可能撞上總線競(jìng)爭(zhēng)實(shí)測(cè)性能能掉到原來(lái)的三分之一。三個(gè)必須檢查的點(diǎn)一是各段的長(zhǎng)度加總別超過對(duì)應(yīng)存儲(chǔ)器的l值鏈接器會(huì)報(bào)錯(cuò)但報(bào)錯(cuò)信息經(jīng)常很含糊自己先算一遍二是棧大小要顯式指定默認(rèn)值往往只有幾KB遞歸或者大局部數(shù)組一用就溢出三是DMA緩沖區(qū)必須按cache line大小對(duì)齊常見是128字節(jié)用#pragma DATA_ALIGN(buffer, 128)聲明這個(gè)后面第5章還會(huì)細(xì)說(shuō)。3.3 第一個(gè)可驗(yàn)證的工程定時(shí)器中斷翻轉(zhuǎn)IO別一上來(lái)就寫FFT。第一個(gè)工程的目標(biāo)是驗(yàn)證工具鏈、時(shí)鐘、中斷三條鏈路都通最經(jīng)濟(jì)的方式就是定時(shí)器中斷里翻轉(zhuǎn)一個(gè)GPIO用示波器看周期對(duì)不對(duì)。初始化順序是有講究的全局關(guān)中斷清IER寄存器清CSR里的全局中斷使能位配置PLL設(shè)置倍頻系數(shù)然后輪詢PLL狀態(tài)位直到鎖定初始化定時(shí)器周期寄存器把中斷服務(wù)函數(shù)地址填入中斷向量表設(shè)置中斷向量表基址寄存器使能定時(shí)器中斷對(duì)應(yīng)的IER位全局開中斷啟動(dòng)定時(shí)器/* 定時(shí)器周期計(jì)算示例目標(biāo) 1kHz 中斷定時(shí)器輸入時(shí)鐘 75MHz */ #define TIMER_INPUT_CLK 75000000UL #define TARGET_FREQ 1000UL /* 周期值 輸入時(shí)鐘 / 目標(biāo)頻率 - 1 */ #define PERIOD_VAL (TIMER_INPUT_CLK / TARGET_FREQ - 1) void timer_init(void) { /* 先停定時(shí)器再改配置運(yùn)行中改寄存器行為不確定 */ TIMER_CTL 0x0; TIMER_PRD PERIOD_VAL; /* 75000 - 1 */ TIMER_CTL 0x3; /* 內(nèi)部時(shí)鐘源 連續(xù)計(jì)數(shù)模式 */ }PERIOD_VAL那個(gè)減一是因?yàn)橛?jì)數(shù)器從0開始計(jì)滿回繞差了這一個(gè)數(shù)頻率就有千分之一的偏差對(duì)音頻這種場(chǎng)景是能聽出來(lái)的。拿到板子先用示波器量IO口理論1kHz實(shí)際量出來(lái)是999.8Hz這個(gè)誤差來(lái)自晶振精度屬于正常范圍如果量出來(lái)是500Hz或者2kHz那就是分頻配置理解錯(cuò)了回去查定時(shí)器的預(yù)分頻字段。4. 存儲(chǔ)架構(gòu)和數(shù)據(jù)搬移DSP開發(fā)真正的分水嶺4.1 三級(jí)存儲(chǔ)該怎么配主流DSP的存儲(chǔ)層次一般是L1程序SRAM、L1數(shù)據(jù)SRAM兩者都可以配置成cache或者純SRAM、L2統(tǒng)一SRAM、外部DDR。配置策略直接決定了程序是「可預(yù)測(cè)地快」還是「平均來(lái)看還行」。我的默認(rèn)配置是L1P和L1D全部配置成SRAM不啟用cache。理由是實(shí)時(shí)系統(tǒng)里「最壞情況執(zhí)行時(shí)間」比「平均執(zhí)行時(shí)間」重要得多。cache看著命中率高但那1%的miss帶來(lái)的抖動(dòng)可能就是你電機(jī)環(huán)路失控或者音頻爆音的原因。代價(jià)是什么你要手動(dòng)把熱數(shù)據(jù)搬進(jìn)L1D。這就引出了DMA——整個(gè)DSP開發(fā)里最重要的一個(gè)外設(shè)。4.2 DMA和乒乓緩沖把搬數(shù)和算數(shù)徹底解耦DMA的價(jià)值在于「CPU在算當(dāng)前幀的時(shí)候下一幀已經(jīng)在搬了」。標(biāo)準(zhǔn)做法是雙緩沖#pragma DATA_ALIGN(buf_a, 128) #pragma DATA_ALIGN(buf_b, 128) static float buf_a[FRAME_LEN]; static float buf_b[FRAME_LEN]; volatile int ready_flag 0; static float * volatile active_buf buf_a; static float * volatile pending_buf buf_b;中斷服務(wù)函數(shù)里只做一件事切換指針、置標(biāo)志、重新掛DMA。真正的算法處理放在主循環(huán)或者低優(yōu)先級(jí)任務(wù)里。這樣中斷服務(wù)函數(shù)的執(zhí)行時(shí)間被壓到幾十個(gè)周期不會(huì)被算法拖長(zhǎng)。4.3 一個(gè)完整的參數(shù)推算過程假設(shè)你要做一個(gè)48kHz采樣率的音頻處理每幀256個(gè)采樣點(diǎn)主頻600MHz。這一串?dāng)?shù)字怎么落到代碼里幀時(shí)長(zhǎng) 256 ÷ 48000 ≈5.33ms這幀內(nèi)可用的CPU周期 5.33ms × 600MHz ≈3.2M 個(gè)周期工程上只敢用掉30%到40%剩下的留給中斷、DMA搶占和其他任務(wù)所以實(shí)際算法預(yù)算約1.0M 到 1.3M 個(gè)周期緩沖區(qū)字節(jié)數(shù) 256 × 4單精度f(wàn)loat1KB雙緩沖就是2KB對(duì)齊要求按128字節(jié)對(duì)齊用ceil(1024/128)×128算下來(lái)正好是1024沒問題但如果幀長(zhǎng)是250點(diǎn)就是1000字節(jié)必須向上補(bǔ)到1024字節(jié)聲明數(shù)組時(shí)直接給1024這個(gè)預(yù)算是你后面做算法選型的硬約束。比如一個(gè)256點(diǎn)的復(fù)數(shù)FFT如果在你的平臺(tái)上需要50萬(wàn)個(gè)周期那一幀只能跑兩次就必須考慮降到128點(diǎn)或者換更高效的實(shí)現(xiàn)。先算預(yù)算再寫代碼比寫完發(fā)現(xiàn)跑不動(dòng)再回來(lái)優(yōu)化省十倍時(shí)間。踩過的坑有一次DMA搬進(jìn)來(lái)的數(shù)據(jù)總是整體偏移幾個(gè)元素查了很久最后發(fā)現(xiàn)是源地址沒按元素大小對(duì)齊。DMA對(duì)地址對(duì)齊的要求比CPU寬松但對(duì)齊沒做好就會(huì)出現(xiàn)這種「錯(cuò)位」而不是「報(bào)錯(cuò)」的靜默故障特別難查。5. 定點(diǎn)還是浮點(diǎn)這道選擇題決定了你的整個(gè)代碼風(fēng)格5.1 Q格式說(shuō)到底就是一個(gè)縮放約定早期DSP沒有浮點(diǎn)單元全用定點(diǎn)。現(xiàn)在中高端DSP基本都帶單精度甚至雙精度浮點(diǎn)但定點(diǎn)仍然在成本和功耗敏感的場(chǎng)合大量使用而且理解定點(diǎn)對(duì)你理解溢出、飽和、精度損失非常有幫助。Q格式的本質(zhì)是「用整數(shù)表示小數(shù)」。Q15的意思是用16位有符號(hào)整數(shù)隱含小數(shù)點(diǎn)在第15位之后。格式位寬表示范圍最小分辨率實(shí)際值換算Q1516位-1.0 ~ 0.999971/32768 ≈ 0.0000305整數(shù) ÷ 32768Q3132位-1.0 ~ 0.99999999951/2^31整數(shù) ÷ 2147483648Q1216位-8.0 ~ 7.99981/4096 ≈ 0.000244整數(shù) ÷ 4096運(yùn)算規(guī)則要記牢加減法直接算但結(jié)果可能溢出乘法要先擴(kuò)位再右移。兩個(gè)Q15相乘結(jié)果是一個(gè)Q30的數(shù)必須右移15位才能回到Q15同時(shí)要做飽和處理防止溢出。static inline short q15_mul(short a, short b) { int prod (int)a * (int)b; /* Q30 */ prod (prod 0x4000) 15; /* 加半個(gè)LSB做四舍五入 */ if (prod 32767) prod 32767; if (prod -32768) prod -32768; return (short)prod; }那個(gè)0x4000是四舍五入用的偏置加不加在單次運(yùn)算里看不出差別但在濾波器這種幾百次累加的場(chǎng)景里會(huì)導(dǎo)致明顯的直流偏移。我第一次做FIR的時(shí)候沒加輸出信號(hào)整體偏了幾十個(gè)LSB還以為是ADC的問題。5.2 飽和處理寧可削頂不要翻轉(zhuǎn)定點(diǎn)運(yùn)算最怕的不是精度不夠是溢出翻轉(zhuǎn)。一個(gè)正的大數(shù)溢出會(huì)變成負(fù)數(shù)在你的音頻輸出里就是一聲刺耳的爆響在電機(jī)控制里可能就是一次危險(xiǎn)的誤動(dòng)作。所以所有可能溢出的地方都要做飽和clamp到最大最小值而不是讓它自然回繞。很多DSP有專門的飽和加法指令編譯器也提供了內(nèi)建函數(shù)能用一個(gè)周期搞定別自己寫if-else。5.3 實(shí)測(cè)對(duì)比同一段算法的兩種實(shí)現(xiàn)我在一塊600MHz的浮點(diǎn)DSP上做過一組對(duì)比32點(diǎn)復(fù)數(shù)FFT輸入數(shù)據(jù)已經(jīng)就位只算內(nèi)核循環(huán)實(shí)現(xiàn)方式平均周期數(shù)相對(duì)精度誤差代碼可讀性單精度浮點(diǎn)開 -O3約 1,050基準(zhǔn)高Q15定點(diǎn)純C約 2,300約 -62dB中需注釋Q15定點(diǎn)內(nèi)建函數(shù) 循環(huán)展開約 780約 -60dB低有意思的地方在于定點(diǎn)用了內(nèi)建函數(shù)之后反而比浮點(diǎn)快因?yàn)楦↑c(diǎn)單元在這個(gè)平臺(tái)上只有兩個(gè)定點(diǎn)單元有八個(gè)。但代價(jià)是代碼幾乎沒法維護(hù)而且每換一個(gè)平臺(tái)就要重寫。所以我的實(shí)踐原則是原型階段全用浮點(diǎn)把算法跑通、把精度驗(yàn)證過等確定性能不夠了再把最內(nèi)層的一兩個(gè)熱點(diǎn)函數(shù)改定點(diǎn)。一上來(lái)就全手寫定點(diǎn)優(yōu)化返工率極高。6. 實(shí)時(shí)性保障中斷延遲和流水線優(yōu)化6.1 中斷延遲怎么測(cè)別靠猜中斷延遲是從外部事件觸發(fā)到中斷服務(wù)函數(shù)第一條指令執(zhí)行的時(shí)間。這個(gè)數(shù)字直接決定你的系統(tǒng)能不能滿足實(shí)時(shí)要求但很多人從沒量過全憑感覺。測(cè)法很樸素在中斷服務(wù)函數(shù)的第一行立刻翻轉(zhuǎn)一個(gè)GPIO不要做任何其他事用示波器雙通道一個(gè)通道接外部觸發(fā)信號(hào)一個(gè)通道接這個(gè)GPIO看兩個(gè)上升沿的時(shí)間差。多測(cè)幾百次取最大值而不是平均值。實(shí)測(cè)經(jīng)驗(yàn)中斷延遲從幾十個(gè)周期到幾百個(gè)周期都可能主要取決于三件事中斷是否被更長(zhǎng)的高優(yōu)先級(jí)中斷屏蔽、中斷向量表是否在片上、以及中斷服務(wù)函數(shù)的第一條指令有沒有撞上取指miss。我測(cè)過一個(gè)案例把中斷向量表從外部存儲(chǔ)搬到片上SRAM之后最壞延遲從400多個(gè)周期降到80個(gè)周期左右效果比優(yōu)化算法本身明顯得多。關(guān)鍵提醒所有中斷服務(wù)函數(shù)都必須能在你算出來(lái)的「最壞延遲 執(zhí)行時(shí)間」預(yù)算內(nèi)完成。超了就要拆分成上半部和下半部上半部只做搬數(shù)據(jù)和置標(biāo)志。6.2 讓編譯器替你展開流水線超長(zhǎng)指令字架構(gòu)的DSP一個(gè)周期能發(fā)射多條指令前提是編譯器能把循環(huán)展開成軟件流水。你要做的是給它創(chuàng)造條件循環(huán)體內(nèi)不要有函數(shù)調(diào)用能內(nèi)聯(lián)就內(nèi)聯(lián)指針參數(shù)加restrict關(guān)鍵字告訴編譯器源和目的不重疊用#pragma MUST_ITERATE告訴編譯器循環(huán)次數(shù)的范圍它才敢放心展開循環(huán)次數(shù)最好是展開因子的整數(shù)倍或者手動(dòng)補(bǔ)零對(duì)齊優(yōu)化有沒有生效別猜去看編譯器生成的匯編文件。開啟「保留匯編輸出」選項(xiàng)找到你的熱點(diǎn)函數(shù)看里面有沒有那個(gè)標(biāo)志性的軟件流水序言一串并行的裝載指令和排空段。如果生成的還是一條一條串行的指令說(shuō)明優(yōu)化沒起來(lái)通常是循環(huán)里有編譯器看不透的依賴。我優(yōu)化過一個(gè)FIR濾波器原始寫法是每個(gè)抽頭一次乘加實(shí)測(cè)約12個(gè)周期一個(gè)抽頭。改成讓編譯器做軟件流水之后穩(wěn)定在1個(gè)周期一個(gè)抽頭同一個(gè)函數(shù)提速12倍代碼一行沒改邏輯只是加了幾個(gè)限定符和pragma。這種優(yōu)化性價(jià)比遠(yuǎn)超你自己手寫匯編。6.3 內(nèi)建函數(shù)的使用邊界編譯器提供的內(nèi)建函數(shù)intrinsics能直接映射到特殊指令比如雙字加載、點(diǎn)積、飽和加法。用它們的前提是你已經(jīng)確認(rèn)這里是熱點(diǎn)。我見過有人一上來(lái)全文用內(nèi)建函數(shù)結(jié)果代碼沒法讀、換平臺(tái)全廢、性能也沒提升多少因?yàn)槠款i根本不在那些地方。7. 常見問題速查這些坑我基本都踩過現(xiàn)象最可能的原因排查動(dòng)作仿真器連不上目標(biāo)板時(shí)鐘沒起振、復(fù)位被拉住、啟動(dòng)模式引腳配置錯(cuò)誤先用示波器量晶振和復(fù)位腳電平再核對(duì)啟動(dòng)模式跳線程序能下載但一運(yùn)行就飛棧溢出、.bss段超出片上存儲(chǔ)、看門狗未喂看鏈接生成的map文件確認(rèn)各段實(shí)際占用把棧臨時(shí)調(diào)大驗(yàn)證DMA搬的數(shù)據(jù)整體錯(cuò)位緩沖區(qū)未按cache line對(duì)齊用#pragma DATA_ALIGN顯式對(duì)齊源和目的都要中斷進(jìn)不去中斷使能位沒開、向量表基址沒設(shè)、外設(shè)中斷標(biāo)志未清逐級(jí)檢查全局使能、IER、外設(shè)使能三層FFT結(jié)果和理論值對(duì)不上輸入順序未做位反轉(zhuǎn)、縮放因子漏乘、輸入實(shí)部虛部搞反先用一組已知的單頻正弦驗(yàn)證再看幅值是否差固定倍數(shù)程序在外部存儲(chǔ)里跑明顯變慢未開cache且代碼段放在了DDR把.text段挪回片上SRAM或?qū)π⊙h(huán)啟用cache浮點(diǎn)累加結(jié)果偏差越來(lái)越大單精度有效位不足長(zhǎng)序列累加誤差累積改用分塊累加或Kahan求和關(guān)鍵處提升到雙精度串口打印卡死輸出函數(shù)依賴的底層驅(qū)動(dòng)沒初始化或阻塞等待標(biāo)志位先用GPIO翻轉(zhuǎn)代替打印確認(rèn)程序走到哪一步上電后偶發(fā)啟動(dòng)失敗電源時(shí)序、PLL鎖定時(shí)間不足、復(fù)位保持時(shí)間不夠逐項(xiàng)對(duì)照手冊(cè)電氣特性用示波器抓上電波形再補(bǔ)兩個(gè)文檔里不會(huì)寫的經(jīng)驗(yàn)。第一個(gè)新工程第一次跑先把所有中斷關(guān)掉只用主循環(huán)點(diǎn)燈確認(rèn)最基本的代碼執(zhí)行沒問題再一層層加中斷和外設(shè)。一次全開出問題你根本不知道是哪一層。第二個(gè)養(yǎng)成看map文件的習(xí)慣鏈接器會(huì)列出每個(gè)段的確切地址和長(zhǎng)度、每個(gè)函數(shù)的地址程序跑飛的時(shí)候用異常返回的地址去map文件里查能直接定位到是哪一行代碼觸發(fā)的比盲猜快得多。8. 學(xué)習(xí)路線的取舍哪些該啃哪些可以先跳過市面上講DSP的書分兩類一類是信號(hào)處理理論一類是芯片手冊(cè)。新手最容易犯的錯(cuò)是抱著理論書從第一章啃起啃到第六章發(fā)現(xiàn)自己連開發(fā)環(huán)境都還沒裝。我的建議路線是這樣的先花兩天把芯片手冊(cè)里「CPU內(nèi)核」和「存儲(chǔ)映射」這兩章讀完不求全懂但要建立起「這個(gè)芯片有哪些資源、它們?cè)谀?、有多快」的空間感。接著裝環(huán)境、跑通廠商提供的定時(shí)器例程自己改參數(shù)、用示波器驗(yàn)證。這一步做完你就算入門了。然后用「三板斧」串起整個(gè)學(xué)習(xí)過程定時(shí)器中斷 → DMA搬數(shù) → FIR濾波。這三個(gè)做完你已經(jīng)覆蓋了實(shí)時(shí)性、數(shù)據(jù)流、算法實(shí)現(xiàn)三大塊后面學(xué)FFT、學(xué)控制算法都是在這個(gè)骨架上掛東西。反過來(lái)先學(xué)FFT蝶形推導(dǎo)你會(huì)發(fā)現(xiàn)推導(dǎo)看懂了代碼還是寫不出來(lái)。至于算法理論夠用就行。做音頻的先把窗函數(shù)和頻譜泄漏搞明白做電機(jī)的先把坐標(biāo)變換和環(huán)路帶寬搞明白別想著一次學(xué)全。DSP這門手藝的核心競(jìng)爭(zhēng)力不在數(shù)學(xué)推導(dǎo)在于你能不能在有限的周期和存儲(chǔ)里把一個(gè)理論上成立的算法可靠地跑起來(lái)。前者查書就有后者只能靠一次次實(shí)測(cè)積累。最后分享一個(gè)我自己的習(xí)慣每做一個(gè)新模塊我都會(huì)在代碼里留一段「性能記錄」注釋寫清楚這段代碼在什么主頻、什么存儲(chǔ)配置下跑了多少周期、占了多少預(yù)算。半年后回頭看這些數(shù)字比任何筆記都有價(jià)值因?yàn)樗鼈兪悄阍谡鎸?shí)硬件上量出來(lái)的換平臺(tái)可以直接按比例推算也能立刻判斷出「這個(gè)需求在當(dāng)前硬件上到底能不能做」。踩過的坑記下來(lái)會(huì)過期量過的數(shù)據(jù)不會(huì)。 DSP芯片開發(fā)實(shí)戰(zhàn)：存儲(chǔ)架構(gòu)、DMA與實(shí)時(shí)性優(yōu)化入門
1. 先別急著寫代碼搞清楚DSP到底指什么剛接觸這塊的時(shí)候我被問得最多的一句話是「學(xué)DSP是不是得先把信號(hào)與系統(tǒng)啃一遍」。問這話的人里有一半連示波器都沒怎么摸過另一半則是從單片機(jī)轉(zhuǎn)過來(lái)、看到「哈佛架構(gòu)」「流水線」「EDMA」這幾個(gè)詞就頭皮發(fā)麻。我特別理解這種感受因?yàn)槲易约寒?dāng)年也是從一份全英文的數(shù)據(jù)手冊(cè)開始翻了三頁(yè)就合上了。先把最容易混淆的概念拆開。DSP這三個(gè)字母在不同語(yǔ)境下指三樣完全不同的東西第一種是芯片一類專門為乘加運(yùn)算和實(shí)時(shí)數(shù)據(jù)流設(shè)計(jì)的處理器第二種是學(xué)科數(shù)字信號(hào)處理這套數(shù)學(xué)理論第三種是開發(fā)范式一種以確定性時(shí)序、數(shù)據(jù)流驅(qū)動(dòng)、資源極度受限為前提的工程方法。小白最容易踩的坑就是把三者當(dāng)成一件事結(jié)果去啃了三個(gè)月的濾波器推導(dǎo)回頭連開發(fā)環(huán)境都裝不利索。這篇文章講的是第一種和第三種——怎么在一塊DSP芯片上把程序跑起來(lái)、跑對(duì)、跑快第二種只講夠用的部分。適合誰(shuí)看如果你有單片機(jī)基礎(chǔ)能看懂C語(yǔ)言和基本的電路圖想在音頻、電機(jī)控制、雷達(dá)前端、工業(yè)檢測(cè)這類場(chǎng)景里做實(shí)時(shí)信號(hào)處理那這篇內(nèi)容就是給你準(zhǔn)備的。零基礎(chǔ)也能跟我會(huì)把每個(gè)術(shù)語(yǔ)都落到具體操作上。2. 選型和硬件準(zhǔn)備別在第一步就把錢浪費(fèi)掉2.1 為什么不用普通單片機(jī)硬扛很多人第一反應(yīng)是「我用一顆主頻兩三百兆的單片機(jī)不行嗎」。短答案看你的算法預(yù)算。長(zhǎng)答案得從架構(gòu)說(shuō)起。通用單片機(jī)的核心思路是「什么活都能干」指令集通用、分支預(yù)測(cè)復(fù)雜、緩存層級(jí)深。而DSP的核心思路是「把乘加做到極致」單周期完成一次乘加MAC、硬件零開銷循環(huán)、位反轉(zhuǎn)尋址專門服務(wù)FFT、飽和運(yùn)算防止溢出、多總線哈佛結(jié)構(gòu)讓取指和取數(shù)互不打架。舉個(gè)實(shí)測(cè)量的對(duì)比同樣做1024點(diǎn)復(fù)數(shù)FFT一顆200MHz左右的通用MCU用CMSIS-DSP庫(kù)跑下來(lái)大概在20萬(wàn)到30萬(wàn)個(gè)時(shí)鐘周期量級(jí)而一顆600MHz的浮點(diǎn)DSP能壓到1萬(wàn)個(gè)周期出頭。差了將近兩個(gè)數(shù)量級(jí)這不只是主頻的差距是架構(gòu)層面的差距。更關(guān)鍵的是確定性。實(shí)時(shí)系統(tǒng)里最怕的不是慢是「有時(shí)候快有時(shí)候慢」。通用MCU的緩存命中與否會(huì)讓同一段代碼的執(zhí)行時(shí)間浮動(dòng)好幾倍而DSP可以把L1配置成SRAM而不是cache讓每次訪問的周期數(shù)固定下來(lái)。這個(gè)特性在電機(jī)環(huán)路、音頻采樣這種硬實(shí)時(shí)場(chǎng)景里是剛需。注意不要因?yàn)椤窪SP聽起來(lái)高級(jí)」就無(wú)腦上。如果你的算法是簡(jiǎn)單的PID加幾個(gè)濾波幾百兆的MCU完全夠用開發(fā)成本還低得多。DSP的價(jià)值在算法復(fù)雜度真正壓到通用核喘不過氣的時(shí)候才體現(xiàn)。2.2 開發(fā)板怎么挑看三個(gè)硬指標(biāo)第一塊板子選錯(cuò)后面全是坑。我總結(jié)下來(lái)就看三樣考察項(xiàng)具體看什么為什么重要外設(shè)閉環(huán)是否同時(shí)有ADC/DAC或音頻編解碼器、至少一個(gè)外部存儲(chǔ)接口、串口、若干GPIO沒有輸入輸出你只能跑純計(jì)算驗(yàn)證不了實(shí)時(shí)性資料完整度是否提供原理圖不是框圖是完整原理圖、BOM、成套例程出問題時(shí)能對(duì)著電路查信號(hào)這是救命的東西調(diào)試鏈路仿真器是板載還是外接、接口定義是否標(biāo)準(zhǔn)便宜的板子經(jīng)常省掉這部分自己配很折騰我并不推薦一上來(lái)就自己畫板。DSP的電源時(shí)序、內(nèi)核電壓、PLL環(huán)路濾波、去耦電容布局任何一個(gè)細(xì)節(jié)出錯(cuò)都可能導(dǎo)致「程序能下載但一跑就飛」。先在成熟評(píng)估板上把軟件跑通等你能獨(dú)立定位「是硬件問題還是軟件問題」了再考慮自制最小系統(tǒng)。2.3 自制最小系統(tǒng)時(shí)的三個(gè)致命細(xì)節(jié)如果你確實(shí)要自己畫這三處必須盯死。上電時(shí)序。多數(shù)DSP要求內(nèi)核電壓先于或同時(shí)于IO電壓建立反過來(lái)的話內(nèi)部ESD二極管可能倒灌長(zhǎng)期會(huì)損傷芯片。穩(wěn)妥做法是用帶時(shí)序控制的電源管理芯片或者至少讓兩路電源的使能信號(hào)有明確的先后。我見過一塊板子因?yàn)閮陕冯娫赐瑫r(shí)上電、上升沿斜率不同導(dǎo)致大約每十次上電有一次啟動(dòng)失敗查了兩周才定位。復(fù)位保持時(shí)間。復(fù)位信號(hào)必須在電源穩(wěn)定之后還要保持一段時(shí)間具體數(shù)值查手冊(cè)的電氣特性表。太短會(huì)導(dǎo)致內(nèi)部PLL沒鎖住就釋放復(fù)位表現(xiàn)是「偶爾能啟動(dòng)大部分時(shí)間卡死」。時(shí)鐘與PLL。外部晶振的負(fù)載電容要按晶振廠家的規(guī)格配不是隨便放兩個(gè)22pF。PLL的環(huán)路濾波電阻電容更要嚴(yán)格按照手冊(cè)給的參數(shù)這顆電容選錯(cuò)輸出時(shí)鐘抖動(dòng)會(huì)大到讓高速外設(shè)工作異常。心得新板子第一次上電先不裝芯片或者先只焊電源部分用示波器把每一路電壓的上升時(shí)間和紋波看一遍。紋波超過手冊(cè)要求的一半就要改別等調(diào)試階段再回頭懷疑電源。3. 開發(fā)環(huán)境落地從裝軟件到點(diǎn)亮第一個(gè)中斷3.1 工具鏈的四個(gè)組成部分一套完整的DSP開發(fā)環(huán)境包含四塊集成開發(fā)環(huán)境編輯器項(xiàng)目管理調(diào)試前端、C編譯器含優(yōu)化器和匯編器、鏈接器配合鏈接命令文件工作、以及芯片廠商提供的外設(shè)驅(qū)動(dòng)庫(kù)和算法庫(kù)。這四塊缺一不可尤其后面兩個(gè)是新手最容易忽略的。我的建議是第一次安裝全部用默認(rèn)路徑、默認(rèn)版本。工具鏈版本不匹配導(dǎo)致的詭異問題比如某個(gè)庫(kù)編譯通過但鏈接報(bào)符號(hào)沖突排查成本極高新手根本扛不住。工程目錄我習(xí)慣這樣組織project/ ├── src/ 源文件 ├── include/ 頭文件 ├── cmd/ 鏈接命令文件 ├── lib/ 第三方庫(kù)和廠商庫(kù) ├── obj/ 編譯中間產(chǎn)物不進(jìn)版本庫(kù) └── tools/ 燒寫腳本、輔助工具3.2 鏈接命令文件DSP開發(fā)的第一道門檻這是從單片機(jī)轉(zhuǎn)過來(lái)的人最先懵的地方。通用MCU通常由IDE自動(dòng)生成分散加載文件你基本不用管而DSP的鏈接命令文件.cmd必須自己寫因?yàn)樗苯記Q定了代碼和數(shù)據(jù)放在哪塊存儲(chǔ)上而這個(gè)選擇對(duì)性能的影響可能有十倍。MEMORY { L2SRAM : o 0x00800000, l 0x00040000 /* 片上 256KB */ DDR : o 0x80000000, l 0x04000000 /* 外部 64MB */ } SECTIONS { .text L2SRAM /* 代碼放片上取指快且確定 */ .cinit L2SRAM .const L2SRAM .switch L2SRAM .bss L2SRAM .stack L2SRAM .data DDR /* 大塊數(shù)據(jù)放外部 */ .far DDR }為什么這么分.text代碼段放在片上SRAM取指延遲固定中斷響應(yīng)時(shí)間才可預(yù)測(cè).data、.far這些放大數(shù)組的段丟到外部DDR容量大但訪問慢正好用DMA去搬運(yùn)。你要是把代碼放DDR里跑每次取指都可能撞上總線競(jìng)爭(zhēng)實(shí)測(cè)性能能掉到原來(lái)的三分之一。三個(gè)必須檢查的點(diǎn)一是各段的長(zhǎng)度加總別超過對(duì)應(yīng)存儲(chǔ)器的l值鏈接器會(huì)報(bào)錯(cuò)但報(bào)錯(cuò)信息經(jīng)常很含糊自己先算一遍二是棧大小要顯式指定默認(rèn)值往往只有幾KB遞歸或者大局部數(shù)組一用就溢出三是DMA緩沖區(qū)必須按cache line大小對(duì)齊常見是128字節(jié)用#pragma DATA_ALIGN(buffer, 128)聲明這個(gè)后面第5章還會(huì)細(xì)說(shuō)。3.3 第一個(gè)可驗(yàn)證的工程定時(shí)器中斷翻轉(zhuǎn)IO別一上來(lái)就寫FFT。第一個(gè)工程的目標(biāo)是驗(yàn)證工具鏈、時(shí)鐘、中斷三條鏈路都通最經(jīng)濟(jì)的方式就是定時(shí)器中斷里翻轉(zhuǎn)一個(gè)GPIO用示波器看周期對(duì)不對(duì)。初始化順序是有講究的全局關(guān)中斷清IER寄存器清CSR里的全局中斷使能位配置PLL設(shè)置倍頻系數(shù)然后輪詢PLL狀態(tài)位直到鎖定初始化定時(shí)器周期寄存器把中斷服務(wù)函數(shù)地址填入中斷向量表設(shè)置中斷向量表基址寄存器使能定時(shí)器中斷對(duì)應(yīng)的IER位全局開中斷啟動(dòng)定時(shí)器/* 定時(shí)器周期計(jì)算示例目標(biāo) 1kHz 中斷定時(shí)器輸入時(shí)鐘 75MHz */ #define TIMER_INPUT_CLK 75000000UL #define TARGET_FREQ 1000UL /* 周期值 輸入時(shí)鐘 / 目標(biāo)頻率 - 1 */ #define PERIOD_VAL (TIMER_INPUT_CLK / TARGET_FREQ - 1) void timer_init(void) { /* 先停定時(shí)器再改配置運(yùn)行中改寄存器行為不確定 */ TIMER_CTL 0x0; TIMER_PRD PERIOD_VAL; /* 75000 - 1 */ TIMER_CTL 0x3; /* 內(nèi)部時(shí)鐘源 連續(xù)計(jì)數(shù)模式 */ }PERIOD_VAL那個(gè)減一是因?yàn)橛?jì)數(shù)器從0開始計(jì)滿回繞差了這一個(gè)數(shù)頻率就有千分之一的偏差對(duì)音頻這種場(chǎng)景是能聽出來(lái)的。拿到板子先用示波器量IO口理論1kHz實(shí)際量出來(lái)是999.8Hz這個(gè)誤差來(lái)自晶振精度屬于正常范圍如果量出來(lái)是500Hz或者2kHz那就是分頻配置理解錯(cuò)了回去查定時(shí)器的預(yù)分頻字段。4. 存儲(chǔ)架構(gòu)和數(shù)據(jù)搬移DSP開發(fā)真正的分水嶺4.1 三級(jí)存儲(chǔ)該怎么配主流DSP的存儲(chǔ)層次一般是L1程序SRAM、L1數(shù)據(jù)SRAM兩者都可以配置成cache或者純SRAM、L2統(tǒng)一SRAM、外部DDR。配置策略直接決定了程序是「可預(yù)測(cè)地快」還是「平均來(lái)看還行」。我的默認(rèn)配置是L1P和L1D全部配置成SRAM不啟用cache。理由是實(shí)時(shí)系統(tǒng)里「最壞情況執(zhí)行時(shí)間」比「平均執(zhí)行時(shí)間」重要得多。cache看著命中率高但那1%的miss帶來(lái)的抖動(dòng)可能就是你電機(jī)環(huán)路失控或者音頻爆音的原因。代價(jià)是什么你要手動(dòng)把熱數(shù)據(jù)搬進(jìn)L1D。這就引出了DMA——整個(gè)DSP開發(fā)里最重要的一個(gè)外設(shè)。4.2 DMA和乒乓緩沖把搬數(shù)和算數(shù)徹底解耦DMA的價(jià)值在于「CPU在算當(dāng)前幀的時(shí)候下一幀已經(jīng)在搬了」。標(biāo)準(zhǔn)做法是雙緩沖#pragma DATA_ALIGN(buf_a, 128) #pragma DATA_ALIGN(buf_b, 128) static float buf_a[FRAME_LEN]; static float buf_b[FRAME_LEN]; volatile int ready_flag 0; static float * volatile active_buf buf_a; static float * volatile pending_buf buf_b;中斷服務(wù)函數(shù)里只做一件事切換指針、置標(biāo)志、重新掛DMA。真正的算法處理放在主循環(huán)或者低優(yōu)先級(jí)任務(wù)里。這樣中斷服務(wù)函數(shù)的執(zhí)行時(shí)間被壓到幾十個(gè)周期不會(huì)被算法拖長(zhǎng)。4.3 一個(gè)完整的參數(shù)推算過程假設(shè)你要做一個(gè)48kHz采樣率的音頻處理每幀256個(gè)采樣點(diǎn)主頻600MHz。這一串?dāng)?shù)字怎么落到代碼里幀時(shí)長(zhǎng) 256 ÷ 48000 ≈5.33ms這幀內(nèi)可用的CPU周期 5.33ms × 600MHz ≈3.2M 個(gè)周期工程上只敢用掉30%到40%剩下的留給中斷、DMA搶占和其他任務(wù)所以實(shí)際算法預(yù)算約1.0M 到 1.3M 個(gè)周期緩沖區(qū)字節(jié)數(shù) 256 × 4單精度f(wàn)loat1KB雙緩沖就是2KB對(duì)齊要求按128字節(jié)對(duì)齊用ceil(1024/128)×128算下來(lái)正好是1024沒問題但如果幀長(zhǎng)是250點(diǎn)就是1000字節(jié)必須向上補(bǔ)到1024字節(jié)聲明數(shù)組時(shí)直接給1024這個(gè)預(yù)算是你后面做算法選型的硬約束。比如一個(gè)256點(diǎn)的復(fù)數(shù)FFT如果在你的平臺(tái)上需要50萬(wàn)個(gè)周期那一幀只能跑兩次就必須考慮降到128點(diǎn)或者換更高效的實(shí)現(xiàn)。先算預(yù)算再寫代碼比寫完發(fā)現(xiàn)跑不動(dòng)再回來(lái)優(yōu)化省十倍時(shí)間。踩過的坑有一次DMA搬進(jìn)來(lái)的數(shù)據(jù)總是整體偏移幾個(gè)元素查了很久最后發(fā)現(xiàn)是源地址沒按元素大小對(duì)齊。DMA對(duì)地址對(duì)齊的要求比CPU寬松但對(duì)齊沒做好就會(huì)出現(xiàn)這種「錯(cuò)位」而不是「報(bào)錯(cuò)」的靜默故障特別難查。5. 定點(diǎn)還是浮點(diǎn)這道選擇題決定了你的整個(gè)代碼風(fēng)格5.1 Q格式說(shuō)到底就是一個(gè)縮放約定早期DSP沒有浮點(diǎn)單元全用定點(diǎn)。現(xiàn)在中高端DSP基本都帶單精度甚至雙精度浮點(diǎn)但定點(diǎn)仍然在成本和功耗敏感的場(chǎng)合大量使用而且理解定點(diǎn)對(duì)你理解溢出、飽和、精度損失非常有幫助。Q格式的本質(zhì)是「用整數(shù)表示小數(shù)」。Q15的意思是用16位有符號(hào)整數(shù)隱含小數(shù)點(diǎn)在第15位之后。格式位寬表示范圍最小分辨率實(shí)際值換算Q1516位-1.0 ~ 0.999971/32768 ≈ 0.0000305整數(shù) ÷ 32768Q3132位-1.0 ~ 0.99999999951/2^31整數(shù) ÷ 2147483648Q1216位-8.0 ~ 7.99981/4096 ≈ 0.000244整數(shù) ÷ 4096運(yùn)算規(guī)則要記牢加減法直接算但結(jié)果可能溢出乘法要先擴(kuò)位再右移。兩個(gè)Q15相乘結(jié)果是一個(gè)Q30的數(shù)必須右移15位才能回到Q15同時(shí)要做飽和處理防止溢出。static inline short q15_mul(short a, short b) { int prod (int)a * (int)b; /* Q30 */ prod (prod 0x4000) 15; /* 加半個(gè)LSB做四舍五入 */ if (prod 32767) prod 32767; if (prod -32768) prod -32768; return (short)prod; }那個(gè)0x4000是四舍五入用的偏置加不加在單次運(yùn)算里看不出差別但在濾波器這種幾百次累加的場(chǎng)景里會(huì)導(dǎo)致明顯的直流偏移。我第一次做FIR的時(shí)候沒加輸出信號(hào)整體偏了幾十個(gè)LSB還以為是ADC的問題。5.2 飽和處理寧可削頂不要翻轉(zhuǎn)定點(diǎn)運(yùn)算最怕的不是精度不夠是溢出翻轉(zhuǎn)。一個(gè)正的大數(shù)溢出會(huì)變成負(fù)數(shù)在你的音頻輸出里就是一聲刺耳的爆響在電機(jī)控制里可能就是一次危險(xiǎn)的誤動(dòng)作。所以所有可能溢出的地方都要做飽和clamp到最大最小值而不是讓它自然回繞。很多DSP有專門的飽和加法指令編譯器也提供了內(nèi)建函數(shù)能用一個(gè)周期搞定別自己寫if-else。5.3 實(shí)測(cè)對(duì)比同一段算法的兩種實(shí)現(xiàn)我在一塊600MHz的浮點(diǎn)DSP上做過一組對(duì)比32點(diǎn)復(fù)數(shù)FFT輸入數(shù)據(jù)已經(jīng)就位只算內(nèi)核循環(huán)實(shí)現(xiàn)方式平均周期數(shù)相對(duì)精度誤差代碼可讀性單精度浮點(diǎn)開 -O3約 1,050基準(zhǔn)高Q15定點(diǎn)純C約 2,300約 -62dB中需注釋Q15定點(diǎn)內(nèi)建函數(shù) 循環(huán)展開約 780約 -60dB低有意思的地方在于定點(diǎn)用了內(nèi)建函數(shù)之后反而比浮點(diǎn)快因?yàn)楦↑c(diǎn)單元在這個(gè)平臺(tái)上只有兩個(gè)定點(diǎn)單元有八個(gè)。但代價(jià)是代碼幾乎沒法維護(hù)而且每換一個(gè)平臺(tái)就要重寫。所以我的實(shí)踐原則是原型階段全用浮點(diǎn)把算法跑通、把精度驗(yàn)證過等確定性能不夠了再把最內(nèi)層的一兩個(gè)熱點(diǎn)函數(shù)改定點(diǎn)。一上來(lái)就全手寫定點(diǎn)優(yōu)化返工率極高。6. 實(shí)時(shí)性保障中斷延遲和流水線優(yōu)化6.1 中斷延遲怎么測(cè)別靠猜中斷延遲是從外部事件觸發(fā)到中斷服務(wù)函數(shù)第一條指令執(zhí)行的時(shí)間。這個(gè)數(shù)字直接決定你的系統(tǒng)能不能滿足實(shí)時(shí)要求但很多人從沒量過全憑感覺。測(cè)法很樸素在中斷服務(wù)函數(shù)的第一行立刻翻轉(zhuǎn)一個(gè)GPIO不要做任何其他事用示波器雙通道一個(gè)通道接外部觸發(fā)信號(hào)一個(gè)通道接這個(gè)GPIO看兩個(gè)上升沿的時(shí)間差。多測(cè)幾百次取最大值而不是平均值。實(shí)測(cè)經(jīng)驗(yàn)中斷延遲從幾十個(gè)周期到幾百個(gè)周期都可能主要取決于三件事中斷是否被更長(zhǎng)的高優(yōu)先級(jí)中斷屏蔽、中斷向量表是否在片上、以及中斷服務(wù)函數(shù)的第一條指令有沒有撞上取指miss。我測(cè)過一個(gè)案例把中斷向量表從外部存儲(chǔ)搬到片上SRAM之后最壞延遲從400多個(gè)周期降到80個(gè)周期左右效果比優(yōu)化算法本身明顯得多。關(guān)鍵提醒所有中斷服務(wù)函數(shù)都必須能在你算出來(lái)的「最壞延遲 執(zhí)行時(shí)間」預(yù)算內(nèi)完成。超了就要拆分成上半部和下半部上半部只做搬數(shù)據(jù)和置標(biāo)志。6.2 讓編譯器替你展開流水線超長(zhǎng)指令字架構(gòu)的DSP一個(gè)周期能發(fā)射多條指令前提是編譯器能把循環(huán)展開成軟件流水。你要做的是給它創(chuàng)造條件循環(huán)體內(nèi)不要有函數(shù)調(diào)用能內(nèi)聯(lián)就內(nèi)聯(lián)指針參數(shù)加restrict關(guān)鍵字告訴編譯器源和目的不重疊用#pragma MUST_ITERATE告訴編譯器循環(huán)次數(shù)的范圍它才敢放心展開循環(huán)次數(shù)最好是展開因子的整數(shù)倍或者手動(dòng)補(bǔ)零對(duì)齊優(yōu)化有沒有生效別猜去看編譯器生成的匯編文件。開啟「保留匯編輸出」選項(xiàng)找到你的熱點(diǎn)函數(shù)看里面有沒有那個(gè)標(biāo)志性的軟件流水序言一串并行的裝載指令和排空段。如果生成的還是一條一條串行的指令說(shuō)明優(yōu)化沒起來(lái)通常是循環(huán)里有編譯器看不透的依賴。我優(yōu)化過一個(gè)FIR濾波器原始寫法是每個(gè)抽頭一次乘加實(shí)測(cè)約12個(gè)周期一個(gè)抽頭。改成讓編譯器做軟件流水之后穩(wěn)定在1個(gè)周期一個(gè)抽頭同一個(gè)函數(shù)提速12倍代碼一行沒改邏輯只是加了幾個(gè)限定符和pragma。這種優(yōu)化性價(jià)比遠(yuǎn)超你自己手寫匯編。6.3 內(nèi)建函數(shù)的使用邊界編譯器提供的內(nèi)建函數(shù)intrinsics能直接映射到特殊指令比如雙字加載、點(diǎn)積、飽和加法。用它們的前提是你已經(jīng)確認(rèn)這里是熱點(diǎn)。我見過有人一上來(lái)全文用內(nèi)建函數(shù)結(jié)果代碼沒法讀、換平臺(tái)全廢、性能也沒提升多少因?yàn)槠款i根本不在那些地方。7. 常見問題速查這些坑我基本都踩過現(xiàn)象最可能的原因排查動(dòng)作仿真器連不上目標(biāo)板時(shí)鐘沒起振、復(fù)位被拉住、啟動(dòng)模式引腳配置錯(cuò)誤先用示波器量晶振和復(fù)位腳電平再核對(duì)啟動(dòng)模式跳線程序能下載但一運(yùn)行就飛棧溢出、.bss段超出片上存儲(chǔ)、看門狗未喂看鏈接生成的map文件確認(rèn)各段實(shí)際占用把棧臨時(shí)調(diào)大驗(yàn)證DMA搬的數(shù)據(jù)整體錯(cuò)位緩沖區(qū)未按cache line對(duì)齊用#pragma DATA_ALIGN顯式對(duì)齊源和目的都要中斷進(jìn)不去中斷使能位沒開、向量表基址沒設(shè)、外設(shè)中斷標(biāo)志未清逐級(jí)檢查全局使能、IER、外設(shè)使能三層FFT結(jié)果和理論值對(duì)不上輸入順序未做位反轉(zhuǎn)、縮放因子漏乘、輸入實(shí)部虛部搞反先用一組已知的單頻正弦驗(yàn)證再看幅值是否差固定倍數(shù)程序在外部存儲(chǔ)里跑明顯變慢未開cache且代碼段放在了DDR把.text段挪回片上SRAM或?qū)π⊙h(huán)啟用cache浮點(diǎn)累加結(jié)果偏差越來(lái)越大單精度有效位不足長(zhǎng)序列累加誤差累積改用分塊累加或Kahan求和關(guān)鍵處提升到雙精度串口打印卡死輸出函數(shù)依賴的底層驅(qū)動(dòng)沒初始化或阻塞等待標(biāo)志位先用GPIO翻轉(zhuǎn)代替打印確認(rèn)程序走到哪一步上電后偶發(fā)啟動(dòng)失敗電源時(shí)序、PLL鎖定時(shí)間不足、復(fù)位保持時(shí)間不夠逐項(xiàng)對(duì)照手冊(cè)電氣特性用示波器抓上電波形再補(bǔ)兩個(gè)文檔里不會(huì)寫的經(jīng)驗(yàn)。第一個(gè)新工程第一次跑先把所有中斷關(guān)掉只用主循環(huán)點(diǎn)燈確認(rèn)最基本的代碼執(zhí)行沒問題再一層層加中斷和外設(shè)。一次全開出問題你根本不知道是哪一層。第二個(gè)養(yǎng)成看map文件的習(xí)慣鏈接器會(huì)列出每個(gè)段的確切地址和長(zhǎng)度、每個(gè)函數(shù)的地址程序跑飛的時(shí)候用異常返回的地址去map文件里查能直接定位到是哪一行代碼觸發(fā)的比盲猜快得多。8. 學(xué)習(xí)路線的取舍哪些該啃哪些可以先跳過市面上講DSP的書分兩類一類是信號(hào)處理理論一類是芯片手冊(cè)。新手最容易犯的錯(cuò)是抱著理論書從第一章啃起啃到第六章發(fā)現(xiàn)自己連開發(fā)環(huán)境都還沒裝。我的建議路線是這樣的先花兩天把芯片手冊(cè)里「CPU內(nèi)核」和「存儲(chǔ)映射」這兩章讀完不求全懂但要建立起「這個(gè)芯片有哪些資源、它們?cè)谀?、有多快」的空間感。接著裝環(huán)境、跑通廠商提供的定時(shí)器例程自己改參數(shù)、用示波器驗(yàn)證。這一步做完你就算入門了。然后用「三板斧」串起整個(gè)學(xué)習(xí)過程定時(shí)器中斷 → DMA搬數(shù) → FIR濾波。這三個(gè)做完你已經(jīng)覆蓋了實(shí)時(shí)性、數(shù)據(jù)流、算法實(shí)現(xiàn)三大塊后面學(xué)FFT、學(xué)控制算法都是在這個(gè)骨架上掛東西。反過來(lái)先學(xué)FFT蝶形推導(dǎo)你會(huì)發(fā)現(xiàn)推導(dǎo)看懂了代碼還是寫不出來(lái)。至于算法理論夠用就行。做音頻的先把窗函數(shù)和頻譜泄漏搞明白做電機(jī)的先把坐標(biāo)變換和環(huán)路帶寬搞明白別想著一次學(xué)全。DSP這門手藝的核心競(jìng)爭(zhēng)力不在數(shù)學(xué)推導(dǎo)在于你能不能在有限的周期和存儲(chǔ)里把一個(gè)理論上成立的算法可靠地跑起來(lái)。前者查書就有后者只能靠一次次實(shí)測(cè)積累。最后分享一個(gè)我自己的習(xí)慣每做一個(gè)新模塊我都會(huì)在代碼里留一段「性能記錄」注釋寫清楚這段代碼在什么主頻、什么存儲(chǔ)配置下跑了多少周期、占了多少預(yù)算。半年后回頭看這些數(shù)字比任何筆記都有價(jià)值因?yàn)樗鼈兪悄阍谡鎸?shí)硬件上量出來(lái)的換平臺(tái)可以直接按比例推算也能立刻判斷出「這個(gè)需求在當(dāng)前硬件上到底能不能做」。踩過的坑記下來(lái)會(huì)過期量過的數(shù)據(jù)不會(huì)。