%ARMA法模態(tài)參數(shù)識別
%%%%%%%%%%%%%%%%%%%%%
clear
clc
close?all?hidden
format?long
[DATAfileDATApath]=uigetfile（‘zzz.txt‘‘請選擇數(shù)據(jù)文件‘）;
FILENAME=[DATApathDATAfile];
if?~ischar（FILENAME）
msgbox（‘請選擇數(shù)據(jù)文件!‘）;
return;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fni=FILENAME;%input（‘input?data?file?name:‘‘s‘）;
fid=fopen（fni‘r‘）;
mn=124;??????%模態(tài)階數(shù)
%定義輸入實(shí)測數(shù)據(jù)類型：
%ig=1?時域數(shù)據(jù)如沖擊響應(yīng)、自由振動、互相關(guān)函數(shù)、隨機(jī)減量法處理結(jié)果
%ig=2?頻域數(shù)據(jù)如頻響函數(shù)實(shí)部虛部數(shù)據(jù)
ig=1;??????%定義輸入實(shí)測數(shù)據(jù)類型，ig=1（時域數(shù)據(jù)如沖擊響應(yīng)、自由振動、互相關(guān)函數(shù)、隨機(jī)減量法處理結(jié)果）ig=2（頻域數(shù)據(jù)如頻響函數(shù)實(shí)部虛部數(shù)據(jù)）
f=200;???????%ig=1時f為采樣頻率sfig=2時f為頻率間隔df
fno=‘jieguo124.txt‘;?????%輸出數(shù)據(jù)文件名
b=fscanf?（fid‘%f‘[iginf]）;?%實(shí)測時域或頻域數(shù)據(jù)
status=fclose?（fid）;
%建立ARMA模型的階數(shù)（為模態(tài)階數(shù)的2倍）
nm=2*mn;
%組織識別計算所用的時域數(shù)據(jù)及參數(shù)
if?ig==1?%實(shí)測時域數(shù)據(jù)
????%取采樣頻率
????sf=f;
????%取時域數(shù)據(jù)1/2的長度
????n=fix?（length?（b）?/2）;
????%將輸入時域數(shù)據(jù)值給列向量h
????h=b?（1?1:?2*n）‘;
????%計算時間間隔
????dt=1/sf;
????%建立離散頻率向量
????t=0:?dt:?（2*n-1）?*dt;
else???%實(shí)測頻域數(shù)據(jù)
????%頻域率間隔
????df=f;
????%取實(shí)測頻響函數(shù)的長度
????n=length?（b?（1:））;
????%建立離散頻率向量
????f=0:?df:?（n-1）?*df;
????%建立對應(yīng)正負(fù)頻率的實(shí)測頻響函數(shù)向量
????H=b?（1:）‘+b?（2:）‘*i;
????H?（n+1）?=real?（H?（n））;
????H?（n+2:?2*n）?=conj?（H?（n:?-1:?2））;
????%頻響函數(shù)經(jīng)IFFT并取實(shí)部變換成脈沖響應(yīng)函數(shù)
????h=real?（ifft?（H））;
????%建立離散時間向量
????t=linspace?（0?1/df?2*n）;
????%計算時間間隔
????dt=t（2）-t（1）;
end
%時間序列響應(yīng)擬合的ARMA參數(shù)建模
%A和B分別為ARMA模型傳遞函數(shù)的分子和分母系數(shù)向量
[A?B]?=p