河流动力学大作业Solutions to Problem 4

SolutiontoProblem4 Problem ConsiderRun20oftheflumeexperimentsofColeman 1986 findthebestvalueof inthemodifiedsuspensionindexforsatisfactoryagreementbetweenprofilesofsuspendedsedimentconcentrationcalculatedfromtheRouseformulaandmeasureddata Reference Coleman N L 1986 Effectsofsuspendedsedimentontheopen channelvelocitydistribution WaterResourcesResearch AGU 22 10 1377 1384 Contents Appendix 6 UseRouseformulatomatchthemeasureddatainRun20 Herethemodifiedsuspensionindexis Problemabstract where MeasuredData AccordingtoColeman sexperiment Run20 in1986 wecangetsomebasicparametersasbelow Table2 Table1 So Hereweget 2 ZhangRJformula 1 empiricalformula Afteranalysisandcalculation nowtheproblemissimplifiedintoamathematicproblem Generalmodel yi与差异最小时 的值最准确 判定yi与差别的指标 x y C Cb y C Cb x x BecauseoftheprecisionlimitationsinCprogram thisresultisnotveryaccurate Figure2 Basicmethod clanguage s 0 0019 ApreferablewaytocalculatewithMatlab HerewesettheGeneralmodel ThenwechoosetheCurveFittingToolinMatlabtofindthebestvalueof Table6 Thecalculationisalsobasedontheleastsquaresolutionandtheresultislistedbelow Asaresult wefind 0 5209 0 5011 0 5407 Figure3 s 0 001858 Conclusion 通过以上方法 我们完成了对Rouse公式中悬浮指标的修正 即得出式中 的最适值 经计算得 0 5209 References 1 张瑞瑾 河流泥沙动力学 北京 中国水利水电出版社 1989 2 ZhixianCao 2014 LecturenotesofMechanicsofSedimentTransport Wuhanuniversity Appendix Programofdichotomy Clanguage include includevoidmain doublex1 x2 x3 doubles1 12 0 023 0 012 0 0082 0 0061 0 0048 0 0033 0 0026 0 0016 0 00076 0 0004 0 0002 0 00011 doubles2 12 6 12 18 24 30 46 69 91 122 137 152 162 doublea1 12 0 b1 12 0 a2 12 0 b2 12 0 inti j k scanf lf lf i doublec1 0 c2 0 for j 0 j 12 j a1 j pow 170 s2 j 1 170 61 0 44906767 x1 b1 j s1 j 0 023 c1 c1 pow a1 j b1 j 2 for k 0 k 12 k a2 k pow 170 s2 k 1 170 6 1 0 44906767 x2 b2 k s1 k 0 023 c2 c2 pow a