# Rのユーティリティ関数を読み込む
attach(DATA+'RUtil.py')
|
|
# Vbeの計算
(1.38065E-23 * (273.15 + 25))/(1.60217E-19) *ln(1E-14/0.25E-3)
-0.615138361567583 -0.615138361567583 |
# トランジスタ本、図2-14のRb2の計算
(Ic, Hfe, Ib, Vcc, Vce, Rc, Vce, Vbe, Rb2) = var('Ic Hfe Ib Vcc Vce Rc Vce Vbe Rb2')
Hfe = 170
Vcc = 5
Vce = 2.6
Vbe = 0.6
Rc = 1.5E3
eq = [ Ic == Hfe * Ib, Vcc == Vce + Rc*(Ib + Ic), Vce == Vbe + Rb2*Ib]
solve(eq, [Ic, Ib, Rb2])
[[Ic == (34/21375), Ib == (1/106875), Rb2 == 213750]] [[Ic == (34/21375), Ib == (1/106875), Rb2 == 213750]] |
#CE01のゲイン
G = -8.86 - ( -35.63)
A = 10^(G/20)
print G, A
26.7700000000000 21.8021839718595 26.7700000000000 21.8021839718595 |
# CE01のゲイン測定値
freq = [20, 30, 40, 60, 100, 200, 400, 600, 1000, 2000, 4000, 6000, 10000, 20000]
orgG = [-34.56, -34.61, -34.64, -34.66, -34.66, -34.67, -34.67, -34.67, -34.67, -34.67, -34.68, -34.68, -34.69, -41.86]
ampG = [-8.51, -8.19, -8.07, -7.98, -7.93, -7.91, -7.90, -7.89, -7.89, -7.87, -7.89, -7.90, -7.90, -15.09]
difG = [ampG[i] - orgG[i] for i in range(len(freq))]
|
|
r_freq = r(freq).name('freq')
r_difG = r(difG).name('difG')
# グラフ表示
graph = preGraph("fig3.pdf")
r('plot(freq, difG, log="x")')
postGraph(graph, fac=0.7)
|
|
(Ic, Is, k, T, q, Vbe)=var('Ic Is k T q Vbe')
Is = 1e-14
k = 1.38065e-23
q = 1.60217e-19
Ic(Vbe, T) = Is*e^(q*Vbe/(k*T))
pT27 = plot(Ic(Vbe, 273.15 + 27) , [Vbe, 0, 1.6], ymax=2e-3)
pT77 = plot(Ic(Vbe, 273.15 + 77) , [Vbe, 0, 1.6], ymax=2e-3)
Vbias = 1.6
Re = 1e3
pVbias = plot((Vbias-Vbe)/Re, [Vbe, 0, 1.6])
(pT27 + pT77 + pVbias).show(figsize=5)
|
|
|
|