base.m

很好用的程序

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);
    
    for i=1:M
        if abs(pmax-Pframe(i))<0.0000000000001
            bottleneck(M)=i;
        end
    end
end

M=1:500;
semilogy(M,bottleneck);