ql.m

2005MCMA题程序

clear,clc;
close all;
tp=2.3125;
a=1.7;
b=2000;
As=2*10^8;
H0=53;
Z0=52.9;
r=0.007;
q=As*r/b;
c=H0-Z0-q;
[T,Z]=ode45('damfail1',[0 tp],Z0);
[t,h]=ode45('HZDY0',[tp,10],H0);
TT=[T',t'];
hh=53*ones(41,1);
h1=[hh',h'];
plot(TT,h1,'linewidth',2);
grid on;
xlabel('Time-----t/hrs');
ylabel('Water Level of Lake Murray-----H/metres');
figure
u=a*h.^0.5;
Ab=b*h;
QZDY0=0.285*u.*Ab;
Q0=r*a*((q+c*exp(-(Z0-Z)/q)).^1.5)*b;
Q=[Q0',QZDY0'];
plot(TT,Q,'linewidth',2);
hold on;
plot(TT,2.6710*10^5*ones(1,length(TT)),'r-.','linewidth',2);
hold off;
grid on;
xlabel('Time-----t/hrs');
ylabel('Flow through the Dam-breach-----Q(t)/ cu.metre/s');