%prony算法模態參數識別
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
format?long
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fni=input（‘prony算法模態識別-輸入數據文件名:‘‘s‘）;
fid=fopen?（fni‘r‘）;???????????????%打開只讀文件
mn=fscanf?（fid‘%d‘1）;????????????%模態階數
%定義輸入實測數據類型：
%ig=1時時域數據如沖擊響應，自由振動，相互關函數，隨機減量法處理結果；
%ig=2頻域數據如頻響函數實部虛部數據
ig=fscanf?（fid‘%d‘1）;
%ig=1時f為采樣頻率sf，ig=2時f為采樣間隔df
f=fscanf?（fid‘%f‘1）;
fno=fscanf?（fid‘%s‘1）;???????????%輸出數據文件名
b=fscanf?（fid‘%f‘[iginf]）;??????%實測時域或頻域數據
status=fclose?（fid）;??????????????
nm=2*mn;???????????????????????????%建立復指數函數的階數（為模態階數的兩倍）

%組織識別計算所用的時域數據及參數
if?ig==1???????????????????????????%實測時域數據
????sf=f;??????????????????????????%取采樣頻率
????n=fix（length（b）/2）;????????????%取時域數據二分之一的長度
????h=b（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3995??2011-04-13?22:43??prony算法\prony.m

?????文件?????242176??2011-04-18?19:10??prony算法\擴展Prony方法_清華大學《現代信號處理》講義_-張賢達.ppt

?????目錄??????????0??2011-04-19?20:34??prony算法

-----------?---------??----------?-----??----

???????????????246171????????????????????3