%最小二乘法模態參數識別（復模態-頻響函數實虛部）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
close?all?hidden
format?long
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%聲明全局變量
global?mn
fni=input（‘最小二乘法模態參數識別-輸入數據文件名:‘‘s‘）;
fid=fopen（fni‘r‘）;?
mn=fscanf（fid‘%d‘1）;???????????%模態階數
df=fscanf（fid‘%f‘1）;???????????%頻率間隔
f0=fscanf（fid‘%f‘mn）;??????????%輸入模態頻率初值數組
d0=fscanf（fid‘%f‘mn）;??????????%輸入模態阻尼比初值數組
fno=fscanf（fid‘%s‘1）;??????????%輸入數據文件名
b=fscanf（fid‘%f‘[2inf]）;??????%實測頻響函數實部虛部數據
status=fclose（fid）;
%建立離散頻率向量
f=0:df:（length（b（1:））?-1）*df;
%建立離散圓頻率向量
w=2*pi*f;
%建立實測頻響函數復數向量
H=b（1:）+b（2:）*i;
%計算模態圓頻率初值向量
w0=2*pi*f0;
%建立模態初參數向量
for?j=1:mn
????l=4*（j-1）;
????x0（l+1:l+4）=[-w0（j）*d0（j）w0（j）*sqrt（1-d0（j）^2）11];
end
%用最小二乘非線性數據擬合法估計復模態參數
x=lsqcurvefit（‘f

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

?????文件????????233??2009-09-22?01:07??LS\0ut.txt

?????文件???????1771??2009-09-21?14:51??LS\frequency.m

?????文件????????383??2009-09-17?23:01??LS\fun82.m

?????文件????????138??2009-09-21?14:53??LS\out.txt

?????目錄??????????0??2010-03-17?21:36??LS

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

?????????????????2525????????????????????5