Rによる三極管モデル(http://www002.upp.so-net.ne.jp/ayumin/audio/pctube/node35.html)を使えば, 以下のようにして,各球の動作点と三定数を求めることができます.
> Ebb <- 240 # 電源電圧 > Ebb2 <- Ebb/2 # 1本分のプレート供給電圧 > Rk <- 1.5e3 # カソード(自己バイアス)抵抗 > Eg0 <- uniroot(function(Eg) Ip(t12AT7, Ebb2 + Eg, Eg) * Rk + Eg, c(-10, 0))$root > Eg0 # 動作点の Eg [1] -1.935316 > Ip(t12AT7, Ebb2+Eg0, Eg0) # 動作点のプレート電流 [1] 0.001290182 > mu(t12AT7, Ebb2+Eg0, Eg0) # 動作点の μ [1] 53.7236 > rp(t12AT7, Ebb2+Eg0, Eg0) # 動作点の rp [1] 25733.88 > gm(t12AT7, Ebb2+Eg0, Eg0) # 動作点の gm [1] 0.002087660
これを何回も繰り返すのはたいへんなので,
プログラムを書いて求めた値は表1のとおりです.