冷杉ナイキスト(L次帯域)フィルター設計

スクリプトを開く

この例では,ローパス冷杉ナイキストフィルターの設計法を示します。また,これらのフィルターを,レイズドコサインフィルターおよびルートレイズドコサインフィルターと比較します。これらのフィルターは,デジタル伝送システムのパルス整形に幅広く使用されます。また,内挿/間引きやフィルターバンクにも使用されます。

振幅応答の比較

このプロットは,等リップルナイキストフィルターとレイズドコサインフィルターの振幅応答を示しています。どちらのフィルターも60の次数と0.5のロールオフファクターをもっています。等リップルフィルターには最適等リップル阻止帯域があるので,同じフィルター次数と遷移幅においてより大きな阻止帯域の減衰をもっています。レイズドコサインフィルターは,解析的インパルス応答を切り捨てることで取得され,どのような見方をしても最適とは言えません。

NBand = 4;N = 60;%过滤器订单R = 0.5;%滚边的因素TW = R / (NBand / 2);%过渡带宽f1 = fdesign.nyquist (NBand,“N, TW”N、TW);情商=设计(f1,“equiripple”，“Zerophase”,真的,“SystemObject”,真正的);多项式系数= rcosdesign (R N / NBand NBand,“正常”）;多项式系数=多项式系数/ max (abs(多项式系数))/ NBand;rc = dsp。FIRFilter (“分子”、多项式系数);fvt = fvtool (eq,钢筋混凝土,“颜色”，“白色”）;传奇(fvt“Equiripple尼奎斯特设计的，提出了余弦设计的）;

実際に,この例で同様の減衰を得るには,レイズドコサイン設計の次数を約1400まで増加しなければなりません。

インパルス応答の比較

ここではインパルス応答を比較します。どちらの場合のインパルス応答も,4番目のサンプルごとにゼロになることに注目してください(中央のサンプルを除く)。ナイキストフィルターは,L次帯域フィルターとも呼ばれます。カットオフ周波数がπ/ Lでインパルス応答がL番目のサンプルごとにゼロになるためです。ここでは4次帯域フィルターを示しています。

f1。FilterOrder = 38;eq1 =设计(f1,“equiripple”，“Zerophase”,真的,“SystemObject”,真正的);多项式系数= rcosdesign (R, f1。FilterOrder / NBand NBand,“正常”）;多项式系数=多项式系数/ max (abs(多项式系数))/ NBand;rc1 = dsp。FIRFilter (“分子”、多项式系数);fvt = fvtool (eq1 rc1,“颜色”，“白色”，“分析”，“冲动”）;传奇(fvt“Equiripple奈奎斯特的，提出了余弦的）;标题('脉冲响应，订单=38,Rolloff = 0.5'）;

勾配のある阻止帯域をもつナイキストフィルター

等リップル設計を使用すると,フィルターの阻止帯域の勾配を制御できます。たとえば,次の設計の減衰の勾配は,0,20日および40 dB /(ラジアン/サンプル)になります。

f1。FilterOrder = 52个;f1。乐队= 8;f1。TransitionWidth = . 05;eq1 =设计(f1,“equiripple”，“SystemObject”,真正的);eq2 =设计(f1,“equiripple”，“StopbandShape”，“线性”，.．.“StopbandDecay”, 20岁,“SystemObject”,真正的);eq3 =设计(f1,“equiripple”，“StopbandShape”，“线性”，.．.“StopbandDecay”现年40岁的“SystemObject”,真正的);fvt = fvtool (eq1 eq2 eq3,“颜色”，“白色”）;传奇(fvt“斜率= 0”，的斜率= 20，的斜率= 40 '）

最小位相設計

全体的なナイキストフィルターの最小位相スペクトル因子(周波数領域内の平方根)を設計できます。このスペクトル因子は,整合フィルターアプリケーションにおけるルートレイズドコサインフィルターと同じように使用できます。フィルターの平方根は送信側に配置され,他方の平方根は受信側に配置されます。

f1。FilterOrder = 30;f1。乐队= NBand;f1。TransitionWidth = TW;eq1 =设计(f1,“equiripple”，“Minphase”,真的,“SystemObject”,真正的);多项式系数= rcosdesign (R, N / NBand, NBand);多项式系数=多项式系数/ max(多项式系数)*(1 /(π* NBand) *(π* (r1) - 4 * R));srrc = dsp。FIRFilter (“分子”、多项式系数);fvt = fvtool (srrc eq1,“颜色”，“白色”）;传奇(fvt“最小相位equiripple设计”，.．.“平方根raised-cosine设计”）;

ロールオフファクターの減少

レイズドコサインフィルターの応答は,ロールオフファクターが減少するにつれて改善されます(ここではロールオフ= 0.2の場合が示されています)。これは,インパルス応答の切り捨てに使用される矩形ウィンドウの周波数応答のメインローブが狭いためです。

f1。FilterOrder = N;f1。TransitionWidth = 1;eq1 =设计(f1,“equiripple”，“Zerophase”,真的,“SystemObject”,真正的);R = 0.2;多项式系数= rcosdesign (R N / NBand NBand,“正常”）;多项式系数=多项式系数/ max (abs(多项式系数))/ NBand;rc1 = dsp。FIRFilter (“分子”、多项式系数);fvt = fvtool (eq1 rc1,“颜色”，“白色”）;传奇(fvt“尼奎斯特equiripple设计”，提出了余弦设计的）;

ウィンドウを適用したインパルス応答ナイキスト設計

ナイキストフィルターは,切り捨てられ,ウィンドウが適用されたインパルス応答方法を使用して設計することもできます。これは,レイズドコサイン設計に代わるもう1つの設計です。たとえば,カイザーウィンドウ方法を使用して,次の初期仕様を満たすフィルターを設計できます。

f1。TransitionWidth = TW;kaiserFilt =设计(f1,“kaiserwin”，“SystemObject”,真正的);

カイザーウィンドウ設計が仕様を満たすには,等リップル設計と同じ次数(60)が必要です(これに反して,阻止帯域の仕様を満たすには,1400もの次数のレイズドコサインフィルターが必要であったことを思い出してください)。

kaiserFilt fvt = fvtool (eq,钢筋混凝土,“颜色”，“白色”）;传奇(fvt“Equiripple设计”，.．.提出了余弦设计的，Kaiser窗设计的）;

内挿用ナイキストフィルター

デジタルデータ伝送に加えて,ナイキストフィルターは内挿の目的にも適しています。その理由は,前述のように,Lサンプルごとにゼロのサンプルが得られるためです(中央のサンプル除く)。これには2つの利点があります。これらの利点はポリフェーズ表現を見ると明らかになります。

调频= fdesign.interpolator (4“奈奎斯特”）;kaiserFilt =设计(fm,“kaiserwin”，“SystemObject”,真正的);fvt = fvtool (kaiserFilt,“颜色”，“白色”）;fvt。PolyphaseView =“上”；

ポリフェーズサブフィルター4号はオールパスフィルターで,純粋な遅延になります(FVToolでインパルス応答を選択するか,FVToolでフィルター係数を見てください)。このため,以下が成り立ちます。1 .その乗数は,そのポリフェーズ分岐の効率的な実装につながる乗数を除き,すべてゼロになります。2 .入力サンプルは,フィルターが理想的でなくても変更せずに内挿フィルターに渡されます。