POPPUR爱换

 找回密码
 注册

QQ登录

只需一步,快速开始

手机号码,快捷登录

搜索
查看: 5333|回复: 18
打印 上一主题 下一主题

惊天大骗局,44.1KHz的采样率真的够吗? -------------- 后续!

[复制链接]
跳转到指定楼层
1#
发表于 2008-1-21 23:24 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
帖子又沉掉了,没办法,只能再开一个! o+ @1 {6 w9 o0 n0 e8 R5 `( w
http://we.pcinlife.com/thread-867632-1-1.html3 [; x' r' k' N
其实到最后也没得到什么结果,就这么不了了之。
( z4 q' k" S- o. l, c/ @- u6 G+ Z# ~) J5 ?% H
HerculesVR老兄推荐我去问一下Dan Lavry。于是我去他论坛上发了个帖子,这些天几乎忘了这回事了,刚回去一看,似乎有些进展。先贴出来大家参考一下。后续有什么进展我会跟踪报道B)  期望我这个无聊的疑问会有一个好的结局。& K  N+ k9 S/ S6 p! u
0 p; `5 C2 C4 ?  i, V* |( ?) c5 I
zifzhu
+ n, _% e4 U1 i
. y5 K( D2 y' T9 L. H/ a" g- J( p+ M$ e* t
, l, K) X1 W& N. g
Joined: 26 Dec 2007
$ O. T! O. S4 I! D6 D' f( h' `, rPosts: 34 Y: h, z5 L" O0 L
Location: Suzhou, China

: s8 o7 O1 K  b! k0 L* w) C1 `) }. l9 S
Posted: Wed Dec 26, 2007 7:37 am
7 i0 s8 B% f2 y# l$ O8 H   Post subject: 44.1Ksps sample rate enough for audio?
We are always be told 2x sample rate is enough. It seems 44.1Ksps for 20KHz audio is enough. ' L7 L, K4 C5 f+ B
- d* [* Q$ s4 t* e5 e# \5 `8 i

7 y& k) Q; z( ~8 L- D: \! W$ QAs this diagram, 2x sample rate regenrate the source signal perfectly. But the question is, it is in the best case. It samples at the peak luckly, but if it sample 90 degree before or after. Then it will simple samples all zero.   u$ M+ W: N- d) H

1 a7 f/ K; D/ P& {( _5 {It seems 8x sample is the minimium reqirement for sample for hi-fi. Thus 192Ksps is the minimum spec to record. % E' l' b' F( ]4 q! p
1 s2 ^9 a7 i) B8 f/ [2 I
Am I wrong at this point?( O; X: c$ Q/ T2 e1 W4 U/ O+ r
_________________& V4 {$ S) f' Q! f; M3 Y& Z

( U/ {/ L4 d  S/ t. m; N- a8 F* G
Back to top

! r$ ?* S0 y( \2 |1 ]# h7 f, g+ B5 D

2 z3 |: \2 r8 P' b
Brad Johnson
0 l! G8 m- U! j  K& ^6 {1 S/ `9 L$ I8 o* d$ v
8 z1 y. ]4 P  n0 t. {3 O
. n+ g! L" }/ T+ W$ T. [- \
Joined: 08 Aug 2006
9 D  v# u4 q& [1 K/ E$ CPosts: 126) j$ f( A* P7 }, e/ c! q

" j, P' G) W: d0 [; ~" x( d
Posted: Wed Dec 26, 2007 6:48 pm0 l' W1 g% j" s7 }
   Post subject:
There is a "white paper" under the SUPPORT tab of the Lavry Engineering website that is a very complete discussion of this subject.
3 M. [5 X- o4 g4 B& e- V& U$ Y  N
, `( B9 ^' z( e2 e) dIt is entitled "Sampling Theory." 0 u2 c9 _4 k7 E0 H2 Q1 R
/ j" r. Y3 o* b" H- z
http://www.lavryengineering.com/documents/Sampling_Theory.pdf
0 T, p) N. X& F# S  i# X/ H9 Y1 f  M" X- r& A. s$ }& t  V+ k0 V5 f
There are many sources of inaccuracy in real-world conversion systems that can cause an apparent loss of high frequency information, and they are not the sample rate. There are also additional trade-offs which vary with the sample rate that occur in Digital Signal Processing; in both the processing that happens in contemporary converter technology and what commonly occurs in mixing or mastering after a signal is digitized . So in the end "higher" does not alway result in "better." $ `; ]& A9 M8 n9 j* e  n( n

( r; ?: M3 s" w* k7 EA recording made accurately at 44.1 kHz will always sound better (or more "analog") than an inaccurate recording made at 96 or 192 kHz. In many cases this would be mostly due to jitter in the AD and DA conversion, as versus the sample rate. But even in these cases, the degradation is almost certain to be caused by more than just this one source of inaccuracy.
- m, l1 j3 \& k
  ]& t3 [: D- q+ ~; p; bBrad Johnson
5 k1 ~" k& }9 ^Lavry Tech Support
  w1 [! @' C6 A$ y9 u& w- Q8 y
Back to top

5 u. d5 g& n, B/ l/ v" Y$ `2 T* c; W! m% _2 w3 ]$ Y8 @! n) ~

2 \6 Z7 N! e) W& _  w4 |5 t: q3 C9 a
cj
" k% p5 N) F# t, v& _8 k% x* S/ U6 ?$ i4 T$ O0 N

( h6 |: G- v! ?2 ]( _" U. `9 H4 H/ E' N; B, m
Joined: 29 Aug 2005
$ E& k) N( q2 c8 X; _6 G$ SPosts: 57! g4 [) }, b1 _# U, E; k
Location: Finland
  t* g4 r( J! L% x
Posted: Fri Jan 04, 2008 7:35 pm: B1 h: {8 {5 O/ X+ |" J
   Post subject: Re: 44.1Ksps sample rate enough for audio?
zifzhu, 8 v/ R' D2 D6 j1 L1 M, w

9 C2 }+ X$ H0 i: JExactly half the fs is kind of a "special case". Check the "critical frequency" part here:
# B; e6 V1 F3 \http://en.wikipedia.org/wiki/Nyquist%E2%80%93Shannon_sampling_theorem
9 W% f- J" W% b/ d% N/ ]
6 C0 t2 s3 F& a& s$ k6 _) rC.J., Finland2 Q3 F  i$ p3 M3 J
_________________; s- h" ]$ J( F" v8 c! E5 w
"Science is a way for us to not fool ourselves." - Richard Feynman
7 H  M0 C; G5 ]( ?( q& {* s+ D" `
Back to top
' e$ T& E5 f4 p! {' B
+ ]  r. a6 O9 J/ q0 L: B2 }( Q! l0 x
: e5 }3 Z% j- A' l0 P3 Z2 D. c
2#
 楼主| 发表于 2008-1-21 23:25 | 只看该作者
Mike Derrick
/ S. q4 J9 o/ q3 F0 D
0 X% }  j! f2 {2 e* \5 A. h$ m7 r  u, k  d% k4 j( _  _
+ l& ]* H: i! F& r  o" {" q
Joined: 29 Jun 2005
3 C! R2 ?( K0 q9 R4 m( nPosts: 16. E2 m$ @0 I# t; b1 }$ W
Location: White Rock, BC

3 c" ?  }5 |5 `! b+ t
Posted: Sat Jan 05, 2008 8:43 pm
$ ^: O# o& T7 V7 s. E8 ^. h1 F$ o' N   Post subject:

! C; v9 k  L! D7 X* v
Brad Johnson wrote:
...There are also additional trade-offs which vary with the sample rate that occur in Digital Signal Processing; in both the processing that happens in contemporary converter technology and what commonly occurs in mixing or mastering after a signal is digitized . So in the end "higher" does not alway result in "better." 0 Y/ H6 `: }5 s  Q* j) W) N
9 ~5 H9 v8 z1 J  f6 l- A4 ^
A recording made accurately at 44.1 kHz will always sound better (or more "analog") than an inaccurate recording made at 96 or 192 kHz. In many cases this would be mostly due to jitter in the AD and DA conversion, as versus the sample rate. But even in these cases, the degradation is almost certain to be caused by more than just this one source of inaccuracy. 5 R) n+ k, A9 z( x
1 R! c6 l- b. T/ Y& I7 R1 \
Brad Johnson
9 E- n- s1 o* z  VLavry Tech Support

& m0 z, Z- U/ O
' R* {/ P$ M% Q. J0 m5 B. _Hi Brad, & g  t, {( u+ ^8 @

3 e, D0 k0 n0 i! Z- @8 K# ^: ]+ ]2 GYou mentioned additional trade-offs in digital signal processing. Specifically could you elaborate on what those trade-offs might be within the context of the mixing?
8 e# D2 ?7 B2 J' [
# k# w; o7 V! oCan recording at 96k be as accurate as recording at 44.1k?
1 q0 k& [8 A. {) v" C* Q5 A+ V9 o: o4 N! f
Mike
+ ]" T; p8 w6 M5 |2 K, x_________________! t" w% ~& P0 g3 @' [
Mike Derrick
: k6 v7 I2 y2 j+ D+ |
Back to top

+ T1 ^0 g- m$ D8 g0 R+ t
3 Z* M7 N+ G4 l/ e3 g' Z; ~; t5 C% a
Brad Johnson; R- c3 p5 e. g4 g# _# J* E* @

: H' b; G" }7 \8 A5 x+ z7 ^# {6 F0 E5 G' e1 G/ R

6 Z/ \: ]6 h# c; G0 F* \Joined: 08 Aug 2006
: a8 K0 w9 h  L" S! p1 o* g( t+ B8 [Posts: 126
8 t5 n- n6 f# ?8 i! A9 {3 ]

& {1 x: u  Y; }/ C( l8 j
Posted: Tue Jan 08, 2008 11:31 pm
, Z* w+ j! e* N   Post subject:
The trade-offs that can occur in processing after a signal is digitized at higher sample rates are mainly in the amount of processing power needed for low audio frequency EQ or filtering. Since all real-world systems have a finite amount of processing power, this can lead to compromise in the accuracy of the processing, or taking much more time for non-real time processing. So this is probably not as much of a consideration as the compromise that happen in the actual conversion. # n4 e. W8 A6 t" y" U! s

- D7 F' [4 ^" m+ V% uIn context of the original posting regarding the accuracy of the digitizing of high audio frequency signal- a 96kHz recording can be as accurate as a 44.1 kHz recording (if not more accurate). The point is that the sample rate is only one of many factors that can affect the accuracy of a digital recording. Present day converter technology can be very accurate at 88.2 or 96 kHz. 4 R/ }" y! J& Q) R* b

2 e6 t1 @- ^8 z" L) }In addition to possible trade-off's in noise-shaping that is typical of present-day over-sampling conversion, there are many possible "analog" sources of inaccuracy in converter circuitry. Effects such as stray capacitance in sample and hold circuitry or settling time of amplifiers become more significant as the sample rate is increased. This does not mean that 96 kHz conversion will alway be less accurate than 44.1, because it is the total effect of all of the factors involved. Careful design can minimize these effects, but if you consider that present-day converters typically are already sampling at many times the output sample rate , doubling or quadrupling these very high sample rates to achieve output sample rates higher than 96 kHz can make these "analog" problems that much more significant. If you don't raise the internal sample rate and output a higher sample rate, you lower the "bit accuracy" of the digital audio output, so you would compromise the quality of the signal across the audio spectrum. 4 I) ?0 B/ C5 n) Z: U* }! K& y

/ w# p! `5 {- g% u. c8 LThe question remains- what is the best over-all performance? Do we really want to make accurate recordings of supersonic frequencies at the expense of sound (music) we can hear?
- R( d8 d" r8 m) h
$ |) C: g4 H! H) X3 xBrad Johnson
4 c: r2 z! o( Q) ^, V
Back to top
1 O- w5 R$ J, `% T2 A) W
, c4 ]( ]" |& L9 f
8 S3 k# z' r4 T- ?
zifzhu
- {' X& l& ^) T) Z2 v: ]% j. w  C# a4 f9 d$ k5 c

2 v: L+ V& @5 N
! D' m+ e+ {+ x+ _, d) vJoined: 26 Dec 2007) s$ q" z) ]& h" F& u
Posts: 3$ \# c! w# D5 X* j* t
Location: Suzhou, China
( `  ^% h$ x7 [2 |, R& }$ M5 _
Posted: Mon Jan 14, 2008 5:21 am
( x, h$ R, D. T& e6 S   Post subject: Re: 44.1Ksps sample rate enough for audio?

: t* g: o- ]% w! g# J, p
cj wrote:
zifzhu,
3 k1 ~- n" x) l# I
; }/ K9 `& ]7 L7 a# w7 q% v: A. {Exactly half the fs is kind of a "special case". Check the "critical frequency" part here:
" u8 ?2 ~, K* n7 ~$ Thttp://en.wikipedia.org/wiki/Nyquist%E2%80%93Shannon_sampling_theorem 7 m" Y! j7 S2 z
1 i0 E0 h& k& A6 ]. l: ]3 M" P
C.J., Finland
" l, G7 l( v9 M) I  M8 p8 F
Unable to open that link. Will try it later. # y0 T) u( b3 S/ r$ z' S4 A
; l4 h$ }: v: P. `2 P5 k% q
If exactly half the fs is not able to be sampled correctly. The next question is 'How much?'. 0.000000001% more? Or 100% more? There must be a value.
9 J: Z: l* y/ L0 ~_________________
5 j" x& K8 {. ]1 ?! \
; Y3 h# a% j$ k+ G" t0 ?
Back to top

2 a: q$ t# c# h6 |" K. s
% E7 b! s6 v# D2 [2 }& J0 k& j) d) Y1 g4 @$ q, I" e
回复 支持 反对

使用道具 举报

3#
 楼主| 发表于 2008-1-21 23:26 | 只看该作者
lavrye% l, i! K9 ~: H8 Q, Z& `' E
Site Admin& L) J( Y2 `! j; v3 S" h/ O# o0 `$ j
9 k* n0 T, o+ k8 w
, }! C% d- C/ [7 w# K- c( h
Joined: 15 May 2005
, l: w. a( m# }% L, l$ ^2 KPosts: 458
3 X& ^9 [+ ]1 T4 T. o$ hLocation: Poulsbo, WA9 {: H* C; Y  r8 x  |
Posted: Tue Jan 15, 2008 7:58 pm
  Y4 Q) {, ^% G+ I( A& I8 _8 y   Post subject: Re: 44.1Ksps sample rate enough for audio?
. p% ~/ O. ?  Z0 h1 ^" Y% Z' ]$ ~, n
zifzhu wrote:
We are always be told 2x sample rate is enough. It seems 44.1Ksps for 20KHz audio is enough.
* ?& N  A  h7 P- Z4 j , f  U5 f- K3 _  L5 |) q. s

# N7 J" H% I% k- w# D" w4 S5 a) GAs this diagram, 2x sample rate regenrate the source signal perfectly. But the question is, it is in the best case. It samples at the peak luckly, but if it sample 90 degree before or after. Then it will simple samples all zero. ( u+ W- h' S( V& Q
& |# x7 h! E) ^$ d0 z' w+ h0 u+ H
It seems 8x sample is the minimium reqirement for sample for hi-fi. Thus 192Ksps is the minimum spec to record.   F7 T- z4 B. V  j

1 Y9 i9 R" r% g. F# F9 XAm I wrong at this point?
9 e( [) ]) x: o$ g) i

. u, a0 g+ a! K; m9 b# y1 kThe sampling theorem states that you have to sample FASTER then twice the highest harmonic. " O, k5 J; h( t

4 g0 j# O% P# L6 ^8 xMany people look at a single cycle near Nyquist (say 20KHz at 44.1KHz), and wonder how it is possible to reconstruct the wave from 2 "dots". That is a wrong view of sampling. 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. $ T- Q# H4 M1 m. C5 x( E# Q
: s; [$ S7 U) N( T8 N
It is correct to note that more recent cycles (nearer in time to the point of interest) contribute more to the reconstruction. It is correct to note that what happened many cycles ago has less influence on the outcome. But one can not ignore previous samples, and the more "history" you take into account, the better the reconstructed wave. Of course there are practical limits to how far back we look.   I( H- e1 `, j. [7 x9 x
8 d% C2 o* \6 a: K; _
At low audio frequencies, there are many samples in each audio cycle. ( E/ |" K* E2 H9 T& z! ]
As you increase the audio frequency, we get to a point where there are only 2 dots per cycle, and things get visually confusing for those that view only one cycle. Further examination shows that the samples "walk" across the audio cycles. Lets take some cycle as a "starting point". Many cycles later, we walk far enough to get back to (or pass by) a picture that looks like the "starting point". I will call that a "beat cycle". That beat cycle is repetitive (for a sine wave or any periodic input). 0 [0 O3 P7 {5 s  n6 j4 X8 W0 j& b

7 j! c8 d0 l( w/ c1 W7 \, m4 d; B' OIf we wanted to reconstruct say 22049Hz with 44.1KHz (Nyquist at 22050Hz), the amount of data we would need to take (to reach even a crude approximation) would be huge. That is because the "distance" between 22049 and Nyquist is only 1HZ. The 1Hz difference has a "beat cycle" time of 1 second. That means that at 44.1KHz, a complete "beat cycle" consists of 44100 samples. That is a lot of data to collect and filter...
/ N- r5 _7 `$ `" I! ]4 }& B& s+ z5 U
But say we want to cleanly reproduce 20KHz with the same Nyquist, then the frequency difference is 2050Hz whias a "beat cycle time" of only 487usec. It takes only 21 samples of data to go through one "beat cycle".
: j2 L) d2 w8 X" m! s8 W/ Q) ^# \! h/ G. M5 }
How many "beat cycle times" does one need for good reconstruction? This needs to be quantified, but clearly reconstruction of 22049Hz with 44.1KHz sampling is impractical! On the other hand, reconstruction of say 20KHz to very high quality is very practical.
' _( ~3 A/ s8 K' }5 p9 T6 I/ e% c' t/ G" I
So instead of trying to do conversion up against Nyquist, we back of slightly and pre filter the signals so that the converter would not need to try for what is virtually impossible. Of course, with newer converters (over sampling front end) the Nyquist is higher by huge factors (in the MHz range) and the issue AD converter anti aliasing filters belong to the past (it was an issue until the early 1990'). The same holds for DA reconstruction which is filtered after a lot of up sampling.
' F1 Q3 f# i0 ]/ ]7 @6 X# G! D7 @8 ^% D; J4 O
The issue does exist for down sampling (decimation) to say 44.1KHz and that is a digital filter FIR's or IIR's. To get the FIRs closer to Nyquist, one needs longer filters, and to decimate to say 22,049KHz we would need probably over a hundred thousands samples and a huge compute engine. It is not worth it. If you are willing to settle for say 1000Hz below Nyquist, the problem becomes vastly easier (a beat cycle time is 1000 times shorter). You can do practicable decimation with reasonable amount of data. Similarly doing IIR has its limitations as well. % l" }) O1 d! C# i  S* {* r

; y6 C2 X! ]: Y; c  cRegards
/ U2 `2 `( w8 B: T. ^/ |! {Dan Lavry
/ x) l3 b+ F; R! s_________________% t, U- ~: x. P/ N) A* \
Dan Lavry
回复 支持 反对

使用道具 举报

4#
 楼主| 发表于 2008-1-21 23:27 | 只看该作者
zifzhu
5 t2 @5 h: Q7 H$ M6 M( A  F3 A+ x! q1 f

  t2 D- R- i- k; A' N+ E9 q0 y, Y9 I- N% I0 w6 p) L+ l( U( K" C& ?
Joined: 26 Dec 2007$ ]) C7 a; d. M
Posts: 37 c, [& @3 V# }- J% I$ |
Location: Suzhou, China
- X* p  N7 a( y- O" d& J3 {
Posted: Mon Jan 21, 2008 3:12 pm
) e  e6 e; y; m/ V8 N/ j1 w4 t   Post subject: Re: 44.1Ksps sample rate enough for audio?

' A' P: D" u2 S! ?# c
lavrye wrote:
zifzhu wrote:
We are always be told 2x sample rate is enough. It seems 44.1Ksps for 20KHz audio is enough. 0 U  ?8 v) }8 Q! N# U% L

: ]5 n3 h! \# o2 D% g. ?" z) m6 j% X0 j1 B$ h$ T, L' u6 Z5 S
As this diagram, 2x sample rate regenrate the source signal perfectly. But the question is, it is in the best case. It samples at the peak luckly, but if it sample 90 degree before or after. Then it will simple samples all zero.
1 M+ y) b/ n; P/ ]* N- v
) i# h: r# h6 L+ p) g: WIt seems 8x sample is the minimium reqirement for sample for hi-fi. Thus 192Ksps is the minimum spec to record.
! D& x2 C- J% G( j5 r2 f( Y0 K2 ?8 M; _1 T* s5 [
Am I wrong at this point?
5 I' A: v4 H* ~$ j: v: `+ h
" G. m& z  s: u7 N8 N
The sampling theorem states that you have to sample FASTER then twice the highest harmonic. . y. n6 {  V# n

5 J5 ~% D+ Y' |1 _Many people look at a single cycle near Nyquist (say 20KHz at 44.1KHz), and wonder how it is possible to reconstruct the wave from 2 "dots". That is a wrong view of sampling. 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. - W# O) c3 f. ?; ~
; M3 n3 K) W' H
It is correct to note that more recent cycles (nearer in time to the point of interest) contribute more to the reconstruction. It is correct to note that what happened many cycles ago has less influence on the outcome. But one can not ignore previous samples, and the more "history" you take into account, the better the reconstructed wave. Of course there are practical limits to how far back we look. , ]" n, v3 q4 P8 W0 F
0 ]% m8 W( F0 t7 A
At low audio frequencies, there are many samples in each audio cycle.
/ m; S5 a( N0 z  Q! N, m+ FAs you increase the audio frequency, we get to a point where there are only 2 dots per cycle, and things get visually confusing for those that view only one cycle. Further examination shows that the samples "walk" across the audio cycles. Lets take some cycle as a "starting point". Many cycles later, we walk far enough to get back to (or pass by) a picture that looks like the "starting point". I will call that a "beat cycle". That beat cycle is repetitive (for a sine wave or any periodic input). ; _% r5 Z- E& g  @* z  {

- q8 Z. {$ {# K# M& u  |If we wanted to reconstruct say 22049Hz with 44.1KHz (Nyquist at 22050Hz), the amount of data we would need to take (to reach even a crude approximation) would be huge. That is because the "distance" between 22049 and Nyquist is only 1HZ. The 1Hz difference has a "beat cycle" time of 1 second. That means that at 44.1KHz, a complete "beat cycle" consists of 44100 samples. That is a lot of data to collect and filter... ) P2 \1 E  J+ C, J: g* q' a
5 ~0 g: T# M+ b0 ~
But say we want to cleanly reproduce 20KHz with the same Nyquist, then the frequency difference is 2050Hz whias a "beat cycle time" of only 487usec. It takes only 21 samples of data to go through one "beat cycle".
/ I* Y  }, d0 l0 [# C
. a+ w) H8 x9 u8 XHow many "beat cycle times" does one need for good reconstruction? This needs to be quantified, but clearly reconstruction of 22049Hz with 44.1KHz sampling is impractical! On the other hand, reconstruction of say 20KHz to very high quality is very practical. 8 k2 l! u  ^) S2 r& c) _& M
; w+ M# |' y9 i2 V* o) O
So instead of trying to do conversion up against Nyquist, we back of slightly and pre filter the signals so that the converter would not need to try for what is virtually impossible. Of course, with newer converters (over sampling front end) the Nyquist is higher by huge factors (in the MHz range) and the issue AD converter anti aliasing filters belong to the past (it was an issue until the early 1990'). The same holds for DA reconstruction which is filtered after a lot of up sampling.
% R) x# z. W3 u) `* |4 x# j# d/ O* m* u. _+ Z' d8 m
The issue does exist for down sampling (decimation) to say 44.1KHz and that is a digital filter FIR's or IIR's. To get the FIRs closer to Nyquist, one needs longer filters, and to decimate to say 22,049KHz we would need probably over a hundred thousands samples and a huge compute engine. It is not worth it. If you are willing to settle for say 1000Hz below Nyquist, the problem becomes vastly easier (a beat cycle time is 1000 times shorter). You can do practicable decimation with reasonable amount of data. Similarly doing IIR has its limitations as well.
! P. ]5 ]% ]. v5 a% S+ c' V% t7 u& S! f# m
Regards
8 K+ g- Y% p7 w% c' WDan Lavry

+ K  R2 ?2 a+ {If I'm not wrong, the 'beat cycle' is another kind of aliasing that Nyquist sampling theorem haven't told us. 9 p% U% M6 W- A1 n! i6 w! m
; v+ p) M* O) O7 `6 |( Z; w3 H
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 o' ^3 T* L0 [4 T  v5 u

9 N: Q/ t  c5 g. Y6 S) CUse 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. ) ]1 S  U. R$ i. |
9 ~/ I+ x8 Z& G/ U4 ?; E
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.
1 X! A9 j) Z/ o; z* h. Y_________________
. Q8 y) A( q/ Q
回复 支持 反对

使用道具 举报

5#
发表于 2008-1-22 08:33 | 只看该作者
:funk: 又来
回复 支持 反对

使用道具 举报

6#
发表于 2008-1-22 08:38 | 只看该作者
翻译一下吧~~~
回复 支持 反对

使用道具 举报

7#
发表于 2008-1-22 09:00 | 只看该作者
楼主其实很无聊...
; u7 @& O4 H9 Y0 l
- s# y: l8 W* Y! s3 r别人写了那么多东西, 楼主就来个标题"惊天大骗局", 说的是自己的这个帖子?
回复 支持 反对

使用道具 举报

8#
发表于 2008-1-22 09:17 | 只看该作者
:sweatingbullets:  太长没耐心看完。。。。
回复 支持 反对

使用道具 举报

9#
发表于 2008-1-22 09:49 | 只看该作者
楼主还是没明白,不是说得很清楚了吗
- D& H* s5 H7 W6 e$ K1 @从你上个帖子里发的图看,你似乎错误的认为采样后重建信号是简单的利用单周期采样结果,但实际上,有限频域信号在使用Nyquist采样所得进行重建时,纯理论的说法是频域有限时域无限,因此可以利用无限多的采样点重建原始输入,当采样频率->2*Fs,采样点数目->无穷,工程上则是在能够接受的范围(数据处理能力、存储器能力、成本等)内,利用有限采样点进行重建(但绝少会出现单周期)# b$ ~: D6 j* D0 J) K! b# D
PS:以上采样实际上均为采样后的量化数据
回复 支持 反对

使用道具 举报

10#
 楼主| 发表于 2008-1-22 12:56 | 只看该作者
boB/ X: {" o& C! _0 y4 k

* ^: F5 ~$ g3 [/ p- C+ |1 L
1 D" I0 N* D' h# q8 K& X: S0 q! D* F, d  {% i
Joined: 28 Aug 2006
" A7 Z. J" p9 O( u- z+ o# xPosts: 34
' O  N5 h3 n1 y. lLocation: Everett, WA/ e) Z0 i  j( W. l- s$ `; j
Posted: Mon Jan 21, 2008 8:45 pm    Post subject: Re: 44.1Ksps sample rate enough for audio?   
& l4 w: T. L4 _, U1 v) f. w: W% b5 ]& u: P
--------------------------------------------------------------------------------& ~2 t' @9 M4 M0 c1 Y# f

  E# T7 u/ a' P! J: N. y( Ylavrye wrote:
$ u) r+ i% x- n4 ?3 {
" i% H+ v' a$ g9 `/ b" v
# G9 i$ b  p7 ]. x6 ^" yYou 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.   q. P, m, ~8 u. ]8 v9 a2 t  g2 v! R
; W: P7 E0 W5 T0 q  k0 r( k6 |
~~~~MORE STUFF~~~~
* h! J7 V4 f# F8 p" X" b1 G( m$ y/ g
Regards # c$ B; j, L  p
Dan Lavry
) V4 E" y+ N- _( J" h
2 Y) o0 w/ C3 v  c8 Z) h  l" M3 Z5 y( S3 u7 m; I3 k, ~+ Q
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 ... ( ^, m% K  a9 U' P

# D' }, z. R% i3 C5 t+ tBut I suppose the history IS because of the same bandwidth constraint ???
7 G) F) d$ t: }5 \
1 x8 Y2 [  l  j1 G4 X. {7 MIs this a valid thought, or another way of saying the same thing ?
! R# `' ^& w- F  X+ ]+ a) Z
4 V; b" [% `  c) g0 `& ^Thanks,
' R0 M* H+ ~/ v6 A" o& YboB
6 f8 `' ?, Y4 H_________________
" B8 X; {/ R! s9 U/ ^+ m/ O! D0 {) VK7IQ
; _: F  J2 D7 }( I+ }# Z) A
7 j% t$ N  m5 f4 D% ^/ B6 f; lBack to top        
8 I; T5 d/ {+ t4 e
& v1 j; O3 r4 e  \% _" C9 _. g: \, U- u% |9 m% m0 ~1 q) u
lavrye* U) D. |" W& q' B) @
Site Admin2 a* @& f2 S& ]' Q
# |0 I. h3 }( E% ?3 }/ r
5 M% ?, f% b" I8 d3 L, r$ a
Joined: 15 May 2005' f3 F" |; h/ ^, r  Z' g& w
Posts: 462
- O$ V: k4 ]% s' l4 l: HLocation: Poulsbo, WA
* w( K' {/ t: J( `; JPosted: Mon Jan 21, 2008 11:53 pm    Post subject: Re: 44.1Ksps sample rate enough for audio?   3 [7 I2 R" \  z- x* R
8 |0 m7 l/ |. ^" u% K( x
--------------------------------------------------------------------------------5 g" G+ e- Q8 D9 R2 b6 l' {
6 v2 p1 v8 D1 j) N% k7 a# f
If I'm not wrong, the 'beat cycle' is another kind of aliasing that Nyquist sampling theorem haven't told us. 0 v, ^7 C3 e0 `4 g; D2 r

4 B. \$ @8 ^' ]8 MWhen 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'.
" C4 u- I& Q9 M( W
0 a' o. y# i( \- O! |1 dUse 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. , c5 J' }) j) v% ]$ o
8 N' l& @3 [3 |: e
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]
6 t5 P. J! g- }% }* `; |6 ^8 e$ b( |
7 n' i' @$ w. ?" Y6 _# NThe term "beat cycle" is one I "generated" to help me explain a concept.
- u2 l* r. ~4 s/ `0 o" ^4 {
7 J. n: q8 \4 I$ o: D. SNyquist did not miss a thing. I think you did not understand what I said. 6 z* R& x! E# o1 u2 n
I think you missed my explanation.
3 H$ D9 V4 [7 }; Q8 }& M: @9 z5 V) r7 @' H
The waves that you are looking at are BEFORE the reconstruction (before the filtering of energy that does not belong to the final outcome).
+ f' o* M1 Y9 D1 q! H, e3 B4 k7 E1 p/ ~# G$ O
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.
  k; |( r2 q. @6 ^- V8 H/ n5 g3 J
That "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. : k9 R! D/ D) |) C5 a5 h7 ?5 W

- q( e) `* \; q( ~( nYou 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." # ]6 J6 b' J, |1 m6 ~# O

! v# Z% s/ U( ]( M. K: C  LIf 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. ; t8 ]% [) }. h- w5 V5 {

. z9 g* X) [8 u5 E  yYou 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. 9 @6 g% k! E' w0 s* K# C
, Z* \+ {  d& Z7 Q% F  D" `
The fact is: you DO NEED to include many samples of the past history. ( K- G, ]6 _# ?
Say you give me 8 points, instead of 2 sample points per cycle. I can draw INFINITE different shapes that go through those 8 points.
! u  |: N# Q0 w( z5 x  z9 }One 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. . @  _- Z$ J0 X& Y5 v7 p( Z$ g
5 Y3 O; B0 ?5 v6 U4 q! |( ~
And 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.
+ d) T) v: N9 o: F) _' c) a. }# v
$ B! H' K, c& b! iThe 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).   p- T" T8 w- \/ ~

. s" a6 ^" o. S# W! ?It 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".
+ }  ~, F- J7 e+ A! \9 b6 ~1 c/ a5 z0 P* y$ h* Y
Regards + e" B7 S5 w; x4 X8 z
Dan Lavry& f3 o: y5 N8 e  L
_________________/ E6 E6 I. A9 z
Dan Lavry / y' u) t$ _$ {2 B
8 I$ t/ Y& p9 d4 M5 O1 y; V
Back to top        ! a5 \& p$ ^4 L' r6 R$ u: E

0 o. E: P) e& |7 B" `% F" X2 w
/ |7 W* d: `0 B8 N/ K$ mlavrye# F$ f2 F. ]2 a; R
Site Admin
& @7 h" d. O' k/ J3 j# ^/ Y
6 z( a1 I, X2 k  i: k% T! a
+ E( l, I. T* X( aJoined: 15 May 2005
, B! N6 {) J% ?6 ~+ pPosts: 4627 g# `- ?: L% L6 }& [* O: ?" `2 \6 d3 B
Location: Poulsbo, WA, |4 K4 z9 z3 F: r" l& g$ x6 f
Posted: Tue Jan 22, 2008 12:12 am    Post subject: Re: 44.1Ksps sample rate enough for audio?   
4 f# ]) ~( ]' b9 {  m. R0 I7 W9 X
" _/ Y. w! p/ H. C# C" }--------------------------------------------------------------------------------
' B. H% H5 Z7 T3 U6 G. N3 I) j& b4 ?; ?2 X: r* `+ E0 d
boB wrote:
3 L& D' ^9 J" E; Y6 X0 Vlavrye wrote:
8 _3 p0 B6 A% f
4 Y/ q6 ~/ a! Y
( ~' g: K0 j! }) e7 J; `8 n" gYou 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. $ X' y% [8 H2 O( x

8 b4 U7 [6 x, l% a~~~~MORE STUFF~~~~
% X1 R* `) ?: A3 a: j( C9 E7 t" |; b7 m$ z; E
Regards
( N8 I. T# m* a. Y0 i* D: nDan Lavry
- s7 y) B' H+ ?* j& o5 a- Y
) K+ L( [9 v$ U( C$ q2 ^$ D! @1 l* ^5 i7 Z
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 ...
; k  [3 o0 }3 m/ _; k( g
; Q; A( X3 g  s6 RBut I suppose the history IS because of the same bandwidth constraint ??? * M, L+ \6 T% o8 C: Q4 y

3 A: r8 V, w2 [" ?' u* PIs this a valid thought, or another way of saying the same thing ? / I. b( E) p  y$ J- e  K

5 \2 P' N. B& G& ]' s6 l, AThanks,
9 T0 i0 C2 p7 s+ v) ^5 sboB
1 T. W' s, c" G8 I2 {" J( s- D
, a1 i! Q" {8 f, j8 i9 [& B6 h! P
3 J. [" {6 l" M8 r2 K6 I* P) FI am not sure I understood what you said. Can you elaborate? & U' Y6 Q! e3 M! t# e& _% R

9 u$ T3 T# P( m4 d0 t: ]; g6 cWhen 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 {1 Y5 h( n$ ]" l0 m

& F7 \9 ]* c5 R/ K: C3 ^$ V% OBut 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. 0 w/ W7 x5 K4 {" ]! E, d5 }' h" J

" ^$ a/ C. X" dI 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.
  z: ^/ f8 ]$ X6 S, Q; [8 }% F! C; a
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. , c/ N' z" _! C: C6 u$ w, [4 X, D

/ u( U4 f5 v$ R2 S  M4 QRegards 2 t: p" h: l. k- ~* M
Dan Lavry
; [/ ]; Z' [2 H0 o1 n_________________: b7 B; P; D2 U
Dan Lavry
% C5 u/ ]+ F; `2 g" T( l4 ?( t; `9 }) V6 X( W
Back to top        
) L/ O& K. m# h$ g
% z; B' ~6 B+ a
7 f, h9 Z/ L/ b3 OboB5 G% Z- {/ a8 d5 B
) i9 f% F5 T* A: i: N: V
3 K5 `7 b: }5 u6 {: K
: ^0 u$ h4 y" {  R2 V/ O" `. `
Joined: 28 Aug 2006
, _. V9 x# J+ d7 t5 N/ X2 {- qPosts: 34
* j; O, d7 p( I) p" a" U3 n0 M. HLocation: Everett, WA
) C; L# O, r4 K/ W0 \Posted: Tue Jan 22, 2008 1:26 am    Post subject: Re: 44.1Ksps sample rate enough for audio?   
+ D, N  b7 y- X% q7 N! R
4 N( y/ u5 t6 q--------------------------------------------------------------------------------& p! t9 @+ S4 ?' U6 M# w

( b" i; e" P2 ?6 C; Clavrye wrote:
2 C, J! U7 O( I7 k! t% F4 C1 U
7 g& [5 `. _7 @" j6 w  g
5 B3 P. u5 A5 o8 V$ F' Y  qI am not sure I understood what you said. Can you elaborate? 9 j8 E$ I, V  @

1 ?9 @4 i" y- Q. JWhen 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 W! ~" P- E$ T; {6 x* ^

7 y. ~" o/ O/ S. u/ b" a" _Regards 6 B  i9 C; q; H7 k: W
Dan Lavry
# c/ S5 E, p0 x- v! I# f7 s8 B0 ~) G1 I- r: E. e3 t6 o+ {% V

5 U6 J; u' U! f" c3 `3 b4 K& FI did indeed mean somewhat below FS/2 and that the filtering was filling in-between the dots. And, Sinewave at this 20kHz frequency % D4 ?( e1 C3 u; ?- 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 $ n& z/ a% U2 [% ^* j4 S1 t
to think about this stuff. . `+ {& ~0 Y% c$ t: d  U

& V& z7 x' v7 a- f4 W  R8 X# l' j% TI've thought that the bandwidth limit (reconstruction) was all that was necessary to "fill in" the rest of the waveform
( `0 x! ~/ F7 n& y# v( V& Q4 Sbecause all of the energy of that sinewave is below FS/2. $ w+ j( [9 C! C8 |- Y7 s

6 Z$ K5 @" S; I7 v7 WIs the "history" coming from the filtering itself ? Does this have to do with the sinc function ? % b; G8 d' n" t! J( o
3 U% ^+ L7 A& d  M- A+ e0 d
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 ? 8 p7 `$ |" e; O& o& i0 ^/ V* I
I can't remember what kind of scope it was but I seem to remember that it was an older one. 8 O' L3 z6 W3 L" S
6 U( b+ f0 T7 m3 p0 Y* z
boB
! b' H# V" L: F1 b* @4 Q' c: I! T3 b$ a_________________
# c; p8 b0 S$ pK7IQ ' s0 d# a7 W- g  D0 m: s

0 o, G4 u8 U+ M- R: SBack to top        
9 {9 _$ l" G3 _4 K% k0 A3 Z, M) v+ C8 q9 T6 W# f- U

% b  {7 A% ^2 V7 S6 G8 J5 Lzifzhu
; l3 f# V. B+ W/ d5 R; ^7 i
& @2 O6 U$ g! A/ [( \& Z  @0 o) Z6 E- H

6 r' s% S/ ]5 L7 {Joined: 26 Dec 2007
: J2 D5 O9 v0 U3 O, h% o9 l* FPosts: 4
; Z6 A/ \3 N2 S1 H3 _Location: Suzhou, China, C5 I2 c9 m) p5 X$ i; V# r
Posted: Tue Jan 22, 2008 4:53 am    Post subject:     
8 A. B6 u1 z  n* [+ k- k" d2 |
- Z4 _1 w# H0 @- b" A% h--------------------------------------------------------------------------------
! O' {3 k0 z( G6 ?9 t8 A6 @* o2 L$ ]' L8 v  T4 H
My understanding about 'history' is just increase the sample count and period, so that will be more accurate and not lost 'something'.
, Z) s  N2 a: |: C- T. `$ B* }0 g5 U9 q- S) p( q$ o
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'.
4 F" W6 f; N7 _+ E4 t9 }* d: a7 k, j1 r& M
9 |0 p- O& g- b! P2 U$ m8 z1 M

5 j$ U$ w% p( L3 r_________________0 w7 [# O' k% ~5 H9 S

+ V  j( i' A+ Y; O% z[ 本帖最后由 zifzhu 于 2008-1-22 12:59 编辑 ]
回复 支持 反对

使用道具 举报

11#
 楼主| 发表于 2008-1-22 13:00 | 只看该作者
原帖由 poposun 于 2008-1-22 09:49 发表 % r$ q( q5 p+ w; b
楼主还是没明白,不是说得很清楚了吗
+ ^8 P' P! B" N2 d% D1 G从你上个帖子里发的图看,你似乎错误的认为采样后重建信号是简单的利用单周期采样结果,但实际上,有限频域信号在使用Nyquist采样所得进行重建时,纯理论的说法是频域有限时域无 ...

+ y; d9 l/ ~& i" i8 N, `3 z7 K
0 x3 v5 j6 t/ U4 z4 _看我上面的回复o:)
回复 支持 反对

使用道具 举报

12#
发表于 2008-1-22 14:56 | 只看该作者
用个PLL不就结了?:huh:
回复 支持 反对

使用道具 举报

13#
 楼主| 发表于 2008-1-23 14:00 | 只看该作者
boB
0 o. K! w5 U! u' E% K6 W# E: d# A5 B- r/ Z/ C+ z4 [
4 e) j# M, ~( ?; W  t

4 @( k5 D# G! \0 \+ A( f$ pJoined: 28 Aug 20067 X* Y( h4 J0 X- h9 w) d5 R
Posts: 35
/ y4 G* h9 L8 M; Q) N8 ~Location: Everett, WA
: c5 {: l/ i: K0 r
Posted: Tue Jan 22, 2008 6:44 am* }. u0 B& B% b( G! j" x, B5 K5 M' |. @
   Post subject:
Hello zifzhu... I hope all is well in Suzhou, China !
# D( J9 [7 P7 X) m+ ^( d. ^6 K+ P" ^. z
I have seen this very same Excel example before. One of the problems is that the beat and triangle waves you see are not real because the bandwidth
7 A+ J; u3 H/ v1 @is not limited to FS/2... I think that if you were to filter the sample points, you will end up with a sinewave again because a triangle wave has harmonics
4 P% l4 s- b# {+ M8 cabove the frequency it has re-constructed. This cannot happen with a filtered reconstruction of a sampled signal close to FS/2. The reconstruction can
* W. k; Y: Z2 ~# j4 b! l9 WONLY be a sinewave, however it could take on a high or low amplitude.   W, y3 P& T: E9 P/ ?4 d. W

) ~, C, P- U% I- N5 R; ODoes this make sense ?? There is a very good web site that talks about this very same example and if I find it again I will definately post the URL here.
$ k" r4 d- Y) n  O& B9 S1 V9 @5 p# _5 S" q

5 W+ d# m3 V) J" G  R& R5 F0 ?boB0 o7 [9 N+ {9 \2 W  r  W
_________________
8 _7 V; I& R9 |+ u' ?K7IQ
5 d( B! Y/ a! o
Back to top

  i4 ^* p. Z3 X: e
) |1 O* w) y; g# _& A  {$ b; i6 A# F- k8 w
cj7 y" o4 Z4 q7 P- C- `
$ W/ _% ^' d8 }. O2 I
' b" \! c2 m. I* L3 h) U2 n: \

/ k( }0 m6 ?7 Z+ e' H$ A- cJoined: 29 Aug 20055 Z5 y# x% G% A, h% i0 H8 }
Posts: 58, a' {. O1 h, H+ N0 z1 Z. F
Location: Finland
; |9 l0 B8 A) ?3 L% q# U
Posted: Tue Jan 22, 2008 6:05 pm! D; W$ I% B, {6 Y, S2 o
   Post subject:

9 }' O0 C: \: k8 ?6 U
zifzhu wrote:
Hard to believe this wave will be reconstructed after filtering. At least, the total amplitude is decreased.
2 T' k# A% m5 j) t4 Z0 X& \The harmonic contains the 'beat cycle'. ; \6 b; F, x0 Q) Y

* C# e% v3 J$ j  E: M0 }2 K' NHard to believe...? OK, here you have it right before your eyes. % a/ @* h9 N" n) q/ D2 \8 |
I've used this same graph many times over the years. It shows the sample points of a 10kHz sine wave. ' u) ^8 y( g5 P1 n+ x
In the upper wave lines are simply drawn between the sample points (Sound Forge). 2 _& E# T; I- }- I% i
In the lower wave, Adobe Audition has correctly reconstructed the samples back to a sine wave.
: D! ?/ V+ ?5 x" CThe sample values are exactly the same in both waves.
% P0 X# i  f; [8 oC-J
3 q- O6 v7 C9 w) a# K- |+ E7 l
/ T8 E* j* V5 }  D# n# }+ o) z4 S. a) e% `4 i
_________________
9 o& g- [' w' [8 a"Science is a way for us to not fool ourselves." - Richard Feynman
回复 支持 反对

使用道具 举报

14#
 楼主| 发表于 2008-1-23 14:00 | 只看该作者
lavrye4 L- n, H  V; }: h6 L
Site Admin
5 s3 w" n+ [5 _+ K  o% H6 b: J& Q1 f- l& L
0 {3 O: a& B& J# w: A3 h, `
Joined: 15 May 2005
, ~0 |1 x5 u% j- K9 ]Posts: 463. H7 X3 ~0 P2 N3 ?) Y
Location: Poulsbo, WA4 |2 o- Z- B- W+ i
Posted: Wed Jan 23, 2008 3:20 am2 o8 R9 F" [+ n% k
   Post subject:

+ \$ `6 @* t; M
cj wrote:
zifzhu wrote:
Hard to believe this wave will be reconstructed after filtering. At least, the total amplitude is decreased. 9 p3 z4 k  O* }) p% n5 E6 W6 i
The harmonic contains the 'beat cycle'. ) }% n0 s3 T: l: N

" p6 V2 j( Y4 EHard to believe...? OK, here you have it right before your eyes. - h. k3 F. A3 ^  l
I've used this same graph many times over the years. It shows the sample points of a 10kHz sine wave. ! ?" m* \6 b. V" T4 G
In the upper wave lines are simply drawn between the sample points (Sound Forge).   D, s/ R" _1 r8 a: a- d" q
In the lower wave, Adobe Audition has correctly reconstructed the samples back to a sine wave. ( {5 t/ C# B3 R: I1 D
The sample values are exactly the same in both waves.
, N, @1 w# ^) N! r8 O9 d. nC-J 7 C7 a! U2 K/ ?5 A

: L9 i3 `9 T  Q. ]' Z, q
9 ?; B* m1 c4 _/ M. O
' C6 q+ _+ P* p% J1 m% L
( @% ?- R2 \3 a7 t
The Excel file shows the wave BEFORE FILTERING. I keep saying that we need to filter. It is very difficult to make a filter that will pass 22049Hz, and reject 22050Hz, but if you could do it, and if your file was much bigger then 44100 points (you would need a lot of "history"), your filtered wave would be a sine wave.
$ }7 J# c- \" h( c% t. O* A3 E: u. Z- e7 e8 t' v
No one seems to have a problem with say 1KHz sine wave. Why? Because it looks intuitivly (visualy) correct. The dots are "tracking" a sine wave shape so the process seems OK. At 1KHz, no one is questioning what happens to the signal between the sample points, and we all accept that the filter will "connect" the dots correctly. At 1KHz, look at one cycle and you see 44 points clearly on a sine wave....We did not need Nyquist to tell us that...
: O( G/ v( w/ j7 `- O! \  e, s1 A) C* Z0 U- n3 i9 i/ M6 V
So lets increase the frequency to say 5KHz. Now we have only about 9 points per cycle, and if you connect them with stright lines, it still looks "sort of" like a sine wave, but with some "sharp corners". You filter it and it becomes a sine wave - the "sharp corners" are smoothed out.... 3 b( m& u3 K; Y9 i9 t

/ n! t7 c; E# U# Y# NCJ shows you the case of 10Kz, with about 4.5 points per cycle. Connect the dots with stright lines, and things are starting to look rather "non sinusoidal". At 15KHz it looks worse, that is when you connect with stright line. But filtering "brings back" the sine wave...
3 x$ z; k. b/ f* U5 S* K& \; |1 l8 p! W& z: V2 q) @' q  `( |- J
Say you start at 1Hz signal. The filtering is very easy - pass signals to 1Hz, reject signals at 22050Hz.
1 F" o; ^; h2 u& a5 G# {' d7 s  a, C
Raise the signal 1Hz at a time and look at the required filter. At 1KHz, the filter needs to pass signals to 1KHz but block at 22050Hz. That is an easy filter, with a pass band of 1KHz and transition band of around 21Khz. At 10Khz, the filter gets tougher, the slope is about twice as steep, but still "doable". At 20KHz, you need a filter that passes 20KHz but blocks 22050Hz, that is a very steep filter - transition band to pass band ratio of 10:1. At 22049Hz the filter is so steep it would take a super computer to simulate... The transition to pass band ratio is 22049:1 ! [* O( y: [+ J! c8 ^

7 Q3 T1 M  G( L1 `- z6 W0 U7 AThere is no "sudden" break at some "magic frequency". Why is 1Khz OK and 22049Hz not OK? Or is there a "sudden break" at say 19KHz? Why not 12KHz? The whole thing is gradual - as you increase the gradualy the tone frequency, you need to make the filter steeper...
# B& L7 N* l% d6 p+ s/ C0 @. n' k) C; b8 R7 N
I said that one needs to look at a lot of history to steer the filter outcome to track the signal appropriatly. Your graph competly violated that, because connecting adjacent samples with a stright line completly discounts all the history except for a single previous sample. That is not how one filters a signal. You do need to filter the signal. # L- d0 P$ p  _& u5 f3 p
+ r* Q( J3 b' r6 ]3 ?6 X/ t
CJ's example is very good. Take a 10KHz sine signal and it looks "terrible", but when you filter it it looks like a perfect sine wave. Those that do not like math can use a scope to show it, and it will match the plot that CJ posted - the "strange" wave, when filtered, will become a clean sine wave.
. z, G1 d7 b: L& W3 ^8 a! K) j( M7 P
Better yet, take a DA converter, a good one (with a filter). Run a 10KHz tone through it, the look at the output with a scope. It does not look like 4-5 dot per cycle connected by a stright line. It looks like a clean sine wave!
6 y0 E- N/ ?" C3 W. l9 v! U0 d/ ]8 f: Y8 Z8 x
And before you point out that the DA is up sampled, I will remind you that up sampling makes a lot more sample points (samples) out of the few original sample point by a process no other FILTERING. 2 g9 e0 _) B% Z. Q, s9 d0 G* K4 W2 }
% b$ x' R$ d2 R8 F: a5 \: i, p
Regards * O( Y! g9 F4 p9 o
Dan Lavry5 M' K- c1 v+ U$ d! M9 n
_________________
' g0 ?/ A! O' ]/ d1 g/ Q) W2 P0 S. ?Dan Lavry
回复 支持 反对

使用道具 举报

15#
 楼主| 发表于 2008-1-23 14:01 | 只看该作者
zifzhu
8 |( T# Z3 b& ^5 h& h7 T# l0 F) N# D
! Q& U/ t: a2 p" h3 C
) r0 O4 K' B  m
Joined: 26 Dec 2007, x) O6 c2 t3 i" E
Posts: 69 Y% p0 C2 K/ o' @; V
Location: Suzhou, China
6 a; |6 o5 }/ h5 n! ~* j
Posted: Wed Jan 23, 2008 5:28 am# v2 d4 i* B% Z  N, `3 S
   Post subject:
Hello boB and Cj. I was also confusion about filtering that triangle wave. The low pass filter should like just smooth the wave. The better idea is to use fourier transform to analyze it's frequency energy.
( s5 Z5 q& B8 \' ~6 h9 R1 X' k8 a6 ?$ E/ j! e4 G9 E  r) m
I've found a excel FFT analyze tool from maxim. http://www.maxim-ic.com.cn/appnotes.cfm/an_pk/3292
& [" ~5 h. w5 K9 d: j8 T% r. B
4 f# z% I9 _8 jAfter input 128 samples of 2.2X sampling wave. The FFT shows it does contain other frequency(not a straight bar in 20Khz). It contains 18Khz at -30db, 11KHz at -40db, 0-1.7KHz at -43.4db. 9 C, A2 M% G/ T! n

7 |/ ^' B0 A. F% hAgain, I input 4.4X sampling wave. The FFT shows a better result. 0-1KHz is only -51db. About 7.6 db better than 2.2X sampling. 6 k& s$ }% o! H5 \2 ?

/ z0 @& O7 e/ _9 z. ?) k1 XSo. Can I put a conclusion '2X sample rate can not do a perfect sample'. To do a better sample. A higher sample rate is prefered." o$ C9 S. S  m8 U+ A- `2 {
_________________3 e1 m* N; ^8 R6 ^

6 o2 {" ?  d' [4 e- t5 x4 ^: x1 B
( G, _" M9 l6 d% z: C7 G. n[ 本帖最后由 zifzhu 于 2008-1-23 14:09 编辑 ]
回复 支持 反对

使用道具 举报

16#
 楼主| 发表于 2008-1-23 14:02 | 只看该作者
zifzhu( ?9 \( c( a' j( D: U
6 L6 A$ C2 N$ o0 a+ Y" q
& ?' L$ Q( F, [0 p# `: W
) u# L; p, [. _# a2 \
Joined: 26 Dec 20071 Y6 Q* M1 c( k" Y; S. A$ I
Posts: 61 P# k. @) k% J: G
Location: Suzhou, China

' W/ T4 \7 |" g7 O, m& c" p
Posted: Wed Jan 23, 2008 5:54 am- t3 V: i* q  d
   Post subject:

; Z2 X% j0 N4 v% s* L. G
lavrye wrote:
cj wrote:
zifzhu wrote:
Hard to believe this wave will be reconstructed after filtering. At least, the total amplitude is decreased.
* }- }7 t5 V  ]8 |$ Z8 s+ g+ Z7 QThe harmonic contains the 'beat cycle'.
1 ]2 l$ A& r4 F& _/ P2 ]* s( n

1 O2 Y1 w" O, B  G, V$ y/ PHard to believe...? OK, here you have it right before your eyes.
+ t. w" o" I6 }2 R" {1 j0 n6 YI've used this same graph many times over the years. It shows the sample points of a 10kHz sine wave. + L/ B/ h+ [9 r# Z. U3 @
In the upper wave lines are simply drawn between the sample points (Sound Forge).
* u' T* [5 o: N7 sIn the lower wave, Adobe Audition has correctly reconstructed the samples back to a sine wave. ; b0 X. J: V1 B% m& E
The sample values are exactly the same in both waves.
8 z1 [( c8 }& Z: `C-J 6 u" X% E' l2 C- O2 @

/ A4 D; W4 C7 w4 h7 y" o! g
0 D. W( S- i, `* `: n- l

- g. c% [+ q; @( \. b8 }& j2 }& G; R+ M4 J4 I$ D
The Excel file shows the wave BEFORE FILTERING. I keep saying that we need to filter. It is very difficult to make a filter that will pass 22049Hz, and reject 22050Hz, but if you could do it, and if your file was much bigger then 44100 points (you would need a lot of "history"), your filtered wave would be a sine wave. 4 z$ e/ O/ x6 \' S; |
* b2 s+ a8 ~/ @& F# {7 s
No one seems to have a problem with say 1KHz sine wave. Why? Because it looks intuitivly (visualy) correct. The dots are "tracking" a sine wave shape so the process seems OK. At 1KHz, no one is questioning what happens to the signal between the sample points, and we all accept that the filter will "connect" the dots correctly. At 1KHz, look at one cycle and you see 44 points clearly on a sine wave....We did not need Nyquist to tell us that...
4 z" t  u: ~: r# ]* X( r! S6 Z! L! Z4 f
So lets increase the frequency to say 5KHz. Now we have only about 9 points per cycle, and if you connect them with stright lines, it still looks "sort of" like a sine wave, but with some "sharp corners". You filter it and it becomes a sine wave - the "sharp corners" are smoothed out.... ! O6 S* _5 Y- J( t) @7 H- |

. h8 D$ X; w( y9 _" J$ n& m- fCJ shows you the case of 10Kz, with about 4.5 points per cycle. Connect the dots with stright lines, and things are starting to look rather "non sinusoidal". At 15KHz it looks worse, that is when you connect with stright line. But filtering "brings back" the sine wave... - R4 n6 V6 `: A. c
1 ^1 Z' ]0 m! n; I+ I5 R
Say you start at 1Hz signal. The filtering is very easy - pass signals to 1Hz, reject signals at 22050Hz. 0 G/ ~  @1 X- m9 j, {6 e5 C
' s0 q3 _( ]% n$ X" B
Raise the signal 1Hz at a time and look at the required filter. At 1KHz, the filter needs to pass signals to 1KHz but block at 22050Hz. That is an easy filter, with a pass band of 1KHz and transition band of around 21Khz. At 10Khz, the filter gets tougher, the slope is about twice as steep, but still "doable". At 20KHz, you need a filter that passes 20KHz but blocks 22050Hz, that is a very steep filter - transition band to pass band ratio of 10:1. At 22049Hz the filter is so steep it would take a super computer to simulate... The transition to pass band ratio is 22049:1   U7 m; f2 J5 _; w! N
" P+ L& x# \/ o* J
There is no "sudden" break at some "magic frequency". Why is 1Khz OK and 22049Hz not OK? Or is there a "sudden break" at say 19KHz? Why not 12KHz? The whole thing is gradual - as you increase the gradualy the tone frequency, you need to make the filter steeper...
; ?6 p" b1 G$ Z9 r( W: l% c1 `  `! F- }" U! `0 I- ?* B
I said that one needs to look at a lot of history to steer the filter outcome to track the signal appropriatly. Your graph competly violated that, because connecting adjacent samples with a stright line completly discounts all the history except for a single previous sample. That is not how one filters a signal. You do need to filter the signal. : @- {- v- ]( f. e4 E& L7 [6 F

4 h7 l! q1 P  ~# a7 h* QCJ's example is very good. Take a 10KHz sine signal and it looks "terrible", but when you filter it it looks like a perfect sine wave. Those that do not like math can use a scope to show it, and it will match the plot that CJ posted - the "strange" wave, when filtered, will become a clean sine wave.
7 G, U- A/ H- z+ x% E2 O* R5 n& z/ d$ W
Better yet, take a DA converter, a good one (with a filter). Run a 10KHz tone through it, the look at the output with a scope. It does not look like 4-5 dot per cycle connected by a stright line. It looks like a clean sine wave!
8 i3 P8 U& x9 t" P6 x0 K- g4 P2 L3 ^
And before you point out that the DA is up sampled, I will remind you that up sampling makes a lot more sample points (samples) out of the few original sample point by a process no other FILTERING. 4 w/ I' |. m8 v& z! j
( `$ e: h. X( u9 D5 R/ q2 {8 g
Regards ; [) y; J1 x5 L; k+ k: \
Dan Lavry
" r& K+ {$ V# F
Hi Dan, 8 K4 H9 y4 F2 j% d6 ]* V; c4 K
I think the problem is the harmonic contains not only higher frequency but also lower frequency, that is not able to be filted by a low pass filter. As my previous post, the higher sample rate can get a cleaner lower frequency.
7 q1 `( x3 Q6 x: X  K+ T7 d; X
& N  o0 w' m/ t4 r) G- ~There is no PERFECT sampling. The higher sample rate loss less information. Nyquist frequency is just a point that avoid aliasing that means samples the freqency correctly, not perfectly.+ ^6 w' L! {. |( E. n0 b; K# n
_________________3 R- K1 q" P, Q4 L' N2 L# J
回复 支持 反对

使用道具 举报

17#
 楼主| 发表于 2008-1-23 14:03 | 只看该作者
又有新的进展了:a) + |. C$ V$ C. E; {  R; [* [) l
' p4 a  R/ _9 P* R, U9 @
估计很快就会有结论了:p
回复 支持 反对

使用道具 举报

18#
发表于 2008-1-23 22:50 | 只看该作者
楼主很执着:lol:
回复 支持 反对

使用道具 举报

19#
发表于 2008-1-24 11:38 | 只看该作者
44.1足够了,不是技术决定的,是市场决定的
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

广告投放或合作|网站地图|处罚通告|

GMT+8, 2025-8-20 02:10

Powered by Discuz! X3.4

© 2001-2017 POPPUR.

快速回复 返回顶部 返回列表