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()