lmh.m

2005MCMA题程序

clear,clc;
k=0.9;
Qm=3.7773*10^5;
L1=0:20:14/5*1000;
L2=14/5*1000:20:14*7/15*1000;
L3=14*7/15*1000:20:14*2/3*1000;
L4=14*2/3*1000:20:14*1000;
u1=5:((3-5)/length(L1)):3;
u23=3:((2-3)/(length(L2)+length(L3)-2)):2;
u4=2:((1-2)/(length(L4)-1)):1;
u=[u1,u23,u4];
L=[L1,L2,L3,L4];
w=3*10^9*ones(1,length(L));
Qcm=w.*(w/Qm+L./(u*k)).^(-1);
plot(L/1000,Qcm,'b','linewidth',2);
grid on;
xlabel('Distance from the dam------L/km');
ylabel('Maximum Flow ------Qcm/cu.metre/s');