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アンプPICCPA強磁性常微分方程式モンテカルロ解析odeトランジスタ状態密度インターフェーススイッチング回路ecaljPDS5022DOS定電流半導体シェルスクリプト乱数レベルシフトHP6632Aブレッドボード分散関係温度解析R6452Aトランジスタ技術I2C可変抵抗反強磁性セミナー数値積分確率論偏微分方程式バンド構造非線形方程式ソルババンドギャップ絶縁熱設計シュミットトリガLEDA/Dコンバータ三端子レギュレータLM358ISO-I2CGW近似カオスフォトカプラマフィンティン半径TL431数値微分PC817Cアナログスイッチ直流動作点解析発振回路USBサーボカレントミラー74HC4053パラメトリック解析LDAbzqltyチョッパアンプ量子力学FFT2ちゃんねるアセンブラBSch開発環境電子負荷ブラべ格子イジング模型補間基本並進ベクトル標準ロジック単振り子キュリー温度繰り返しMaxima状態方程式失敗談相対論スピン軌道相互作用FETランダムウォーク熱伝導六方最密充填構造コバルトewidthTLP621GGAQSGW不規則合金位相図抵抗SMPcygwinラプラス方程式スレーターポーリング曲線gfortranスイッチト・キャパシタ詰め回路TLP552三角波格子比熱TLP521条件分岐LM555MCUNE555QNAPマントルテスタ過渡解析FXA-7020ZRダイヤモンドデータロガーガイガー管自動計測Writer509UPSシュレディンガー方程式ブラウン運動awk差し込みグラフ熱力学平均場近似仮想結晶近似VCAfsolve井戸型ポテンシャルVESTA起電力スーパーセルOpenMP第一原理計算ubuntu固有値問題L10構造OPA2277interp12SC1815fccウィグナーザイツ胞面心立方構造フィルタジバニャン方程式ヒストグラム確率論マテリアルデザインspecx.f等高線正規分布PGAフェルミ面非線型方程式ソルバ初期値固定スピンモーメントスワップ領域ルチル構造リジッドバンド模型edeltquantumESPRESSO岩塩構造BaOSIC二相共存ZnOウルツ鉱構造フォノンデバイ模型c/aノコギリ波全エネルギーFSMTeXgnuplotmultiplotハーフメタルCapSense半金属合金結晶磁気異方性Ubuntu文字列入出力TS-110TS-112疎行列Excel直流解析ヒストグラム円周率不規則局所モーメントトラックボールPC等価回路モデルパラメータ・モデルキーボードRealforce三次元マンデルブロ集合フラクタル化学反応重積分縮退日本語最小二乗法関数フィッティングGimpMAS830LHiLAPW熱拡散方程式両対数グラフナイキスト線図負帰還安定性陰解法Crank-Nicolson法P-10クーロン散乱境界条件連立一次方程式片対数グラフEAGLEPIC16F785LMC662トランスシンボルCK1026線種凡例MBEAACircuitグラフの分割軸ラベルifort

最新コメント
リンク

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