POPPUR爱换

标题: 讨论观点：整数倍提高采样率可以在不失真的前提下改善听感 [打印本页]
作者: 归德居士    时间: 2010-4-1 21:53
标题: 讨论观点：整数倍提高采样率可以在不失真的前提下改善听感
本帖最后由 归德居士 于 2010-4-1 21:56 编辑 - B8 ]; S1 B/ e8 w

2 u( L$ @0 I( O7 m2 B( ~+ `: F. j和大家讨论个观点啊，我分析认为整数倍重采样提高采样率可以在不失真的前提下改善听感。$ f% D/ k' N3 o: v$ H) v
, C0 T& W! G$ Z" N$ K! _, c
因为我想，整数倍提高重采样不需要破坏原来的采样点，仅仅是在原来的采样点每两个之间增加几个采样点来过渡。
& E. V" K7 X* D( t- F. @3 j' N, m( j0 v
比如第一个采样点数值为a，第二个数值为b，如果把采样率提高整数倍，就只需要在两个采样点之间增加一个数值在前后两个采样点之间的采样点作为过渡，比如说采用直线算法升频两倍，就只需要在两个采样点之间的中间点增加一个数值为 二分之一（a+b）的采样点，升频三倍，只需在两个采样点之间的三分之一处和三分之二处各增加一个采样点，共增加两个采样点，数值分别为 三分之二a+三分之一b 和 三分之一a+三分之二b。
- j1 ?; X. N1 r$ {# U6 a7 ?
* F( s3 f" j* C这样，原来的采样点都没有失真，新增加的过渡采样点使声音更细腻，原来是从一个采样点猛跳到另一个，现在是在中间加几个，逐步过渡过去。
$ w" f+ f) u. e$ c: K1 E0 m- o( F) X, w$ \) l$ Y
鄙人才学疏浅，经验匮乏，还请大家多多指教和批评。
作者: 深水的鱼    时间: 2010-4-2 00:31
整数倍转换是可以改善听感的，但非整数倍转换不行，如图
& t# g# e. C5 t1 f[attach]1250028[/attach]
作者: chenxy    时间: 2010-4-2 00:52
讨论下假如40Hz采样率提高到60Hz会怎么样还比较有意思。
作者: bull    时间: 2010-4-2 01:26
你的结论正确，但是原理不对
作者: sd-iori    时间: 2010-4-2 13:21
N久前，，就流行88Khz采样率去播放的啦，，，
. ~. {/ ^9 ]  n/ L  O也确实感觉好一些，，，9 F) Z# a6 w* z% H
不过原理，，不懂，，，
作者: caileili    时间: 2010-4-2 13:29
本帖最后由 caileili 于 2010-4-2 13:31 编辑
. \" L7 U0 V9 t8 V9 b; d
; S0 S. {& J/ A/ e+ F4 Y- Q9 g1 S单纯的整数倍提高采样率，然后中间采用简单的线性插值算法填充的话，只能说可以减少失真，但对提升音质意义不大！
  K% u0 k2 p+ n$ }如果能配合扰动算法，按音频信号曲线或趋向计算插值，这样可使波形接近采样前的原始声音，从而一定程度上提高音质！  I6 q& [7 n( y/ c3 f
但这个毕竟类似于把低像素的照片放大，最先进的放大算法画质也提升不了多少的，所以没什么人钻这个牛角尖！
作者: bull    时间: 2010-4-2 19:59
简单说一下，这个做法意义很大。8 ^1 z4 H# G; x, E2 ]6 \
因为可以大幅度提高SNR  r: G$ @8 e$ t2 P4 i" s
大到了示波器上都能看出区别来
作者: Timme    时间: 2010-4-2 21:26
在PC上软件实现这个转换没有任何意义6 T) L4 M$ p/ a* s& p- R( Q

5 }1 W7 b& m8 k0 u$ ~# L  S对于固定的DAC，采样率越低非线性失真越小（就是频率范围满足的前提下采样率越低越好）
作者: bull    时间: 2010-4-2 23:11
在PC上软件实现这个转换没有任何意义
& b7 j4 f0 K& W0 }- G0 l8 @6 ^% A9 R+ |
对于固定的DAC，采样率越低非线性失真越小（就是频率范围满足的前提 ...1 O. S9 H3 C$ S3 q* L; C
Timme 发表于 2010-4-2 21:26
2 L4 V* O; ~& f; u( _8 h- y
抱歉，并不是越低越好
: f3 Z; i& k4 }6 R就现在的声卡 2倍频是个很好的选择
作者: Leading    时间: 2010-4-3 00:33
支持6楼         。。。
作者: Metaverse    时间: 2010-4-3 10:29
如果能配合扰动算法，按音频信号曲线或趋向计算插值，这样可使波形接近采样前的原始声音，从而一定程度上提高音质！+ i  {' P, v$ Y6 p8 t ..., B% J9 z4 v7 b, z
caileili 发表于 2010-4-2 13:29
* b; B' t& R5 k: `2 W

1 L! x* U; ~' VFOOBAR有没有这样功能的插件？
作者: Timme    时间: 2010-4-3 11:03
本帖最后由 Timme 于 2010-4-3 11:11 编辑
( ?5 B% G5 i: I
* T2 R- c* M9 `8 C6 N& x无论整数倍还是非整数倍，算法都是一样的。2楼的图中波形没有经过低通，这种图在N年前就被我批过了
9 [$ Y! _! k& K, p. B4 b) r" i  N2 k8 F) t
引用我自己设计的一个SRC算法（只是自己独立想出来的而已，不敢说第一个发明啊。后来我好像看到IEEE 70年代就有人发表这种算法了）
. u+ o; X! i+ }7 \0 S+ V9 W6 @$ \( h, w+ k/ U/ g. A
（1）采样频率2n->4n的算法，整数倍" I' ]) ?/ _0 H/ T8 q& P4 ^

' i5 `8 k' ]6 Y! |3 M原始数据：a,b,c,d...7 x) {! [8 F3 X+ F) ?
先机械填充为a,a,b,b,c,c,d,d...（同样时间采样点数目翻倍，多出来的采样点数值复制原采样点）5 c9 q* ~; L$ J  @7 K
通过一个低通，变成a0,a1,b0,b1,c0,c1,d0,d1...1 B* E6 c' u6 [2 l
SRC计算完成。同样时间内，4个采样点变成了8个采样点。
1 n: I- w- E2 z7 l% X8 |（当低通性能足够好并且相差接近零时，有 a约等于a0，b约等于b0...）/ l+ S$ L0 L  t* A' X

1 r6 w6 V  o! ^/ {4 q: o（2）采样频率2n->3n的算法，非整数倍
$ Z" B- p( R3 |/ c: x0 v2 i8 f5 @/ h8 O/ J
原始数据：a,b,c,d...
: K6 H# ^* l: M! x9 X先机械填充为2，3的最小公倍数的采样率：6，即a,a,a,b,b,b,c,c,c,d,d,d...! {' d& m" h+ O; l" p) M2 E4 `, ^
通过一个低通，变成a0,a1,a2,b0,b1,b2,c0,c1,c2,d0,d1,d2...
! C( F' c3 c6 [4 b: B7 m我们跳一个取数：a0,a2,b1,c0,c2,d1...
7 N: A5 x. q, Q* }- dSRC完成。同样时间内，4个采样点变成了6个采样点。" J% g  b. O0 @. B
5 o. q2 \: y8 a! h& U
**************************************************6 r4 b" t6 q+ q9 m2 b% B" ~
. n2 `) M  `0 a+ ]+ H+ e
你看到整数倍SRC和非整数倍有什么不同吗？没有的。其实非整数倍升频 = 整数倍升频 + 整数倍降频 而已，而且降频时还根本不需要计算，直接跳着取数就可以了。
作者: Timme    时间: 2010-4-3 11:21
另外bull，所谓硬件升频提高SNR是DAC内部电路而已，但如果你把DAC工作模式从48k改为96k，那输出频段24k-48k这范围的噪音DAC内部就不过滤了，SNR反而低了。也可以理解为外部选择的采样频率越高，DAC内部升频的倍数就越低。如果说96k比48k好，我唯一能想到的就是低通在拐点附近性能劣化而已，但这个点离听觉极限已经有几khz的余量我觉得足够了。
4 a# r3 q* J) y3 F
& s; O8 `) U0 t我觉得这是个很简单的问题，不知我有没有理解错
作者: bull    时间: 2010-4-3 14:55
另外bull，所谓硬件升频提高SNR是DAC内部电路而已，但如果你把DAC工作模式从48k改为96k，那输出频段24k-48k ...0 V! Z+ s( C) w. h+ J6 a; e. m2 g
Timme 发表于 2010-4-3 11:21
( g9 @* |- n2 M" B3 f. J8 ?
你对现在的DAC芯片怎么工作的并不清楚。
作者: Timme    时间: 2010-4-3 16:10
比教科书上写的Multibit + OverSampling + Noise Shaping还多了啥？
作者: caileili    时间: 2010-4-3 22:36
本帖最后由 caileili 于 2010-4-4 02:05 编辑
( Y4 z  C+ D- `, F
（2）采样频率2n->3n的算法，非整数倍…………
  f/ P* W" R4 v/ z: s整数倍SRC和非整数倍有什么不同吗？没有的。其实非整数倍升频 = 整数倍升频 + 整数倍降频 而已，而且降频时还根本不需要计算，直接跳着取数就可以了。+ l# K  {4 b, h0 p

7 M6 _  Z" g/ d* C% ^+ n引用 ...
5 L, f: p: G" ^1 oTimme 发表于 2010-4-3 11:03
3 W! s2 m; O8 w6 Z, C- J
% v! a1 L, q1 M1 p, z3 I

! K0 H8 w! e8 D! f" C5 Y, [, f. T' g- ?+ ^4 h
那就按这位兄台的高见对图例中非整数倍取样重新进行演算：（没时间画图，大家看步骤理解吧！）
1 J& u: e( k/ ?% F5 Y) o1 e' e原始：0,1,0,-1,0,  0,1,0 ……. v! l& h6 T4 b, ]9 v. I0 {
整数倍升频：[0,0,0],[1,1,1],[0,0,0],[-1-1,-1],[0,0,0],  [0,0,0],[1,1,1],[0,0,0] ……
/ u( ]: o0 }4 F7 c! k$ m低通(简单)：[0,0,0.33],[0.66,1,0.66],[0.33,0,-0.33],[-0.66,-1,-0.66],[-0.33,0,0],  [0,0,0.33],[0.66,1,0.66],[0.33,0,0] ……
; C) o  R/ W1 Q9 _1 u( b% r. o1 A整数倍降频：[0,0.33],[1],[0.33,-0.33],[-1],[-0.33,0],  [0],[0.66,0.66],[0] ……
+ Q2 z4 c/ \2 \  M% @1 a; w输出：0,0.33,1,0.33,-0.33,-1,-0.33,0,  0,0.66,0.66,0 ……
3 T# n+ M4 |; f9 B4 f# b* h+ s3 G  V3 Y, [+ ^& t  F( y5 k) R( i
原始的第二三四取样点之间是条直线波形[1,0,-1]，经非整数倍取样后输出波形变成了[1,0.33,-0.33,-1]仍是条直线，此时没有失真；
8 N3 y, b. j" W. W原始的第四五取样点之间是条直线波形[-1,0]，经非整数倍取样后输出波形变成了[-1,-0.33,0]是条折线，即对原始波形间存在变形/失真  \% S7 z5 F/ h% ]0 L  J2 Y
原始的第六七八取样点之间是条上三角波形[0,1,0]，经非整数倍取样后输出波形变成了[0,0.66,0.66,0]是个被削波的梯形，即对原始波形间存在较大变形/失真
& J7 t: U" ]7 {( u# m7 c看来这种非整数倍SRC算法的失真确实是存在的（演示图片描述是恰当的），而且可以看到是在整数倍降频（直接跳着取数）过程中产生了失真，在输出时的时间位移上也会产生失真，这些不是改进低通平滑算法能解决的，所以不能说是我演示时低通算法的问题哦！
9 F1 i9 H: R2 _* C  Z; G- D% X4 s) X
# F4 [. m7 V; {换个角度来说吧，如果这种简单的非整数倍SRC算法不存在失真，那AC97的SRC问题也就不会出现了，很多声卡也不会搞双晶振了，Foobar里面也不会开发这么多重采样插件了……！
  T" c5 S# e; I+ v" iFoobar里的SSRC X重采样插件的保真度挺高，所需CPU运算也挺高的，所以它的算法原理应该相当复杂的，不然也不会达到这么高的转换音质了！
作者: Timme    时间: 2010-4-3 23:20
无语，低通这一步直接无视了么。。。
作者: bull    时间: 2010-4-3 23:36
无语，低通这一步直接无视了么。。。
1 l( I+ B3 {2 q- h7 ?Timme 发表于 2010-4-3 23:20

( k2 _$ _* b5 R$ ~哈哈哈，一句话直接击中要害
作者: Timme    时间: 2010-4-4 01:10
本帖最后由 Timme 于 2010-4-4 01:23 编辑 " [0 E! ]  k- f' |" ]# ^& ~7 _
* l- l  E) |/ o) }% {
在此给出一个例子，是用我楼上的算法把6khz采样率升为8khz采样率。
6 q( D! ^7 }+ k! d% t0 M% c5 e
1 B4 a: T! A- ~' x6 K对6khz的原始数据，我们取12个数，循环：“0 1w 2w 3w 2w 1w 1w 0 -2w -3w -3w -2w”: ~7 @$ F' p# T+ _: S2 |
7 \: c1 \7 n' ]
【设计原始数据是个技术活，你不能产生边界频率分量（在此是3khz），此外不要有直流分量。】' _8 L7 F* z/ O1 f% M9 j! D6 w

  Z. C  T% U' K; o4 f简单扩展为24khz：“0 0 0 0 1w 1w 1w 1w 2w 2w 2w 2w 3w 3w 3w 3w 2w 2w 2w 2w 1w 1w 1w 1w 1w 1w 1w 1w 0 0 0 0 -2w -2w -2w -2w -3w -3w -3w -3w -3w -3w -3w -3w -2w -2w -2w -2w”  `( ~" r9 k+ m' q7 t( S- {% K

6 r7 Y3 _8 ]) R& `* u! F1 }% V1 E+ V$ e低通后，取一个数跳两个数，得“1402 9151 15496 24899 28865 21793 12681 9932 9094 1594 -12561 -24898 -30031 -30037 -24944 -12433”: w+ A& m, E8 v
: A7 Z6 L0 ~4 c# E& l+ [! t$ w, L/ F
至此SRC计算完成。
9 c4 \8 R" D, H3 O1 O; s; q
: @1 ^; [  I' w附SRC前后波形，自己看看有没有改变吧。另外附上每个过程产生的WAV文件供下载，可以拿去慢慢研究，看频谱，等等。。。
4 B+ P' X8 X+ ]
/ Z& `0 m, a7 h2 E% Q, X6khz下的：
, p7 \, m! W8 G
, M3 {/ q3 W1 d0 Z[attach]1251881[/attach]0 Z. i3 q! y, E; @% r7 w! H# m1 E
; t4 J  A( O+ [4 t

" j, W( z( ?% H) f: c3 o+ H24khz 低通后的：
* B6 B- h. I- ~' `4 S2 ~6 L1 z; q5 `1 K% C. a& e7 J* G! p( s
[attach]1251884[/attach]
2 H( P4 V2 H' R3 |8khz下的：
9 T$ ?/ [- H' I& M( d/ o
5 N, z! X- a/ _% K[attach]1251882[/attach]+ J- J6 G( X6 a( y$ a
; l8 ?1 M( Y+ I) K1 v; Z  V
为什么24khz-lowpass.wav频谱那么好看，而8khz.wav多了一些很小的谐波，这是因为由于低通不理想，低通后的波形并不是严格的以48个采样点为循环的。而我只对其中48个采样点进行了采样并拿来当做降频的循环初始数，所以降频后的波形。。。“欠Dither”
, t, x+ V8 J% c- H7 t+ A6 @. I6 r; z* g
[attach]1251883[/attach]



欢迎光临 POPPUR爱换 (https://we.poppur.com/) Powered by Discuz! X3.4