カーブフィットを行うプログラムについて

質問

カーブフィットを行うプログラムについて質問があります。

下記は、周波数応答関数を入力し、
カーブフィットを行うプログラムですが、
プログラムの下４行の意味が全く分かりません。
（下４行とは、■マークが先頭に書いてあるものです）

どうか何卒ご教授お願い致します。
わかりにくい質問の仕方をしてしまい、申し訳ないです。
----------------------
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.4

noocyte

「+MATLAB +"Structural Dynamics Toolbox"」で Google 検索
http://www.google.co.jp/search?sourceid=navclient-ff&ie=UTF-8&rls=GGGL,GGGL:2006-34,GGGL:ja&q=%2BMATLAB+%2B%22Structural+Dynamics+Toolbox%22

してみましたが，日本語ではほとんど情報がないですね．
次のものが多少は参考になるかも．

ANSYSとMATLABを用いた構造系と制御系の同時最適化 ‐スマート構造への適用‐
(2006/03/03 日本機械学会「設計研究会」発表スライド)
http://www.jsme.or.jp/dsd/A-TS12-05/minutes/18/Ishizuka_2006-03-03ANSYS_MATLAB.pdf

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/applications/t_m/application1.shtml

お礼

noocyteさんへ

本当に、色々な情報をありがとうございます。
凄くありがたいです。
教えていただいた情報を無駄にしないよう、
一生懸命勉強いたします。

ありがとうございました。
失礼します。

回答No.3

noocyte

#1 です．

「+MATLAB +id_rm」で Google 検索
http://www.google.co.jp/search?sourceid=navclient-ff&ie=UTF-8&rls=GGGL,GGGL:2006-34,GGGL:ja&q=%2BMATLAB+%2Bid_rm
してみると，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文字のテキスト×５行)．

その他の組み込み関数は，

・res2xf：留数モデルに関連付けられた多項式表現を作成する．
　留数と極の集合に対し，対応する FRF (周波数ポイントｗで評価される) を生成する．

・res2nor：複素留数を，ノーマルモード留数または比例減衰ノーマルモード形式に近似変換する．

・XF：Database wrapper object のためのユーザ・インタフェース．

・iicom：FRF データ可視化のための UI (ユーザ・インタフェース) コマンド関数．

　yasuyasu19 さんが，振動モード解析 (または制御工学) と MATLAB について
どの程度の知識をお持ちなのかわかりませんが，前者について充分な知識がおありなら，
上記の Structural Dynamics Toolbox のマニュアルを読めばわかると思います．
その場合，MATLAB についてわからないことは検索するなり，ここで質問すれば
いいと思います．

　しかしそうでない場合 (例えば上の訳を見てもピンとこない，「留数」や「極の多重度」
の意味がわからない) には，ここで誰かに (私は無理です orz) 教えてもらったり，検索
したりしたとしても，失礼ながら (この分野を勉強しない限り) どだい無理です．職場に
この専門分野に詳しい人がいるはずですので，その人に手取り足取り教えてもらうか，
さもなくば「制御工学がわからないのでできません」宣言をして他の人に代わってもらう
しかないと思います．

●参考 (になるかもしれない)

「"モーダルパラメータ" "カーブフィット"」で Google 検索
http://www.google.co.jp/search?sourceid=navclient-ff&ie=UTF-8&rls=GGGL,GGGL:2006-34,GGGL:ja&q=%22%E3%83%A2%E3%83%BC%E3%83%80%E3%83%AB%E3%83%91%E3%83%A9%E3%83%A1%E3%83%BC%E3%82%BF%22+%22%E3%82%AB%E3%83%BC%E3%83%96%E3%83%95%E3%82%A3%E3%83%83%E3%83%88%22

モーダル解析技術 (教材情報)
http://www.tetras.uitec.ehdo.go.jp/db/kyouzai/kyouzai_detail.php?stat=3&id=177

振動モード解析の理論と実際 (教材情報)
http://www.tetras.uitec.ehdo.go.jp/db/kyouzai/kyouzai_detail.php?stat=2&id=52

古典制御へのMATLABの利用
http://ichiro.maruta.googlepages.com/matlabtutorial

お礼

こんにちは！
noocyteさん、色々と調べて頂き、本当に感謝します。
教えて頂いたＨＰをもとに、調べてみようと思います。

再度質問をするかもしれません。
その際は、どうぞよろしくお願い致します。

本当にありがとうございました。
失礼します。

回答No.2

rabbit_cat

matlabのスクリプトでしょうか。
見慣れない関数は、Structural Dyanamics Toolboxのもののようですが、私は使ったことがないのでよくわかりません。

回答No.1

noocyte

> どうか何卒ご教授お願い致します。

それは無理です．理由は (少なくとも) ２つ．

(1) この言語が何という言語なのか示されていません．
　　だからこの言語を知らない人には構文の意味さえわからないし，
　　調べようもありません．
　　(たまたまこの言語を使っている人がここを見ていれば，
　　 何かコメントがつくかもしれませんが．)

(2) (特に■の行の) どれがユーザ定義関数で，どれが組み込み関数
　　なのかわかりません．ユーザ定義関数ならば，その定義も示して
　　いただかないと，それが何をしているかなんて誰にもわかりません．
　　組み込み関数ならば，(1) の理由で調べようがありません．

お礼

説明不足でした。
折角見ていただいたのに、本当に申し訳ございません。
noocyteさんの問いに対する回答を書きましたので
どうかご覧下さい。

----------

(1) この言語が何という言語なのか示されていません．
　　だからこの言語を知らない人には構文の意味さえわからないし，
　　調べようもありません．
↑↑
言語はＭＡＴＬＡＢです。

(2) (特に■の行の) どれがユーザ定義関数で，どれが組み込み関数
　　なのかわかりません．ユーザ定義関数ならば，その定義も示して
　　いただかないと，それが何をしているかなんて誰にもわかりません．
↑↑
組み込み関数と入力データを以下に示します。それ以外は、組み込み関数です。

'wagawara_ryou10_'←入力データ
strcat←組み込み関数
textread←組み込み関数
complex←組み込み関数
'wagawara_ryou10_mif_pk_dm.txt'←入力データ
textread←組み込み関数
idcom←組み込み関数
■id_rm←組み込み関数
■res2xf←組み込み関数
■iicom←組み込み関数
■res2nor←組み込み関数

----------

以上です。
返信が遅くなってしまい、大変申し訳ございませんでした。
まだまだ説明不足の箇所があるかと思いますが、
これに回答していただけると本当にありがたいです。
どうかよろしくお願い致します。

失礼します。