直角反射器CST仿真實(shí)例（2）- 距離分辨率，ICZT與逆傅里葉IFT

上期寫了用Field sources模式和ICZT算法來計(jì)算一維距離像，物體只是一個(gè)，所以很容易分辨出來。

這期我們用同樣的案例，物體是兩個(gè)，看看如何分辨距離，還有ICZT與IFT的對比。 

先通過transform復(fù)制直角反射器，使徑向的距離相差20mm：     

場源距離兩個(gè)物體是2m和2.02m。

A求解器還是計(jì)算的寬頻，頻域采樣間距Fstep=0.0125還是滿足奈奎斯特條件和我們Tmax=20ns的假設(shè)： 

  計(jì)算F21 (antenna coupling)，開始仿真：

同樣使用逆ChirpZ變換的后處理，我們將求解器計(jì)算的全部帶寬14GHz都用上，結(jié)果也只是一個(gè)波峰，分辨不出兩個(gè)物體：

這里橫坐標(biāo)是時(shí)域結(jié)果，我就不換成距離像了。可以簡單驗(yàn)證：峰值13.4ns是往返時(shí)間，所以測得的物體距離是13.4/2*1e-9*3e8=2.01m左右。

那么我們現(xiàn)在的距離分辨率是多少呢？dR=CLight/Bandwidth/2, 所以距離分辨率目前是3e8/14e9/2=0.01m，也就是說，2.01+/-0.01m的范圍內(nèi)有物體，所以分辨不出來2m和2.02m的兩個(gè)物體了（兩個(gè)波峰幾乎重合了）。

那么如何分辨出這0.02m的差別呢？當(dāng)然是增加帶寬，可我們這里用的的遠(yuǎn)場源當(dāng)初定義的時(shí)候就只包括70-84GHz的帶寬，所以增加帶寬是不行了。為了演示分辨兩個(gè)物體，我們可以將徑向間距擴(kuò)大到0.04m，這樣14GHz帶寬是可以區(qū)分兩個(gè)物體的：

重新開始仿真，自動(dòng)運(yùn)行后處理，用的還是14GHz帶寬和Tmax，這種情況我們就可以用兩個(gè)波峰分辨兩個(gè)距離了：

可以驗(yàn)證一下兩個(gè)時(shí)域波峰時(shí)差和0.04m距離差的關(guān)系：0.268/2*1e-9*3e8=0.04m，沒錯(cuò)。 

下面我們用逆傅里葉IFT驗(yàn)證一下為什么我們設(shè)置的Fstep=0.0125GHz在這里滿足奈奎斯特條件。上期提到過，假如我們的時(shí)域信號最大看到20ns，那么最小采樣密度df=1/2/20ns，最大頻率的最少采樣數(shù)為2*84G/df= 6720個(gè)，進(jìn)入傅里葉變換后處理：

對F21進(jìn)行IFT計(jì)算，最大時(shí)間同樣是20ns，選擇Nyquist重新采樣，如果我們用同樣的采樣數(shù)6720, 則計(jì)算順利：

這里IFT雖然也能看出兩個(gè)波峰，但功能上不如ICZT，這個(gè)結(jié)果也不能完全和ICZT相比，所以我們就不放在一起了，算法還是有差別的。

若少一點(diǎn)IFT采樣，則提示采樣不夠：

小結(jié)：

1）場源帶寬和求解器帶寬可以設(shè)置相同，這個(gè)是仿真項(xiàng)目能用到的最大帶寬；后處理進(jìn)行ICZT時(shí)，可用最大帶寬內(nèi)的任意帶寬。

2）距離分辨率只受帶寬影響：dR=CLight/Bandwidth/2。有的表達(dá)式也寫成Rres=c/BW之類的，若不除以2，可能是往返距離。

3）IFT只用來驗(yàn)證ICZT的采樣，不能替代ICZT計(jì)算距離像。ICZT模板有更多功能，比如設(shè)置帶寬、濾波器以及補(bǔ)零功能等等，ICZT配合A求解器中的Nyquist采樣，更適合用來手動(dòng)計(jì)算距離像。

4）A求解器中的頻域采樣數(shù)太多會減慢計(jì)算速度，收益也不大，所以要按Nyquist設(shè)置比較好。