def testDiodeParameterExtraction():
"""
пробуем экстрагировать коэффициенты из модели диода
коэффициенты модели: Ток утечки Is, коэффициент неидеальности N, омическое сопротивление, параллельное диоду R
входные параметры: напряжение, приложенное источником к системе резистор-диод
+-----------|||||---------->|--------- -
Резистор подключен до диода
:return:
"""
jacf=jacdiode
funcf=diode
#теперь попробуем сделать эксперимент.
c={}
Ve=np.asmatrix( [0.1] )
btrue=[1.238e-14, 1.8]
bstart=np.array(btrue)-np.array(btrue)*0.2
bend=np.array(btrue)+np.array(btrue)*0.2
binit=[1e-10,1.1]
xstart=[0.01]
xend=[2]
N=30
print("performing normal research:")
startplan = o_p.makeUniformExpPlan(xstart, xend, N)
print(len(startplan))
measdata = o_p.makeMeasAccToPlan(funcf, startplan, btrue, c, Ve)
planplot1=[x[0] for x in startplan]
measplot1=[x['y'][0] for x in measdata]
plt.plot(planplot1, measplot1, 'bo')
plt.show()
gknu=o_e.grandCountGN_UltraX1 (funcf, jacf, measdata, binit, c, NSIG=6, sign=1)
#как мы помним, в случае неявных функций должно ставить sign=0
print (gknu[0])
print (o_q.getQualitat(measdata, gknu[0], Ve, funcf, c))
#plotting Sk graph
#TODO better organize: this code to estimation or somewhere
rng=np.arange(0,len(gknu[3]))
plt.plot(rng , gknu[3], label='Sk drop')
plt.legend(loc='upper left')
plt.ylabel('Sk')
plt.xlabel('Interation')
plt.grid()
plt.show()
N=20
print("performing aprior plan:")
oplan=o_ap.grandApriornPlanning (xstart, xend, N, bstart, bend, c, Ve, jacf, funcf, Ntries=6)[1]
o_p.writePlanToFile(oplan)
measdata = o_p.makeMeasAccToPlan(funcf, oplan, btrue, c,Ve )
planplot1=[x[0] for x in oplan]
measplot1=[x['y'][0] for x in measdata]
plt.plot(planplot1, measplot1, 'bo')
plt.show()
gknu=o_e.grandCountGN_UltraX1 (funcf, jacf, measdata, binit, c, NSIG=6, sign=1)
print (gknu[0])
print (o_q.getQualitat(measdata, gknu[0], Ve, funcf, c))
#plotting Sk graph
#TODO better organize: this code to estimation or somewhere
rng=np.arange(0,len(gknu[3]))
plt.plot(rng , gknu[3], label='Sk drop')
plt.legend(loc='upper left')
plt.ylabel('Sk')
plt.xlabel('Interation')
plt.grid()
plt.show()
