%[A,B,C,D,E]=textread('tmp.txt','%d%s%f%f')

%[A]=textread('111.txt','%d%s%f%f')
%load 111.txt          %讀入111.txt數(shù)據(jù)文件，該文件為頻率文件或者周期數(shù)據(jù)文件，

syms freq                 %定義變量
syms alfa
syms time
freq = load ('111.txt');     %讀入數(shù)據(jù)文件
n_freq=length(freq);
alfa=freq                    %定義alfa與freq的數(shù)組維數(shù)和數(shù)值形式一致
time=freq
double alfa
for i=1:n_freq                  %如果111.txt文件輸入的是頻率文件，那么要求出對應(yīng)的周期alfa，然后再進(jìn)行，如果是周期文件則不用執(zhí)行該循環(huán)轉(zhuǎn)換
    time(i)=1/freq(n_freq-i+1);
end
damp=0.05                                   %阻尼比
alf_max=0.08                                %地震影響系數(shù)最大值
tg=0.45                                     %卓越周期
r=0.9+(0.05-damp)/(0.5+5*damp);   #p#分頁標(biāo)題#e#        
cof_1=0.02+(0.05-damp)/8;                   %系數(shù)1
cof_2=1+(0.05-damp)/(0.06+1.7*damp);        %系數(shù)2
for i=1:n_freq                              %開始循環(huán)求解標(biāo)準(zhǔn)反應(yīng)譜譜值
    if (time(i)>=0)&(time(i)<0.1)
        alfa(i)=0.45*alf_max+(cof_2-0.45)*alf_max*time(i)/0.1;
    elseif (time(i)>=0.1)&(time(i)<tg)
        alfa(i)=cof_2*alf_max;
    elseif (time(i)>=tg)&(time(i)<(5*tg))
        alfa(i)=((tg/time(i))^cof_2)*alf_max;
    else
        alfa(i)=(cof_2*0.2^r-cof_1*(time(i)-5*tg))*alf_max;
    end
end
alfa                   %顯示最終結(jié)果
out=[time alfa]