function [ww,zhx1,A,m]=plzhx(a,b,dn)
%求解自振頻率ww，振型zhx1---輸入變量：質(zhì)量陣a
%                                     層間剛度矩陣b
%                                     結(jié)構(gòu)層數(shù)dn
%                           返回變量：自振頻率ww，振型zhx1
m=diag(a);
[ik]=kzong(b,dn);

[A,x]=eig(ik,m);
x=diag(sqrt(x));

for i=1:dn
    [x1(i),j]=min(x);
    Agd(:,i)=A(:,j);
    x(j)=max(x)+1;
end
x=x1;
A=Agd;
for i=1:dn
        ww(i)=x(i);
        zhx1(:,i)=A(:,i)/A(dn,i);
end
if abs(min(min(zhx1)))>abs(max(max(zhx1)))
        T=abs(min(min(zhx1)));
    else
        T=abs(max(max(zhx1)));
    end
%繪圖
for i=1:dn
        yi=0:0.25:dn;
        xi=interp1(0:dn,[0,zhx1(:,i)'],yi,'spline');
        subplot(2,dn,i)
        plot([0,zhx1(:,i)'],0:dn,'.',xi,yi);
        ylabel('X');
        xlabel('Z');
        title(['振型圖\omega=',num2str(ww(i)),'']);
        text(zhx1(1,i),dn,'1','HorizontalAlignment','right','FontSize',16)
        axis([-1*T*1.5,T*1.5,0,dn]);
    end