基于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(i,j)=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(i,j)=tline22(3); j=j+1; continue end end fclose(fid22); end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %计算频响函数的实频、虚频、幅频、相频 fs=4096; %采样频率 N=4096; %采样点数。一般FFT程序中所用数据点数与原含有信号数据点数相同， %这样的频谱图具有较高的质量，可减小因补零或截断而产生的影响。 for i=1:3 x=xf(i,:); xfft(i,:)=fft(x,N)*2/N; %FFT变换得到的真实值 y=ya(i,:); yfft(i,:)=fft(y,N)*2/N; end Hfft=yfft./xfft; H=mean(Hfft); H_real=real(H); %频响函数的实频 H_imag=imag(H); %频响函数的虚频 H_f=abs(H); %频响函数的幅频 H_j=atan(H_imag./H_real); %频响函数的相频 f=(0:N-1)*fs/N; %真实频率 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %绘制频响函数的实频、虚频、幅频、相频及Nyqiust图 figure(1); subplot(2,2,1); plot(f(1:N/12.8),H_real(1:N/12.8),'b'); xlabel('频率（Hz）'); ylabel('实部'); title('频响函数的实频曲线'); grid on; subplot(2,2,2); plot(f(1:N/12.8),H_imag(1:N/12.8),'g'); xlabel('频率（Hz）'); ylabel('虚部'); title('频响函数的虚频曲线'); grid on; subplot(2,2,3); plot(f(1:N/12.8),H_f(1:N/12.8),'c'); xlabel('频率（Hz）'); ylabel('幅值'); title('频响函数的幅频曲线'); grid on; subplot(2,2,4); plot(f(1:N/12.8),H_j(1:N/12.8),'r'); xlabel('频率（Hz）'); ylabel('相位'); title('频响函数的相频曲线'); grid on; figure(2); plot(H_real(1:N),H_imag(1:N),'k'); xlabel('实部'); ylabel('虚部'); title('Nyquist图'); grid on; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %以写的方式建立一个新文件，存储频响函数的实部和虚部 fid=fopen(['pinxianghanshu',num2str(n),'.txt'],'w'); %输出频响函数数据 for k=1:N/2; %每行输出2个实型数据，频响函数的实部、虚部 fprintf(fid,'%f %f\n',H_real(k),H_imag(k)); end fclose(fid);