John Smith
Guest
Sat Jun 21, 2008 1:01 am
Quote:
Uzyskanie takich wyników jest nawet teoretycznie niemożliwe (dla ADC 16bit)
z uwagi na szum kwantyzacji.
Możliwe, możliwe - szum kwantyzacji jest rozłożony w widmie.
Wyniki symulacji:
http://www.zind.ikem.pwr.wroc.pl/tmp/elektronika/spectrum1.pdf
Wykres niebieski sama kwantyzacja obliczona wg: floor(sin()*2^15)
Tutaj:
www.analog.com/library/analogdialogue/archives/39-06/Chapter%202%20Sampled%20Data%20Systems%20F.pdf
na stronie 2.42 dół strony od akapitu zaczynającego się od "Figure 2.44 shows the FFT output for an
ideal 12-bit ADC" jest przedstawiony sposób liczenia noise floor dla przetwornika idealnego.
Po podstawieniu do odpowiednich wzorów, N=16 bit i M=4096 punktów FFT otrzymuję:
RMS QUANTIZATION NOISE LEVEL = 98.08dB zaś
FFT NOISE FLOOR = 98.08 + 33.11 = 131.19dB
a to jest znacznie odległe od uzyskanych przez Ciebie 300dB.
Za chiny ludowe nie chce wyjść inaczej. Może robię coś źle?
K.
J.F.
Guest
Sat Jun 21, 2008 11:15 am
On Fri, 20 Jun 2008 11:41:02 +0000 (UTC),
Quote:
Mnie wyszło (dla tego fs), że sygnał ma 4973.7 Hz (odchylenie std 0.5Hz)
z wyjątkiem ...
... no właśnie - w przebiegu masz 6 chwil czasowych (widać je ładnie na
wykresie plot(diff(X))), w których są duże zakłócenia impulsowe (sieją po
całym paśmmie) i one sporo zaniżają SNR.
I tak na oko odbiegaja o blisko 4096 i 8192.
Ciekawe .. przeklamuje transmisja cyfrowa, czy czasem przetwornikowi
przeklamuje SAR ?
Quote:
Jak weźmiesz fragment przebiegu y=x(2500+(1:2*1024)); to na widmie ładnie
widać 2 i 3 harmoniczną oraz coś przy 60 kHz
A tak swoja droga - usiluje dopasowac sinusoide ... i nie wychodzi :-(
0.1064546 rad/probke [ 4235,69 Hz przy 250ksps].
Najlepiej pasujaca amplituda - 30502
A jednak roznica siega ok 50 i jest pierwsza harmoniczna
[zmodulowana kolejna sinusoida].
No i nie wiem - nie trafilem z czestotliwoscia, czy amplituda
minimalnie plywa ?
Ale zawartosc harmonicznych i tak swietna.
J.
Guest
Sat Jun 21, 2008 7:19 pm
J.F. <jfox_xnospamx@poczta.onet.pl> wrote:
Quote:
Jak weźmiesz fragment przebiegu y=x(2500+(1:2*1024)); to na widmie ładnie
widać 2 i 3 harmoniczną oraz coś przy 60 kHz
A tak swoja droga - usiluje dopasowac sinusoide ... i nie wychodzi :-(
0.1064546 rad/probke [ 4235,69 Hz przy 250ksps].
Greg pisał, że fs = 75MHz/256 ~ 292.97 kHz - czyli masz 4963.7 Hz
Quote:
Najlepiej pasujaca amplituda - 30502
A jednak roznica siega ok 50 i jest pierwsza harmoniczna
[zmodulowana kolejna sinusoida].
Przesadzasz

20*log10(50/30502) = -55.70 dB
A to moje próby estymacji:
1) funkcja expfit (Prony) z octave:
f=0.1064541 rad/probke
A=30534.14 amplituda
różnica przekracza 33,
ale rms różnicy to tylko SNR=63.740 dB (w odniesieniu do rms sygnału)
A Prony jak to Prony

w takim sygnale nic więcej nie znajdzie.
2) mój faworyt - metoda pęku macierzowego Sarkar, a konkretnie to
implementacja:
http://groups.google.pl/group/comp.soft-sys.matlab/msg/c71f1639d5253b1c
2a) Dla standardowych ustawień (tol=1e-3):
f=0.10645458 rad/probke
A=30521.52
SNR=65.370 dB, max różnica 21.4
2b) dla tol=1e-5 mpencil znalazł 2 i 3 harmoniczną oraz wspomniane
zakłocenie 60kHz, ale także jakieś fluktuacje f1 (nie wiem czy to błędy
numeryki czy własność sygnału ;)
SNR=84.303 dB, max różnica 4.48
dB Hz rad/sample
-59.52220 0.00000 0.00000
-81.22014 4.89945 0.10508
0.00029 4.96370 0.10645
-99.80767 5.20857 0.11171
-84.77337 9.93031 0.21297
-86.62847 14.89060 0.31935
-92.12940 60.04270 1.28771
Quote:
No i nie wiem - nie trafilem z czestotliwoscia, czy amplituda
minimalnie plywa ?
Wszystko możliwe
Quote:
Ale zawartosc harmonicznych i tak swietna.
Potwierdzam
Guest
Sat Jun 21, 2008 7:21 pm
John Smith <dam9723@buziaczek.pl> wrote:
Quote:
Dziękuję za śliczny link, jak znajdę czas to odpowiem na twoje obserwacje.