畢業(yè)設(shè)計(jì)---基于matlab的dpll設(shè)計(jì)與仿真

1、<p><b>　　摘要</b></p><p>　　現(xiàn)代通信系統(tǒng)中，同步問(wèn)題是決定系統(tǒng)性能和應(yīng)用的根本問(wèn)題。鎖相環(huán)路作為同步技術(shù)的核心部件,已在模擬和數(shù)字通信及無(wú)線電電子學(xué)等各個(gè)領(lǐng)域中得到了極為廣泛的應(yīng)用，特別是在數(shù)字通信的調(diào)制解調(diào)和位同步中常常要用到各種各樣的鎖相環(huán)。</p><p>　　本文分析了全數(shù)字鎖相環(huán)(NR—DPLL) 基本組成和工作原理，利用

2、MATLAB對(duì)鎖相環(huán)的進(jìn)行設(shè)計(jì)，了解其各個(gè)部分功能具體結(jié)構(gòu)，利用MATLAB工具箱對(duì)鎖相環(huán)的進(jìn)行仿真，改變鎖相環(huán)參數(shù)最后進(jìn)行了參數(shù)分析。同時(shí)重點(diǎn)針對(duì)DDS做詳細(xì)了解，對(duì)其做出仿真并進(jìn)行分析</p><p>　　關(guān)鍵詞：奈奎斯特型全數(shù)字鎖相環(huán)；模塊設(shè)計(jì)；仿真；參數(shù)分析</p><p><b>　　Abstract</b></p><p>　　Mo

3、dern communication systems, the synchronization problem is to determine system performance and application of the fundamental problems. PLL synchronization technology as the core components in analog and digital communic

4、ation and radio electronics and other fields has been very widely used, especially in the digital communications modem and bit synchronization often use a variety of all kinds of PLL.</p><p>　　This paper ana

5、lyzes the all-digital phase-locked loop (NR-DPLL) basic composition and working principle of the PLL using MATLAB to design, understand the function of each part of the specific structure of the PLL using MATLAB simulati

6、on toolbox change the PLL parameters Finally, a parametric analysis. At the same time focusing on a detailed understanding of DDS, making simulation and analysis of its</p><p>　　Keywords: Nyquist-type all-d

7、igital phase locked loop, module design, Simulation, Parameter</p><p><b>　　目 錄</b></p><p><b>　　摘要I</b></p><p>　　AbstractII</p><p><b>　　目

8、 錄III</b></p><p><b>　　第一章 緒論2</b></p><p>　　1.1 選題的背景與意義2</p><p>　　1.2 數(shù)字鎖相環(huán)的分類2</p><p>　　1.3 論文的主要研究工作3</p><p>　　第二章 鎖相環(huán)基本原理4<

9、;/p><p>　　2.1 鑒相器(PD)5</p><p>　　2.2 環(huán)路濾波器(LPF)6</p><p>　　2.3 壓控振蕩器(VCO)6</p><p>　　第三章 數(shù)字鎖相環(huán)的經(jīng)典結(jié)構(gòu)設(shè)計(jì)7</p><p>　　3.1 奈奎斯特采樣鑒相器7</p><p>　　3.

10、2 數(shù)字環(huán)路濾波器7</p><p>　　3.3 數(shù)字控制振蕩器（NCO）8</p><p>　　3.4 NR-DPLL的動(dòng)態(tài)方程與相位模型13</p><p>　　3.5 NR-DPLL的性能分析14</p><p>　　第四章 奈奎斯特?cái)?shù)字鎖相環(huán)（NR-DPLL)的仿真16</p><p>　　4

11、.1 仿真原理與過(guò)程16</p><p>　　4.2 仿真結(jié)論22</p><p>　　第五章 論文以后的研究工作24</p><p>　　5.1 異或鑒相器24</p><p>　　5.2 一階環(huán)路濾波器25</p><p>　　5.3 壓控振蕩器25</p><p>　

12、　第六章 結(jié)束語(yǔ)26</p><p><b>　　參考文獻(xiàn)27</b></p><p><b>　　謝辭28</b></p><p>　　附錄一 MATLAB簡(jiǎn)介29</p><p>　　附錄二 H( )的幅頻響應(yīng)30</p><p>　　附錄三 誤差電壓

13、和控制電壓 31</p><p><b>　　第一章 緒論</b></p><p>　　1.1 選題的背景與意義</p><p>　　鎖相環(huán)路已在模擬和數(shù)字通信及無(wú)線電子電子學(xué)的各個(gè)領(lǐng)域中得到了幾位廣泛的應(yīng)用。隨著大規(guī)模，超高速數(shù)字2集成電路的發(fā)展以及計(jì)算機(jī)的普遍應(yīng)用，在傳統(tǒng)的模擬鎖相環(huán)路（APLL）應(yīng)用領(lǐng)域中，一大部分已被數(shù)字鎖相環(huán)路（D

14、PLL）所取代。全數(shù)字鎖相環(huán)ADPLL(All Digital Phase-Locked Loop)，顧名思義，其環(huán)路中的所有部件都是用數(shù)字電路來(lái)實(shí)現(xiàn)的，到20世紀(jì)70年代開始出現(xiàn)的數(shù)字鎖相環(huán)，在現(xiàn)代科技的很多領(lǐng)域都有鎖相環(huán)的應(yīng)用。數(shù)字鎖相技術(shù)在數(shù)字通信的調(diào)制解調(diào)、位同步、頻率合成中常常要用到各種各樣的鎖相環(huán)，本文主要討論的奈奎斯特全數(shù)字鎖相環(huán)模塊設(shè)計(jì)與仿真，對(duì)于加深對(duì)數(shù)字鎖相環(huán)的理解，以及對(duì)其他種類數(shù)字鎖相環(huán)的理解應(yīng)用都有知道很好的指

15、導(dǎo)作用。</p><p>　　1.2 數(shù)字鎖相環(huán)的分類</p><p>　　所謂全數(shù)字鎖相環(huán)，就是全數(shù)字化的鎖相環(huán)路，即各個(gè)環(huán)路各個(gè)部件全部數(shù)字化，采用數(shù)字鑒相器，數(shù)字環(huán)路濾波器，數(shù)字壓控振蕩器構(gòu)成的鎖相環(huán)路。按照數(shù)字鑒相器的的形式把數(shù)字鎖相環(huán)分成四類，分別是：</p><p><b>　　過(guò)零型數(shù)字鎖相環(huán)路</b></p>&

16、lt;p>　　這種數(shù)字鎖相環(huán)路采用過(guò)零采樣數(shù)字鑒相器，即本地估算信號(hào)在輸入信號(hào)的過(guò)零點(diǎn)上采樣進(jìn)行A/D變換，得到數(shù)字相位誤差信號(hào)輸出</p><p>　　觸發(fā)器型數(shù)字鎖相環(huán)路</p><p>　　這類全數(shù)字鎖相環(huán)路使用的數(shù)字鑒相器是觸發(fā)器型數(shù)字鑒相器。其特點(diǎn)是利用輸入信號(hào)和本地估算信號(hào)的正向或者負(fù)向過(guò)零點(diǎn)對(duì)觸發(fā)器進(jìn)行觸發(fā)，在觸發(fā)器的置0或者置1的時(shí)間間隔內(nèi)得到相位誤差信號(hào)。置0或置1

17、的時(shí)間間隔寬度就表征了輸入信號(hào)和本地估算信號(hào)的指尖的相位誤差大小。 </p><p>　　導(dǎo)前-滯后型數(shù)字鎖相環(huán)路</p><p>　　這種數(shù)字鎖相環(huán)路采用的鑒相器是導(dǎo)前-滯后型數(shù)字鑒相器，導(dǎo)前-滯后型數(shù)字鑒相器在每一個(gè)周期內(nèi)得到輸入信號(hào)的相位比本地估算信號(hào)相位超前或滯后的信息，因此，這種鑒相器的相位誤差輸出只有道歉和滯后兩種狀態(tài)。</p><p>　　奈奎斯特速率

18、采樣型數(shù)字鎖相環(huán)路</p><p>　　在這種數(shù)字鎖相環(huán)路中，對(duì)輸入信號(hào)的采樣按奈奎斯特速率進(jìn)行。對(duì)于輸入信號(hào)進(jìn)行A/D變換的采樣速率必須按奈奎斯特速率進(jìn)行，以使輸入信號(hào)能夠按奈奎斯特取樣定理再現(xiàn)。A/D變換后的輸入信號(hào)與本地估算信號(hào)進(jìn)行數(shù)字相乘，得到需要的相位誤差信號(hào)，以完成鑒相功能。只有這一類為均勻采樣DPLL，前三中均不是。</p><p>　　1.3 論文的主要研究工作<

19、/p><p>　　（1）本次研究主要是奈奎斯特?cái)?shù)字鎖相環(huán)，前面討論的三種數(shù)字鑒相器，即過(guò)零型數(shù)字鎖相環(huán)路、觸發(fā)器型數(shù)字鎖相環(huán)路、導(dǎo)前-滯后型數(shù)字鎖相環(huán)路，它們都是以本地估算信號(hào)為基準(zhǔn)相位對(duì)輸入信號(hào)進(jìn)行比較，得到相位誤差信息。這種采樣方式通常稱為非均勻采樣。而奈奎斯特速率采樣數(shù)字鑒相器的采樣方式不同，它是以某一固定的頻率基準(zhǔn)作為輸入信號(hào)的采樣脈沖。這時(shí)，采樣頻率必須足夠高，以使采樣后的采樣后的輸入信號(hào)能夠依據(jù)奈奎斯特采

20、樣定理再現(xiàn)輸入信號(hào)。 </p><p>　　（2）闡述奈奎斯特?cái)?shù)字鎖相環(huán)的典型結(jié)構(gòu)與相位模型，并分析一階和二階的跟蹤性能。</p><p>　　（3）對(duì)DDS利用Simulink進(jìn)行仿真，并對(duì)其性能進(jìn)行分析</p><p>　　第二章 鎖相環(huán)基本原理</p><p>　　先了解鎖相環(huán)的基本概念，所謂鎖相，就是相位同步

21、的自動(dòng)控制。完成兩個(gè)信號(hào)間相位同步的自動(dòng)控制系統(tǒng)的環(huán)路叫做鎖相環(huán)，也稱PLL(Phase Locked Loop)。最典型的鎖相環(huán)由鑒相器(Phase Detector)、環(huán)路濾波器(Loop Filter)、壓控振蕩器(Voltage Controlled Oscillator)三部分組成，如圖所示：</p><p>　　圖2.1 鎖相環(huán)典型結(jié)構(gòu)</p><p>　　當(dāng)壓控振蕩器的頻率

22、由于某種原因而發(fā)生變化時(shí)，必然引起相位的變化，該相位變化在鑒相器中與參考晶體的穩(wěn)定相位(對(duì)應(yīng)于頻率 )相比較，使鑒相器輸出一個(gè)與相位誤差信號(hào)成比例的誤差電壓 (t)，經(jīng)過(guò)低通濾波器，取出其中緩慢變動(dòng)數(shù)值，將壓控振蕩器的輸出頻率拉回到穩(wěn)定的值上來(lái)，從而實(shí)現(xiàn)了相位負(fù)反饋控制。</p><p>　　當(dāng) 與 相等時(shí)，兩矢量以相同的角速度旋轉(zhuǎn)，相對(duì)位置固定，即夾角維持不變，通常數(shù)值又較小，這就是環(huán)路的鎖定狀態(tài)。<

23、/p><p>　　從輸入信號(hào)加到鎖相環(huán)路的輸入端開始，一直到環(huán)路達(dá)到鎖定的全過(guò)程，稱為捕獲過(guò)程。設(shè)系統(tǒng)最初進(jìn)入同步狀態(tài) 的時(shí)間為 ，那么從 的起始狀態(tài)到達(dá)進(jìn)入同步狀態(tài)的全部過(guò)程就稱為鎖相環(huán)路的捕獲過(guò)程。捕獲過(guò)程所需的時(shí)間 。稱為捕獲時(shí)間。顯然，捕獲時(shí)間幾的大小不但與環(huán)路的參數(shù)有關(guān)，而且與起始狀態(tài)有關(guān)。</p><p>　　對(duì)一定的環(huán)路來(lái)說(shuō)，是否能通過(guò)捕獲進(jìn)入同步取決于起始頻差 。若 超過(guò)某一范

24、圍，環(huán)路就不能捕獲了。這個(gè)范圍的大小是鎖相環(huán)路的一個(gè)重要性能指標(biāo)，稱為環(huán)路的捕獲帶 。</p><p>　　捕獲狀態(tài)終了，環(huán)路的狀態(tài)穩(wěn)定在 </p><p>　　這就是同步狀態(tài)的定義。只要在整個(gè)變化過(guò)程中一直滿足以上兩式，那么仍稱環(huán)路處于同步狀態(tài)。由上可知，在輸入固定頻率信號(hào)的條件之下，環(huán)路進(jìn)入同步狀態(tài)后，輸出信號(hào)與輸入信號(hào)之間頻差等于零，相差等于常數(shù)，即</p><p

25、><b>　　, </b></p><p>　　這種狀態(tài)就稱為鎖定狀態(tài)。</p><p>　　實(shí)際應(yīng)用中有各種形式的環(huán)路，但它們都是由圖1這個(gè)基本環(huán)路演變而來(lái)的。下面逐個(gè)介紹基本部件在環(huán)路中的作用。</p><p>　　2.1 鑒相器(PD)</p><p>　　鑒相器是一個(gè)相位比較裝置，用來(lái)檢測(cè)輸入信號(hào)相位與反

26、饋信號(hào)相位之間的相位差。輸出的誤差信號(hào)是相差的函數(shù)，即鑒相特性可以是多種多樣的，有正弦形特性、三角形特性、鋸齒形特性等等。常用的正弦鑒相器可用模擬相乘器與低通濾波器的串接作為模型。如圖2.2所示。</p><p>　　圖2.2 正弦鑒相器模型</p><p>　　設(shè)乘法器的相乘系數(shù)為Km，單位為 ，則</p><p><b>　　= </b>&

27、lt;/p><p><b>　　= </b></p><p>　　在經(jīng)過(guò)低通濾波器(LPF)慮除2 。成分后，得到誤差電壓</p><p>　　令 為鑒相器輸出的最大電壓。則正弦鑒相器特性為</p><p>　　上述鑒相器的功能可分解為兩個(gè)作用，首先是相位相減，即取得兩個(gè)輸入信號(hào)之間的相位差;其次是將相位差轉(zhuǎn)換為誤差電壓輸出

28、，所以它是相位轉(zhuǎn)換為電壓的裝置。</p><p>　　隨著半導(dǎo)體技術(shù)的不斷發(fā)展，近十年出現(xiàn)的數(shù)字鎖相環(huán)(DPLL)的鑒頻鑒相器核心是由鑒相器和電荷泵(CP)構(gòu)成的。這種結(jié)構(gòu)有利于實(shí)現(xiàn)數(shù)字統(tǒng)一控改善鎖相環(huán)的捕捉帶。</p><p>　　2.2 環(huán)路濾波器(LPF)</p><p>　　環(huán)路濾波器是線性電路，由線性元件電阻、電容或運(yùn)算放大器組成。它的輸入是鑒相器的輸出

29、電壓 (t)，具有低通特性，它可以起到低通濾波器的作用，慮除誤差電壓 (t)的高頻分量，取出平均電壓 (t)去控制壓控振蕩器。它可以改善控制電壓的頻譜純度，提高系統(tǒng)穩(wěn)定度，更重要的是它對(duì)環(huán)路參數(shù)的調(diào)整起到了決定性的作用，直接影響到輸出信號(hào)的穩(wěn)定、頻譜純度、鎖定時(shí)間等。</p><p>　　表示環(huán)路濾波的輸出輸入電壓關(guān)系的是濾波器的傳遞函數(shù)Z(S),表示為 (s)=Z(s)× (s

30、)</p><p>　　式中s是復(fù)頻率。當(dāng)將s= 代入上式時(shí)，就得到了濾波器的頻率響應(yīng)特性。</p><p>　　比較常用的環(huán)路濾波器有簡(jiǎn)單RC濾波器、無(wú)源比例積分濾波器、有源比例積分濾波器，其它濾波器都是在上述三種基礎(chǔ)上的改進(jìn)。</p><p>　　環(huán)路濾波器的主要指標(biāo)是帶寬、直流增益、高頻增益，它由濾波器的時(shí)間常數(shù)和濾波器的類型決定。比例積分濾波器在高頻時(shí)有一定

31、增益，這對(duì)鎖相環(huán)的捕捉特性有利，而且比例積分引入了一個(gè)零點(diǎn)，有利于增加環(huán)路的穩(wěn)定性。有源濾波器還要考慮它的線性動(dòng)態(tài)范圍，要求輸出的電壓能夠供給VCO在鎖定頻率所需的控制電壓。</p><p>　　2.3 壓控振蕩器(VCO)</p><p>　　壓控振蕩器是一個(gè)電壓-頻率變換裝置，在環(huán)路中作為被控振蕩器，它的振蕩頻率應(yīng)隨輸入控制電壓 (t)線性地變化。VCO理想的頻率受控特性應(yīng)為線性的，

32、即</p><p>　　當(dāng)控制電壓 (t)=0時(shí)，VCO的輸出為自由振蕩頻率 ， 為壓控靈敏度。VCO反饋到鑒相器上，由于鑒相器輸出誤差電壓起作用的是其相位 ，而相位是頻率的積分，即</p><p><b>　　= </b></p><p><b>　　其中 </b></p><p>　　實(shí)際應(yīng)用中

33、的壓控振蕩器的控制特性只有有限的線性控制范圍，超出這個(gè)范圍之后控制靈敏度將會(huì)下降。</p><p>　　第三章 數(shù)字鎖相環(huán)的經(jīng)典結(jié)構(gòu)設(shè)計(jì)</p><p>　　數(shù)字鎖相環(huán)全部部件都是數(shù)字電路，所以信號(hào)都是二進(jìn)制的或多進(jìn)制的數(shù)字信號(hào)。數(shù)字鎖相環(huán)常用來(lái)提取位同步信號(hào)。</p><p>　　3.1 奈奎斯特采樣鑒相器</p><p>　　數(shù)字鑒相

34、器也稱采樣鑒相器，用于檢測(cè)出兩個(gè)輸入信號(hào)的相位差，它在很大程度上決定著鎖相環(huán)的性能，本文重點(diǎn)講解奈奎斯特采樣鑒相器，框圖如圖所示：</p><p>　　圖3.1 奈奎斯特采樣鑒相器</p><p>　　= 3.1-1</p><p><b>　　輸入信號(hào)</b></p><p>　　=

35、 3.1-2</p><p><b>　　其中</b></p><p>　　= , = 3.1-3</p><p>　　輸入信號(hào)經(jīng)A/D采樣后，第k個(gè)采樣時(shí)刻量化后的數(shù)字信號(hào)為</p><p>　　= 3.1-

36、4</p><p>　　對(duì)輸入信號(hào)進(jìn)行A/D變換的采樣速率有帶通信號(hào)奈奎斯特采樣定理確定，但為防止信號(hào)頻譜混疊并保證信號(hào)相位信息的有效抽取，采樣速率一般選取前置帶通濾波器的兩倍帶寬以上。</p><p>　　令 = (k), = (k),即 (k)= (k)相乘后，經(jīng)低通濾波器得到數(shù)字誤差信號(hào)</p><p>　　=

37、 3.1-5</p><p>　　- 3.1-6</p><p>　　3.2 數(shù)字環(huán)路濾波器</p><p>　　數(shù)字環(huán)中使用的數(shù)字環(huán)路濾波器與模擬環(huán)中使用的環(huán)路濾波器作用一樣，都對(duì)噪聲及高頻分量起抑制作用，并且控制著環(huán)路相位校正的速度與精度。適當(dāng)選擇濾波器參數(shù)，可以改善環(huán)路的性能。一

38、般構(gòu)成形式由A／D、數(shù)字計(jì)數(shù)器和D／A三部分組成。</p><p>　　數(shù)字濾波器在實(shí)際應(yīng)用中一階數(shù)字環(huán)路濾波器的實(shí)現(xiàn)形式如圖：</p><p>　　圖3.2 一階數(shù)字環(huán)路濾波器</p><p><b>　　其Z域傳遞函數(shù)</b></p><p>　　F(z)= = +

39、 3.2-1</p><p>　　按照?qǐng)D3.2所實(shí)現(xiàn)的數(shù)字濾波器，其頻率特性與式F( )= 所示理想積分濾波器的頻率特性一致，兩種濾波器參數(shù)之間也有著一定的對(duì)應(yīng)關(guān)系。對(duì)理想積分濾波器的傳遞函數(shù)式F(s)= 采用雙線性變換,即令s= ，代入式F(s)= ，得到</p><p>　　F（z）= + = + </p><p>　　式中， 為采樣周期，</p>

40、<p>　　= 3.2-2</p><p>　　= 3.2-3</p><p>　　式3.2-2和3.2-3表明了兩種濾波器參數(shù)之間的對(duì)應(yīng)關(guān)系。</p><p>　　二階數(shù)字環(huán)路濾波器可由兩個(gè)一階數(shù)字濾波器串聯(lián)得到

41、。</p><p>　　3.3 數(shù)字控制振蕩器（NCO）</p><p><b>　　(一)基本工作原理</b></p><p>　　NCO采用直接數(shù)字相位綜合技術(shù)(DDS）,該技術(shù)主要是由時(shí)鐘驅(qū)動(dòng)讀取三角函數(shù)表，功能框圖如圖3.3所示</p><p><b>　　(k) </b></p&g

42、t;<p><b>　　取樣時(shí)鐘</b></p><p>　　圖3.3 DDS功能框圖</p><p>　　DDS是根據(jù)正弦函數(shù)的產(chǎn)生，從相位出發(fā)，由不同的相位給出的電壓幅度，即相位—正弦幅度變換，最后濾波，拼花輸出所需要的頻率。為了更好地理解DDS，有必要簡(jiǎn)單地回顧正弦函數(shù)的產(chǎn)生，在此基礎(chǔ)上建立DDS的基本概念。</p><p&g

43、t;　　圖3.4 表示了半徑R長(zhǎng)度為1的單位圓。R以原點(diǎn)為中心旋轉(zhuǎn)與 軸的正方向形成夾角 ,即相位角。設(shè)R在 軸上的投影為 。當(dāng)R的斷點(diǎn)連續(xù)不斷的繞圓旋轉(zhuǎn)時(shí)，S將取＋1～－1之間的任何值，而 則從 變化。</p><p>　　圖3.4 單位圓表示的正弦函數(shù) </p><p>　　S的長(zhǎng)度就是正弦函數(shù)的幅度，也就是 ，如圖3.5所示。如果R的斷點(diǎn)不是連續(xù)不斷地繞圓旋轉(zhuǎn)，

44、而是以等步長(zhǎng)的相位增量節(jié)約是旋轉(zhuǎn)，比如7步旋轉(zhuǎn)一周，那么S值將形成階梯式的近似正弦函數(shù)，如下圖所示是表示了一個(gè)64步長(zhǎng)的近似函數(shù)。</p><p>　　圖3.5 64步長(zhǎng)的階梯狀的正弦函數(shù)</p><p>　　從這個(gè)簡(jiǎn)單的討論過(guò)程可以看出，通過(guò)改變相位增量的步長(zhǎng)大小，就可以改變繞圓旋轉(zhuǎn)的步長(zhǎng)數(shù)目，從而產(chǎn)生變換的方法。步長(zhǎng)越小，旋轉(zhuǎn)一周就越快，因此近似正弦波的頻率就越高。相反，當(dāng)步長(zhǎng)增加時(shí)

45、，頻率就下降，最后變成方波輸出。</p><p>　　典型的DDS原理方框圖如圖3.6 所示。主要包括以下部件：相位累加器；相位—幅度轉(zhuǎn)換器，即正弦查找表ROM；D/A變換器；低通濾波器。</p><p>　　圖3.6 典型DDS原理方框圖</p><p>　　相位累加器相當(dāng)于一個(gè)簡(jiǎn)單的計(jì)數(shù)器，它是由L位存儲(chǔ)數(shù)字相位增量字的頻率寄存器、后接L位全加器和相位寄存器組成

46、，后二者常常是合并在一起的。數(shù)字輸入相位增量字進(jìn)入頻率寄存器后，在每個(gè)參考時(shí)鐘周期或者脈沖期間，表示相位增量的數(shù)值就加到容量正比于輸出信號(hào)相位的數(shù)字累加器中。在同樣的時(shí)間內(nèi)，累加器的數(shù)值傳遞給正弦查找表ROM，然后ROM輸出給D/A變化器。D/A變換器產(chǎn)生一系列的表示以時(shí)間脈沖速率抽樣的電壓階躍，最后經(jīng)低通濾波器平滑輸出。當(dāng)相位累加器由于重復(fù)相加而溢出時(shí)，它的最高有效位(MSB)就從1變成0，又開始一個(gè)新的周期。只要檢測(cè)到MSB有輸出，

47、DDS就產(chǎn)生一個(gè)改變符號(hào)的方波，改變的符號(hào)與抽樣正弦波交叉時(shí)的符號(hào)相一致。因?yàn)榉讲ㄖ辉趨⒖紩r(shí)鐘周期時(shí)改變符號(hào)，所以它的相位一般相對(duì)于被產(chǎn)生信號(hào)的理想相位要延遲一個(gè)變化量，相位延遲正比于累加器溢出好后剩余的數(shù)值。對(duì)于L位累加器，使用模 溢出的性質(zhì)，那么，在 個(gè)參考時(shí)鐘周期后，累加器中存儲(chǔ)的余數(shù)就是零。</p><p>　　下面建立DDS輸出頻率與其他一些參數(shù)之間的基本關(guān)系：

48、

49、 ——參考時(shí)鐘頻率， ； 3.3-1</p><p>　　——輸出頻率， 3.3-2</p><p>　　——相位增量字(非定值)，由二進(jìn)制碼表示的頻率建立字。</p><p>　　設(shè)累加器的長(zhǎng)度是L位。通常，實(shí)際中使用W位(L位中的最高W位，W<L,即相位截?cái)?

50、。其中 就是最小的相位增量。顯然，F(xiàn)SW對(duì)應(yīng)的相位增量是 。這樣，完成一整周期的正弦波輸出需要經(jīng)過(guò) 個(gè)參考時(shí)鐘周期。</p><p>　　或 3.3-3 其中 </p><p>　　相位累加器的輸出就是被合成信號(hào)線性增加的相位值。用合成信號(hào)頻率 和參考時(shí)鐘頻率 表示的<

51、/p><p><b>　　3.3-4</b></p><p><b>　　其基本結(jié)構(gòu)如下：</b></p><p>　　圖3.7 DDS結(jié)構(gòu)</p><p>　　取累加器長(zhǎng)度為32位，正弦ROM為256*sin( [0:2*pi/(2^10):2*pi] )+256，當(dāng)控制字為134217728是其輸出

52、波形為如圖</p><p>　　圖3.8 32點(diǎn)正弦波</p><p>　　每周期為32點(diǎn)，功率譜為</p><p>　　圖3.9 功率譜密度</p><p><b>　　NCO的數(shù)字模型</b></p><p>　　在相位累加器字長(zhǎng)N,采樣頻率為 確定的情況下，設(shè)NCO自由振蕩頻率為 ， (0)

53、=0 并由式（3-1-7）和式（3-1-6)確定頻率控制字初始值為 </p><p><b>　　及初始值為 </b></p><p>　　當(dāng)數(shù)字環(huán)路濾波器輸出數(shù)字控制電壓為 (k)時(shí)，加注于NCO控制端，相應(yīng)的頻率控制字變化量為 = (k)，NCO輸出頻率為</p><p><b>　　輸出相位為</b><

54、/p><p>　　其中 ,定義 rad/(sV))為NCO的頻率控制增益， rad/V)為相位控制增益， </p><p>　　NCO有相位累加器的功能</p><p>　　設(shè)初始狀態(tài) (0)=0時(shí)，利用Z變換，有 (k+1)=z (k)，則NCO輸出相位與控制電壓的關(guān)系為</p><p><b>　　（3.3-5）</b&

55、gt;</p><p>　　即NCO的數(shù)學(xué)模型。</p><p>　　3.4 NR-DPLL的動(dòng)態(tài)方程與相位模型</p><p>　　根據(jù)奈奎斯特?cái)?shù)字鎖相環(huán)的結(jié)構(gòu)，綜合分析數(shù)字誤差信號(hào)，環(huán)路濾波器Z域傳遞函數(shù)式以及NCO數(shù)學(xué)模型可得到ND-DPLL的動(dòng)態(tài)方程為</p><p>　　(k)= F(z) </p><p>

56、;<b>　　= </b></p><p>　　其中 (rad)，是數(shù)字環(huán)路增益。</p><p>　　NR-DPLL的Z域相位模型如圖3-7。</p><p>　　圖3-7 NR-DPLL的z域相位模型</p><p>　　3.5 NR-DPLL的性能分析</p><p>　?。ㄒ唬?一

57、階環(huán)的跟蹤性能</p><p>　　當(dāng)F(z)=1時(shí)，得到一階環(huán)的動(dòng)態(tài)方程</p><p><b>　　(3-5-1)</b></p><p>　　式（3-3-1）可等效為 </p><p><b>　?。?-5-2）</b></p><p><b&g

58、t;　　1。 </b></p><p>　　此時(shí) ，k>=0，再由式（3-5-2）可得</p><p><b>　　(3-5-3)</b></p><p>　　設(shè) ，當(dāng)k趨于無(wú)窮大時(shí)環(huán)路進(jìn)入同步狀態(tài)時(shí)，穩(wěn)態(tài)相位誤差 為一固定值 ，再由式（3-5-3）得 ，即一階數(shù)字環(huán)可無(wú)凈差的跟蹤相位階躍信號(hào)。</p>&l

59、t;p>　　2. =? ×k 得 </p><p><b>　?。?-5-5)</b></p><p>　　當(dāng)K趨于無(wú)窮大時(shí)環(huán)路進(jìn)入同步狀態(tài)時(shí)，希望有 (k+1)= ，則</p><p><b>　?。?-5-6）</b></p><p>　　令 ，得到一階數(shù)字環(huán)的同步

60、帶</p><p>　　超出環(huán)路同步帶，環(huán)路將不能跟蹤。</p><p>　　（二）二階數(shù)字環(huán)路跟蹤性能</p><p>　　在二階數(shù)字環(huán)路中，數(shù)字環(huán)路濾波器的域傳遞函數(shù)即</p><p>　　F(z)=G1+G2 (3-5-7)</p><p>　　其中A=G1=

61、G2,a= .</p><p>　　環(huán)路在同步狀態(tài)下，相位誤差較小時(shí)，可將鑒相器線性化，令</p><p>　　即將環(huán)路動(dòng)態(tài)方程式線性化，可得到線性化動(dòng)態(tài)方程</p><p><b>　　（3-5-8）</b></p><p>　　將式（3-5-7）代入式（3-5-8），得到</p><p>&l

62、t;b>　?。?-5-9）</b></p><p><b>　　1. </b></p><p>　　此時(shí) ，則式（3-5-9）中</p><p>　　2. 式（3-5-9）中</p><p>　　上述表明了，在輸入信號(hào)為相位階躍和頻率階躍時(shí)，環(huán)路式（3-5-7）的左邊都為0，可整理為</p>

63、;<p><b>　　（3-5-10）</b></p><p>　　令 , ,則式（3-5-10）可簡(jiǎn)寫為</p><p>　　若希望無(wú)靜差跟蹤頻率階躍，則需滿足</p><p><b>　　(3-5-11)</b></p><p>　　此條件下，環(huán)路對(duì)輸入相位階躍和頻率階躍信號(hào)都

64、可以無(wú)靜差地跟蹤。 有上述可見，在輸入頻率階躍信號(hào)時(shí)，環(huán)路方程與頻率階躍量無(wú)關(guān)，但其環(huán)路捕獲帶仍與參數(shù)有關(guān)。由此可見在數(shù)字環(huán)路中，環(huán)路參數(shù)的設(shè)計(jì)是十分重要的。</p><p>　　奈奎斯特?cái)?shù)字鎖相環(huán)（NR-DPLL)的仿真</p><p>　　4.1 仿真原理與過(guò)程</p><p>　　Z域相位模型如圖4-1</p><p>　　圖

65、4-1 NR-DPLL線性化相位模型</p><p>　　由圖4-1可得環(huán)路Z域傳遞函數(shù)為</p><p>　　H (z) = （4-1-1）</p><p>　　選擇式（3-1-5）描述的數(shù)字環(huán)路濾波器，代入式（4-1-1）得</p><p>　　H (z) = (4-1-2)</

66、p><p>　　式（4-1-2）對(duì)應(yīng)的模擬環(huán)路是理想二階環(huán)。模擬理想二階環(huán)的捕獲帶與環(huán)路參數(shù)無(wú)關(guān)，主要受控制電壓的范圍和壓控振蕩器的頻率變化范圍限制；而上述數(shù)字環(huán)路的捕獲帶不僅取決于環(huán)路參數(shù)還與采樣速率密切相關(guān)。當(dāng)環(huán)路參數(shù)設(shè)計(jì)合適時(shí)，可得到最大的捕獲帶約為采樣速率一半</p><p>　　設(shè)計(jì)一個(gè)捕獲帶 可取采樣速率 =6kHz，設(shè)計(jì)方法如下：</p><p>　　先設(shè)

67、計(jì)相應(yīng)的模擬理想二階環(huán)的參數(shù)，為得到好的時(shí)域響應(yīng)，一般取 =0.707且取 =10 。</p><p>　　令NCO的相位控制增益 ，鑒相靈敏度 ，得NCO的頻率控制增益 = ，則數(shù)字環(huán)路增益K= = 再由</p><p><b>　　Z= , </b></p><p><b>　　得到</b></p&

68、gt;<p><b>　　=0.707</b></p><p><b>　　解得 </b></p><p><b>　　， </b></p><p>　　時(shí)，環(huán)路的3dB帶寬 </p><p>　?。?）確定數(shù)字環(huán)路濾波器的參數(shù)</p><p

69、><b>　　G2= </b></p><p><b>　　G1= </b></p><p>　　所設(shè)計(jì)的環(huán)路參數(shù)滿足式（3-5-9)</p><p>　　(3)將上述設(shè)計(jì)的參數(shù)代入式（4-1-2）得到所設(shè)計(jì)的數(shù)字環(huán)路z域傳遞函數(shù)</p><p>　　H(z)=

70、 (4-1-3)</p><p>　　將Z= 代入式（4-1-3）得到數(shù)字環(huán)路的頻率響應(yīng)</p><p>　　圖4-2是H( )的幅頻響應(yīng)曲線。</p><p><b>　　圖4-2</b></p><p>　　由圖可見，-3dB截止頻率點(diǎn)在130Hz附近，與所設(shè)計(jì)的環(huán)路帶寬值相近，仿真結(jié)果符合設(shè)計(jì)結(jié)果。</p

71、><p>　　圖4-3是用所設(shè)計(jì)的數(shù)字環(huán)路，在參數(shù)G1=0.045,G2=0.0025,，輸入載波頻率偏移 時(shí)，環(huán)路入所過(guò)程中，誤差電壓 和控制電壓 的波形。此波形清楚地表明，環(huán)路可以捕獲入鎖且穩(wěn)態(tài)相位誤差為0。</p><p>　　令 G1=g1，G2=g2其中，g1g2數(shù)值。</p><p>　　用所設(shè)計(jì)的數(shù)字環(huán)路參數(shù)仿真得到</p><p&g

72、t;　?。?）環(huán)路捕獲帶 。</p><p>　　(2)將數(shù)字環(huán)路濾波器參數(shù) G1,G2的值同時(shí)減小為 。</p><p>　　此時(shí)G2= = =0.04875</p><p><b>　　推得 則</b></p><p><b>　　求得3dB帶寬</b></p>

73、<p><b>　　圖4-3</b></p><p>　　圖4-4為改變參數(shù)后的H( )的幅頻響應(yīng)曲線。</p><p><b>　　圖4-4</b></p><p><b>　　圖4-6</b></p><p>　　由圖可見-3dB截止頻率點(diǎn)在100Hz附近，與更

74、改參數(shù)后的環(huán)路帶寬值相近，仿真結(jié)果符合設(shè)計(jì)結(jié)果</p><p>　　圖4-6為改變參數(shù)后的H( )的幅頻響應(yīng)曲線。</p><p>　　由圖可見，-3dB截止頻率點(diǎn)在200Hz附近，與更改參數(shù)后的環(huán)路帶寬值相近，仿真結(jié)果符合設(shè)計(jì)結(jié)果。</p><p>　　圖4-7為改變參數(shù)后的誤差電壓 和控制電壓 的波形。</p><p><b>

75、　　圖4-7</b></p><p><b>　　圖4-8</b></p><p>　　可見，當(dāng)數(shù)字環(huán)路濾波器參數(shù)G1,G2的值同時(shí)增大為2g1,2g2時(shí)， ，若G1,G2的值同時(shí)再增加時(shí)，環(huán)路捕獲帶不再變化，只是捕獲時(shí)間短了一些。因?yàn)镚1,G2的值同時(shí)再增加時(shí)，相當(dāng)于增大了 ,使得捕獲時(shí)間減小，但因采樣頻率的限制，捕獲帶不可能再增加。</p>

76、;<p>　　(4)在環(huán)路參數(shù)值G1,G2同的條件下，增大、減小采樣頻率。</p><p>　　圖4-8為采樣頻率為增大為7000Hz時(shí)的環(huán)路幅頻響應(yīng)。</p><p>　　圖4-9為采樣頻率為減小為5000HZ時(shí)的環(huán)路幅頻響應(yīng)。</p><p><b>　　圖4-9</b></p><p>　　由圖可見，

77、當(dāng)采樣頻率發(fā)生變化時(shí)，環(huán)路的幅頻響應(yīng)曲線也有相應(yīng)變化。圖4-8中的-3dB截止頻率點(diǎn)與圖4-9中的明顯不同。說(shuō)明采樣頻率參數(shù)對(duì)環(huán)路也有一定影響。</p><p><b>　　4.2 仿真結(jié)論</b></p><p>　　數(shù)字環(huán)路的捕獲帶不僅取決于環(huán)路參數(shù)，還與采樣速率密切相關(guān)。環(huán)路參數(shù)值G1,G2的增大減小使得捕獲帶相應(yīng)變化，但增大有上限；采樣速率增大減小也使得捕獲

78、帶相應(yīng)變化。所以設(shè)計(jì)環(huán)路參數(shù)G1,G2時(shí)，在滿足捕獲帶的基礎(chǔ)上盡量使G1,G2的值小，以便得到較小的 ，以增加環(huán)路抑制輸入噪聲的能力</p><p>　　第五章 論文以后的研究工作</p><p>　　使用Simulink工具箱對(duì)DPLL進(jìn)行仿真</p><p>　　5.1 異或鑒相器</p><p>　　構(gòu)成鑒相器的電路形式很多，按期電路

79、性質(zhì)可分為模擬鑒相器，數(shù)字鑒相器和取樣保持鑒相器，這里只簡(jiǎn)單的介紹其中的一種。異或門鑒相器 </p><p>　　異或門的邏輯真值表示于表5.1，邏輯符號(hào)圖如下。</p><p>　　從表1可知，如果輸入端A和B分別送入占空比為50%的信號(hào)波形， 則當(dāng)兩者存在相位差 時(shí)，輸出端F的波形的占空比與 有關(guān)。將F輸出波形通過(guò)積分器平滑，則積分器輸出波形的平均值，它同樣與 有關(guān)，這樣，我們就可以利

80、用異或門來(lái)進(jìn)行相位到電壓的轉(zhuǎn)換，構(gòu)成相位檢出電路。于是經(jīng)積分器積分后的平均值（直流分量）為：</p><p>　　U = * Δ θ/π (1)</p><p>　　不同的 ，有不同的直流分量 。 與U的關(guān)系可用邏輯圖來(lái)描述。從圖中可知，兩者呈簡(jiǎn)單線形關(guān)系：</p><p>　　= * (2)； 為鑒相靈敏度。</p><p><

81、;b>　　、 關(guān)系圖</b></p><p>　　5.2 一階環(huán)路濾波器</p><p>　　數(shù)字濾波器在實(shí)際應(yīng)用中一階數(shù)字環(huán)路濾波器的實(shí)現(xiàn)形式如圖：</p><p>　　圖5.1 一階數(shù)字環(huán)路濾波器</p><p>　　利用Simulink實(shí)現(xiàn)形式為</p><p>　　5.3 壓控振蕩器&

82、lt;/p><p>　　具體思路由以下仿真圖</p><p><b>　　第六章 結(jié)束語(yǔ)</b></p><p>　　本文介紹了DDS、PLL的原理，并對(duì)DDS行了建模、仿真，介紹了MATLAB的Simulink的使用，分析了DPLL的入鎖過(guò)程，以及環(huán)路參數(shù)對(duì)其影響，為系統(tǒng)設(shè)計(jì)提供了理論依據(jù)。DDS利用Simulink工具箱給于了實(shí)現(xiàn)。對(duì)信號(hào)源的

83、輸出方法進(jìn)行了仿真實(shí)驗(yàn)，利用不同的采樣點(diǎn)數(shù)仿真不同的輸出頻率，本文只是對(duì)DPLL進(jìn)行了M文件的仿真，做了硬件初步設(shè)計(jì)，要達(dá)到DPLL全部Simulink實(shí)現(xiàn)，還有很多工作要做，如仿真調(diào)試，算法的改進(jìn)等等，這些問(wèn)題的解決還需要付出很多的努力與汗水。</p><p><b>　　參考文獻(xiàn)</b></p><p>　　季仲梅、楊洪生、王大鳴、劉正軍 編著 通信中的同步技

84、術(shù)及應(yīng)用. 北京：清華大學(xué)出版社，2008</p><p>　　張磊、畢靖、郭蓮英 編著 MATLAB實(shí)用教程. 北京：人民郵電出版社，2008</p><p>　　胡華春、石玉 編著 數(shù)字鎖相環(huán)路原理與應(yīng)用. 上海：上?？茖W(xué)技術(shù)出版社，1990</p><p>　　羅偉雄、韓力、丁志杰 編 鎖相技術(shù)及其應(yīng)用. 北京：北京理工大學(xué)出版社，1990

85、</p><p>　　王家禮、孫璐 編 頻率合成技術(shù). 西安：西安電子科技大學(xué)出版社，2009</p><p><b>　　謝辭</b></p><p>　　首先，我謹(jǐn)向我的導(dǎo)師季仲梅教授致以衷心的謝意。在季教授的悉心指導(dǎo)和言傳身教下，我才能得以順利完成本科生畢業(yè)論文工作。本文的選題和研究工作傾注了季老師的大量心血和諄諄教誨。在學(xué)術(shù)上，

86、季教授具有高瞻遠(yuǎn)矚的目光，又具有堅(jiān)實(shí)廣博的知識(shí)。在工作上，季教授既具旺盛持久的工作熱情，又具嚴(yán)謹(jǐn)認(rèn)真的治學(xué)態(tài)度。在我大三、大四階段的學(xué)習(xí)期間，無(wú)論學(xué)習(xí)科研，工作態(tài)度，以及為人處事，都受益匪淺并將影響我終身。</p><p>　　感謝與我同期畢業(yè)的李現(xiàn)慧同學(xué)。在我的論文準(zhǔn)備和寫作期間，他們自始至終給予我巨大的幫助和指導(dǎo)。</p><p>　　我要衷心地感謝我的父母和家人，沒有他們的關(guān)心和愛護(hù)

87、，我不可能度過(guò)這漫長(zhǎng)而艱難的求學(xué)生涯，并最終完成學(xué)業(yè)。</p><p>　　最后，我要感謝評(píng)審本論文的老師、教授和專家:，感謝他們抽出寶貴的時(shí)間來(lái)閱讀本文，并提出寶貴的意見和建議。</p><p>　　附錄一 MATLAB簡(jiǎn)介</p><p>　　MATLAB是Matrix laboratory(矩陣實(shí)驗(yàn)室)的縮寫，它是以線性代數(shù)軟件包LINPACK和特征值計(jì)算軟

88、件包EISPACK的子程序?yàn)榛A(chǔ)發(fā)展起來(lái)的一種開放型程序設(shè)計(jì)語(yǔ)言。MATLAB將計(jì)算、可視化和編程等功能集于一個(gè)易于使用的環(huán)境，具體功能強(qiáng)大、簡(jiǎn)單易學(xué)、編程效率高的特點(diǎn)，是目前世界上最流行的仿真計(jì)算軟件之一。</p><p><b>　　MATLAB的組成</b></p><p>　　1.MATLAB是這一系列軟件的基石，提供了高級(jí)科學(xué)計(jì)算語(yǔ)言，是進(jìn)行數(shù)據(jù)分析算法開發(fā)

89、的集成開發(fā)環(huán)境，具有數(shù)值計(jì)算，圖形繪制，程序設(shè)計(jì)等功能。</p><p>　　2.MATLAB擴(kuò)展包含MATLAB編譯器，WEB服務(wù)器，數(shù)據(jù)庫(kù)工具箱,報(bào)表生成器等，用來(lái)支持在MATLAB環(huán)境中對(duì)系統(tǒng)的實(shí)施與開發(fā)</p><p>　　3.工具箱是針對(duì)解決特定種類問(wèn)題而特別制作的一系列MATLAB函數(shù)庫(kù)。它具有開放性和可擴(kuò)展性。</p><p>　　4. Simulin

90、k是交互式動(dòng)態(tài)系統(tǒng)建模，仿真和分析的圖形環(huán)境，也是進(jìn)行基于模型的嵌入式系統(tǒng)開發(fā)的基礎(chǔ)開發(fā)環(huán)境。它把模塊圖形界面和MATLAB主要數(shù)值，圖形和語(yǔ)言函數(shù)有效地組合起來(lái)，從而具有生動(dòng)的模擬能力。他對(duì)線性非線性，連續(xù)，離散系統(tǒng)均適用，可以針對(duì)控制系統(tǒng)，信號(hào)處理以及通信系統(tǒng)等進(jìn)行系統(tǒng)的建模，仿真，分析等工作。</p><p>　　附錄二 H( )的幅頻響應(yīng)</p><p>　　prompt ={&

91、#39;鍵入載波頻率G1','鍵入載波頻率G2：'};</p><p>　　tittle='設(shè)置基本參數(shù)';line=1;</p><p>　　def={'0.0975','0.005'};</p><p>　　XY=inputdlg(prompt,tittle,line,def);</p

92、><p>　　g1=str2num(XY{1});g2=str2num(XY{2});</p><p>　　w=0:0.001:(0.5*pi);</p><p><b>　　% g1=g1*8</b></p><p><b>　　%g2=g2*8</b></p><p>　　h

93、z=((g1+g2)*exp(-j*w)-g1*exp(-2*j*w))./(1+(g1+g2-2)*exp(-j*w)+(1-g1)*exp(-2*j*w));</p><p>　　%(0.45*exp(-j*w)-0.18*exp(-2*j*w))</p><p>　　m1=abs(hz);</p><p>　　m2=angle(hz);</p>

94、<p>　　mang=20*log10(abs(hz));</p><p>　　semilogx((w/(2*pi))*6000,mang,'linewidth',2);</p><p>　　xlabel('f(hz)');ylabel('20lg(abs(H(j2\pif))');</p><p>　　TI

95、TLE('閉環(huán)頻率特性')</p><p><b>　　grid on</b></p><p>　　附錄三 誤差電壓 和控制電壓 </p><p>　　prompt={'鍵入環(huán)路濾波器參數(shù)G1:','鍵入環(huán)路濾波器參數(shù)G2:'};</p><p>　　title='

96、;設(shè)置環(huán)路濾波器參數(shù)';line=1;</p><p>　　def={'0.045','0.00225'};</p><p>　　Y=inputdlg(prompt,title,line,def)</p><p>　　G1=str2num(Y{1});G2=str2num(Y{2});</p><p>

97、　　original_seq=ones(1,50);</p><p>　　original_seq=[original_seq,original_seq];</p><p>　　original_seq=[original_seq,original_seq,original_seq,original_seq,original_seq];</p><p>　　orig

98、inal_seq=2*original_seq-1;</p><p>　　sample_impluse=ones(1,32);</p><p>　　sampled_seq=kron(original_seq, sample_impluse);</p><p>　　sample_num=length(sampled_seq);</p><p>

99、　　m1=1:sample_num;</p><p>　　fc=1500000;</p><p>　　delta_f=300;</p><p><b>　　fs=48000;</b></p><p>　　sampled_sig=cos(((fc+delta_f)/fs)*2*pi*m1).*sampled_seq;<

100、;/p><p>　　length_sampled=length(m1);</p><p>　　w=blackman(60);</p><p>　　lpf_coff=fir1(59,0.5,'low',w);</p><p>　　Ud_i=zeros(1,length_sampled);</p><p>　　

101、Ud_i_filtered=zeros(1,length_sampled);</p><p><b>　　phase=0;</b></p><p><b>　　acc_1f=0</b></p><p>　　1=length(sample_impluse);</p><p>　　l=length(sa

102、mple_impluse);</p><p>　　Uc=zeros(1,length(original_seq));</p><p>　　for m1=1:(sample_num/1)</p><p>　　for m2=1:1</p><p>　　phase=phase+Uc(m1);</p><p>　　Ud_i_t

103、emp(m2)=sampled_sig(1*(m1-1)+m2)*sin(((m1-1)*1+m2)*pi/2+phase);</p><p><b>　　end;</b></p><p>　　Ud_i_filtered(m1)=sum(Ud_i_temp)/1;</p><p>　　Ud(m1)=Ud_i_filtered(m1);</

104、p><p>　　acc_1f=acc_1f+Ud(m1);</p><p>　　Uc(m1+1)=Ud(m1)*G1+acc_1f*G2;</p><p>　　tp1(m1)=sin(phase);</p><p>　　tp2=cos(phase);</p><p><b>　　end;</b><

105、;/p><p>　　figure(1);</p><p>　　subplot(2,1,1)</p><p><b>　　plot(Ud);</b></p><p>　　xlabel('K');ylabel('ud(k)')</p><p>　　subplot(2,1,2