ЭвЭ - Лаб 3 / raschet

clear, clc, close all

Qp=[25 25 37.5 37.5 100 100 75 75 100 100 25 25 0];
t=[0 4 4 6 6 10 10 16 16 22 22 24 24];
Qmax=76;
a=1.1;
T=24;
Hg=16;
Hm=26;
kp=0.93;
kpd=63;

Q=Qp*Qmax/100;
Qsum=0;
for N=1.0: 2.0: 11.0
Qsum=Q(1,N)*(t(1,(N+1))-t(1,N))+Qsum;
disp(N)
end
Qch=a*Qsum/T;
[m,~]=min(Q(Q>0));
Qchmin=m*a;
[m,i]=max(Q);
Qchmax=m*a;

Htr=Hg+(Hm-Hg-1)+1;
k=(Hm-Hg)/10000;
m=0:10:100;
Hp=k*m.*m+Hg;
H=k*Qp.*Qp+Hg;
Hmax=Hg+max(H)+1;
Hmin=Hg+min(H)+1;

N=100000*9.81*Qchmax*Hmax/3600/0.95/kp/kpd;