PDS5022SとExcelで高速フーリエ変換

オシロスコープPDS5022SとExcelを用いて、回路の測定結果に対してフーリエ変換を行い周波数スペクトルを求めることができます。
本エントリでは、単電源三角波/方形波発振回路で生成した方形波の周波数スペクトルをもとめました。

001_20090908015052.png 001_20090918013724.png 002_20090918013723.png


Excelの分析ツール


PDS5022SとExcelでヒステリシス曲線では、PDS5022Sで取得したデータをPC上のExcel編集してグラフを描画する方法を書きました。

今回は同様にして、取得データに対してFFTを行います。
Excelを用いた具体的なFFTの方法は、(有)ゴッドフット企画さんエクセルを用いたフーリエ変換(FFT)の記事で解説しています。窓関数のかかっていない(矩形窓の)フーリエ変換です。

オシロスコープPDS5022Sで取得したデータに対してFFTをする際には、データ点数を4096個とします。
タイムスパンをどのようにとっても、使えるサンプリング数は変わらないので、見たい周波数近辺の分解能を上げるためには少しコツが必要です。
オシロスコープでデータを取得する段階で、長時間データを取るほど低周波側の分解能があがり高周波側の分解能が下がります。

方形波のFFT


実際にPDS5022SとExcel(とgnuplot)で描いた周波数スペクトルをfig.1に示します。


001_20090918013724.png
fig.1: 実測データから作成した周波数スペクトル


測定した方形波は、単電源三角波/方形波発振回路で生成しました。
周波数スペクトルの比較用にLTspiceによるシミュレーション結果をfig.2に示します。


002_20090918013723.png
fig.2: LTspiceによる周波数スペクトルのシミュレーション


シミュレーションモデルは、LTspiceで三角波/方形波発振回路のタイムスパンを一部変更したものをつかいました。


001_20090908015052.png
fig.3: 実測した回路の回路図

003_20090907230309.png
fig.4: LTspiceシミュレーションのスケマティック


関連エントリ




参考URL




参考文献




フィードバック



にほんブログ村 その他趣味ブログ 電子工作へ

 ↑ 電子工作ブログランキング参加中です。1クリックお願いします。


コメント・トラックバックも歓迎です。 ↓      


 ↓ この記事が面白かった方は「拍手」をお願いします。


tag: LTspice PDS5022 FFT 

comment

Secret

ウェブ参照ありがとうございます

Excelを用いたフーリエ変換を参照しただきありがとうございます。すごい計算なので驚いています。
Excelを用いた科学技術計算を書籍化しました。
署名:Excelで操る!ここまでできる科学技術計算
丸善 出版事業部
こちらもよろしくお願いします。

Re: ウェブ参照ありがとうございます

神足史人さん、こんにちは。

測定データに対する数値処理を、私もやってみたいな、と思いつつ二の足を踏んでいます。ほとんどプログラム言語と変わらない環境は敷居が高いですし、GUIで操作できる(らしいという)Originとかは、私には少々お値段がお高いです。
Excelでこういったことができてしまうなら万々歳です。機会があれば読ませていただきます。

PDS5022

今さらながら、DCにも使えるオシロスコープが個人的に欲しくなってきました。2チャンネル、トリガ入力別というのが魅力的です。帯域25MHzというのは問題になりませんか?

Re: PDS5022

のりたんさん、こんばんは。

> 帯域25MHzというのは問題になりませんか?
PSoCのアナログブロックで扱える程度の周波数の信号なら問題ないと思います。
とはいえ、デジタル信号、たとえばマイコンのクロックに重畳しているジッタを探すとかいった用途には、帯域不足を感じることがあります。

PDS5022がこの価格帯では最高のオシロであることは確かです。ワンランク上の性能を求めるとなると、新品ならお値段がもう一桁上がってしまいそうです。多少ギャンブルになりますが、使用目的によってはヤフオクで中古品を探すといった手もありだと思います。
FC2カウンター
カテゴリ
ユーザータグ

LTspiceAkaiKKRmachikaneyamaScilabKKRPSoCOPアンプCPA強磁性PICモンテカルロ解析常微分方程式odeトランジスタecalj状態密度DOSインターフェース定電流スイッチング回路PDS5022半導体シェルスクリプト乱数レベルシフトHP6632A温度解析分散関係I2Cトランジスタ技術R6452A可変抵抗ブレッドボードセミナーバンドギャップ数値積分確率論反強磁性偏微分方程式バンド構造絶縁熱設計非線形方程式ソルバフォトカプラシュミットトリガLEDLM358カオスISO-I2C三端子レギュレータGW近似A/Dコンバータカレントミラーアナログスイッチ数値微分マフィンティン半径TL431発振回路サーボPC817CUSB直流動作点解析74HC4053補間FFTBSch開発環境パラメトリック解析2ちゃんねるチョッパアンプ量子力学bzqlty電子負荷イジング模型LDA標準ロジックアセンブラ基本並進ベクトルブラべ格子単振り子熱伝導位相図TLP621キュリー温度繰り返し状態方程式MaximaVESTAスイッチト・キャパシタ相対論FETランダムウォークスピン軌道相互作用SMP六方最密充填構造抵抗不規則合金ewidthスレーターポーリング曲線GGAラプラス方程式cygwingfortranQSGW失敗談コバルト条件分岐TLP521テスタLM555Writer509TLP552格子比熱マントルデータロガー自動計測詰め回路ガイガー管ダイヤモンドQNAPMCUFXA-7020ZR過渡解析三角波UPSNE555固有値問題熱力学ブラウン運動フェルミ面awk起電力第一原理計算OpenMPfsolveubuntu最大値xcrysden最小値最適化仮想結晶近似VCA差し込みグラフスーパーセル井戸型ポテンシャル平均場近似シュレディンガー方程式FSMフラクタルOPA2277固定スピンモーメント2SC1815全エネルギー合金multiplotgnuplotc/aTeX結晶磁気異方性interp1ウィグナーザイツ胞初期値マンデルブロ集合疎行列面心立方構造fcc不純物問題非線型方程式ソルバフィルタL10構造PGA半金属二相共存SICZnOウルツ鉱構造BaO重積分クーロン散乱磁気モーメント電荷密度三次元CIF岩塩構造CapSenseノコギリ波デバイ模型ハーフメタル正規分布フォノンquantumESPRESSOルチル構造スワップ領域リジッドバンド模型edelt縮退キーボード軸ラベルグラフの分割凡例トラックボールPC不規則局所モーメント片対数グラフトランス両対数グラフCK1026MAS830L直流解析Excel円周率パラメータ・モデルヒストグラム日本語最小二乗法等価回路モデルGimp線種シンボルTS-110TS-112PIC16F785LMC662化学反応文字列specx.f入出力ifortマテリアルデザインヒストグラム確率論Realforce等高線ジバニャン方程式P-10Ubuntuナイキスト線図Crank-Nicolson法陰解法熱拡散方程式HiLAPWAACircuit連立一次方程式負帰還安定性境界条件EAGLEMBE関数フィッティング

最新コメント
リンク

にほんブログ村 その他趣味ブログ 電子工作へ