基于matlab实现模态分析的代码

function motaifenxi_1（n）
%求频响函数的实频、虚频、幅频、相频及Nyqiust图
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%循环读入激励和响应数据
for i=1:3                           %读入激励数据
    fid11=fopen（[‘D:\matlab2009\work\测点7\‘num2str（n）‘\Time - InputX‘num2str（0）num2str（0）num2str（i）‘.txt‘]‘r‘）;
    j=1;
    while ~feof（fid11）
        tline1=fgetl（fid11）;
        if double（tline1（1））>=48&&double（tline1（1）<=57）
            tline11=str2num（tline1）;
            xf（ij）=tline11（3）;
            j=j+1;
            continue
        end
    end
    fclose（fid11）;
end
for i=1:3                            %读入响应数据
    fid22=fopen（[‘D:\matlab2009\work\测点7\‘num2str（n）‘\Time - InputY‘num2str（0）num2str（0）num2str（i）‘.txt‘]‘r‘）;
    j=1;
    while ~feof（fid22）
        tline2=fgetl（fid22）;
        if double（tline2（1））>=48&&double（tline2（1）<=57）
            tline22=str2num（tline2）;
            ya（ij）=tline22（3）;
            j=j+1;
            continue
        end
    end
    fclose（fid22）;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%计算频响函数的实频、虚频、幅频、相频
fs=4096;                          %采样频率
N=4096;                           %采样点数。一般FFT程序中所用数据点数与原含有信号数据点数相同，
                                  %这样的频谱图具有较高的质量，可减小因补零或截断而产生的影响。
for i=1:3
    x=xf（i:）;