function [qx,a]=zhxcyxsh(B,b,ww,dn)
%求解振型參與系數(shù)qx，地震影響系數(shù)a---輸入變量：振型矩陣B，
%                                              質(zhì)量b
%                                              結(jié)構(gòu)層數(shù)dn
%                                    返回變量：振型參與系數(shù)qx，地震影響系數(shù)a

Tg=0.4;amax=0.16;%求解地震影響系數(shù)所需參數(shù),查規(guī)范

r=0.9;n1=0.02;n2=1.0; %一般結(jié)構(gòu)的阻尼參數(shù)

l=diag(ones(dn));b=diag(b);

ww=diag(ww);T=2*pi*inv(ww);T=diag(T);%求解自振周期T

for j=1:dn
    qx(j)=(B(:,j))'*diag(b)/((B(:,j))'*diag(B(:,j))*diag(b));
    if T(j)<=0.1
        a(j)=(n2*amax-0.45*amax)*10*T(j)+0.45*amax;
    elseif T(j)>0.1&T(j)<=Tg
        a(j)=n2*amax;
    elseif T(j)>Tg&T(j)<=(5*Tg)
        a(j)=(Tg/T(j))^r*n2*amax;
    else
        a(j)=(n2*0.2^r-n1*(T(j)-5*Tg))*amax;
    end
end