clear all

R=20*60;         %the update period
H=10;            %the water depth
Et=32.4;
d=1000;
f=50;

alpha=0.504*f-11.2;
        
TL=10*log10(d)+alpha*d*10^(-3);
        
SL=TL+85;
        
It=10^(SL/10)*0.67*10^(-18);
        
Pt=2*pi*H*It;
        
for M=1:500
   
    for i=1:M
        
        Pframe(i)=1.23*Pt*(M-i+1)*(1+1/5+1/5+1/25);
        
    end
    pmax=max(Pframe);
    
    Tactive=Et/pmax;
        
    Ttotal=1.23*(M+M/5);
        
    Tlifetime(M)=Tactive/Ttotal*R/24;
  
end

M=1:500;
semilogy(M,Tlifetime);
ylim([1,10^8])