求助高手指点怎样编写绘制幅频特性的程序，谢谢

罗鹏3121

求助高手指点怎样编写绘制幅频特性的程序，谢谢

ChaChing

help bode ???

dingxinran

%Calculate the natural frequencies,the normalized mode shapes and the FRFS
%六自由度系统振动分析的实例
clear
close
a=0.001;b=0.001;
M=[687 0 0 0 0 0;0 687 0 0 0 0;0 0 687 0 0 0;0 0 0 43.5309 -2.9 6.3;0 0 0 -2.9 152.32 -3.7;0 0 0 6.3 -3.7 135.2691];
K=[4171600 33 33 33 33 333;11 2358900 11 111 1111 333;1111 1111 3750000 0 0 333;0 0 0 330300 0 0;0 0 0 0 363700 0;55 0 55 0 55 569600];
C=a*M+b*K;           %比例阻尼
[V,D]=eig(K,M);
for j=1:1:6
omega(j)=sqrt(D(j,j));
for k=1:1:6
Y(k,j)=V(k,j)/V(6,j);
end
end

w=0:0.1:100;               % 频率取值范围和步长
H=zeros(1001,6,6);
for j=1:1:801              %此处801对应的是80，如果为701，则对于的w为70
Z=K-w(j)^2*M+w(j)*C*i;    % 动刚度矩阵 2001 李德葆 实验模态分析及其应用
HH=inv(Z);               % 频响函数矩阵        
H(j,:,:)=HH;
end;
Hr=real(H);Hi=imag(H);    %频响函数的实部和虚部   
Habs=abs(H);Hang=(180/pi)*angle(H); % 频响函数的模为幅频特性
figure(1)
subplot(2,1,1),plot(w,Hr(:,1,1)),grid on
title('实频曲线'),xlabel('\omega'),ylabel('Hr')
subplot(2,1,2),plot(w,Hi(:,1,1)),grid on
title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
figure(2)
subplot(2,1,1),plot(w,Habs(:,1,1)),grid on
title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
subplot(2,1,2),plot(w,Hang(:,1,1)),grid on
title('相频曲线'),xlabel('\omega'),ylabel('\psi')

%画跨点 的频响
figure(4)
subplot(2,1,1),plot(w,Hr(:,2,1)),grid on
title('实频曲线'),xlabel('\omega'),ylabel('Hr')
subplot(2,1,2),plot(w,Hi(:,2,1)),grid on
title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
figure(5)
subplot(2,1,1),plot(w,Habs(:,2,1)),grid on
title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
subplot(2,1,2),plot(w,Hang(:,2,1)),grid on
title('相频曲线'),xlabel('\omega'),ylabel('\psi')

%画跨点 的频响
figure(6)
subplot(2,1,1),plot(w,Hr(:,3,1)),grid on
title('实频曲线'),xlabel('\omega'),ylabel('Hr')
subplot(2,1,2),plot(w,Hi(:,3,1)),grid on
title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
figure(7)
subplot(2,1,1),plot(w,Habs(:,3,1)),grid on
title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
subplot(2,1,2),plot(w,Hang(:,3,1)),grid on
title('相频曲线'),xlabel('\omega'),ylabel('\psi')

%画跨点 的频响
figure(8)
subplot(2,1,1),plot(w,Hr(:,4,1)),grid on
title('实频曲线'),xlabel('\omega'),ylabel('Hr')
subplot(2,1,2),plot(w,Hi(:,4,1)),grid on
title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
figure(9)
subplot(2,1,1),plot(w,Habs(:,4,1)),grid on
title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
subplot(2,1,2),plot(w,Hang(:,4,1)),grid on
title('相频曲线'),xlabel('\omega'),ylabel('\psi')

%画跨点 的频响
figure(10)
subplot(2,1,1),plot(w,Hr(:,5,1)),grid on
title('实频曲线'),xlabel('\omega'),ylabel('Hr')
subplot(2,1,2),plot(w,Hi(:,5,1)),grid on
title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
figure(11)
subplot(2,1,1),plot(w,Habs(:,5,1)),grid on
title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
subplot(2,1,2),plot(w,Hang(:,5,1)),grid on
title('相频曲线'),xlabel('\omega'),ylabel('\psi')

%画跨点 的频响
figure(12)
subplot(2,1,1),plot(w,Hr(:,6,1)),grid on
title('实频曲线'),xlabel('\omega'),ylabel('Hr')
subplot(2,1,2),plot(w,Hi(:,6,1)),grid on
title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
figure(13)
subplot(2,1,1),plot(w,Habs(:,6,1)),grid on
title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
subplot(2,1,2),plot(w,Hang(:,6,1)),grid on
title('相频曲线'),xlabel('\omega'),ylabel('\psi')