カーブフィットを行うプログラムについて質問があります。
下記は、周波数応答関数を入力し、
カーブフィットを行うプログラムですが、
プログラムの下4行の意味が全く分かりません。
(下4行とは、■マークが先頭に書いてあるものです)
どうか何卒ご教授お願い致します。
わかりにくい質問の仕方をしてしまい、申し訳ないです。
----------------------
clear variables global;comgui('close all')
fm='wagawara_ryou10_';
n1=101;
n2=142;
n=0;
for i=n1:n2
fnm=strcat(fm,int2str(i),'_h.txt')
[IIw,Mag,Phi,Q,Q]=textread(fnm,'%f %f %f %f %f','headerlines',3);
n=n+1;
IIxf(:,n)=complex(Mag.*cos(Phi*pi/180),Mag.*sin(Phi*pi/180));
end
%---初期固有振動数データと減衰比データの読み込み.2005.9.13---
fn='wagawara_ryou10_mif_pk_dm.txt'
[peak,zeta]=textread(fn,'%f %f','headerlines',2);
pkze=[peak,zeta];
%---処理開始---
iiplot %伝達関数のプロット
XF(5).po=pkze %初期固有振動数と減衰比を代入。
idcom('e 15 371'); %idcom('e i w') i:帯域幅,w:中心周波数
idcom('est');
idcom('eup .05 .002 -10'); %idcom('eup dstep fstep num i')
■out=id_rm(XF(5),[1 1 1 1]);
■XF(3)=res2xf(out,IIw);
■iicom('IIxhOn');
■[som2,ga2,pbs2,cps2]=res2nor(IIres,IIpo,IDopt);
----------------------
よろしくお願い致します。
失礼します。
No.1ベストアンサー
- 回答日時:
> どうか何卒ご教授お願い致します。
それは無理です.理由は (少なくとも) 2つ.
(1) この言語が何という言語なのか示されていません.
だからこの言語を知らない人には構文の意味さえわからないし,
調べようもありません.
(たまたまこの言語を使っている人がここを見ていれば,
何かコメントがつくかもしれませんが.)
(2) (特に■の行の) どれがユーザ定義関数で,どれが組み込み関数
なのかわかりません.ユーザ定義関数ならば,その定義も示して
いただかないと,それが何をしているかなんて誰にもわかりません.
組み込み関数ならば,(1) の理由で調べようがありません.
説明不足でした。
折角見ていただいたのに、本当に申し訳ございません。
noocyteさんの問いに対する回答を書きましたので
どうかご覧下さい。
----------
(1) この言語が何という言語なのか示されていません.
だからこの言語を知らない人には構文の意味さえわからないし,
調べようもありません.
↑↑
言語はMATLABです。
(2) (特に■の行の) どれがユーザ定義関数で,どれが組み込み関数
なのかわかりません.ユーザ定義関数ならば,その定義も示して
いただかないと,それが何をしているかなんて誰にもわかりません.
↑↑
組み込み関数と入力データを以下に示します。それ以外は、組み込み関数です。
'wagawara_ryou10_'←入力データ
strcat←組み込み関数
textread←組み込み関数
complex←組み込み関数
'wagawara_ryou10_mif_pk_dm.txt'←入力データ
textread←組み込み関数
idcom←組み込み関数
■id_rm←組み込み関数
■res2xf←組み込み関数
■iicom←組み込み関数
■res2nor←組み込み関数
----------
以上です。
返信が遅くなってしまい、大変申し訳ございませんでした。
まだまだ説明不足の箇所があるかと思いますが、
これに回答していただけると本当にありがたいです。
どうかよろしくお願い致します。
失礼します。
No.4
- 回答日時:
「+MATLAB +"Structural Dynamics Toolbox"」で Google 検索
http://www.google.co.jp/search?sourceid=navclien …
してみましたが,日本語ではほとんど情報がないですね.
次のものが多少は参考になるかも.
ANSYSとMATLABを用いた構造系と制御系の同時最適化 ‐スマート構造への適用‐
(2006/03/03 日本機械学会「設計研究会」発表スライド)
http://www.jsme.or.jp/dsd/A-TS12-05/minutes/18/I …
Structural Dynamics Toolbox 5 等のカタログ
http://www.sdtools.com/pdf/datasheet_japan.pdf
MATLAB:適用分野:Test & Measurement:アプリケーション例 (サイバネットシステム)
Structural Dynamics Toolboxを使用した振動解析/制御
http://www.cybernet.co.jp/matlab/solution/applic …
noocyteさんへ
本当に、色々な情報をありがとうございます。
凄くありがたいです。
教えていただいた情報を無駄にしないよう、
一生懸命勉強いたします。
ありがとうございました。
失礼します。
No.3
- 回答日時:
#1 です.
「+MATLAB +id_rm」で Google 検索
http://www.google.co.jp/search?sourceid=navclien …
してみると,id_rm のマニュアル (後述) と,#2 さんのおっしゃっていた
Structural Dynamics Toolbox のマニュアルが見つかりました.
Structural Dynamics Toolbox FEMLink (For Use with MATLAB(R))
http://aertia.com/docs/sdtools/sdt_manual.pdf
んで,色々検索しまくって,■の行にある id_rm のマニュアルの一部だけ訳してみました.
(合ってるかどうかわかりません.制御工学なんて20数年ぶりですし,こんなに検索しまくったのは初めてです.)
●OUT = id_rm(IN, multi)
http://www.sdtools.com/help/id_rm.html
用途:MIMO (多入力多出力) 制御系の最小モデルを作成し,
scaled modal 入出力を得るために reciprocity (相互?) 制約を適用する.
入力:(1) IN:必要なフィールドは下記.
IN.res:留数.
IN.po:極.
IN.idopt:同定オプション.
.FittingModel:Posit,Complex,Normal モード.
.NSNA:センサ/アクチュエータの数.
.Reciprocity:未使用,1 FRF (周波数応答関数?) または真の MIMO.
.Collocated:Reciprocity を使用する場合,collocated FRF の
indices (指数?).
(2) multi:IN.po 内の各極の (重根の) 多重度を表すベクトル (オプション).
出力:OUT:下記のフィールドを持つ.(将来変更されそう.)
.po:(適切な多重度を持つ) 極.
.def:出力 shape (シェープ?) 行列 (CPSI).
.DOF:センサの DOF (自由度?)
.psib:入力 shape 行列 (PSIB).
.CDOF:collocated FRF の指数.
.header:ヘッダ (最大72文字のテキスト×5行).
その他の組み込み関数は,
・res2xf:留数モデルに関連付けられた多項式表現を作成する.
留数と極の集合に対し,対応する FRF (周波数ポイントwで評価される) を生成する.
・res2nor:複素留数を,ノーマルモード留数または比例減衰ノーマルモード形式に近似変換する.
・XF:Database wrapper object のためのユーザ・インタフェース.
・iicom:FRF データ可視化のための UI (ユーザ・インタフェース) コマンド関数.
yasuyasu19 さんが,振動モード解析 (または制御工学) と MATLAB について
どの程度の知識をお持ちなのかわかりませんが,前者について充分な知識がおありなら,
上記の Structural Dynamics Toolbox のマニュアルを読めばわかると思います.
その場合,MATLAB についてわからないことは検索するなり,ここで質問すれば
いいと思います.
しかしそうでない場合 (例えば上の訳を見てもピンとこない,「留数」や「極の多重度」
の意味がわからない) には,ここで誰かに (私は無理です orz) 教えてもらったり,検索
したりしたとしても,失礼ながら (この分野を勉強しない限り) どだい無理です.職場に
この専門分野に詳しい人がいるはずですので,その人に手取り足取り教えてもらうか,
さもなくば「制御工学がわからないのでできません」宣言をして他の人に代わってもらう
しかないと思います.
●参考 (になるかもしれない)
「"モーダルパラメータ" "カーブフィット"」で Google 検索
http://www.google.co.jp/search?sourceid=navclien …
モーダル解析技術 (教材情報)
http://www.tetras.uitec.ehdo.go.jp/db/kyouzai/ky …
振動モード解析の理論と実際 (教材情報)
http://www.tetras.uitec.ehdo.go.jp/db/kyouzai/ky …
古典制御へのMATLABの利用
http://ichiro.maruta.googlepages.com/matlabtutor …
こんにちは!
noocyteさん、色々と調べて頂き、本当に感謝します。
教えて頂いたHPをもとに、調べてみようと思います。
再度質問をするかもしれません。
その際は、どうぞよろしくお願い致します。
本当にありがとうございました。
失礼します。
No.2
- 回答日時:
matlabのスクリプトでしょうか。
見慣れない関数は、Structural Dyanamics Toolboxのもののようですが、私は使ったことがないのでよくわかりません。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# C言語プログラム変更 2 2022/12/21 15:03
- Visual Basic(VBA) ファイル全てを .xlsm に変更したところ、プログラムが途中で落ちてしまっています 17 2022/12/07 12:03
- Visual Basic(VBA) 列と行の名前(重複あり)が交差するセルに、データを入力したい 3 2022/06/12 11:17
- Visual Basic(VBA) 列と行の名前(重複あり)が交差するセルに、データを入力したい 2 2022/06/25 22:42
- C言語・C++・C# [C言語] コメント文字列を無視して、数値データを読み込むプログラム部分について 5 2022/10/05 11:03
- Visual Basic(VBA) VBAが止まります。 2 2022/09/02 14:02
- クラシック 楽譜の読み方についての質問 10 2022/09/07 15:00
- 統計学 t値の計算方法 1 2022/11/29 18:37
- C言語・C++・C# このプログラミング誰か教えてくれませんか 1 2022/06/02 15:27
- 物理学 なぜ、問題文からこの式が出てくるのでしょうか? 5 2023/08/15 15:07
関連するカテゴリからQ&Aを探す
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
組み込みマイコンでのソフトウ...
-
PLCとPCの通信
-
H8S2368マイコンのTCSR_1レジス...
-
PICマイコンのCLRWDT動作について
-
RTCで取得した値をLCDに表示したい
-
PICにおけるパルス波形の検...
-
Macターミナルで実行中のプログ...
-
緯度、経度の 10進法と 60進法...
-
sil_rew_memのメモリアクセスに...
-
VBA kernel32 の意味
-
VB6のwinsockでconnectできない
-
VBSの処理中一旦処理を止めて再...
-
Excelでのセル内容の高速消去方法
-
バックグラウンドのプロセスの...
-
メッセージボックスのボタン名変更
-
win10で、正確な待ち時間の作り方
-
VBAの配列サイズとメモリに関して
-
C言語で、メモリを解放しないで...
-
explorer.exeが異様にメモリを食う
-
サスペンド(休止やスタンバイ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
PLCとPCの通信
-
レジスタマップって何ですか?
-
組み込みマイコンでのソフトウ...
-
アセンブラプログラム
-
SH7144でパルス幅を測定したい。
-
マイクロ秒で精度が出るタイマー
-
機械語のレジスタペアについて
-
シリアル通信でのデータ取りこぼし
-
H8マイコンのプロテクト方法
-
マイコンからSDカード内のテキ...
-
フラッシュメモリへの書き込み
-
TCNTとTIFRの意味
-
RTCで取得した値をLCDに表示したい
-
2つのアナログ信号が一致して...
-
USBハブの自作について
-
グラフィックLCD SG12...
-
古い計量機器からのrs232cを使...
-
マイコンSPI入門
-
printf文で変数のアドレスを表...
-
PICマイコンのCLRWDT動作について
おすすめ情報