在本文中，我們將重點介紹與傳統(tǒng)編程語言的數值方案相比，符號方案在二極管電路仿真方面的優(yōu)勢。

什么是限幅電路？

限幅電路(也稱為限幅器或削波器)與簡單的半波整流器相似，不同之處在于其輸出信號是二極管兩端的壓降，而不是電阻器兩端的壓降。第二個特別之處則在于其二極管上串聯了一個電池。在分析這樣的電路之前，有必要先回想一下在開發(fā)仿真軟件時遇到的困難。在本期中，我們將向大家展示，傳統(tǒng)編程語言需要用多行代碼才能實現的一個例程，使用Mathematica軟件則只需一行代碼就能解決問題。

數值方案

讓我們從最簡單的電路開始：一個電阻R，輸入信號為Vin(t)=VM sin ωt。我們將這個量視為輸出信號，而電流按系數R?1縮放。在這種情況下，仿真可以說是完全相同的，因為它再現了輸入信號。然后將二極管D與電阻器R串聯，如圖1所示。

圖1：串聯電阻和二極管

通過應用基爾霍夫第二定律并考慮電壓-電流特性，經過簡單的步驟，我們就得出以下方程：

VT是以伏特為單位的熱當量，在室溫下等于26mV。因此，串聯RD中流動的電流就是方程(1)的解。這就是所謂的“函數方程”，因為未知數是函數而不是數字。很明顯方程(1)不能以封閉形式求解，因為我們無法獲得電路中流動的電流i(t)的解析表達式。之后，我們可以使用Mathematica進行數值計算，對連續(xù)變量t在周期區(qū)間[0,2π/ω]內進行離散采樣，即方程(1)的第二個因素。結果得到一個由N個解耦方程組成的系統(tǒng)：

而這可以通過Mathematica的Solve語句輕松解決。例如，對于：

解列表{x0,x1,…,xN}被繪制為tn的函數，如圖2所示。通過連接連續(xù)線段，我們獲得如圖3所示的趨勢。

圖2：以數值方式找到的解xn=in的趨勢

圖3：使用ListLinePlot指令獲得的解xn=in的趨勢

符號方案

將函數方程(1)改寫為以下形式更為方便：

其中x(t)是反向飽和電流歸一化的電流強度，a=Ri0/VT，并且f(t)=vin(t)。在符號模式下，我們強制Wolfram語言內核通過Solve指令找到解x(t)。相應的輸出通過Lambert W Function(朗伯W函數)表示，該函數在Wolfram語言中是內置的，并由指令ProductLog調用。對于正弦輸入，最好假設電阻R和峰值VM作為參數；對于角頻率，我們照常假設ω=20rad/s，因為結果與該參數無關。而R和VM才是決定性的，尤其是峰值對二極管整流效果的影響最大。例如，如果VM=1mV，則無論R的值是多少，二極管都不會進行整流，如圖4-5所示。

圖4：R=1Ω時以安培為單位的電流趨勢

圖5：R=1MΩ時以安培為單位的電流趨勢

基本上，這種表現是由于輸入信號歸一化為熱當量(以伏特為單位)，即VT=26mV。相對于VT，輸入信號的峰值必須足夠高。圖6的圖形證實了這一點，從圖中我們可以看出，整流器實際上是理想的，但由于R值較低，它幾乎在所有情況下都是正向偏置短路。當R=1MΩ并保持輸入峰值時，我們就可以得到圖7的圖形。由于歐姆元件的線性，R上的壓降將具有與電流相同的趨勢。否則，我們就會發(fā)現二極管兩端的壓降趨勢是扭曲的。如果我們使用Mathematica計算數值，就會發(fā)現圖8的奇怪趨勢，即輸入信號的負半波的最小峰值被截斷。這是Mathematica內核在搜索方程(4)的解時發(fā)出警告的結果。為了驗證上述內容，只要使用基爾霍夫第二定律計算電壓vD并繪制解的圖形即可(見圖9)。

圖6：R=1Ω且VM=10V時的電流強度趨勢

圖7：R=1MΩ且VM=10V時的電流強度趨勢

圖8：R=10Ω且VM=1V時二極管兩端壓降的趨勢

圖9：R=10Ω且VM=1V時二極管兩端壓降的正確趨勢