-
大小: 1KB文件類型: .rar金幣: 2下載: 4 次發(fā)布日期: 2021-03-13
- 語言: Matlab
- 標簽: 模態(tài)參數(shù)??MATLAB??polymax法??模態(tài)識別??
資源簡介
頻域內(nèi)的模態(tài)參數(shù)識別方法之多參考點最小二乘復(fù)頻域法,完全自己編寫,希望有同樣做模態(tài)參數(shù)識別的前輩給予指教與幫助
代碼片段和文件信息
%polymax
clc
clear
%%?導(dǎo)入各通道加速度數(shù)據(jù)并繪制時域圖
[NUM]=xlsread(‘19-a30‘3‘A:C‘);
f=2048;%采樣頻率
a=[NUM];
l=3;%通道數(shù)
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)
%%?數(shù)據(jù)分段
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
評論
共有 條評論