資源簡介
頻域內的模態參數識別方法之多參考點最小二乘復頻域法,完全自己編寫,希望有同樣做模態參數識別的前輩給予指教與幫助
代碼片段和文件信息
%polymax
clc
clear
%%?導入各通道加速度數據并繪制時域圖
[NUM]=xlsread(‘19-a30‘3‘A:C‘);
f=2048;%采樣頻率
a=[NUM];
l=3;%通道數
nt=size(a(:1));
t=?[0:1/f:(nt-1)/f];
plot(ta(:1)‘r‘‘linewidth‘2);grid?on;
set(gca‘FontName‘‘Times?New?Roman‘‘FontSize‘20);
xlabel([‘\fontname{宋體}時間\fontname{Times?New?Roman}/s?‘]‘FontSize‘20)
ylabel([‘\fontname{宋體}加速度\fontname{Times?New?Roman}/m/s^2?‘]‘FontSize‘20)
%%?數據分段
p=867;q=877;%截取位置
an=a(f*p+1:f*q:);
ntn=size(an(:1));
tn=?[0:1/f:(ntn-1)/f];
plot(tnan(:1)‘r‘‘linewidth‘2);grid?on;
set(gca‘FontName‘‘Times?New?Roman‘‘FontSize‘20);
xlabel([‘\fontname{宋體}時間\fontname{Times?New?Roman}/s?‘]‘FontSize‘20)
ylabel([‘\fontname{宋體}加速度\fontname{Times?New?Roman}/m/s^2?‘]‘FontSize‘20)
%%?建立功率譜矩陣
n=f*(q-p);
noverlap=0;
window1=hamming(100);%漢明窗
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????2374??2016-05-09?23:46??polymax.m
-----------?---------??----------?-----??----
?????????????????2374????????????????????1
評論
共有 條評論