IIRフィルタで疑問があり、過去のページを調べてもわからなかった
ので新たに質問を作りました。
今検討しているのは下記HPの図1のIIRフィルタの1次版(a1,b1,b0のみ)
で伝達関数(Z領域)、サンプリング周波数、各係数が分かっています。
(逆に言うとこれしか分かってません)
ここで上記の資産を流用して新たにサンプリング周波数
だけを変更して使いまわしたいのですが係数をどのように
変更したらよいのか分からずに困っています。
(アナログの伝達特性などはわかりません)
どなたか分かる方、教えていただけないでしょうか。
よろしくお願いします。
http://adsp2191.hp.infoseek.co.jp/misc/009_diffi …
もしくは
http://www.ics.es.yamanashi.ac.jp/~hanawa/asp/pp …(39枚目スライドの右図)
No.4ベストアンサー
- 回答日時:
>係数をカットアンドトライで決めながら元の周波数特性のグラフのような形になるまでひたすら試行するということでしょうか?
一般には「カットアンドトライ」になるのでしょう。
今の例は次数が低いので、三点だけ強引にフィットさせるよう未定係数{a1, b0, b1}を求める、という手もありそう.... 。
EXCEL 使用のシミュレーション例を一つだけ。
H(Z) = (b0+b1*Z)/(1-a1*Z)
を想定。< z^(-1) を Z と略記 >
a0=2.0, b1=b2=0.5
のとき、100kHz サンプリングでカットオフ(3dB ダウン)がほぼ 10kHz 。このとき、25kHz にてほぼ 10dB ダウン。
50kHz サンプリングで 25kHz までほぼ合わせようとすると、
a0=5.0, b1=3.0, b2=1.0
あたりになりそうです。お試しのほどを。
178tallさん、早速の回答ありがとうございます。
ExcelでもIIRのシミュレーションができるんですね!
係数も調べていただいてありがとうございます。
私はMATLABを使用してシミュレーションをやっていて
今ちょっとできない状況なのであとで試してみます。
数回にわたりご解答いただきありがとうございました。
No.5
- 回答日時:
蛇足です。
EXCEL 使用のシミュレーション例をもう一つだけ。
H(Z) = (b0+b1*Z)/(a0+a1*Z) ただし、Z = exp(-jθ)
として Butterworth 一次特性の H(Z) を想定し、カットオフ(3dB ダウン)だけ合わせた場合です。
(1) サンプリングが 100kHz(θ/π=2)、カットオフ(3dB ダウン)が 10kHz(θc/π=0.2) ならば、
a0=-2.08, a1=4.08, b0=b1=1.00
(2) サンプリングが 50kHz(θ/π=2)、カットオフ(3dB ダウン)が 10kHz(θc/π=0.4) ならば、
a0=-0.38, a1=2.38, b0=b1=1.00
-------------------------------
いずれも、ナイキスト周波数(θ/π=1) にて減衰が無限大。
[a0, a1 の算定式]
K=SQRT{(1+cosθc)/(1-cosθc)}
として、
a0=1-K
a1=1+K
178tallさん、お世話になっております。
追加情報ありがとうございました。
a0,a1の算定式に当てはめて計算してみると
確かに178tallさんの係数を求めることが
できました。SIMに関しては後日やって
見ようと思います。
また調べてみるとExcellを使ったIIRやFIR
の設計について書いてある参考書もあるみたい
ですね。今度見てみようと思います。
色々と教えていただきありがとうございました。
No.3
- 回答日時:
>同じ回路を違うサンプリング周波数でも流用したい場合は計算式中のサンプリングの項に代入する値を変更し、再計算することでその特性を得られるIIRの係数が求まると思います。
......>伝達関数が分からず、計算済みのIIRの係数が分かっているのですがここからサンプリング周波数を変えて同じ動作をさせるというのは不可能なのでしょうか?
サンプリング周波数 = 96kHzで使っていた LPF は、ナイキスト周波数(48kHz)を軸として対称な周波数特性ですね。
また、サンプリング周波数 = 48kHz だとそのナイキスト周波数(24kHz)に対して対称な周波数特性になります。
つまり「サンプリング周波数を変えて同じ動作をさせるというのは不可能」じゃありませんか ?
入力スペクトルを 24kHz 以下に制限できるのなら、もとの LPF の 24kHz 以下の周波数特性に近似するという手はありそうですけど。
(もとの IIR 係数に周波数変換を施せば良い、といものではなさそうです。いかがでしょうか)
178tallさん、いつもご解答ありがとうございます。
なるほど。確かにサンプリングを変えると折り返し点が変わって
きてもともとの特性そのままと言うわけにはいかないですね。
今回は96kHzサンプリングで設計されたフィルタ(例えばカットオフ10kHz)
を48kHzサンプリングにして同じく10kHzのカットオフを持たせるという
ことで
入力は24kHz以内に制限されているので問題ないと思います。
>LPF の 24kHz 以下の周波数特性に近似するという手はありそうですけど。
これは係数をカットアンドトライで決めながら元の周波数特性の
グラフのような形になるまでひたすら試行するということでしょうか?
No.2
- 回答日時:
>例えば、サンプリング周波数96kHzで使っていたLPFをカットオフ周波数などの特性を変えずにサンプリング周波数48kHzで使いたい .....
ディジタル・フィルタの周波数特性は周期関数ですから、「全特性を変えずに」サンプリング周波数を下げるのは不可能じゃありませんか ?
たとえば、もとの周波数特性の半周期分に近似するとか.... 。
この回答への補足
通常(かどうかは分かりませんが・・)、デジタルフィルタを設計する際、実現したいアナログでの回路での伝達特性が分かると、そこから双一次変換を用いてデジタルフィルタでの係数を算出します。
ここで計算式の中にサンプリング周波数のパラメータがあり、もし同じ回路を違うサンプリング周波数でも流用したい場合は計算式中のサンプリングの項に代入する値を変更し、再計算することでその特性を得られるIIRの係数が求まると思います。
今回はアナログ回路の伝達関数が分からず、計算済みのIIRの係数が分かっているのですがここからサンプリング周波数を変えて同じ動作をさせるというのは不可能なのでしょうか?
No.1
- 回答日時:
>... サンプリング周波数だけを変更して使いまわしたいのですが係数をどのように変更したらよいのか分からずに困っています。
サンプリング周波数だけ変更し係数をそのままにしておけば、規準化周波数特性は不変です。
係数を変更するとその特性も変わりますから、再計算せねばなりません。
.... って、回答になってないのかな。
ご回答ありがとうございます。
説明が分かりづらくて申し訳ありませんでした。
今回は例えば、サンプリング周波数96kHzで使っていたLPFを
カットオフ周波数などの特性を変えずにサンプリング周波数
48kHzで使いたいということです。
単にサンプリング周波数を変えてしまうと上記例ですとカットオフも
当然低域にずれてきてしまうので係数の変更をしようと考えました。
アナログの伝達関数など分かっていればそこから計算する方法は
分かったのですが、アナログ特性が未知なのでどうしてよいものか・・・。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Visual Basic(VBA) エクセルのマクロで対象ごとにシート分けしてその内容をセルに書き込みたい 9 2022/08/24 13:23
- その他(AV機器・カメラ) デジタル音楽にかけるエフェクトの仕組みについて教えて下さい。 たとえば『イコライザー』のように、周波 3 2023/03/06 18:30
- その他(自然科学) 風車音の測定 3 2023/04/28 07:12
- 物理学 風車から出る音(その②) 8 2023/04/17 12:25
- その他(コンピューター・テクノロジー) PIC16F1シリーズマイコンのNCO機能について 1 2023/04/18 08:41
- クラシック 楽譜の読み方についての質問 10 2022/09/07 15:00
- 物理学 風力発電での音 1 2023/04/16 08:55
- 工学 フィルタの遮断域の損失あるいは、ケーブルの漏話減衰量などのように60〜100dB程度の高損失を測定す 1 2023/07/07 20:50
- 新卒・第二新卒 就活 新卒 自己PR 強み 添削お願いします。 以下の文章で自己PRを書いてみたのですが、どこが競争 3 2023/02/25 21:27
- 数学 フーリエ変換後の負の周波数成分の扱いについて 4 2022/09/03 10:18
このQ&Aを見た人はこんなQ&Aも見ています
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
5元連立方程式を解きたい
-
|(x-a1)/b1|>|(x-a2)/b2|を満た...
-
最小二乗法での指数関数の計算
-
大好きなソフト(めも電卓ver0.9...
-
EXCELで2つの数値のうち大きい...
-
エクセルで二つの数字の小さい...
-
Excelで隣のセルと同じ内容に列...
-
エクセルで最初のスペースまで...
-
PowerPointで表の1つの列だけ...
-
エクセルでオートフィルタのボ...
-
エクセル 文字数 多い順 並...
-
2つのエクセルのデータを同じよ...
-
SUMIFとCOUNTIFを合わせたよう...
-
エクセルのオートフィルタで最...
-
エクセルで、2種類のデータを...
-
エクセルでの複数条件下での標...
-
エクセルで時刻(8:00~20:00)...
-
エクセルの表から正の数、負の...
-
VBAで文字列を数値に変換したい
-
エクセル 同じ値を探して隣の...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
エクセルで1と入力すると〇倍 2...
-
ベクトルの内積を考えるとき、θ...
-
なんでa₁b₁を分けて考えないと...
-
最小二乗法での指数関数の計算
-
5元連立方程式を解きたい
-
Maximaで非可換代数計算を行う方法
-
箱桁の重心位置の公式
-
再投稿ですいません。前の質問...
-
この問題の解き方を教えて欲し...
-
数II・Bの数列の問題です a1+a3+...
-
大好きなソフト(めも電卓ver0.9...
-
エクセルでA1のセルの文字に...
-
統計ソフトRの重回帰分析につい...
-
|(x-a1)/b1|>|(x-a2)/b2|を満た...
-
同時分布関数について。 Pr{a1...
-
IIRフィルタについて
-
エクセルのVBA 条件分岐を...
-
vlookup関数でエラー表示を空白...
-
エクセルのVBAで範囲指定に変数...
-
どれだけ似てるかを統計的に解析
おすすめ情報