ノイズをフィルタに通した後の減衰についてこんにちはゲストさん。noiseを作ろうとしています。作成法としては、まず、rundom関数を用いて、平均値0、分散1、サンプル時間を0.00002(サンプリング周波数50kHz)と設定し雑音を作っています。この時の標準偏差の値は1となっています。それを、5.3倍に増幅し、一次のローパスフィルタを通しました。カットオフ周波数は10kHzです。その後、標準偏差を調べてみたところ、約4.0となりました。分散を2にした時、つまり初期値の標準偏差を約1.4にした時、同じ条件でローパスフィルタ後の標準偏差は約5.6となりました。シミュレーション上こうなったのでこの値を信用しようと思っています。今これを計算で求めたいと思っています。何か公式のようなものはあるのでしょうか?ネットで色々調べてみたのですが、分からなかったのでこちらに投稿させていただきました。どなたかご存知の方いらっしゃいましたら、教えていただけると幸いです。19:08:44通報するQNo.5268082noname#137276すぐに回答ほしいです質問者が選んだベストアンサー(1)基本帯域が0~25kHzのわけシャノンのサンプリング定理と離散時間信号のフーリエ変換に関係します。サンプリングするとき、サンプリング周波数を信号の上限周波数の2倍以上にします。逆にいうと、サンプリング周波数の2分の1が信号の上限周波数となります。離散時間信号x(n)のフーリエ変換は、X(exp(jω))=Σ[-∞~∞]x(n)exp(-jωn)となります。式からすぐに、これはωの周期関数でその周期は2πです。そこで、基本帯域は、[-π,π]となりますが、通常は周波数は正の側だけ考えます。ωはサンプリング周波数で規格化した角周波数ですから、2πがサンプリング周波数で、πがその2分の1に対応します。従って、基本帯域は0~25kHzとなるのです。関係を式で表せそうですね。減衰域の傾斜は-6dB/オクターブです。フィルター通過後の雑音の分散が理想ローパスフィルタの場合より大きい理由がわかりました。係数の求め方の方針だけ書きましょう。1次のローパスフィルタの周波数特性はH(exp(jω))=b/{1-a・exp(-jω)}ω=0のときH(exp(jω))=1とすると、b=1-a係数aはカットオフ周波数から求められます。このようにして求められた周波数特性の絶対値の2乗を0からπで積分し、πで割ると、係数が求められます。回答(5)ANo.5sinisorsa補足します。試しに、#4での双一次変換による方法で1次のディジタルフィルタを作って、このフィルタの出力ノイズの標準偏差を求めてみました。結果は、あなたの4.0にかなり近い値になりました。(3.989)多分、matlabで作ったローパスフィルタはアナログフィルタに基づいて双一次変換によりディジタルフィルタに変換する方法のようです。(もっと正確には、カットオフ周波数のプレワープも入れて)この推測に基づいて、補正係数を計算してみました。返事が送れてすみません。sinisorsaさんの言うとおり双一次変換によるディジタルフィルタに変換する方法のようですね。シミュレーション結果も、求めていただいた補正係数でかけた値と同じになりました。フィルタに関しては無知ですが、これから勉強し理解していきたいと思います。ありがとうございました。b1・z^(-1))/(1-a1・z^(-1))となります。分子のzの多項式が0となるzのことを零といいます。分母のzの多項式が0となるzのことを極といいます。零においては、分子が0になりますから、H(z)は0になります。上の式は、極も零も存在する形ですので、極零型といいます。b1=0のときには、零が存在しませんので、全極型といいます。この辺の話は、もっと次数の高い場合について、ディジタル信号処理の大抵の教科書にありますので、調べてみてください。フィルタについて、あまりご存じないとすると、matlabで使ったローパスフィルタがどのようなものだったのかちょっと心配です。noiseの周波数特性がどのようなものなのか。1次フィルタでよいのかなど、心配になります。減衰の傾斜が緩いものでよいのでしょうか。(Fcはカットオフ周波数)を双一次変換でディジタルフィルタに変換する方法があります。この辺の話がわかればよいのですが。(0)通報するANo.3sinisorsa#1,2です。ローパスフィルタを全極型1次フィルタでは、カットオフ周波数が少し大きくなると、高域の減衰量が少なくなります。(ちょっとやってみました)あなたの使ったローパスフィルタは、極零型ですか。14:51:52お礼とても分かりやすい説明ありがとうございます。申し訳ないのですが、全極型と極零型の違いが分かりません。フィルタに関してはあまり知識がなく、ネットで探してみましたがよく分かりませんでした。信号の平均エネルギーは、10/25倍になります。元の信号の平均エネルギーは、振幅を5.3倍しているので、5.3^2=28.09(これは信号の分散に等しい)従って、フィルタ通過後の平均エネルギーは、28.09×10/25=11.236(フィルタ通過後の信号の分散)標準偏差はこれの平方根で、約3.352となります。あなたの使ったローパスフィルタは理想フィルタではないので、カットオフ以上の周波数成分が一部通過していますので、この数字より大きくなると思います。標準偏差は変わると思います。23:17:45お礼お早い返事ありがとうございます。基本帯域が25kHzとなっておりますが、今回ランダム信号のサンプリング周波数は50kHzですが、なぜ基本帯域は25kHzとなるのでしょうか?もしよろしければ教えていただいてもよろしいでしょうか?また、ローパスフィルタのタイプは一次なのですが、どのように傾斜量を考えればよろしいのでしょうか?カットオフ周波数fc×A(指標のようなもの)='通過域のエネルギー'みたいな式があると便利なのですが...。デジタルフィルタでローパスフィルタを作っています。技術翻訳をやっております。電子回路の知識はありませんので、わかりやすく説明してください。今、ローパスフィルタを作っています。そこでわからないことがあります。行き詰ってしまったので質問させてください。(1/(Ts+1))で除去しようと思っています。外乱オブザーバの設計について質問です。Reserved.掲載情報の著作権は提供元企業等に帰属します。