簡(jiǎn)介三種數(shù)據(jù)輔助型定時(shí)同步算法的MATLAB仿真實(shí)現(xiàn)專為通信、信號(hào)處理領(lǐng)域的學(xué)生和研究人員準(zhǔn)備用于理解不同定時(shí)同步策略的實(shí)際效果。SC算法利用自相關(guān)特性估計(jì)符號(hào)邊界Minn算法基于分段滑窗尋找相關(guān)峰值Park算法則以匹配濾波預(yù)處理抑制噪聲三者覆蓋了數(shù)據(jù)輔助同步中的主要技術(shù)路線適合在不同信道條件和噪聲環(huán)境下對(duì)比同步性能。壓縮包共有3個(gè)文件均為.m格式的MATLAB源碼整體大小僅3KB輕量簡(jiǎn)明可在MATLAB中直接運(yùn)行并觀察同步曲線。三個(gè)腳本相互獨(dú)立可分別運(yùn)行或改參重測(cè)適合動(dòng)手驗(yàn)證算法原理也可作為課程設(shè)計(jì)與科研編碼的起始模板。已有1821人學(xué)習(xí)下載足見其在定時(shí)同步入門和工程驗(yàn)證中的實(shí)用價(jià)值。1. 數(shù)據(jù)輔助型三種定時(shí)同步算法SC、Minn、Park到底在同步什么數(shù)據(jù)輔助型三種定時(shí)同步算法——SC、Minn、Park——是OFDM接收機(jī)里繞不開的符號(hào)定時(shí)方案。我最初在某跨平臺(tái)系統(tǒng)的基帶仿真里被問“FFT窗到底從哪個(gè)采樣點(diǎn)開始切”才真正意識(shí)到這三個(gè)算法不是擺設(shè)它們都靠收發(fā)雙方約定的訓(xùn)練序列在接收端用滑動(dòng)相關(guān)構(gòu)造定時(shí)度量再從度量里找符號(hào)起始點(diǎn)。選錯(cuò)算法或調(diào)錯(cuò)參數(shù)時(shí)偏誤差會(huì)直接變成子載波干擾星座圖整個(gè)散掉。這篇筆記適合正在做OFDM鏈路仿真、或者在寫同步模塊的開發(fā)者目標(biāo)是讓你能動(dòng)手復(fù)現(xiàn)、能比較、能避開常見的坑。2. 三種算法的原理分水嶺從重復(fù)半段到共軛對(duì)稱的訓(xùn)練序列設(shè)計(jì)這一章不講復(fù)雜推導(dǎo)只講清楚三套前導(dǎo)結(jié)構(gòu)各自長(zhǎng)什么樣、度量公式為什么長(zhǎng)那樣以及它們最關(guān)鍵的差異點(diǎn)。后面所有仿真和排錯(cuò)都依賴這章的基本判斷。2.1 SC算法兩個(gè)重復(fù)半段和那個(gè)讓人又愛又恨的平頂SC算法的訓(xùn)練符號(hào)由兩個(gè)等長(zhǎng)的半段組成前半段記為A后半段也是A兩段完全一樣每段長(zhǎng)度是N/2。接收端利用這個(gè)重復(fù)性用一個(gè)長(zhǎng)度為N/2的滑動(dòng)窗把當(dāng)前位置d之后的第一個(gè)半段和間隔N/2之后的第二個(gè)半段做相關(guān)。度量表達(dá)式可以寫成P(d) Σ_{m0}^{N/2-1} r(dm) · r*(dmN/2)R(d) Σ_{m0}^{N/2-1} |r(dmN/2)|2M(d) |P(d)|2 / (R(d))2這里r是接收序列*表示共軛。P(d)度量的是兩個(gè)重復(fù)半段的相關(guān)能量R(d)做歸一化把度量約束在0到1之間方便和門限比較。這套設(shè)計(jì)最直觀的優(yōu)點(diǎn)是對(duì)載波頻偏不敏感。因?yàn)轭l偏只會(huì)在相鄰半段的樣本上引入一個(gè)恒定的相位差P(d)的模長(zhǎng)不受影響。但缺點(diǎn)也很明顯當(dāng)滑動(dòng)窗的起點(diǎn)落在循環(huán)前綴和第一個(gè)半段之間時(shí)窗內(nèi)后半段和前半段仍然保持著重復(fù)結(jié)構(gòu)M(d)會(huì)連續(xù)保持一個(gè)接近1的峰值形成一個(gè)“平頂”。平臺(tái)長(zhǎng)度差不多等于循環(huán)前綴長(zhǎng)度這就讓峰值位置變得模糊argmax之后每次估計(jì)都會(huì)在平臺(tái)內(nèi)隨機(jī)跳。SC的平頂在單徑AWGN信道里還勉強(qiáng)能忍入了多徑信道以后平臺(tái)會(huì)被延遲副本抬得更加復(fù)雜定時(shí)誤差方差會(huì)明顯變大。所以在實(shí)際仿真里很少有人直接把a(bǔ)rgmax(M)當(dāng)最終結(jié)果多數(shù)會(huì)先取平臺(tái)再取平臺(tái)中點(diǎn)或者重心。2.2 Minn算法四段符號(hào)與取負(fù)改造如何削掉平臺(tái)Minn算法的出發(fā)點(diǎn)是干掉SC的平臺(tái)。它把訓(xùn)練符號(hào)分成四段每段長(zhǎng)度是N/4符號(hào)模式為A、A、-A、-A。也就是說前兩段相同后兩段是前兩段的相反數(shù)。接收端把整個(gè)窗分成四段分別計(jì)算第一段和第二段的正相關(guān)、第三段和第四段的正相關(guān)再把兩個(gè)相關(guān)項(xiàng)相加P(d) Σ_{m0}^{N/4-1} r(dm) · r*(dmN/4) Σ_{m0}^{N/4-1} r(dmN/2) · r*(dm3N/4)R(d) Σ_{m0}^{N/4-1} |r(dmN/4)|2 Σ_{m0}^{N/4-1} |r(dm3N/4)|2M(d) |P(d)|2 / (R(d))2因?yàn)楹髢啥稳×素?fù)號(hào)當(dāng)滑動(dòng)窗錯(cuò)開一小段時(shí)相關(guān)項(xiàng)里會(huì)出現(xiàn)正負(fù)抵消度量會(huì)快速掉下來所以不會(huì)再出現(xiàn)一個(gè)寬平臺(tái)。這是Minn算法最核心的改進(jìn)。不過Minn算法也付出了代價(jià)它在主峰兩側(cè)會(huì)留下副峰副峰高度大約是主峰的一半。如果信噪比低或者搜索窗內(nèi)同時(shí)有多個(gè)OFDM符號(hào)副峰很容易被閾值判決誤選成主峰。這是Minn算法在工程落地時(shí)最主要的一個(gè)坑后面第五章會(huì)專門排查。2.3 Park算法用共軛對(duì)稱換一個(gè)沖激式的峰值Park算法換了一條路不再做前后段重復(fù)而是讓訓(xùn)練符號(hào)具備共軛對(duì)稱性。設(shè)前半段為A長(zhǎng)度為N/2后半段取前半段的反向共軛即fliplr(conj(A))。整個(gè)序列滿足s(n) conj(s(N1-n))這種結(jié)構(gòu)在時(shí)域上以中點(diǎn)為中心呈共軛對(duì)稱。接收度量也不再使用共軛相關(guān)而是直接把滑動(dòng)窗前半段和后半段反向求乘積P(d) Σ_{k0}^{N/2-1} r(dk) · r(dN-k)R(d) Σ_{k0}^{N/2-1} |r(dk)|2M(d) |P(d)|2 / (R(d))2注意這里的第二項(xiàng)r(dN-k)沒有取共軛。在正確的符號(hào)起點(diǎn)位置s(k)和s(N-k)互為共軛乘積正好等于|s(k)|2所有項(xiàng)同相疊加M(d)會(huì)形成一個(gè)接近單樣本寬度的尖峰比Minn更加尖銳。這個(gè)“不要取共軛”的細(xì)節(jié)是Park算法最容易寫錯(cuò)的地方。如果順手寫成conj(r(dN-k))在頻偏不為零時(shí)乘積相位會(huì)隨k變化相干積累會(huì)被破壞峰值會(huì)塌下來。第5.4節(jié)會(huì)單獨(dú)說這個(gè)翻車現(xiàn)象。3. 用 MATLAB 跑通三種定時(shí)同步算法完整仿真鏈路與核心代碼這一章的代碼可以直接復(fù)制到MATLAB里跑只保留最核心的鏈路發(fā)射端構(gòu)造前導(dǎo)、接收端算度量、定時(shí)誤差統(tǒng)計(jì)。仿真目的不是追求極限性能而是讓你自己看到三種算法的峰形曲線并得到一個(gè)可以橫向比較的MSE趨勢(shì)。3.1 仿真參數(shù)與發(fā)射端前導(dǎo)構(gòu)造我建議用下面的參數(shù)作為基線FFT點(diǎn)數(shù)不要太大否則蒙特卡洛跑起來太慢也不要太小否則循環(huán)前綴和信道時(shí)延的比例不好模擬。參數(shù)取值說明FFT點(diǎn)數(shù) N256OFDM有效符號(hào)長(zhǎng)度循環(huán)前綴 CP16常規(guī)比例約1/16數(shù)據(jù)調(diào)制QPSK便于觀察星座圖也足夠驗(yàn)證定時(shí)數(shù)據(jù)OFDM符號(hào)數(shù)10幀里前導(dǎo)后跟10個(gè)數(shù)據(jù)符號(hào)信道3徑多徑延遲[0,3,7]個(gè)采樣增益[1,0.5,0.3]歸一化頻偏0.05子載波間隔的0.05倍蒙特卡洛次數(shù)1000統(tǒng)計(jì)MSE需要足夠樣本SNR范圍0到30 dB步進(jìn)2 dB看低信噪比下的趨勢(shì)發(fā)射端前導(dǎo)生成代碼N 256; cp 16; nSym 10; snrList 0:2:30; nMonte 1000; % QPSK符號(hào)生成函數(shù) qpskSym (M) exp(1j*pi/4 1j*pi/2 * randi([0 3], M, 1)); switch alg case SC A qpskSym(N/2); preamble [A; A]; case Minn A qpskSym(N/4); preamble [A; A; -A; -A]; case Park A qpskSym(N/2); preamble [A; conj(fliplr(A))]; end % 給前導(dǎo)加上循環(huán)前綴 txPreamble [preamble(end-cp1:end); preamble];這段代碼的邏輯是三種算法共用同一個(gè)QPSK符號(hào)生成函數(shù)只是把前導(dǎo)組織成不同結(jié)構(gòu)。最后一句把前導(dǎo)最后CP個(gè)樣本復(fù)制到開頭這是OFDM符號(hào)進(jìn)入信道之前的標(biāo)準(zhǔn)操作。參數(shù)上A的長(zhǎng)度分別為N/2或N/4對(duì)應(yīng)不同分段數(shù)不能搞混。數(shù)據(jù)符號(hào)可以另外生成txData zeros(nSym * (N cp), 1); for k 0:nSym-1 sym qpskSym(N); txData(k*(Ncp)1 : k*(Ncp)cp) sym(end-cp1:end); txData(k*(Ncp)cp1 : (k1)*(Ncp)) sym; end txBurst [txPreamble; txData];這里每個(gè)數(shù)據(jù)OFDM符號(hào)的CP來自當(dāng)前符號(hào)尾部而前導(dǎo)的CP來自前導(dǎo)尾部?jī)烧叻珠_處理。仿真時(shí)把txBurst當(dāng)成一幀發(fā)送即可。3.2 三種定時(shí)度量的 MATLAB 核心實(shí)現(xiàn)寫一個(gè)獨(dú)立的函數(shù)輸入接收向量、滑動(dòng)起點(diǎn)d、算法標(biāo)識(shí)和FFT點(diǎn)數(shù)輸出當(dāng)前d處的定時(shí)度量M。這個(gè)函數(shù)只關(guān)心訓(xùn)練符號(hào)的有效部分不把CP納入相關(guān)窗所以d的搜索范圍要保證窗內(nèi)取到的是完整N個(gè)采樣。function M timingMetric(rx, d, alg, N) switch alg case SC L N / 2; a rx(d1 : dL); % 前半段 b rx(dL1 : dN); % 后半段 P sum(conj(a) .* b); R sum(abs(b).^2); case Minn L N / 4; a rx(d1 : dL); b rx(dL1 : d2*L); c rx(d2*L1 : d3*L); e rx(d3*L1 : d4*L); P sum(conj(a) .* b) sum(conj(c) .* e); R sum(abs(b).^2) sum(abs(e).^2); case Park L N / 2; c rx(d1 : dL); e rx(dN : -1 : dL1); % 后半段反向取 P sum(c .* e); % 注意這里沒有 conj R sum(abs(c).^2); end M abs(P).^2 / R.^2; end這段代碼里SC和Minn都是標(biāo)準(zhǔn)的共軛相關(guān)Park則是共軛對(duì)稱序列的直接乘積。邏輯說明SC和Minn的P都是一段能量與另一段能量的相關(guān)積累歸一化R必須用參與相關(guān)的那半段能量Park的R用前半段能量即可后半段不需要再用因?yàn)楣曹棇?duì)稱保證了兩半段能量相等。參數(shù)說明d是MATLAB下標(biāo)從1開始rx中至少要能取出dN個(gè)樣本不然取后半段反向時(shí)會(huì)越界搜索窗口通常取真實(shí)起點(diǎn)前后半個(gè)符號(hào)長(zhǎng)度比如真實(shí)起點(diǎn)是start則d從start-32到start32這樣可以看到完整的平臺(tái)或峰形。3.3 仿真結(jié)果怎么讀定時(shí)度量波形與 MSE 對(duì)比表拿到M(d)曲線后先畫出來看一眼。SC會(huì)看到一個(gè)寬平臺(tái)Minn是帶副峰的尖峰Park是幾乎只有一個(gè)采樣寬的脈沖。這是三種算法最直觀的區(qū)別也是后續(xù)所有參數(shù)判斷的基礎(chǔ)。[~, imax] max(M); if strcmp(alg, SC) % 把平臺(tái)區(qū)域找出來取中點(diǎn) idxRange find(M 0.9 * M(imax)); dEst round(mean(idxRange)); % 更穩(wěn)妥用平臺(tái)中點(diǎn)的均值 else dEst imax; end timingErr(mc) dSearch(dEst) - trueStart;這段是定時(shí)估計(jì)的核心處理。SC如果直接取imax誤差會(huì)在平臺(tái)內(nèi)部隨機(jī)抖動(dòng)先圈出接近最大值的平臺(tái)區(qū)間再取中點(diǎn)或重心能顯著降低方差。Minn和Park則直接用最大峰位置即可。三種算法的典型性能趨勢(shì)可以先用一張表記住具體數(shù)值會(huì)隨信道和隨機(jī)種子變化但相對(duì)關(guān)系基本穩(wěn)定對(duì)比項(xiàng)SCMinnPark度量峰形寬平臺(tái)尖峰帶副峰單樣本尖峰低信噪比下估計(jì)方差大中等較小對(duì)循環(huán)前綴長(zhǎng)度的敏感度高低低歸一化復(fù)雜度低略高低抗頻偏能力強(qiáng)強(qiáng)強(qiáng)前提是P不加conj4. 真實(shí)信道下三種算法的表現(xiàn)多徑、頻偏和門限設(shè)定理想AWGN信道上三個(gè)算法都能工作但實(shí)際系統(tǒng)躲不開多徑和頻偏。這章說清楚它們各自在什么條件下會(huì)退化以及門限該按什么邏輯去設(shè)。4.1 多徑信道下平臺(tái)偏移與峰值退化多徑信道里接收信號(hào)是多個(gè)延遲副本的疊加。SC的度量平臺(tái)會(huì)被延遲副本“抬高”或“填平”最大峰的位置往往偏向最強(qiáng)徑而不是第一徑。OFDM允許的定時(shí)窗口是“最大時(shí)延擴(kuò)展落在CP內(nèi)”也就是說定時(shí)點(diǎn)不能太靠后否則上一符號(hào)的拖尾會(huì)污染當(dāng)前符號(hào)。如果你直接把SC的最大峰位置當(dāng)符號(hào)起點(diǎn)而最強(qiáng)徑又晚于第一徑幾個(gè)采樣FFT窗起點(diǎn)就會(huì)偏晚等效時(shí)延擴(kuò)展變大星座圖出現(xiàn)不可收斂的旋轉(zhuǎn)和彌散。Minn和Park的主峰寬度窄對(duì)多徑的畸變相對(duì)不敏感但延遲副本也會(huì)在主峰附近產(chǎn)生毛刺。尤其是Park晚到一兩個(gè)采樣的副本會(huì)在主峰旁邊形成旁瓣如果旁瓣高度接近主峰也會(huì)被插值算法誤判。一個(gè)常見做法是先用能量粗同步把搜索范圍縮小到前導(dǎo)附近再做精細(xì)定時(shí)這樣旁瓣出現(xiàn)在搜索窗外的概率就低很多。4.2 載波頻偏對(duì)三種度量的影響以及頻偏估計(jì)從哪來只要訓(xùn)練序列的結(jié)構(gòu)不因?yàn)轭l偏而破壞定時(shí)度量的模長(zhǎng)就基本不受頻偏影響。SC的P(d)在頻偏下會(huì)額外乘一個(gè)e^{-j2πLε/N}的常數(shù)相位取模后消失Minn的兩個(gè)相關(guān)項(xiàng)也各有一個(gè)常數(shù)相位取模后同樣消失Park的P(d)雖然不帶共軛但每個(gè)乘積的相位也是常數(shù)取模后不受頻偏影響。不過有一個(gè)前提Park的代碼里不能對(duì)后半段取共軛。如果取了共軛乘積相位就會(huì)變成隨k變化的函數(shù)頻偏越大主峰衰減越明顯。這一點(diǎn)在第五章會(huì)再踩一遍。頻偏信息本身可以從P(d)的相位里估出來這也是數(shù)據(jù)輔助同步的額外福利。以SC為例P(d)的相位大約是-2πLε/N-πε所以ε估計(jì)值就是-angle(P)/π。這個(gè)估計(jì)范圍大致是±1個(gè)子載波間隔適合做小數(shù)倍頻偏粗估計(jì)整數(shù)倍頻偏還需要另外用導(dǎo)頻或序列相關(guān)性去解。4.3 幀檢測(cè)閾值怎么設(shè)SC 看平臺(tái)中點(diǎn)Minn/Park 看最大峰突發(fā)幀接收時(shí)通常先在接收流上滑動(dòng)計(jì)算M(d)再和閾值比較判定有沒有前導(dǎo)到達(dá)。閾值設(shè)太高會(huì)漏檢太低則噪聲尖峰造成虛警。我一般這樣設(shè)SC先用0.8到0.9倍的最大值圈出平臺(tái)范圍然后取平臺(tái)中點(diǎn)作為幀起點(diǎn)。Minn閾值取0.5到0.7倍的最大值找到第一越過閾值的位置再往后找局部最大峰。Park閾值取0.5倍最大值以上直接取全局最大峰。低信噪比下絕對(duì)值固定門限不可靠建議用相對(duì)門限也就是和當(dāng)前搜索窗內(nèi)的最大M做比值。更講究一點(diǎn)可以在前導(dǎo)前預(yù)留一小段空白噪聲估計(jì)噪聲功率用“最大值×0.5 噪聲功率開方×3”作為門限后面最后一章會(huì)給出具體代碼。5. 定時(shí)同步仿真避坑指南五個(gè)典型翻車現(xiàn)場(chǎng)與排查方法這章是我自己反復(fù)踩過的坑每個(gè)都按“現(xiàn)象、原因、解決”拆開。癥狀都比較典型排查時(shí)可以直接對(duì)照。5.1 SC的平頂導(dǎo)致定時(shí)點(diǎn)隨機(jī)漂移現(xiàn)象每次蒙特卡洛仿真估計(jì)的符號(hào)起點(diǎn)都在跳星座圖EVM抖動(dòng)大低信噪比尤其明顯。原因SC度量的最大值在平臺(tái)區(qū)間內(nèi)都成立argmax受噪聲影響每次都落在平臺(tái)內(nèi)不同位置。解決不要直接取argmax。先找到所有大于0.9倍最大值的索引再取索引區(qū)間的中點(diǎn)或者按M值做重心加權(quán)。平臺(tái)中點(diǎn)的估計(jì)方差比隨機(jī)取點(diǎn)小一個(gè)數(shù)量級(jí)。這一招對(duì)任何帶平臺(tái)的度量都適用。5.2 Minn的副峰被誤判為主峰現(xiàn)象定時(shí)估計(jì)結(jié)果總是差N/4或N/2個(gè)采樣低信噪比下錯(cuò)誤概率很高。原因Minn訓(xùn)練符號(hào)中有符號(hào)翻轉(zhuǎn)窗錯(cuò)開一個(gè)段時(shí)第二個(gè)相關(guān)項(xiàng)會(huì)貢獻(xiàn)部分能量形成副峰。副峰高度約為主峰一半閾值若設(shè)在0.6倍最大值以下就會(huì)把副峰也納入候選。解決閾值設(shè)高一點(diǎn)同時(shí)要求主峰寬度不超過3個(gè)采樣。更穩(wěn)妥的是先做一次粗同步把范圍限制在真實(shí)起點(diǎn)附近再在局部范圍內(nèi)找最大峰副峰落在范圍外的概率就低很多。還可以用定時(shí)度量的一階差分來判斷副峰是否出現(xiàn)在固定間隔上。5.3 Park的循環(huán)前綴帶來周期性“鬼峰”現(xiàn)象Park度量的主峰附近每隔一個(gè)CP長(zhǎng)度或一個(gè)FFT長(zhǎng)度出現(xiàn)一個(gè)小峰小峰有時(shí)候比主峰還高。原因OFDM符號(hào)的循環(huán)前綴是符號(hào)尾部的復(fù)制在Park的共軛對(duì)稱窗里CP和數(shù)據(jù)部分的邊界會(huì)產(chǎn)生局部的對(duì)稱相關(guān)形成周期性偽峰。前導(dǎo)本身也帶CP所以鬼峰幾乎是必然存在的。解決不要把CP納入對(duì)稱窗的計(jì)算范圍滑動(dòng)搜索時(shí)確保d之后正好是N個(gè)訓(xùn)練符號(hào)本體。另外搜索窗不要覆蓋多個(gè)OFDM符號(hào)只在前導(dǎo)粗同步到的位置附近搜索。如果鬼峰仍在可以對(duì)度量曲線做長(zhǎng)度約束比如主峰位置附近的M曲線應(yīng)該是一個(gè)單峰而不是連續(xù)多個(gè)峰。5.4 Park度量里多寫一個(gè)conj導(dǎo)致峰值塌了現(xiàn)象Park算法在無噪聲信道下M(d)最大值也到不了0.3主峰模糊看起來像沒同步上。原因把P(d)寫成了conj(r(dN-k))。共軛對(duì)稱序列要求的是s(k)與s(N-k)互為共軛接收端直接乘積即可。多寫conj以后信號(hào)本身的對(duì)稱性被抵消相關(guān)項(xiàng)相位隨k變化無法相干疊加。解決把conj去掉改成P sum(rx(d1:dL) .* rx(dN:-1:dL1))。驗(yàn)證方法也很簡(jiǎn)單在無噪聲、無頻偏、定時(shí)準(zhǔn)確的條件下P的實(shí)部應(yīng)該接近前半段能量虛部接近0。如果不滿足說明訓(xùn)練序列或方向取反了。5.5 歸一化分母用錯(cuò)導(dǎo)致門限失靈現(xiàn)象M(d)整體偏小最大值只有0.2或者整體偏大噪聲里也能超過0.9門限怎么調(diào)都不對(duì)。原因R(d)應(yīng)該用參與相關(guān)的那段能量。SC用后半段Minn用第二段和第四段Park用前半段。有人圖省事用整個(gè)FFT窗能量做歸一化甚至把循環(huán)前綴都算進(jìn)去分母多出一塊無關(guān)能量度量自然被壓低。如果分母里只有部分能量度量又會(huì)偏高。解決嚴(yán)格按3.2的代碼寫。如果確實(shí)想用全窗能量做簡(jiǎn)化必須重新標(biāo)定門限不要沿用按定義推導(dǎo)出來的0.7、0.9這類經(jīng)驗(yàn)值。6. 把仿真算法推向工程峰值細(xì)化與門限自適應(yīng)的實(shí)用技巧6.1 用拋物線插值把定時(shí)偏移做到分?jǐn)?shù)采樣級(jí)仿真驗(yàn)證的是整采樣點(diǎn)定時(shí)但工程上有時(shí)希望把起點(diǎn)精度細(xì)化到亞采樣級(jí)。Park和Minn的峰形接近拋物線可以做局部二次插值[~, k] max(M); if k 1 k length(M) num M(k-1) - M(k1); den 2 * (M(k-1) - 2*M(k) M(k1)); fracOffset -num / den; % 相對(duì)整數(shù)峰值位置的偏移 dFine k fracOffset; end代碼里的num和den來自三點(diǎn)二次插值公式fracOffset通常落在-0.5到0.5之間。需要注意如果系統(tǒng)沒有過采樣插值出來的小數(shù)采樣不代表真實(shí)采樣時(shí)鐘能對(duì)準(zhǔn)它更多是用于統(tǒng)計(jì)分析和定時(shí)誤差跟蹤。SC的平臺(tái)不能這么插值平臺(tái)內(nèi)插值得到的位置是假的必須先取平臺(tái)中點(diǎn)再考慮細(xì)調(diào)。6.2 門限自適應(yīng)用前導(dǎo)前空白段估計(jì)噪聲底在幀設(shè)計(jì)中前導(dǎo)前面留一段空白不出信號(hào)接收端就能估計(jì)噪聲功率然后用自適應(yīng)門限代替固定值noisePower mean(abs(rx(noiseWin)).^2); th 0.5 * max(M) 3 * sqrt(noisePower);這個(gè)門限有兩個(gè)分量一半來自當(dāng)前搜索窗的峰值另一半來自噪聲底。3是經(jīng)驗(yàn)系數(shù)相當(dāng)于抑制虛警的啟發(fā)式設(shè)置。低信噪比時(shí)噪聲項(xiàng)會(huì)把門限抬住防止噪聲尖峰被當(dāng)信號(hào)高信噪比時(shí)峰值項(xiàng)占主導(dǎo)門限接近最大值的一半照顧到了Minn副峰問題。我現(xiàn)在拿到一個(gè)同步算法第一件事永遠(yuǎn)是畫出M(d)曲線看峰形看它到底是平臺(tái)、尖峰還是帶毛刺的脈沖然后再談門限和插值。這個(gè)習(xí)慣來自在某跨平臺(tái)系統(tǒng)上被SC平臺(tái)坑了一整晚的教訓(xùn)——當(dāng)時(shí)EVM怎么調(diào)都下不來后來發(fā)現(xiàn)問題不在頻偏補(bǔ)償而在于定時(shí)點(diǎn)一直在平臺(tái)里隨機(jī)游走。希望這篇筆記能幫你少走幾步同樣的彎路。本文還有配套的精品資源點(diǎn)擊獲取