请问，BJT的beta曲线如何仿真呢？

hspice examples
beta.sp
* beta.sp
* level 2beta optimization
*
*.. optimize ib vs ie, (practically the same as optimizing beta vs ic)
*..bf,ikf,nkf,ise,ne
*
*.. print: betam, betad vs icm or icd
*
*.. note: bf turns out to be much higher than max beta.
*..bf is actually a weak parameter is this example, because
*..high-current and low-current roll-off regions have merged
.option kcltest=0
.option nomodvntol=1e-12 reli=1e-4 relv=1e-4 gmindc=1e-14
+ post
.model optmod opt relin=1e-4 grad=1e-5 close=10 cut=2
+ cendif=1e-6 relout=1e-4 max=1e6
+ itropt=15
* common base circuit, forcing ie and vbc
vc1 0 c1 vbc1
vs1 s1 0 vsb1
iie1 e1 0 "ib1+ic1"
q1 c1 0 e1 s1 npn1
.model npn1 npn level=2
+ eg = 1.1af=1xcjc=0.95subs=1
+ cjs=0.5p tf=1uitf=1m vtf=20 xtf=1.0e-01
+ tr=1u cje=1.2e-12cjc=0.3p fc=0.9 mje=0.5 mjc=0.5
+ rb=rbrbm=rbmirb=irb
+ bf=bf ikf=ikf nkf=nkf
+ vaf=vaf
+ is=is ise=ise
+ rc=rc re=re
+ nf=nf ne=ne
+ br=1 ikr=1m var=10
+ iss=1e-17 brs=20
.param
+ bf=opt1(450,10,800)
+ ise=opt1(1e-17,1e-18,1e-16)
+ ne=opt1(1.15,1.10,1.5)
+ ikf=opt1(1m,10u,10m)
+ nkf=opt1(0.5,0.3,0.7)
+ vaf=100
+ re=5 rc=70 rb=100 rbm=10 irb=100u
+ is=3.0928e-16
+ nf=1.0184
* optimization using circuits 1,3
.dc data=dat1 sweep optimize=opt1
+ results=errb model=optmod
* beta error, ie forced
.meas dc errb err1 par(ib1) i2(q1)
* print-back, same as opt
.dc data=dat1
.print icm=par('i1(q1)') icd=par(ic1)
+ betam=par('i1(q1)/i2(q1)') betad=par('par(ic1)/par(ib1)')
.data dat1
vbe1 vbc1 ic1 ib1 vsb1 is1
0.48-2.001.7338e-069.9448e-09 -5.009.9448e-09
0.50-2.003.3684e-061.8730e-08 -5.001.8730e-08
0.52-2.006.5428e-063.5076e-08 -5.003.5076e-08
0.54-2.001.2652e-056.5198e-08 -5.006.5198e-08
0.56-2.002.4450e-051.2070e-07 -5.001.2070e-07
0.58-2.004.7040e-052.2264e-07 -5.002.2264e-07
0.60-2.008.9598e-054.0884e-07 -5.004.0884e-07
0.62-2.001.6920e-047.5054e-07 -5.007.5054e-07
0.64-2.003.1084e-041.3622e-06 -5.001.3622e-06
0.66-2.005.4988e-042.4304e-06 -5.002.4304e-06
0.68-2.009.2706e-044.2354e-06 -5.004.2354e-06
0.70-2.001.4734e-037.1378e-06 -5.007.1378e-06
0.72-2.002.1940e-031.1550e-05 -5.001.1550e-05
0.74-2.003.0646e-031.7992e-05 -5.001.7992e-05
0.76-2.004.0358e-032.7142e-05 -5.002.7142e-05
0.78-2.005.0692e-033.9746e-05 -5.003.9746e-05
0.80-2.006.1632e-035.6590e-05 -5.005.6590e-05
0.82-2.007.3108e-037.8304e-05 -5.007.8304e-05
0.84-2.008.5092e-031.0542e-04 -5.001.0542e-04
0.86-2.009.7526e-031.3844e-04 -5.001.3844e-04
0.88-2.001.1040e-021.7808e-04 -5.001.7808e-04
0.90-2.001.2362e-022.2452e-04 -5.002.2452e-04
.enddata
.end

对于BJT来说，应该是接成 “共射放大”形式。
1.在基极上串联一个idc=ib（参数）的理想电流源，变量名为ib；
2.在Virtuoso Analog Design Environment中的choosing Analysis中选择DC Analysis
并选择sweep variable 选变量ib
sweep range：1nA~1uA
3.仿真
4.仿出基极和集电极的电流波形来
如：IC随变量ib变化有一条曲线；
IB随变量ib变化应是一条斜率为1通过原点的直线
5.调用计算器计算出IC/IB随IB的变化曲线；
6.最后双击波形坐标轴的横坐标，将横坐标改为IC
最后就出来了IC/IB~IC即Belta随IC的变化曲线

Hspice的demo里有BJT beta仿真文件，是和semi_bamboo说的一样，仿真不同IB下，IC随VCE的变化曲线