Interesante tutorial sobre matlabDescripción completa
Ejercicios de matlabDescripción completa
Descripción: Basic Matlab Programming
1. La presión atmosférica (p) varía en función de la altura (h) según la siguiente expresión: p=1035*e-0.12h, donde la altura se mide en kilómetros y la presión en milibares. a) Escribir un…Descripción completa
Matlab Book
HW-Chapter HW-Chapter 9
Naoki Mizukami
9-1. Detention basin with 700ft x 700ft with !1 si"e s#ope. $op $op of the berm is at e#e%ation 9 ft. &nf#ow h'"ro(raph is a trian(u#ar shape with a peak "is)har(e of 1*+ )fs at a time of , hr an" with a base time of * hr. $he "etention basin is initia##' empt'. oute the inf#ow h'"ro(raph throu(h the "etention basin an" "etermine the peak outf#ow an" sta(e. ta(e-stora(e re#ationship is )ompute" as fo##ows z
z
0
0
V z / A z / dz 700 z * dz
, z 700 700 *
ta(e-outf#ow re#ationship is (i%en an" summarize" a#on( with stora(e in the tab#e be#ow. Table. Stage-storage-outflow relationships based on given outlet structures and detention basin geometry
Table. reservoir routing using storage-indication method
HW-Chapter 9
Naoki Mizukami
$he peak outf#ow is 2.7 )fs at 19 hr. $he maximum stora(e is , ft at 19 hr. 140
120
100
] 80 s f c [ Q
60
40
20
0 0
5
10
15
20
25
20
25
Time hr
Figure. Inflow and outflow hydrograph 10 9 8 7 6 ] t f [ 5 Z
4 3 2 1 0 0
5
10
15 Time [hr]
Figure. Temporal change in stagee !
HW-Chapter 9
Naoki Mizukami
9-. out the inf#ow h'"ro(raph (i%en usin( Muskin(um metho" with 3t 4 * hr5 3x 4 *+000 ft5 6 4 .2 hr5 4 0.*. 8#ot the inf#ow an" outf#ow h'"ro(raphs an" "etermine the per)ent re"u)tion in the inf#ow peak as we## as the tra%e# time of the peak
HW-Chapter 9
Naoki Mizukami
Figure. Inflow and outflow hydrograph
HW-Chapter 9
Naoki Mizukami
9-12. sin( Muskin(um-Cun(e metho" to route the inf#ow h'"ro(raph an" ri%er in examp#e 9.2 for 0 4 0.001.
Autf#ow is )ompute" usin( the Muskin(um routin( e?uations written b' * 1 * 1 * * "i 1 ( 0 "i ( 1 "i ( * "i 1 * where the "i refers to the f#ow rate at spa)e * an" time i. $he resu#t of routin( is as fo##ows
Table. +us%ingum-(unge routing with constant parameter.
HW-Chapter 9
Naoki Mizukami
Figure. Inflow and outflow hydrographs from +us%ingum-(unge routing.
HW-Chapter 9
Naoki Mizukami
9-1,. Write the )o"e for Muskin(um-Cun(e routin( usin( the point %ariab#e parameter metho" an" app#' it to Bxamp#e 9.2. Mat#ab Co"e Muskin(um-Cun(e with %ariab#e )oeffi)ient metho" Coeffi)ients 5 an" )k wa%e )e#erit'/ are "'nami) in time an" spa)e "ue to %ar'in( ;5 resu#tin( in %ariab#e C05 C15 an" C*. brief "es)ription of point %ariab#e parameter metho" the point metho" )omputes an a%era(e ; usin( inf#ow an" outf#ow at the first step an" next step an" use the a%era(e ; to estimate M-C parameters. $his )omputation is performe" e%er' time outf#ow at the next time step is )ompute". in)e outf#ow at the next time step is unknown5 three point ; are use" for first (uess for a%era(e ;. &teration is ne)essar' to make a%era(e ; )on%er(e" at a )ertain error to#eran)e. )#ear a## )#ose a## ea"in( inf#ow h'"ro(raph "ata1 4 #oa"inh'"ro1,."at/E t 4 "ata1!51/ time FhrG ;!51/ 4 "ata1!5*/ inf#ow F)fsG "t 4 t*/-t1//=20=20E time inter%a# FhrG i%er (eometr' 4 1,000E ri%er rea)h #en(th FftG 0 4 0.000+E be" s#ope Fft
HW-Chapter 9
Naoki Mizukami
for i 4 *!#en(tht/ for K 4 *!#en(thx/ initia# (uess for ; usin( three points inf#ow at 1st time step an" *n" time step an" outf#ow at 1st time step. ;a%e 4 ;i-15K-1/J;i5K-1/J;i-15K//<E iteration ti## point a%era(e ; are )on%er(e" at 0.1 of error to#eran)e whi#e 1 Define mannin( e?uation fun) 4 L'/ Ih
Compute outf#ow at next time step ;i5K/4C0=;i5K-1/JC1=;i-15K-1/JC*=;i-15K/E
$ake a%era(e of points ;a%e14;i-15K-1/J;i5K-1/J;i-15K/J;i5K//<E
&f "ifferen)e between )urrent point a%era(e an" pre%ious one is #ess than 0.1 #ea%e iteration #oop5 otherwise use new a%era(e to )ompute new estimate of outf#ow at next time step if abs;a%e-;a%e1/ 0.1E breakE en" ;a%e 4 ;a%e1E en" en" en" en" of )o"e
HW-Chapter 9
Naoki Mizukami
esu#ts app#ie" this )o"e to examp#e 9.2/ Table. +us%ingum-(unge routing ,inflow and routed hydrograph at ' /// and ' 01///2