5.フーリエ解析

本ページでは、Excelアドイン機能で行えるフーリエ変換の方法について説明します。

なお本サイトでは、Excelで行える フーリエ変換ソフトも無料で公開しています。そちらもご活用ください。


Excelでは、フーリエ変換による時刻歴データの周波数分析が可能です。 また、逆変換も用意されていますので、周波数応答から時刻歴応答への変換も可能です。

Excelのフーリエ変換はFFTのため、データ数は2n個でなければなりません。 また、データ数は最大4096(212)個までとなります。


では、実際にやってみましょう。


フーリエ変換のサンプルとして、次のファイルを使用します。

サンプルファイル:フーリエ変換サンプル.xls


サンプルファイルのフーリエ変換を行う時刻歴応答関数(サンプル関数):

サンプル関数

サンプルファイルの振幅、角速度(次数)のデフォルト設定値:

関数グラフ

サンプリング間隔は0.01sec、サンプリング数は1024個で解析します。

時刻歴データは、下図のように計算しています。


分析ツールメニュー

Excelのフーリエ解析では周波数を計算しませんので、 あらかじめ周波数列を以下の計算に基づき作成します (周波数列はExcelのI列になります)。


周波数の計算

次に、データ/データ分析/フーリエ解析を選択します。


分析ツールメニュー

すると、フーリエ解析ダイアログが開きます。

ここではフーリエ変換を行いますので、入力範囲に「フーリエ変換する数値」の入ったセルを選択し、 出力先に「変換データ」を出力する際の列の先頭セルを選択します。


データ選択

フーリエ変換するセルを選択します。その際、データの数は必ず2n個(n=1~12の整数)とします (最大4096個まで)。今回は1024個のデータで解析をしますので、 “合成”列(ExcelのG列=関数f(x))の数式入力セルを1024個選択します。


データ選択

次に、変換データを出力するための先頭セルを選択します。


データ選択

両方を選択後、“OK”ボタンをクリックすると、計算が始まります。

変換されたデータは複素形式で表示されます。


変換データ

振幅の複素数表示を関数:IMABS(対象セル)によって、絶対値にします。 さらに、Excelではこの絶対値を(サンプル数÷2)で割る必要があります。 この複素数の絶対値に後処理をしたものが、各周波数帯における振幅になります。


解析結果

さて、変換された複素数の絶対値と周波数のセルを用いてグラフを作成します。

変換したデータすべて(1024個)をグラフに取り込むと次のようになります。

今回100Hzまで計算していますが、50Hzを境に対称形となってきます。

フーリエ変換する際に、波形を認識するのに必要な最小の周波数は、 解析上の基本周波数の2倍であるためで、これをサンプリング定理と呼びます。


解析結果

つまり解析可能な最大周波数は、“基本周波数×サンプル数/2”となります。


解析結果

よって、解析結果としては0~50Hzの周波数帯を採用すれば良く、 それ以上のデータについては不要になります。


解析結果

さて、このグラフは、前述した「期待される周波数応答グラフ」と酷似したものとなりました。

このデータと、最初に設定した各正弦波の振幅(ai)と、次数×ω(niω)を比較してみます。


周波数10Hz、40Hz、48Hzで振幅は2と設定していました。 それに対して、上記結果は周波数はおおむね一致しているものの、振幅は若干小さめにでています。

また、周波数20Hzで振幅は5と設定していました。これも前述同様、少し振幅は小さめに出ています。


残念ながら、フーリエ変換したデータが元の関数と完全に一致することはまずありません。 その最大の理由は、Excel等で計算(数値計算)する場合、データが離散化、量子化されているためです。

フーリエ変換で得られたデータは「あくまで近似計算結果である」と、とらえておく必要があります。


また、フーリエ変換で現れる各種エラー (エリアジング、打切り誤差等) には対応できません。 各自でこれらエラーに対する予防を行う必要があります。


フーリエ解析に関する詳細な説明は、機械設計ページ/データ解析/ フーリエ解析を、 参照ください。

また、フーリエ解析に用いられるフーリエ級数展開の数学的な説明は、 基礎科学ページ/数学/フーリエ級数展開 を参照ください。


コメント


役に立った 役に立たなかった

広告

広告募集中
サイズ:150×150

広告募集中
サイズ:150×150

広告に関するお問い合せ

イベント・学会・展示会


学べる探せる設計技術-沐エンジニアリング

Ads by Google

Ads by Google

Ads by Google


技術計算製作所のソフトウエア購入は


Vectorソフトライブラリ/ビジネス

Amazon

技術計算製作所

技術計算製作所

画像をクリックするとpdfが開きます。