作者：周盛 王曉春 計建軍 楊軍 王延群 單位：中國醫學科學院生物醫學工程研究所

1方法

1．1發射編碼

發射編碼按其內容分為相位編碼和頻率編碼，也可以按照發射次數分為單次發射碼和多次發射碼，如表1所示。Chirp碼脈沖壓縮效果受衰減和波束合成等非線性因素的影響相對較小，而且只需要一次發射就可以完成脈沖壓縮，受組織運動影響較小，因此是超聲成像系統編碼信號的理想選擇之一。但是由于Chirp碼需要多幅值發射電壓，對發射電壓和硬件的要求較高。因此在數字化高頻眼科超聲診斷設備中，基于發射編碼、系統復雜程度等考慮，一般采用二相編碼信號作為發射編碼。Golay互補序列雖然可以完全消除旁瓣達到最佳脈沖壓縮效果，但兩次發射的機制，容易受到組織運動的影響，引入誤差。雖然眼睛是基本靜止的生物組織，但由于我們所選用的是機械扇形掃描探頭，所以還是會在一定程度上影響成像效果。M序列在碼長較短時的脈沖壓縮效果較Barker碼差;又根據文獻［10－11］給出的最優編碼序列的概念、判斷標準和計算方法，通過窮舉搜索可以得到不同碼長的最優編碼序列［8］，Barker碼除了碼長N=11時不是最優編碼序列，其余都是。而且文獻［12］中也對13位Barker碼應用于高頻超聲成像中做了一定的實驗室研究。綜合比較，最終采用13位Barker碼作為本項目數字超聲眼科診斷系統中編碼激勵的發射碼型。

1．2解碼壓縮方式

目前醫學超聲編碼激勵系統中常用的脈沖壓縮方法分為匹配濾波和非匹配濾波兩種。對于Golay碼和Chirp碼，匹配濾波都能取得很好的脈沖壓縮效果，并且設計及實現方法相對簡單。而對于Barker碼，采用匹配濾波脈沖壓縮方法具有最低的距離旁瓣，但主瓣幅度也只是旁瓣幅度的N倍，其中N是Barker碼的長度。對于最長的13位Barker碼，匹配脈沖壓縮得到的信噪比增益為11．1dB，而峰值旁瓣水平(peaksidelobelevel，PSL)也只有(22dB，對距離旁瓣的抑制不夠好［9］。為此設計出了逆濾波器［13－14］、維納濾波器［15］、尖峰濾波器［16］以及失配濾波器等非匹配濾波脈沖壓縮方法［12，17］，以損失很小的信噪比增益為代價來換取更好的距離旁瓣脈沖壓縮效果。對比各種解碼壓縮方式的優缺點，本研究對匹配濾波與非匹配濾波分別進行了仿真實驗，結果表明。失配濾波方法可以達到最好的脈沖壓縮效果。通過此方法，可以使平均距離旁瓣水平(integratedsidelobelevel，ISL)和PSL都盡可能的低。

1．3仿真實驗

二進制13位Barker碼表示為:［+1，+1，+1，+1，+1，－1，－1，+1，+1，－1，+1，－1，+1］，由于采樣頻率一般選用信號回波頻率的8倍，這里用以下序列來構造每周期8個采樣點的正弦波，分別用［0，－7，－10，－7，0，7，10，7］表示－1，用［0，7，10，7，0，－7，－10，－7］表示+1，圖1～圖3為Barker碼的Matlab仿真波形。其中，圖1為13位Barker碼波形。圖2為經過匹配濾波后得到的結果，主瓣峰峰值達到了±312，而旁瓣峰峰值為±24，主副比為13。圖3為13位Barker碼匹配濾波結果的峰值包絡波形，可見旁瓣衰減達到了－45dB。13位Barker碼匹配濾波的FPGA實現相對簡單，需要開辟一個深度13×8，14bit的存儲空間進行卷積解碼。將解碼濾波系數序列中為+1的系數對應的回波數據保持原數，將解碼濾波系數序列中為(1的系數對應的回波數據進行補碼取反，然后進行對應數據并行累加，即可得到匹配濾波的結果。圖4是用Modelsim仿真得到的結果圖，上端波形為13位Barker碼編碼激勵信號，下端波形為經過匹配濾波解碼壓縮后的信號。可見，Barker碼匹配濾波解碼壓縮后存在較大的旁瓣。圖5為13位Barker碼經過失配濾波后得到的結果，實驗中我們設計的失配濾波器的階數為P=3N=39。主瓣峰峰值達到了±127，而旁瓣峰峰值為±1，主副比為127，是匹配濾波主副比的9．7倍。由此說明失配濾波雖然降低了主瓣峰值，但有效地降低了旁瓣水平，提高了信噪比。圖6為13位Barker碼失配濾波結果的峰值包絡波形，可見旁瓣衰減達到了－60dB。在FPGA的設計中，具體的實現步驟如圖7所示。首先將失配濾波器系數進行量化，乘以27=28，并取整存入FPGA的RAM中，然后構建實時延時回波矩陣，將系數分別與回波矩陣各行相乘，隨后將矩陣39行相乘的結果累加，并除以量化倍數，即可得到失配濾波解碼壓縮的結果。圖8是用Modelsim仿真的結果圖，上端波形為13位Barker碼編碼激勵信號，下端波形為經過失配濾波解碼壓縮后的信號。可見，13位Barker碼失配濾波解碼壓縮后有效地壓縮了旁瓣，在保證分辨率的前提下可以有效地提升信噪比。

1．4系統方案與設計

編碼激勵系統和傳統的脈沖回波成像系統不同之處在于:(1)發射電路采用編碼發射激勵，激勵脈沖的長度遠大于單脈沖;(2)接收電路需要對回波信號進行脈沖壓縮，脈沖壓縮算法的計算量很大。因此，編碼發射的實驗系統設計有以下特殊要求:(1)發射電路能夠產生各種不同的編碼激勵;(2)該系統成像幀頻是10幀/s，每幀由512條線組成，每條線包含768個點，要求FPGA能夠實時地對編碼激勵的回波進行脈沖壓縮并且成像。根據以上系統設計要求，圖9為系統的整體框架結構圖。主要包括上位機主機、硬件電路板、3個超聲探頭、LCD液晶顯示、腳閘、USB外設等。系統中的上位機選用嵌入式的工控主板，這樣可以實現儀器的小型化和便攜性。上位機主要用來實現對底層數據的接收和處理、以及軟件的實時控制。實驗系統采用單陣元機械扇形掃描探頭，3個超聲探頭分別為中心頻率為10MHz的B超探頭、中心頻率為20MHz的B超探頭和中心頻率為10MHz的A超探頭，其分別用來作為眼球全景成像、眼部組織高頻成像和眼球生物測量等模式的檢測探頭。LCD液晶顯示用來對圖像實時成像。雙鍵腳閘實現探頭的啟動、凍結。USB外設包括存儲器、打印機、鼠標和鍵盤。USB外設的使用，使得整機更加人性化和提高了可操作性。

硬件電路板是整個系統的核心部件，其結構如圖10所示。整個硬件電路以FPGA及其外圍電路為核心進行設計。FPGA采用Altera公司的CycloneIII系列EP3C55型產品。機械扇掃探頭的偏轉控制、圖像顯示的同步、USB接口邏輯等系統控制信號均由FPGA產生。FPGA產生的4組編碼脈沖序列通過MOSFET觸發電路將脈沖幅度由3.3V轉化為12V，以激勵發射電路，如圖11所示。發射電路由兩組雙極性脈沖發射芯片構成，并帶有歸零功能。編碼序列通過發射電路后，成為雙極性脈沖序列，激勵超聲換能器晶片。回波信號經高壓隔離電路后進入前置放大與可變增益控制電路。前置放大電路增益范圍為10～15dB，FPGA產生TGC曲線數據，控制可變增益放大器，可以使增益變化50dB左右。實驗系統時鐘為80MHz，高速A/D芯片按照不同模式工作在80MHz和120MHz，采樣精度為14bits，其采樣模式由FPGA通過SPI接口完成設定。FPGA完成的超聲回波數字信號處理包括:高頻數據采集與存儲、數字濾波、數字檢波、對數變換、二次采樣等，如圖12所示。數字信號處理之后的回波數據通過SRAM緩存后，經USB2．0接口電路上傳至計算機，并由人機界面完成整個系統的操作與控制。