帯域制限された信号のサンプリング中に、サンプリングインスタンス間で「失われる」信号値


9

ナイキストシャノンサンプリング定理によれば、帯域幅をもつ任意の連続時間信号 B ナイキスト周波数よりも小さい fN=fs/2 (と fs サンプリング周波数)、これはサンプリング周波数でサンプリングされます fssinc補間(つまり、Whittaker–Shannon補間公式)によって完全に再構築できます。

一定のサンプリング時間で、大きさが制限された連続時間信号である未知のサンプルをサンプリングすると仮定します。 T=1/fs サンプルインスタンスで kT、(kZ)、サンプリングジッタまたは量子化なし。次の制約を追加しますB=αfN0α1

私が理解したいのは以下です: サンプルの瞬間k、それぞれについて決定したい αサンプル間の任意の連続時間信号の最悪の場合の部分的な「オーバーシュート」k1 そして k、私が持っていたかもしれないこと。つまり、連続時間信号が、サンプルの瞬間の最高(絶対)サンプリング値よりどれだけ高かったかk そして k1。サンプリングによって「失われた」連続信号または再構築(sinc補間は完璧なので!!)。

例: 設定α=1 離散時間信号[1,0,1,0,1,1,0,1,0,1]と仮定します(中央付近のdouble 1に注意してください。この信号は α=1?)。サンプル(黒のインパルス)からのsincの再構成(青い線)は次のようになります(各サンプルに属するsincを灰色でプロットしました): シーケンスのSinc Reconstruction サンプル間の「オーバーシュート」k=0 そして k=10.7 または 70%。したがって、元の帯域制限された連続時間、または「完全に帯域制限された再構築された」信号で、値1.7のピークを逃しました。私が3つ以上の連続した1を置いた場合、オーバーシュートは少なくなります(ギブス現象は結局はるかに小さくなります)。したがって、このような2つの連続したサンプルは「最悪のケース」です。

信号を両方向に拡張すると、オーバーシュートが大きくなります。 ここに画像の説明を入力してください これは、1.1 ほぼ2.1の値に。

シーケンスの長さ 2m、この「オーバーシュート」 o(m) 無期限に成長し、 o(m)ln(m)に行く いつ m。これは、sincの各サンプルが建設的な「干渉」を生み出し、1/πn (単位sincのすべてのエンベロープの貢献) n 収束しません。

これは(次のように)似ています:値0を絶えずサンプリングする場合、値0のノードでのみサンプリングされる無限振幅の連続時間信号を再構築することもできます。 sinπfst。これは私にも同じことを教えてくれます。信号をナイキスト周波数にすることを許可した場合、「見逃す」ことができる最悪のオーバーシュートは無限です。

これで、 o(m)|α=1=。そして、私たちはそれを推論することができますo(m)|α=0=0 (帯域制限されていることがわかっている一定の信号をサンプリングすると、一意の一定の再構成が行われます)。

仮に α<1

ここで、これと同じsinc補間を行うと仮定した場合でも、 α<1、 お気に入り α=0.5。次に、(私の直感は言う)この効果は低下し、有限であるはずです(m)!。用するので、任意の信号煉瓦壁は、帯域幅に限定されるものでαfs/2、のフィルターインパルス応答を取得します h(t)sinc(tkTαT)(正しい?)。したがって、信号の遷移は上記の変化するインパルス列の例ほど速くはできません。したがって、再構築中の各sinc関数の寄与は、無限の建設的干渉を生み出すことはできません。

私の問題: ここから先に進む方法がわかりません。2つの連続するサンプル間で見つけることができた最悪の場合のオーバーシュートの「証拠」を形成する方法α<1、 ために any信号(必ずしもこれらの単位インパルス列の例とは限りません)。の指定された値α 坂道をくれ dh(t)dt 帯域制限畳み込みカーネルの h(t)、これはどれだけ連続したサンプル異なる必要があるかについて何かを教えてくれるはずですが、そこから一般的な結論に到達するための手順はわかりません。


このような病理学的シーケンスについては、2002年にcomp.dspで議論しました。件名:サンプルデータ補間Poser、groups.google.com / d / msg / comp.dsp / EQ31d-2SS2o / wT5HXbjQpogJおよび2003、件名:再構築のための最悪の場合の信号、groups.google.com/d/msg/comp.dsp/xwb9p3awrOg/zl20Wl2EiesJ
Olli Niemitalo

関数の帯域幅をゼロ交差の平均密度の上限に関連付ける定理があると思います。さて、ほとんどどこでも無限関数の場合、これらの有限値サンプルはおそらく、ゼロクロッシングが有限値関数に似ているようです。それらの平均密度には上限があります。
Olli Niemitalo 16年

おかげで、時間があれば、両方のGoogleグループディスカッションを詳しく読みます(図での回答はどこに行きましたか?)。それでも、MBazの答えは、最悪の場合の最大絶対微分が存在することを示唆しているようです。sup|x|有限です、有限になります。したがって、帯域制限された信号の場合、無限の値にはなりません。それはあなたが提案しているものとどのように関係していますか?
2016年

離散シーケンスがローパスフィルタリングをそのまま存続させるようなものでなければならないことを考慮しなかったので、私の回答を削除しました。したがって、で何が起こるかについては正しいかもしれませんが、上の私のコメントはそれに同意します。サンプル間で見つかるかもしれないので、あまり言いません。α<1sup|x|
Olli Niemitalo 16年

期間が無限長の周期列の場合、数学はより単純になるのではないかと思います
Olli Niemitalo

回答:


7

私には本当の答えはありませんが、この結果があなたを助けるだろうと感じています。バーンスタインの不等式は、信号がに帯域制限されている場合、ここで、は「最小上限」を表します。x(t)|f|B

|dx(t)dt|4πBsupτR|x(τ)|,tR
sup

この不平等については、Amos Lapidothの優れた(そしてPDF形式で無料の)本「A Foundation in Digital Communication」で見つけました。証明はMA Pinskyの「Introduction to Fourier Analysis and Wavelets」にあります。


ありがとう!それは本当に便利です。連続時間信号からの推論。これは、直線的サンプルで「前進」外挿したときにことを意味すべきであるとサンプルの「後方」、我々は、連続時間信号は、その三角形の下になければならないことを知っている三角形になるだろう。そうでない場合は、より高い周波数のコンテンツが含まれます。次に、最大値とは言えませんでした私が許可する最大振幅で実際に最高周波数()(co)sine によって制限されていますか?(それがどのように関連しているかを確認するためにそれを見つけることができる場合、私はピンスキーの証明を読まなければなりません)k1k|dx(t)/dt|αfN
Retinite

私が理解している証拠はまだ見つかりません。1つの証拠を取得するためだけにPinksyの本に100米ドル以上を費やしたくありません。私の直感では、我々は特定のかもしれないと言うことの代わりに、(対付き)、最大許容信号値。私はここでいくつかの一般的な証拠を見つけましが、でのノルムの使用を理解していません。また、答えがが関数になることを示唆しているかどうかはわかりません(周波数領域での近似)。|dx(t)dt|2πBAmax2π4πAmaxL1||g||1g(ξ)rect
2016年

1
気にしないでください。前述の投稿の説明を入力するだけで、証明を作成しました。
Retinite 2016年

5

観察

私はあなたの1と0の代わりに+1と-1をシーケンスで使用しました。を使用すると、最初の2つの図のバンド制限連続関数は(上記の変更を加えて)次のようになります。α=1fm(T)

(1)fm(T)=k=1mmsign(sinc(πkπ/2))sinc(πTπk),

どこ:

sinc(T)={sin(T)/Tif T01if T=0sign(x)={1if x<00if x=01if x>0.

fm(1/2)はの対数まで線形に増加します。m

$ m $の関数としてのピークの成長
図1.の関数としてプロットされた対数の横軸は、成長をとして線形化します。fm(1/2)log2(m)m

Wolfram Alphaの助けを借りてを簡略化することができます:fm(1/2)

(2)fm(1/2)=2k=1m|sin(π(k0.5))π(k0.5)|={4/πif n=12(ψ(0)(1/2)ψ(0)(m+1/2))πotherwise,

ここで、はディガンマ関数です。に関するシリーズの主要な用語は次のとおりです。ψ(0) (2)m=

2log(m)π,

これは、図1に見られる線形化を説明1。我々は今、正規化されたバージョンを構築することができる。関数のそのbandlimitednessを継承するが、として爆破ない:gm(T)fm(T)m

gm(T)=πfm(T)2log(m)

、の零点でサンプリングされたナイキスト周波数の正弦波に近づくように見えます。mgm(T)

$ g_ {100000}(T)$
図2.は爆発しません。g100000(T)

元のナイキストシャノンサンプリング定理では、最高周波数がサンプリング周波数の半分を下回る必要があるため、境界線のケースではカバーされないようです。ただし、任意に大きい有限と、その結果として任意に大きい有限は、まだカバーされています。mfm(1/2)

証明の概要

元のステートメントの証明の概要は次のとおりです。サンプリング周期を1とします周波数以下に帯域制限します。ここで、は周期2と周波数を表します。してみましょう整数すべてのための有限で。すべてのについて、自明なケースを除外します。してみましょう。そのため、一部のについてはとなります。どちらか:f(T)αππα<1f(T)Tf(T)=0Tg(T)=f(T)/supTf(T)g(T)0T

ケース1. ある整数場合は。はすべてのに対して有限です。g(T)0TsupTf(T)T

ケース2. すべての整数に対して。は、いくつかのに対して無限大です。スケール係数までは、はそのゼロの小数によって決定されます。残りのゼロをもう1つ使用して、関数を消滅させます。すべてのに対してです。これは矛盾しています。以前に、一部のについてためです。ケース2は当てはまりません。g(T)=0TsupTf(T)Tg(T)αg(T)=0Tg(T)0T

したがって、ケース1はtrueであり、すべてのに対しては有限です。f(T)T

均一に分布したゼロの一部を使用して、それらのゼロの平均密度と比較して帯域幅が比較的低い場合に関数を再構築できるという明確な証拠を見つけるとよいでしょう。場合、をにはサンプリング定理で十分だと思います。文献の中で、私はいくつかの興味深い声明を見つけました:α<1g(T)

定理4.1のパート2の証明は、信号が点消失するという性質を持つ帯域制限信号が同じように消失する必要があることを示しました。=πx=nZ

Jeffrey Rauch、「フーリエ級数、積分、および基本的な複雑な分析からのサンプリング」。

おおまかに言えば、は同じように消失することなくcosよりも多くのゼロを持つことができないことはよく知られています。gλt

BF Logan、Jr.「帯域通過信号のゼロクロッシングの情報」、Bell System Technical Journal、vol。56、pp.487-510、1977年4月

1次元信号の一意の仕様に関するほとんどの結果は、帯域制限された関数が全体(どこでも分析)であり、そのためゼロ(実数および複素数)によって定数および指数係数内で一意に指定されるという事実に基づいています。すべてのゼロが実数であることが保証されている場合、任意の帯域制限関数は、その(実数)ゼロクロッシングによって一意に指定されます。
...
追加の作業には、それらが複素数のゼロも含んでいるという事実にもかかわらず、それらの(実際の)ゼロクロッシングによって一意に指定される信号の識別が含まれます。これは、ゼロクロスレートが何らかの意味で情報レートまたは信号の帯域幅よりも高い場合に可能です。

SR Curtis、「ゼロクロッシングからの多次元信号の再構築」、論文、MIT、1985。


面白い。他の投稿からの最大微分アプローチは、非常に最悪のケースの推定を与えることが証明されています。私の最初の(そしてあなたの)側からもう一度この問題に取り組みたいと思います。基本的に、信号は2つのコサイン(1つは順方向、もう1つは逆方向)がサンプル0と1でつなぎ合わされていると言えます。およびも依存するこのg(m)またはf(1/2)の関数またはいくつかの推定値を作成するには、など。α=0.50.25α
2016年

@Retiniteは、サンプルが実際にアドバタイズされた帯域制限された関数を実際にエンコードすることを確認する必要があるため、おそらくそれほど簡単ではありません。
Olli Niemitalo 2016年

証明をありがとう!場合:。これはシーケンス[... 1 0 -1 0 1 1 0 -1 0 1 ...]を与えます。(これは実際にはへのBL ですか?!)これは、周りの級数展開を取得できるように(自動的に)単純化できませんでした。私が得るものは、合計が有限であるかどうか(そしてそれが有限であるかどうか)とその値が何であるかを確認できる、かなり明白な幾何学的系列です。しかし、これはまだ部分的な総当たり法です。α=0.5fm(T)=k=1mmsign(cos(kπαπ4(sign(k12)+1))))sinc(kT)α=0.5mm
2016年

シーケンスがハーフバンド帯域制限関数を表すかどうかをテストできます。フルバンドsincカーネルとハーフバンドsincカーネルによって与えられる「補間」を比較します。あるで2つがとして収束しない場合、答えはノーです。(サンプルポイントでもテストできるため、引用符。)Tm
Olli Niemitalo

1

サンプルに中央のみが含まれている場合でも、秒間隔でサンプルから完全に(補間以降)回復できるフーリエ変換を使用した帯域制限関数を考えてみましょう。sinc関数の他のすべての極大値および極小値をピークおよびミスします。サンプラーが中央のピークを完全に見逃すようにsinc関数を秒遅延させますが、代わりに同じ値の隣接するサンプルを取得します したがって、最大値のオーバーシュートはsinc(t)rect(f)112

sinc(12)=sinπ/2π/2=2π.
12π。証明はありませんが、これは場合の最大のオーバーシュートになると思います。値が小さいほど、サンプルはピーク値近くなり、オーバーシュートもそれに応じて小さくなります。α=1α1
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.