水灾害课程设计.doc

《水灾害课程设计.doc》由会员分享，可在线阅读，更多相关《水灾害课程设计.doc（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、水灾害课程设计任务：1、 分析1961-2009年，站点每年发生干旱总日数、轻旱日数、中旱日数、重旱日数、特旱日数的年均值，要求用柱状图表示各种干旱类型的年均值，在图表中标注，并对图表进行分析说明。2、 分析1961-2009年，站点每年发生干旱总日数、轻旱日数、中旱日数、重旱日数、特旱日数及其逐年变化情况，要求制作表格来表示各种干旱类型发生日数及各种干旱类型的平均强度，用折线图表示每种干旱类型逐年变化情况，做出趋势线并对图表进行分析说明（变化过程、规律、极值）。3、 分析1961-2009年，站点发生干旱过程数量、每个干旱过程持续时间、起止日期及干旱过程强度，要求用表格表示，并分析出干旱过程

2、持续时间及干旱过程强度的最大最小值。二、课程设计具体数据处理步骤1、对2000-2009年的数据CI值*0.1折算为真实值，然后通过VB编程对数据处理，计算结果输出在excel中。具体编程代码如下：Private Sub Command1_Click()Dim i, j As IntegerDim n As IntegerDim gr, qr, zr, zhr, tr As IntegerDim wgrsum, grsum, qrsum, zrsum, zhrsum, trsum, CIsum As Single gr，qr，zr，zhr，tr分别代表各种干旱类型日数Dim gpj, qpj,

3、 zpj, zhpj, tpj As SingleDim CI(1 To 3653) As SingleDim kaishi, tianshu As IntegerOpen I:Users暮光雨泽Desktop本科生课程CI真实值.txt For Input As #1For i = 1 To 3653Input #1, CI(i)Next iClose #1Dim ExlSet Exl = CreateObject(Excel.Application)Set myworkbook = Exl.Workbooks.Open(I:Users暮光雨泽Desktop本科生课程计算结果输出.xlsx)

4、 打开excel文件等待输入数据Exl.Visible = TrueSet xlsheet = myworkbook.Worksheets(sheet1)If Option1.Value = True Then 判断年数，同时闰年非闰年 kaishi = 1: tianshu = 366 ElseIf Option2.Value = True Then kaishi = 367: tianshu = 365 ElseIf Option3.Value = True Then kaishi = 732: tianshu = 365 ElseIf Option4.Value = True Then

5、kaishi = 1097: tianshu = 365 ElseIf Option5.Value = True Then kaishi = 1462: tianshu = 366 ElseIf Option6.Value = True Then kaishi = 1828: tianshu = 365 ElseIf Option7.Value = True Then kaishi = 2193: tianshu = 365 ElseIf Option8.Value = True Then kaishi = 2558: tianshu = 365 ElseIf Option9.Value =

6、True Then kaishi = 2923: tianshu = 366 ElseIf Option10.Value = True Then kaishi = 3288: tianshu = 365End IfFor i = 1 To tianshu CI(i) = CI(kaishi + i - 1) CIsum = CIsum + CI(i)Next in = 0: gr = 0: qr = 0: zr = 0: zhr = 0: tr = 0For i = 1 To tianshu If CI(i) -0.6 Then 判断干旱类型 n = n + 1: wgrsum = wgrsu

7、m + CI(i) gr = tianshu - n: grsum = CIsum - wgrsum ElseIf CI(i) -1.2 Then qr = qr + 1: qrsum = qrsum + CI(i) ElseIf CI(i) -1.8 Then zr = zr + 1: zrsum = zrsum + CI(i) ElseIf CI(i) -2.4 Then zhr = zhr + 1: zhrsum = zhrsum + CI(i) ElseIf CI(i) = -2.4 Then tr = tr + 1: trsum = trsum + CI(i) End IfNext

8、igpj = grsum / grqpj = qrsum / qrzpj = zrsum / zrIf zhrsum = 0 Or zhr = 0 Then zhpj = 0 Else zhpj = zhrsum / zhrEnd IfIf trsum = 0 Or tr = 0 Then tpj = 0 Else tpj = trsum / trEnd If xlsheet.cells(1, 1) = 日期: xlsheet.cells(1, 2) = 干旱总日数: xlsheet.cells(1, 3) = 轻旱日数: xlsheet.cells(1, 4) = 中旱日数 xlsheet.

9、cells(1, 5) = 重旱日数: xlsheet.cells(1, 6) = 特旱日数: xlsheet.cells(1, 7) = 平均干旱强度: xlsheet.cells(1, 8) = 平均轻旱强度 xlsheet.cells(1, 9) = 平均中旱强度: xlsheet.cells(1, 10) = 平均重旱强度: xlsheet.cells(1, 11) = 平均特旱强度If Option1.Value = True Then 输出结果到excel中 xlsheet.cells(2, 1) = 2000: xlsheet.cells(2, 2) = gr: xlsheet.

10、cells(2, 3) = qr: xlsheet.cells(2, 4) = zr: xlsheet.cells(2, 5) = zhr xlsheet.cells(2, 6) = tr: xlsheet.cells(2, 7) = gpj: xlsheet.cells(2, 8) = qpj: xlsheet.cells(2, 9) = zpj: xlsheet.cells(2, 10) = zhpj xlsheet.cells(2, 11) = tpj ElseIf Option2.Value = True Then xlsheet.cells(3, 1) = 2001: xlsheet

11、.cells(3, 2) = gr: xlsheet.cells(3, 3) = qr: xlsheet.cells(3, 4) = zr: xlsheet.cells(3, 5) = zhr xlsheet.cells(3, 6) = tr: xlsheet.cells(3, 7) = gpj: xlsheet.cells(3, 8) = qpj: xlsheet.cells(3, 9) = zpj: xlsheet.cells(3, 10) = zhpj xlsheet.cells(3, 11) = tpj ElseIf Option3.Value = True Then xlsheet.

12、cells(4, 1) = 2002: xlsheet.cells(4, 2) = gr: xlsheet.cells(4, 3) = qr: xlsheet.cells(4, 4) = zr: xlsheet.cells(4, 5) = zhr xlsheet.cells(4, 6) = tr: xlsheet.cells(4, 7) = gpj: xlsheet.cells(4, 8) = qpj: xlsheet.cells(4, 9) = zpj: xlsheet.cells(4, 10) = zhpj xlsheet.cells(4, 11) = tpj ElseIf Option4

13、.Value = True Then xlsheet.cells(5, 1) = 2003: xlsheet.cells(5, 2) = gr: xlsheet.cells(5, 3) = qr: xlsheet.cells(5, 4) = zr: xlsheet.cells(5, 5) = zhr xlsheet.cells(5, 6) = tr: xlsheet.cells(5, 7) = gpj: xlsheet.cells(5, 8) = qpj: xlsheet.cells(5, 9) = zpj: xlsheet.cells(5, 10) = zhpj xlsheet.cells(

14、5, 11) = tpj ElseIf Option5.Value = True Then xlsheet.cells(6, 1) = 2004: xlsheet.cells(6, 2) = gr: xlsheet.cells(6, 3) = qr: xlsheet.cells(6, 4) = zr: xlsheet.cells(6, 5) = zhr xlsheet.cells(6, 6) = tr: xlsheet.cells(6, 7) = gpj: xlsheet.cells(6, 8) = qpj: xlsheet.cells(6, 9) = zpj: xlsheet.cells(6

15、, 10) = zhpj xlsheet.cells(6, 11) = tpj ElseIf Option6.Value = True Then xlsheet.cells(7, 1) = 2005: xlsheet.cells(7, 2) = gr: xlsheet.cells(7, 3) = qr: xlsheet.cells(7, 4) = zr: xlsheet.cells(7, 5) = zhr xlsheet.cells(7, 6) = tr: xlsheet.cells(7, 7) = gpj: xlsheet.cells(7, 8) = qpj: xlsheet.cells(7

16、, 9) = zpj: xlsheet.cells(7, 10) = zhpj xlsheet.cells(7, 11) = tpj ElseIf Option7.Value = True Then xlsheet.cells(8, 1) = 2006: xlsheet.cells(8, 2) = gr: xlsheet.cells(8, 3) = qr: xlsheet.cells(8, 4) = zr: xlsheet.cells(8, 5) = zhr xlsheet.cells(8, 6) = tr: xlsheet.cells(8, 7) = gpj: xlsheet.cells(8

17、, 8) = qpj: xlsheet.cells(8, 9) = zpj: xlsheet.cells(8, 10) = zhpj xlsheet.cells(8, 11) = tpj ElseIf Option8.Value = True Then xlsheet.cells(9, 1) = 2007: xlsheet.cells(9, 2) = gr: xlsheet.cells(9, 3) = qr: xlsheet.cells(9, 4) = zr: xlsheet.cells(9, 5) = zhr xlsheet.cells(9, 6) = tr: xlsheet.cells(9

18、, 7) = gpj: xlsheet.cells(9, 8) = qpj: xlsheet.cells(9, 9) = zpj: xlsheet.cells(9, 10) = zhpj xlsheet.cells(9, 11) = tpj ElseIf Option9.Value = True Then xlsheet.cells(10, 1) = 2008: xlsheet.cells(10, 2) = gr: xlsheet.cells(10, 3) = qr: xlsheet.cells(10, 4) = zr: xlsheet.cells(10, 5) = zhr xlsheet.c

19、ells(10, 6) = tr: xlsheet.cells(10, 7) = gpj: xlsheet.cells(10, 8) = qpj: xlsheet.cells(10, 9) = zpj: xlsheet.cells(10, 10) = zhpj xlsheet.cells(10, 11) = tpj ElseIf Option10.Value = True Then xlsheet.cells(11, 1) = 2009: xlsheet.cells(11, 2) = gr: xlsheet.cells(11, 3) = qr: xlsheet.cells(11, 4) = z

20、r: xlsheet.cells(11, 5) = zhr xlsheet.cells(11, 6) = tr: xlsheet.cells(11, 7) = gpj: xlsheet.cells(11, 8) = qpj: xlsheet.cells(11, 9) = zpj: xlsheet.cells(11, 10) = zhpj xlsheet.cells(11, 11) = tpjEnd IfEnd Sub如果处理更多的数据，就需要用控件数组了。2、用excel进行绘图分析三、课程设计具体结果及分析：1、任务一的结果及分析：（所绘制图表如图所示）分析：如图所示，南阳站所测地区在这10

21、年内，平均干旱总日数为126.5天，平均轻旱日数为69.4天，平均中旱日数为35.7天，平均重旱日数为14.7天，平均特旱日数为6.7天。年均值代表了这一地区干旱的平均情况，同时也代表了这一地区的平均干旱的发展趋势，由长系列的数据可以对未来的一定时间的干旱情况进行预测。2、任务二对应各种表格：（如下所示）南阳站对应干旱类型发生日数和干旱强度表如上所示，上表列出了从2000-2009年这10年间的干旱状况和强度，为下面绘趋势图做准备。图表一：分析：上图为每年的干旱总日数和对应干旱强度变化趋势，对应过程如图中虚线代表了此组数据的发展趋势。整个过程在线性趋势上下波动，存在干旱日数极大值213（代表此

22、年干旱很持久），极小值40（代表此年干旱很少，大部分时间属于正常情况）；平均干旱强度极大值-0.95，极小值-1.93，干旱强度值越小，干旱程度越严重。图表二：分析：此图为轻旱日数及对应年干旱强度年变化过程，整体过程在虚线上下波动，存在轻旱日数极小值15，极大值151，轻旱日数所占比重越多，干旱程度越弱；存在轻旱强度极小值-0.94，极大值-0.80，值越小对应2004年轻旱比较严重。图表三：分析：此图为中旱日数及对应强度变化过程，整体在虚线上下波动，有中旱日数极小值13，和极大值84；存在中旱强度极大值-1.32，极小值-1.48，说明2002年中旱程度重，2003年有所减轻。图表四： 分析

23、：此图为重旱日数及对应强度变化过程，重旱日数整体范围比较小，在0-30天之内，变化趋势线接近水平，代表未来一段时间重旱有可能稍微增加。存在重旱日期极小值0，极大值28；存在强度最小值0，最大值-2.14，说明2009年干旱程度较大，对应时间段内造成损失较大。图表五：分析：此图为特旱日数及对应强度变化过程，整体趋势呈下降态势，基本维持在0-10天内，除非遇见特别干旱的情况（例如2000年的29天），存在特旱日数极小值0，极大值29，说明2000年遭遇大旱；强度最小值也出现在2000年-2.97，最大值为0。（以上所有图表均认为如果没有重旱和特旱情况，其对应的强度为0。）3、 任务三对应图表：年份

24、发生次数起始日期结束日期天数强度200022000/1/12000/1/2020-13.42000/3/72000/6/1298-178200122001/3/252001/8/6135-180.22001/8/292001/12/12106-117.3200222002/1/102002/3/858-53.42002/8/42002/11/190-65200322003/6/52003/7/127-27.42003/7/302003/8/2426-19.8200422004/4/102004/7/25107-131.32004/10/252004/11/1825-15.4200522005/

25、4/202005/7/577-48.12005/12/132006/1/2645-38.9200632006/3/172006/6/278-702006/8/32006/9/1342-18.52006/9/142006/12/482-74.8200732007/4/112007/6/2778-892007/8/132007/9/1332-16.72007/10/282008/1/267-42.8200832008/6/82008/7/3154-60.92008/10/102008/10/3021-11.12008/11/212009/2/2799-137.8200922009/4/62009/4/2823-13.72009/9/162009/11/2066-65.4统计结果如图所示，包括对应的极大值，极小值。