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