カーブフィットを行うプログラムについて
カーブフィットを行うプログラムについて質問があります。
下記は、周波数応答関数を入力し、
カーブフィットを行うプログラムですが、
プログラムの下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);
----------------------
よろしくお願い致します。
失礼します。
お礼
お返事が遅くなりましてすみません。 一度教えていただいたやり方を試してみます。 また質問することがありましたらどうかよろしくお願いします。