% SPECTRA - Бсчейп кейменпх (script file) - Хрплпгйумпт цбумбфйкзт ерйфбчхнузт, вбуей фпх E.A.K. 2000
% To бсчейп кейменпх ухнфбчизке брп фпн Кбфубнп Вбггелз. Бнбрфхчизке уфйт 10/1/2008
% Брп фп бсчейп кейменпх spectra.m хлпрпйей фб еозт:
% б)хрплпгйжей фйт фймет фщн цбумбфйкщн ерйфбчхнуещн гйб кбие мйб брп фйт феууесейт
%   едбцйкет кбфзгпсйет, фйт прпйет рспфейней п Е.Б.К. 2003
% в)еобгей фб брпфелеумбфб уе енб еощфесйкп бсчейп кейменпх
% г)брпфхрщней уе феууесейт гсбцйкет рбсбуфбуейт фб цбумбфб ерйфбчхнуещн 
% ====================================================================================================
% Енбсоз
% ====================================================================================================

% Бсчйкпрпйзуз фщн мефбвлзфщн ейубгщгзт
% ====================================================================================================
im=1;       % ухнфелеуфзт урпхдбйпфзфбт
a=2.3544;   % мегйуфз псйжпнфйб уейумйкз ерйфбчхнуз едбцпхт,A=0.24g (уе m/sec2)
h=1;        % дйпсищфйкпт ухнфелеуфзт гйб брпувеуз дйбцпсз фпх 5%
theta=1;    % ухнфелеуфзт ерйсспзт иемелйщузт
bita=2.5;   % ухнфелеуфзт цбумбфйкзт енйучхузт
q=3;        % ухнфелеуфзт ухмресйцпсбт фзт кбфбукехзт

T1A=0.10;,T2A=0.40; % чбсбкфзсйуфйкет фймет йдйпресйпдщн гйб фп едбцпт Б    
T1B=0.15;,T2B=0.60; % чбсбкфзсйуфйкет фймет йдйпресйпдщн гйб фп едбцпт В
T1C=0.20;,T2C=0.80; % чбсбкфзсйуфйкет фймет йдйпресйпдщн гйб фп едбцпт Г
T1D=0.20;,T2D=1.20; % чбсбкфзсйуфйкет фймет йдйпресйпдщн гйб фп едбцпт Д

% Хрплпгйумпт фймщн цбумбфйкзт ерйфбчхнузт
% ====================================================================================================
jA=0;
jB=0;
jC=0;
jD=0;
for T=0:0.01:4
    if T<T1A 
        jA=jA+1;
        Acc_A(jA)=im*a*(1+((T/T1A)*((h*theta*bita/q))-1)); % Цбумбфйкз ерйфбчхнуз,Едбцпт Б,0<=Ф1
    elseif T>=T1A & T<=T2A 
        jA=jA+1;
        Acc_A(jA)=im*a*h*theta*bita/q; % Цбумбфйкз ерйфбчхнуз,Едбцпт Б,Ф>=Ф1&Ф<=Ф2
    elseif T>T2A
        jA=jA+1;
        Acc_A(jA)=im*a*(h*theta*bita/q)*((T2A/T)^(2/3)); % Цбумбфйкз ерйфбчхнуз,Едбцпт Б,Ф>Ф2
    end
        
    if T<T1B 
        jB=jB+1;
        Acc_B(jB)=im*a*(1+((T/T1B)*((h*theta*bita/q))-1)); % Цбумбфйкз ерйфбчхнуз,Едбцпт В,0<=Ф1
    elseif T>=T1B & T<=T2B 
        jB=jB+1;
        Acc_B(jB)=im*a*h*theta*bita/q; % Цбумбфйкз ерйфбчхнуз,Едбцпт В,Ф>=Ф1&Ф<=Ф2
    elseif T>T2B
        jB=jB+1;
        Acc_B(jB)=im*a*(h*theta*bita/q)*((T2B/T)^(2/3)); % Цбумбфйкз ерйфбчхнуз,Едбцпт В,Ф>Ф2
    end
    
    if T<T1C 
        jC=jC+1;
        Acc_C(jC)=im*a*(1+((T/T1C)*((h*theta*bita/q))-1)); % Цбумбфйкз ерйфбчхнуз,Едбцпт Г,0<=Ф1
    elseif T>=T1C & T<=T2C 
        jC=jC+1;
        Acc_C(jC)=im*a*h*theta*bita/q; % Цбумбфйкз ерйфбчхнуз,Едбцпт Г,Ф>=Ф1&Ф<=Ф2
    elseif T>T2C
        jC=jC+1;
        Acc_C(jC)=im*a*(h*theta*bita/q)*((T2C/T)^(2/3)); % Цбумбфйкз ерйфбчхнуз,Едбцпт Г,Ф>Ф2
    end
    
    if T<T1D 
        jD=jD+1;
        Acc_D(jD)=im*a*(1+((T/T1D)*((h*theta*bita/q))-1)); % Цбумбфйкз ерйфбчхнуз,Едбцпт Д,0<=Ф1
    elseif T>=T1D & T<=T2D 
        jD=jD+1;
        Acc_D(jD)=im*a*h*theta*bita/q; % Цбумбфйкз ерйфбчхнуз,Едбцпт Д,Ф>=Ф1&Ф<=Ф2
    elseif T>T2D
        jD=jD+1;
        Acc_D(jD)=im*a*(h*theta*bita/q)*((T2D/T)^(2/3)); % Цбумбфйкз ерйфбчхнуз,Едбцпт Д,Ф>Ф2
    end
    
end

% Еобгщгз фщн хрплпгйуменщн фймщн цбумбфйкзт ерйфбчхнузт уе бсчейп .table
% H еобгщгз хлпрпйейфбй меущ фзт чсзузт фщн ухнбсфзуещн фщн file I/0 function
% ====================================================================================================
T=0:0.01:4;
z=[T;Acc_A;Acc_B;Acc_C;Acc_D];

fid=fopen('Notes.table','w');
fprintf(fid,'Рйнбкбт фймщн Цбумбфйкзт ерйфбчхнузт\n');
fprintf(fid,'=======================================================================\n');
fprintf(fid,'   T           Rd(T)        Rd(T)        Rd(T)        Rd(T)    \n');
fprintf(fid,'             Едбцпт Б     Eдбцпт В     Едбцпт Г     Едбцпт Д   \n');
fprintf(fid,'  sec        (m/sec2)     (m/sec2)     (m/sec2)     (m/sec2)   \n');
fprintf(fid,'=======================================================================\n');
fprintf(fid,'%6.2f   %10.2f   %10.2f   %10.2f   %10.2f\n',z);
fclose(fid);

% Брпфхрщуз фщн хрплпгйуиенфщн фймщн цбумбфйкзт ерйфбчхнузт уе гсбцйкет рбсбуфбуейт уфп йдйп рбсбихсп 
% ====================================================================================================
subplot(2,2,1)
plot(T,Acc_A)
grid on
title('Spectral values - Soil A','Fontweight','bold')
ylabel('Rd(T) (m/sec2)','FontSize',9,'FontWeight','bold')
xlabel('T (sec)','FontSize',9,'FontWeight','bold')
subplot(2,2,2)
plot(T,Acc_B)
grid on
title('Spectral values - Soil B','Fontweight','bold')
ylabel('Rd(T) (m/sec2)','FontSize',9,'FontWeight','bold')
xlabel('T (sec)','FontSize',9,'FontWeight','bold')
subplot(2,2,3)
plot(T,Acc_C)
grid on
title('Spectral values - Soil C','Fontweight','bold')
ylabel('Rd(T) (m/sec2)','FontSize',9,'FontWeight','bold')
xlabel('T (sec)','FontSize',9,'FontWeight','bold')
subplot(2,2,4)
plot(T,Acc_D)
grid on
title('Spectral values - Soil D','Fontweight','bold')
ylabel('Rd(T) (m/sec2)','FontSize',9,'FontWeight','bold')
xlabel('T (sec)','FontSize',9,'FontWeight','bold')