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

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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

?????文件?????242176??2011-04-18?19:10??prony算法\擴(kuò)展Prony方法_清華大學(xué)《現(xiàn)代信號(hào)處理》講義_-張賢達(dá).ppt

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

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

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