|
boB
1 z% q" M+ @( V- K
4 _3 N% m3 f/ _
+ \# v5 u' B# E. n1 A$ V( r# L6 y9 j* {0 J
Joined: 28 Aug 20068 }0 H+ }0 H$ ~3 l5 M' b* y5 U
Posts: 34
9 |0 ]7 |. X7 K! J8 u2 pLocation: Everett, WA9 N- j/ X( b; P% j
Posted: Mon Jan 21, 2008 8:45 pm Post subject: Re: 44.1Ksps sample rate enough for audio?
8 w- q9 N8 |/ B- M# ^
6 C4 g8 E3 n2 I t+ c--------------------------------------------------------------------------------
% o# U; ` z/ t3 w: {. l: l% Z6 A, _! n) j$ ?- Y1 J
lavrye wrote: / O, d5 x% H1 V
3 l" |- [; ~3 d6 X: a: S- V! b# Y( A/ D' B q' h0 J3 C
You can not tell the wave shape of a cycle from 2 dots. To understand sampling, you need to take into account very many cycles prior to that "cycle with "2 dots". You will soon realize that the samples (dots) hit different values on each of the previous cycles. With data from many prior cycles, we can reconstruct the original sampled wave. ; v3 x0 j! ^0 R" d1 @1 y: E
: q* r" ]1 \1 @
~~~~MORE STUFF~~~~
0 Y+ @" A% @. q" c6 t6 d2 y8 ~" L1 [7 U! Z7 t n* G: P
Regards
/ V' B3 H; E5 SDan Lavry ' I0 Q* w4 K# [# T3 ^- F
$ ?" n5 A! |3 N
1 `, O0 v7 @! X- M& ]$ j/ DMaybe it's exactly the same thing, Dan, but I look at the re-make of the sine wave at close to FS/2 from two "dots" as reconstructing the sine wave because of the bandwidth limitation at FS/2 itself, rather than the "history" of the sampled waveform ... 7 G, I! J3 n, X+ e* |5 e7 b6 @
; a3 `1 |1 L3 \
But I suppose the history IS because of the same bandwidth constraint ???
& B" [' J' @1 C5 T
8 V+ Y# H1 h; I3 r" |" }Is this a valid thought, or another way of saying the same thing ?
$ c9 @6 @8 J# Y) t5 b, m: _" p" d
/ g) W2 \8 j0 a) V0 gThanks,
' x2 ^ B% K0 e& sboB( e r" g% ]$ Z
_________________
( `7 M! M9 l8 M* ~* P8 w. WK7IQ
- L; ]% E- R0 r, { F% L, Y
8 m! {( i0 ~7 e$ K; L! bBack to top 5 I' Z" b* Y: m- H! R L
W- [- o( X9 m$ f( ]( _1 i# e
3 j1 L& b! y+ n t
lavrye8 X" d( ^( c: x. a% B8 L
Site Admin
3 @% t4 w0 ~0 B: c* _0 N/ M' @% I4 d8 k
5 u' G/ ?$ g/ w2 S* k' ?
Joined: 15 May 2005
F9 g O! x- u7 y/ S" YPosts: 4627 P2 `( p4 J. ^) e3 j& J
Location: Poulsbo, WA
9 u }: g( t/ V! N) pPosted: Mon Jan 21, 2008 11:53 pm Post subject: Re: 44.1Ksps sample rate enough for audio?
7 o- \1 a0 H, j5 {5 j% b. M1 @9 e) _. f
$ P9 i+ B/ \5 C9 y6 p--------------------------------------------------------------------------------
5 k! v: I2 X4 q7 M
4 o& S6 p# h- {& _/ M/ V. p. nIf I'm not wrong, the 'beat cycle' is another kind of aliasing that Nyquist sampling theorem haven't told us.
8 |# j' C R+ A* l1 |- E A. x, U8 {* u' Q p/ g
When the sample rate is close to 2x at the frequency. It will generate 'beat cycle'(or aliasing) no matter it's over or less than 2x. I will call it 'forward aliasing' and 'backward aliasing'. 9 q0 D6 T# C7 b/ K
& Q2 B8 x5 u$ N5 c* W
Use 44.1KHz sample rate to sample 20KHz freqency will generate a 2050 Hz aliasing. No matter if it is 'forward aliasing' or 'backward aliasing', that 2050Hz will be reconstructed. ) M) ~; N7 n% [
& O7 j( E+ X. l1 }" S
Maybe the solution is sampling the frequency over 4X sample rate, then the 'beat cycle' will be over 20Khz, that will be easy to use a low pass filter.[/quote] % p- ]- T" C+ p* k" I
$ c2 F0 p; ?3 n- o) @0 |% m; D1 S+ _The term "beat cycle" is one I "generated" to help me explain a concept. $ ? ^" b1 t- Z& `1 g& ?2 ]
& J- \9 A- n% x) A, m$ [Nyquist did not miss a thing. I think you did not understand what I said.
" j0 K0 j7 c0 j- L0 W: R& _I think you missed my explanation.
6 Q% E$ X3 T3 Z- j! T' S- t6 [
( |: N8 w% |1 vThe waves that you are looking at are BEFORE the reconstruction (before the filtering of energy that does not belong to the final outcome). ' `! F, _ |# j
% R; ~' R" x# H- L, X, }What you see as aliasing is not aliasing. The "modulating envelope" of the high frequency is not a signal, and it does not carry energy. & ~# l( Z2 P3 d0 ^
# G2 M$ L- I. K. W; AThat "modulating envelope" will disappear when you filter the signal with an anti imaging filter. The problem is that an anti imaging filter gets to be impractical when you are too close to Nyquist. That is why, as I already stated, we used to cut the bandwidth a little, to increase the margins and make filtering easier. , R2 K+ j0 Z0 [
. u3 h/ T0 ^, w' d5 u4 f" G3 [' E0 n
You said: "Maybe the solution is sampling the frequency over 4X sample rate, then the 'beat cycle' will be over 20Khz, that will be easy to use a low pass filter."
$ Y- y& S( k' @" k6 T+ U0 u- n9 H) [0 i. d/ L" U, M0 E
If you read what I said again, you will note that I stated that these days (and for the last 15 or more years) we do use up sampling, to make filtering easier.
% r& W ?. x7 O5 V9 y( Y* s/ Z& [: E- B
You suggested to go X4 up sampling, and these days, it is common place to go X64, X128, X256 and even X1024! The comment about X4 is confusing and may be misleading, because it suggests 192KHZ, whias nothing to do with the question here.
* P d7 U9 N1 M
^. u+ z( o% c6 |* m$ m9 mThe fact is: you DO NEED to include many samples of the past history.
5 s2 [5 p+ g# u! \/ E% FSay you give me 8 points, instead of 2 sample points per cycle. I can draw INFINITE different shapes that go through those 8 points.
* i( S) z$ K; P0 f$ I. T+ HOne can not figure the original wave shape of a cycle from 8 points. The 8 points only define that the wave goes through those 8 points. What happens between the 8 points? You need a lot more information, thus you need to relay on a lot of previous history. 6 Z' ]3 d0 x, S5 s( Y
_0 _/ D7 B' T; X2 y( h- XAnd more importantly, music is not made out of repetitive identical cycles. With real music, the wave shape is changing all the time, making the signal itself much more complex, yet by filtering away all the energy above Nyquist (both at the AD and DA), you can reconstruct the signal perfectly.
9 H2 e) V4 n5 f. T ~6 w$ d8 i. s: P% a" q/ d
The concept was discovered by Dr. Nyquist, a man with a Phd. in Physics and a lot of solid math background. He was working at Bell Labs, with other greats such as Shannon (the father of information theory). The concept is not a theory it is a theorem. There is a PROOF that if you filter bellow the sample rate at both ends (AD and DA), you will end up with a perfect replica of the original waveform. In addition, Nyquist theorem stood the test of time (over 80 years in many areas of technology).
. G% i9 U5 W. @9 f# |
. F, c' a- ?9 U j) iIt is difficult to construct a filter that will keep say 24049Hz intact but reject 22050Hz (1 cycle away). But when the PRACTICAL margin (from a filter standpoint) is met, there is no "beat cycle" and no "another kind of aliasing".
6 H2 |0 @/ S# q9 X( x E
7 k' B# X; I+ h: \3 z- T- V, {Regards # N+ E% ~7 S) b( e
Dan Lavry
6 X1 d5 v6 U* }) A5 __________________; U/ ]3 P6 e) G% m
Dan Lavry
; ^9 ^7 A" @0 w% u( h) O4 Z7 j! @- z0 t0 B' f
Back to top
( C* q9 V5 v3 }7 `6 U2 H/ c, e+ ~! R
$ z6 @* }+ Y, R5 z
lavrye
) b9 a+ G: n! ?& h$ T1 y% p8 ?" aSite Admin# w; J5 } k0 f3 g& c5 `$ ~
. I7 c! N. m7 i& V
! p1 M2 n8 W+ V; |' }- D* n! O- _Joined: 15 May 2005
7 `. z1 B J1 l7 SPosts: 462
6 }3 D- M8 [6 ^5 C: O% G) \7 }Location: Poulsbo, WA9 C5 d; f0 a: [ B3 t. W4 |6 m' }
Posted: Tue Jan 22, 2008 12:12 am Post subject: Re: 44.1Ksps sample rate enough for audio? 6 V: Q5 I, v8 C( l0 V. v
X* D1 C' h. w2 k1 G: r
--------------------------------------------------------------------------------: z3 z) K5 u$ n- V2 B
c8 [7 F ?9 D5 L* k
boB wrote:
9 n' ?$ v2 m Q6 r/ i$ H; {$ ?lavrye wrote:
" I: p [, O% k4 \7 l6 G, G, k% I% |3 r% }* i+ D# O) @
; V$ {/ \, l% M% }
You can not tell the wave shape of a cycle from 2 dots. To understand sampling, you need to take into account very many cycles prior to that "cycle with "2 dots". You will soon realize that the samples (dots) hit different values on each of the previous cycles. With data from many prior cycles, we can reconstruct the original sampled wave.
; `1 X; M, p: t1 W( e2 M
" Q5 @9 ]5 F5 X+ _3 A~~~~MORE STUFF~~~~ ~7 i) V3 ~! e; m' C" V
6 x- a1 y5 `4 t( ARegards 5 g# j6 k5 [; E7 K
Dan Lavry 4 h" M6 I/ }2 q: W) G3 \
, c2 a4 \: I/ \3 P t( h+ V
) `- S3 p s% {( \4 y. W5 _Maybe it's exactly the same thing, Dan, but I look at the re-make of the sine wave at close to FS/2 from two "dots" as reconstructing the sine wave because of the bandwidth limitation at FS/2 itself, rather than the "history" of the sampled waveform ...
# T2 P3 o( E/ A+ K% K* U) f+ Z& b5 K) A5 j/ h% Y
But I suppose the history IS because of the same bandwidth constraint ??? + u4 u7 H8 o @9 R5 m2 S
1 t9 @) L5 X/ |: K/ X, d3 ~Is this a valid thought, or another way of saying the same thing ?
, Q3 I" c e9 s0 _' Z" S' s O4 U) o# x2 M
Thanks, 9 y( R/ S L9 C- q( I
boB
$ I- g- E" h) f0 J; d5 x3 u7 U8 D7 t- E( `8 @9 ` m2 f$ x7 v( H
, x) R- v. S D3 U I3 qI am not sure I understood what you said. Can you elaborate? . j9 X7 g1 ?+ e8 U% C" {- T) L
6 c6 Z1 U0 G! \7 _0 YWhen you sample at exactly fs/2, you violate Nyquist, because the points at each cycle are at the same location relative to the cycle. 4 o K5 A. k, A$ _0 N, E' z: \
- a+ I: |7 |9 b; }$ n) K7 R
But when you sample a wave that is slightly below fs/2, look at previous cycles, and you can see that the dots are at a different points of each of the previous cycles. That is where the information is to enable the reconstruction of the waveform.
2 g i, Y9 G/ D7 n9 A% ~) a5 O, f% r+ |
I do not like talking about sine waves, because it simplifies matters too much. But OK, lets stay simple. Say you sample a sine wave at 22049Hz. Say the first sample is on a positive going zero crossing. The next one is very slightly earlier then the first negative zero crossing. The next sample is slightly earlier then the second positive zero crossing... and so on. If you take around 44100 cycles, and "lay them on top of each other" (think of a scope triggering on a positive zero crossing with a time base set to show one cycle time), you end up with a "solid looking sine wave" made out of 44100 points. In other words, the history of many previous cycles is what "fills in the gaps" between a single cycle made of only 2 dots.
! m5 s2 Q0 F; a* a, U# |! x) i% n1 J: K) p( @
Of course we do not keep 44100 cycles and "lay them on top of each other". But we do need the information from many previous cycles to "steer" the signal to fill in the values between the samples. Filtering is what does it for us. ' e( b0 T+ i% Y) \5 r# E6 V1 o
5 m" I2 X) s3 `" p
Regards # P4 K4 Q9 O5 R0 L, I$ h
Dan Lavry- {% F: V( ?7 [) | t/ m- q2 |: U
_________________7 B3 u$ |' Z( u" l$ A) n" C
Dan Lavry
5 F# U* ^7 s( G$ S5 I9 L: q, p, q- e9 `" y3 \4 [
Back to top 3 X2 R1 Q2 o% h" ^8 Z! `6 d
, q1 R q0 K E- M6 T3 [6 z
* j K4 r( x" w6 Y
boB/ y5 K0 T8 l& H: z9 ?: V$ D
2 o) w. N2 b4 x& C- M3 [% Z6 w8 W7 W
' t- U, F! r( d3 l2 `9 PJoined: 28 Aug 2006
% J" j. g- b2 A; _" q N1 tPosts: 34
( x$ u% e* Q9 C7 CLocation: Everett, WA
* z; o& R2 O% a. y" uPosted: Tue Jan 22, 2008 1:26 am Post subject: Re: 44.1Ksps sample rate enough for audio? - t+ u8 ` W$ x& S
" M8 J, _5 N: W& K: \$ R--------------------------------------------------------------------------------
3 J0 c8 r3 z# p: b( `+ } R* H
lavrye wrote:
5 `' n( p* c3 T: K6 g& V" _
2 t! a8 _ L# Y
2 H: O. O. G: V$ ^( x' T+ AI am not sure I understood what you said. Can you elaborate? $ S! F, u& P( b! t; h, n4 i2 _, T
p( x: _: n4 i. EWhen you sample at exactly fs/2, you violate Nyquist, because the points at each cycle are at the same location relative to the cycle. 3 Y) p$ F) y2 r) C: X
# W. ]- U5 p# y) l" _ T
Regards
p1 i% _9 l3 }9 C {+ M; iDan Lavry
, z. A. J) B! r7 Z
$ G$ P) F' E- x" b8 p* C
' E2 l0 q2 B- h! z1 dI did indeed mean somewhat below FS/2 and that the filtering was filling in-between the dots. And, Sinewave at this 20kHz frequency ( G3 Z' \2 ]! J4 f& m
area because that is all that can really be represented at this fraction of sample rate.. These 2 requirements make it much easier for me
: }. F& n' }! j# T G# ]: `. jto think about this stuff.
! u7 c9 z: |2 z( s) o
. I4 [. I+ W$ S7 ?I've thought that the bandwidth limit (reconstruction) was all that was necessary to "fill in" the rest of the waveform
1 e2 _. {+ ^6 s6 ubecause all of the energy of that sinewave is below FS/2. 5 b: W c, _) y( K8 t
' Q' A. \' w- m
Is the "history" coming from the filtering itself ? Does this have to do with the sinc function ?
& r" o5 W: W9 B5 t; Z1 q M4 I* \6 K* X' r2 Y! I
In the past, I have seen unfiltered digital samples suddenly turn into sin-waves when the horizontal time base was just right. Is this related ?
2 l2 S& g: K' F/ \( I! BI can't remember what kind of scope it was but I seem to remember that it was an older one.
2 `8 X D- b2 w) H* R& k9 G& s# M$ X! I: X- G/ I# ~& N
boB( s$ v r7 x; n
_________________
4 s% b1 o1 o, ]4 yK7IQ ! u+ ]% X$ ]3 z+ \* ]+ A; n
$ I8 W; t5 H. B0 S
Back to top
( k6 U, M* F$ S. Z2 P9 {+ @+ a6 s% {& C* y ?
1 T5 S5 K6 a/ P, ^. u
zifzhu
5 U/ V4 p3 Z- _% A [
" N3 P/ @! d# |7 G
3 F$ `2 ]$ t$ G% J8 l4 n, Y* R
3 j, M; ?* h AJoined: 26 Dec 2007: j! k& ~9 ^5 B
Posts: 4
) x6 e4 B* m6 p: j: I) W" aLocation: Suzhou, China
; Y2 A' u0 W/ Z% d4 gPosted: Tue Jan 22, 2008 4:53 am Post subject: * r+ w u1 u% n" T j" u
; a: C$ D% g, |2 \
--------------------------------------------------------------------------------
& m2 u/ J" D# ]4 q! a' ?: g( ~3 \9 W7 e% _4 C; Y
My understanding about 'history' is just increase the sample count and period, so that will be more accurate and not lost 'something'. 1 r; L$ {9 U4 S( _" B
/ N5 g" H8 B: i# ]; q
I've create a excel file that demo the 44.1Ksps sample 20KHz sine wave. The amplitude of the frequency is just increasing and decreasing at a 'beat cycle'. It's not a simple sine wave not, but a harmonic. Hard to believe this wave will be reconstructed after filtering. At least, the total amplitude is decreased. The harmonic contains the 'beat cycle'.
6 y& N* ~ X+ q4 X$ W2 W! ~. s
+ s( C, |7 n5 ~: V. y4 N+ w![]()
& t Z* ]& C( o& ^5 s6 s% I6 W/ f7 v+ i) S( h
_________________
" y# ^, ~2 I+ o% \: f
3 I+ n$ U' z3 p8 ?* J! G( q( j[ 本帖最后由 zifzhu 于 2008-1-22 12:59 编辑 ] |
|