diff --git a/meritev-020321/oeo.pdf b/meritev-020321/oeo.pdf new file mode 100644 index 0000000..1e0cab0 Binary files /dev/null and b/meritev-020321/oeo.pdf differ diff --git a/meritev-020321/plot b/meritev-020321/plot index 1021f14..b13897e 100755 --- a/meritev-020321/plot +++ b/meritev-020321/plot @@ -12,8 +12,9 @@ set key off #bottom left set grid set grid mxtics set logscale x +set format x '10^{%T}' set output n.".pdf" -set title "Meritev faznega suma" +set title "Meritev faznega šuma" plot n.".csv" u 1:2 w lines diff --git a/meritev-020321/rezultat.pdf b/meritev-020321/rezultat.pdf deleted file mode 100644 index 07fe759..0000000 Binary files a/meritev-020321/rezultat.pdf and /dev/null differ diff --git a/meritev-050321/fazni-sum b/meritev-050321/fazni-sum new file mode 100755 index 0000000..4e15f7a --- /dev/null +++ b/meritev-050321/fazni-sum @@ -0,0 +1,81 @@ +#!/usr/bin/python3 + +print('*** Preprost zvočni spektralni analizator - S53MV 28.09.2020 ***') +print('*** popravek odziva za EliteBook 8540w - Andrej 11.2.2012 ***') +print('*** meritev faznega suma (interpolacija in povprecenje v frekvenci-sai) - Andrej 15.2.2012 ***') + +fs=192000 #Frekvenca vzročenja zvočnega signala (max 192000Hz) +pts=2048 #Število točk spektra 0...fs/2 (potenca 2 za učinkovit FFT) +avg=128 #Število povprečenj spektra + +import sounddevice as sd #Poiskati na spletu in naložiti s pip3! +import numpy as np #Uporaba učinkovitih funkcij numpy +import matplotlib.pyplot as plt #Risanje rezultata z matplotlib +import matplotlib.animation as animation + +sd.default.device = 0 +fskala=np.linspace(0,fs/2,2*pts+1)[:-1] #Izračunaj frekvenčno skalo [Hz] + +### podatki za kalibracijo sistema +dSB = -51.68+3.36 #dB +L = 1256 #metrov +tau = L*1.4676/3e8 + +### popravek odziva +odziv = np.loadtxt("5mhz-polinom-2048tock.csv", delimiter=',') +a = odziv[:,1] +popravek = np.max(a[4:-1])-a[4:-1] + + +def spekter(fs,pts,avg,fskala): #Zajem signala z zvočno kartico in povprečenje spektra [dB] + s=sd.rec(pts*(1+avg),samplerate=fs,channels=1) #Vzorčenje ADC zvočne kartice, 2D polje! + sd.wait() #Počakaj do konca vzorčenja + nz=1.e-12 #Dodatek za neničelni argument logaritma + s2=np.square(s) #Izpis jakosti/max [dB] + print('Jakost:',"%.1f"%(10*np.log10(nz+np.mean(s2))),'dB Vrh:',"%.1f"%(10*np.log10(nz+np.amax(s2))),'dB',end=' ') + o=1+np.cos(np.linspace(-np.pi,np.pi,2*pts)) #Okno dvignjeni kosinus + k=0 + s2=np.zeros(2*pts) + while k