import matplotlib.pyplot as plt # %% Muda o current working directory do Python para a pasta onde este script # se encontra cwd = os.path.dirname(__file__) # Pega a pasta de trabalho atual os.chdir(cwd) # %% # Cria objeto para stream de dados com o LabJack U3 com o sensor # de temperatura e umidade EI1050 # tempHumid = lju3ei1050.main() tempHumid = None # Para testes com LabJack offline # %% Caso já tenha uma medição em curso, carregue o MeasurementSetup e # MeasurementData MS, D = rmr.med_load('med-teste') # %% Carrega sinais de excitação e cria dicionário para o setup da medição excitationSignals = {} excitationSignals['varredura18'] = pytta.generate.sweep( # Geração do sweep (também pode ser carregado projeto prévio) freqMin=20, freqMax=20000, fftDegree=18, startMargin=0.05, stopMargin=3, method='logarithmic', windowing='hann', samplingRate=44100) excitationSignals['varredura17'] = pytta.generate.sweep( # Geração do sweep (também pode ser carregado projeto prévio)
analysisFileName = '11_05_20_revCh_sem-nada.hdf5' V_revCh = 207 MPP = pytta.roomir.MeasurementPostProcess(nthOct=3, minFreq=100, maxFreq=10000) IREndManualCut = 7 skipInCompensation = True skipOutCompensation = True skipBypCalibration = True skipIndCalibration = True skipRegularization = False skipSave = True plots = False # %% MeasurementData MSrevCh, DrevCh = rmr.med_load('med-final_revCh-Lpe') analysisFileNameTR = '11_05_20_revCh-TR.hdf5' # %% Calcula channelcalibir getChCalib = DrevCh.get('channelcalibration', 'Mic1') ChCalibir = DrevCh.calculate_ir(getChCalib, calibrationTake=1, skipSave=skipSave) for name, msdThng in ChCalibir.items(): print(name) if plots: #p1 = [IR.plot_freq(ylim=[-10, 120]) for IR in msdThng.measuredSignals] p2 = msdThng.measuredSignals[0].plot_freq(xLim=None) p3 = msdThng.measuredSignals[0].plot_time_dB(xLim=None) #%% Calcula sourcerecalibir