polymax 频域内的模态参数识别方法之多参考点最小二乘复频域法

频域内的模态参数识别方法之多参考点最小二乘复频域法，完全自己编写，希望有同样做模态参数识别的前辈给予指教与帮助

%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