• 締切済み

エクセルを用いた移動平均法による高周波除去について

以下のエクセルデータから得た信号波形から、移動平均法を用いて高周波成分(3kHz以上)を除去した波形を得たいのですが、具体的にどのようにすればいいのか全くわかりません。 VBAを用いたプログラム等でも構いませんので、高周波を除去するやり方をどなたか優しく教えて頂けると嬉しいです。 よろしくお願いします。

みんなの回答

noname#221368
noname#221368
回答No.1

 例えば目標が10Hzだったとします。  波形をフーリエ分解した時を想像すると、10Hzの成分は周期0.1secの正弦(余弦)波なので、0.1secの長さで積分(平均)すれば0になる、と考えます。よって、周期0.1secをサンプリング周波数で割ってやった値が移動平均の区間数だ、という感覚でやると思います。同様な感覚で、10Hz以上の成分もハイカット出来てるはずだと・・・(^^;)。  もっとも上記は目安です。正確にはフーリエ分解してハイカットするべきでしょう。Excelには、フーリエ分解のアドインが付いてるはずです。あまり使いやすくはないようですが。

全文を見る
すると、全ての回答が全文表示されます。

関連するQ&A

  • 移動平均とローパスフィルタについて(matlab)

    入力信号10Hz(サンプリング周波数1000Hz)に対して、 1周期分(100点)のデータを取り出して平均値を算出し、1サンプル移動して同様に100点のデータの平均値を算出する、というような移動平均のプログラムを書きました(下に載せました)。 この移動平均プログラムと同様の結果が得られるようなローパスフィルタを設計したいのですが、方法がよくわかりません。 ご存じの方、教えて頂けると幸いです。 宜しくお願いします。 以下、移動平均プログラムになります。 clear all fs = 1000; % サンプリング周波数 t = 1/fs:1/fs:3; f = 10; % 入力信号の周波数 y = sin(2*pi*f*t); % 入力信号の生成 M = round(fs/f); % 一周期分のサンプル数(10Hzの場合、100点) for i = 1:1:length(y)-M+1 y2(i) = mean(y(i:1:i-1+M)); % 移動平均 end plot(t(M:end),y2);

  • 信号処理による高調波の除去のついて

    現在,生体信号の解析を周波数領域(FFT)で行っているのですがある周波数の高調波により 所望信号の周波数と被ってしまうことがあり推定が行えないという課題に直面しています. (第4高調波まで確認できます) 高長波は基本周波数のおおよそ等倍の位置に立つのですが,高長波は非線形成分ですので 簡単に等倍した周波数の除去を周波数領域で行ってもうまくいかず, また周波数が被っている場合,除去することはできません. そのため何らかの信号処理で時間領域で高調波の除去を行おうと考えているのですが, 具体的な解決方法が思いつきません. 何か時間領域で高調波の除去ができる,または出来そうな信号処理を知っている方がいれば 是非教えていただきたいです. よろしくお願いいたします.

  • FFT分析、成分の最大値

    FFTで周波数成分分析を作っているところなのですが疑問があります 例えば、以下の信号をFFTした場合 振幅0-3V 矩形波 *かなり理想形 周波数 10Khz 10KHz 2.0V  30KHz 0.7V  50KHz 0.4V  ... とかの周波数成分が分析できました 10Kが基本周波数になるので最大振幅になっています その振幅値なのですが、 今、信号振幅が3.0Vだと2.0Vと算出は正しいのでしょうか 最大成分だから。3.0V以上になるような気がします 理想波形なら、何Vが正解でしょうか? もしかして、2.0vではなく倍の4.0Vが正解でしょうか?

  • デジタル信号処理における収録周波数

    不規則な変動の信号波形を収録周波数1kHz(1000Hz)で記録したものがあります。  1)このデータを10個おきに抽出する処理を行って得た波形と,  2)元波形を100Hzでローパスフィルタ処理した波形とは, どのような性質(周波数特性)の差が生じるでしょうか?  また,1)にできる限り近い性質を得るためには,逆に2)のフィルタ周波数を100Hzより変える方がよいのでしょうか? エリアジング現象を考えると50Hzという気もするのですが。

  • 20kHz正弦波をCDで再生するとどうなりますか

    デジタル音響の開発者または詳しい方への質問です。 今更ですが、CDにおける高周波成分の再現性について疑問があります。 サンプリング周波数の1/2以下の周波数しか記録・再生できないというサンプリング定理については、ある程度理解しているつもりです。 また、サンプリング周波数の1/2以上の信号をAD変換するとエリアシングが発生するということも、原理として理解しています。 疑問は44.1kHzで20kHzないしそれに近い周波数の信号をサンプリングして記録し、記録したデータをDA変換した場合、どの程度忠実に原波形が再現できるかということです。 20kHzを44.1kHzでサンプリングすると、1周期に2.205回しかサンプリングできません。サンプリングデータを直線でつないでも、原周波数成分は残るものの、正弦波とはかけ離れた波形になってしまいます。 20kHzまで聞こえる人は稀ですし、聞こえても相当減衰しているので、音楽再生では大きな問題にならないかもしれませんが、13kHzだと大抵の人は聞こえるし、感度もある程度あります。13kHzの場合、1周期のサンプル数は約3.4個になりますが、これでも正弦波とはかなり違います。 仮に、まずまず正弦波に近いと言えるのが1周期5サンプル以上だとすると、44.1kHzで正弦波として記録再生できるのは、せいぜい9kHz程度ということになってしまいますが、なにか考え違いをしているでしょうか? あるいは、高度なデジタル信号処理技術によって、20kHzないしそれに近い周波数も正弦波に近い形で再現できるのでしょうか? サンプリング周波数の1/2に近い高周波信号の記録・再生について図解しているサイト等があればURLをご提示いただけると幸いです。 当方、↓このようなサイトは見つけました。 http://www2.oninet.ne.jp/ts0905/deeg/deeg11.htm

  • 周波数スペクトルについて

    時系列データを周波数スペクトルにして解析すると、同じ波形(同じ物質が通過したときの波形)を比較するとスペクトルの形(中心周波数および大きさ)が違っています。これは何か別の周波数成分が波形に影響させていると考えてよいのでしょうか?

  • 高周波増幅器の出力波形

    任意の入力波形に対し、高周波増幅器の出力波形がどうなるかを 計算で予測したいのですが、どうやって計算したらいいのか、 基本的な考え方のわかる参考書か文献をご紹介いただけませんでしょうか。 いくつか要点を書きますと 1)入力はDCから50 GHz程度までの周波数成分を含む、いわゆる通信信号 2)トータル入力強度は、増幅器の利得が抑圧を受けるレベルである 3)単一周波数入力に対しての利得抑圧出力レベルは典型データがある。 高周波領域でのレベル低下は無視できないものとする。 4)小信号利得の周波数特性は典型データがある。 群遅延もとりあえずデータはあるが無視する。 5)増幅器の内部構造は全くわからない。 6)雑音は無視する また、予測した波形を実測結果と比較する際の注意点がありましたら、お教えいただけるとうれしいです。

  • 隣接平均と移動平均

    計測機器のデータ解析を行っています。 あるデータ解析に隣接平均を行ったとの記述がありました。データにスムージングするために、隣接平均は行ったようです。 そこで質問なのですが、隣接平均とは具体的にどのようなことを行ったのでしょうか?移動平均を行ったのでしょうか? 3点移動平均や5点移動平均ならば定義式は分かるのですが、隣接平均の定義がいまいちわかりません。 どなたかご存知の方がおられましたら、お教えください。 よろしくお願い致します。

  • 三角波の除去

    三角波にのった信号成分だけをうまく取り出す方法は無いでしょうか? 信号成分は、正弦波に近く、三角波の100~1000倍程度の周波数で、振幅が、三角波振幅の1%程度しかありません。 信号である正弦波の数を(FFTを使わず)カウンタなどでカウントしたいのですが、 OPアンプを用いて2回微分した信号をみると、三角波の頂点付近の三角波成分が除去しきれていません。 よろしくお願いします。

  • 筋電の周波数成分について教えて下さい。

    表面電極から差動アンプを通して検出した生の筋電波形をFFT演算してみたところ 全周波数成分は数Hz~数kHzでした。しかし、純粋な筋電成分は何Hz付近なのでしょうか? 負荷と筋電レベルの関係を調べたく、どの周波数成分の電圧レベルを見たら良いか分からないのです。 ちなみに電極固定箇所は前腕部の撓側手根屈筋付近です。 よろしくお願い致します。