简单的拟合电子管实测曲线生成仿真模型的教程。

壳壳

这几年多款曲线图示测胆仪面市，不少胆机爱好者都添置了“核”武器，有了曲线，画图设计电路也方便了，有些朋友上升到想用仿真电路验证一下设计，仿真需要模型，以前的有些前辈已经结合多种办法建立和分享了好多电子管的模型（功德无量！点赞）。


三极管的模型不难，难的是多极管模型，需要很多次的反复操作和学习，多极管的帘栅电流很关键，没有帘栅流的模型是废的，仿真时跑不出帘栅电流。所以，当我学了一丁点皮毛的时候，我就沾沾自喜，迫不及待想做一个分享。水平有限，如有错漏，请批判并指正以免因愚学的差错误导别人。拟合软件在这里下载（有可能被屏蔽打不开）：http://www.dmitrynizh.com/tubeparams_image.htm

话不多，以图为例，不明之处不及解释，请包涵谅解。

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

完！





**** 6J8P ******************************************
* Created on 04/22/2026 20:43 using paint_kip.jar
* www.dmitrynizh.com/tubeparams_image.htm
* Plate Curves image file: 6J8P.jpg
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 6J8P P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 6J8P P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=20 KG1=3695 KP=150 KVB=12 VCT=0.115 EX=1.4 KG2=4200 KNEE=4.2 KVC=2.082
+ KLAM=2.209E-8 KLAMG=6.9E-5   KNK=-0.044 KNG=0.006 KNPL=50 KNSL=11 KNPR=120 KNSR=29
+ CCG=7P CGP=0.007P CCP=7P RGI=2000.0
* Vp_MAX=400 Ip_MAX=12 Vg_step=0.5 Vg_start=0 Vg_count=11
* X_MIN=98 Y_MIN=147 X_SIZE=958 Y_SIZE=455 FSZ_X=1616 FSZ_Y=876 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=2.8 Vg_qui=-2.5 Vp_qui=300
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y
* showIg1=n isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n
* UL=0.43 EG2=100 gridLevel2=n addKink=y isTanhKnee=n advSigmoid=n
*----------------------------------------------------------------------------------
RE1  7 0  1G    ; DUMMY SO NODE 7 HAS 2 CONNECTIONS
E1   7 0  VALUE=  ; E1 BREAKS UP LONG EQUATION FOR G1.
+{V(G2,K)/KP*LOG(1+EXP((1/MU+(VCT+V(G,K))/SQRT(KVB+V(G2,K)*V(G2,K)))*KP))}
RE2  6 0  1G    ; DUMMY SO NODE 6 HAS 2 CONNECTIONS
E2  6 0  VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))} ; Kg1 times KIT current
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*ATAN(V(P,K)/KNEE)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))}
G1   P  K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
* Alexander Gurskii screen current, see audioXpress 2/2011, with slope and kink added
RE43 43 K 1G ; Dummy
E43  43 G2 VALUE={0} ; Dummy
G2   43 K  VALUE={V(6)/KG2*(KVC-ATAN(V(P,K)/KNEE))/(1+KLAMG*V(P,K))-V(22)}
RCP  P K  1G     ; FOR CONVERGENCE
C1   K G  {CCG}  ; CATHODE-GRID 1
C2   G P  {CGP}  ; GRID 1-PLATE
C3   K P  {CCP}  ; CATHODE-PLATE
R1   G 5  {RGI}  ; FOR GRID CURRENT   
D3   5 K  DX     ; FOR GRID CURRENT }
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N)
.ENDS
*$

闭关锁国

多级管要勾选KINK选项，否则仿真时没有帘栅极电流。

壳壳

另一个讨论仿真拟合的帖子：diyAudio -- Model Paint Tools Tutorial | diyAudio  https://www.diyaudio.com/communi ... ls-tutorial.310199/

128032949qq

这个软件还有负载线功能，建模做好后勾选Loadline for load Rp可以设置负载线计算谐波失真和增益等

P9M39

太好了，要认真实践！！！

P9M39

6J8P这几个字怎么打上去的？

初烧小子

感谢壳爷的分享精神

壳壳

2P29S




**** 2P29S ******************************************
* Created on 04/23/2026 22:36 using paint_kip.jar
* www.dmitrynizh.com/tubeparams_image.htm
* Plate Curves image file: 2P29S.jpg
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 2P29S P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 2P29S P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=10 KG1=4840 KP=156 KVB=12 VCT=0.2 EX=1.47 KG2=4200 KNEE=13.08 KVC=1.67
+ KLAM=2E-6 KLAMG=3E-4 KNEE2=20 KNEX=30  KNK=-0.044 KNG=0.006 KNPL=50 KNSL=11 KNPR=120 KNSR=29
+ CCG=5P CGP=0.06P CCP=6P RGI=2000.0
* Vp_MAX=300 Ip_MAX=30 Vg_step=1 Vg_start=0 Vg_count=13
* X_MIN=61 Y_MIN=44 X_SIZE=950 Y_SIZE=454 FSZ_X=1616 FSZ_Y=876 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=2 Vg_qui=-6 Vp_qui=300
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y
* showIg1=n isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n
* UL=0.43 EG2=120 gridLevel2=n addKink=y isTanhKnee=y advSigmoid=n
*----------------------------------------------------------------------------------
RE1  7 0  1G    ; DUMMY SO NODE 7 HAS 2 CONNECTIONS
E1   7 0  VALUE=  ; E1 BREAKS UP LONG EQUATION FOR G1.
+{V(G2,K)/KP*LOG(1+EXP((1/MU+(VCT+V(G,K))/SQRT(KVB+V(G2,K)*V(G2,K)))*KP))}
RE2  6 0  1G    ; DUMMY SO NODE 6 HAS 2 CONNECTIONS
E2  6 0  VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))} ; Kg1 times KIT current
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))}
G1   P  K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
* Alexander Gurskii screen current, see audioXpress 2/2011, with slope and kink added
RE43 43 K 1G ; Dummy
E43  43 G2 VALUE={0} ; Dummy
G2   43 K  VALUE={V(6)/KG2*(KVC-ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2))/(1+KLAMG*V(P,K))-V(22)}
RCP  P K  1G     ; FOR CONVERGENCE
C1   K G  {CCG}  ; CATHODE-GRID 1
C2   G P  {CGP}  ; GRID 1-PLATE
C3   K P  {CCP}  ; CATHODE-PLATE
R1   G 5  {RGI}  ; FOR GRID CURRENT   
D3   5 K  DX     ; FOR GRID CURRENT }
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N)
.ENDS
*$

壳壳

4J1S



**** 4J1S ******************************************
* Created on 04/23/2026 22:47 using paint_kip.jar
* www.dmitrynizh.com/tubeparams_image.htm
* Plate Curves image file: 4J1S.jpg
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 4J1S P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 4J1S P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=19.6 KG1=3256 KP=67.57 KVB=12 VCT=0.2 EX=1.512 KG2=4200 KNEE=11.76 KVC=2.056
+ KLAM=2.5E-8 KLAMG=6E-5 KNEE2=20 KNEX=30  KNK=-0.044 KNG=0.006 KNPL=50 KNSL=11 KNPR=120 KNSR=29
+ CCG=4P CGP=0.007P CCP=5P RGI=2000.0
* Vp_MAX=300 Ip_MAX=10 Vg_step=0.5 Vg_start=-1 Vg_count=11
* X_MIN=86 Y_MIN=139 X_SIZE=956 Y_SIZE=451 FSZ_X=1616 FSZ_Y=876 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=2 Vg_qui=-3.5 Vp_qui=300
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y
* showIg1=n isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n
* UL=0.43 EG2=100 gridLevel2=n addKink=y isTanhKnee=y advSigmoid=n
*----------------------------------------------------------------------------------
RE1  7 0  1G    ; DUMMY SO NODE 7 HAS 2 CONNECTIONS
E1   7 0  VALUE=  ; E1 BREAKS UP LONG EQUATION FOR G1.
+{V(G2,K)/KP*LOG(1+EXP((1/MU+(VCT+V(G,K))/SQRT(KVB+V(G2,K)*V(G2,K)))*KP))}
RE2  6 0  1G    ; DUMMY SO NODE 6 HAS 2 CONNECTIONS
E2  6 0  VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))} ; Kg1 times KIT current
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))}
G1   P  K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
* Alexander Gurskii screen current, see audioXpress 2/2011, with slope and kink added
RE43 43 K 1G ; Dummy
E43  43 G2 VALUE={0} ; Dummy
G2   43 K  VALUE={V(6)/KG2*(KVC-ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2))/(1+KLAMG*V(P,K))-V(22)}
RCP  P K  1G     ; FOR CONVERGENCE
C1   K G  {CCG}  ; CATHODE-GRID 1
C2   G P  {CGP}  ; GRID 1-PLATE
C3   K P  {CCP}  ; CATHODE-PLATE
R1   G 5  {RGI}  ; FOR GRID CURRENT   
D3   5 K  DX     ; FOR GRID CURRENT }
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N)
.ENDS
*$

tipware

感谢壳爷的分享精神