台风论文最终成果数学建模免积分.docx

1、台风论文最终成果数学建模免积分 论文题目：台风路径的预测 姓名：陈子浙 专业：电力工程及其自动化0601姓名：王天福 专业：农业工程及其自动化06K1 姓名：吝洪涛 专业：机械工程06k3 摘要台风给广大的地区带来了充足的雨水，成为与人类生活和生产关系密切的降雨系统。但是，台风也总是带来各种破坏，它具有突发性强、破坏力大的特点，是世界上最严重的自然灾害之一。准确及时的台风预报，可以尽大可能减少生命和财产的损失。台风预测是在原有数据的基础上预测将来的数据。台风是由于气压差的作用而形成和行走的，在海平面上影响台风行走路径和风速的因素有很多，例如陆地的影响和海平面的湿度、温度等因素的影响。通过对不同

2、插值方法的比较，结合假设，考虑到台风的运动轨迹应该是时间的连续光滑函数，下面采用三次样条插值模型，对初始点及预测点之间的缺失点进行补充。本文利用分型分布模型预测台风的未来位置，用首尾相接的分段折线将全部数据点一次相连，成为分段变维分型。利用这种方法，算得分维数D，从而预测未来二十四小时的经纬度。其次，本文利用灰色预测模型预测台风的未来风速，利用灰色预测方法建立风速的GM（1，1）模型，从而预测了未来24小时的风速。本文利用Matlab计算得经度的分维数D和纬度的分维数D，并根据模型分别预测出未来24小时的经纬度：时间0:001:002:003:004:005:006:00经度116.79116

3、.582116.374116.166115.959115.752115.545纬度21.977521.935821.894121.852421.810821.769321.7278时间7:008:009:0010:0011:0012:0013:00经度115.339115.133114.928114.722114.517114.313114.108纬度21.686421.64521.603721.562421.521221.4821.4389时间14:0015:0016:0017:0018:0019:0020:00经度113.904113.701113.497113.294113.091112

4、.889112.687纬度21.397921.356921.315921.27521.234221.193421.1527时间21:0022:0023:00经度112.485112.283112.082纬度21.11221.071421.0309根据灰色预测预测出未来24小时的风速，根据Matlab编程计算获得。关键词：分维分布模型 ；灰色预测；插值预测； 一、 问题的提出台风给广大的地区带来了充足的雨水，成为与人类生活和生产关系密切的降雨系统。但是，台风也总是带来各种破坏，它具有突发性强、破坏力大的特点，是世界上最严重的自然灾害之一。准确及时的台风预报，可以尽大可能减少生命和财产的损失。请利

5、用附件中的数据，预测未来24小时内的台风的位置和风速。 二、问题的分析根据某一台风的已知路径，建立两个独立的分段变维分型模型，一个代表已知路径的经度与时间的关系，另一个代表已知路径的纬度与时间的关系。根据这两个模型，即可分别独立地预测未来台风路径的经度和纬度。灰色系统是指既含有己知信息，又含有未知信息的系统，它是介于自色系统和黑色系统之间的一种系统。灰色系统理论是由华中理工大学邓聚龙教授首先提出并加以发展的它研究的是信息不完全的对像，内涵不确定的概念，关系不明确的机制其任务就是要开拓一条新的途径，克服概率统汁的弱点，从杂乱无章、有限的、离散的数据中找出规律.建立灰色系统模型.灰色预测是一种对含

6、有不确定因素的系统进行预测的方法，是根据过去的及现在已知的或非确定的信息.建立一个从过去延伸到将来的灰色模型( GreyModel.简称GM)，从而确定系统未来发展变化的趋势，并为规划、决策提供依据。 三、模型假设1、假设台风的行走路径没有障碍。2、气流涡旋不发生突变性的变化。3、风速与压强相关，但在压强没有太大变化时随时间的变化而变。4、台风的运动路径近似为台风中心点的运动路径。四、 符号的说明t-时间-经度-纬度n-纬度累计和n-经度累计和D-纬度分维数D-经度分维数 五、模型建立与分析一、台风路径的预测1、分形分布的相关公式分形分布可用如下幂指数分布定义 （1）式中：r为特征线度，如时间

7、长度等：N 为与r有关的物体数目，如温度，应力等；c为侍定常数，D 为分维数。在目前一般应用的分形方法中，D为常数这种分形可称为常维分形。它在双对数坐标上是一条直线但是对于非直线的函数关系常准分形就无法处理。为克服这一困难，笔者引入变维分形的概念即分维数D 是特征线度r的函数：D g(r) (2)对于N 与r之间的任一函数关系N=f(r)均可转化为变维分形的形式，为此只须令 (3)即可解得 (4)亦即得到变维分形的形式。2、经纬度的预测如何构造分形模型，是一个需要认真研究和反复探索的问题。经过多次试验和对比，发现可以根据一台风的已知路径，建立两个独立的连续变维分形模型，一个代表已知路径的经度与

8、时间的关系，另个代表已知路径的纬度与时间的关系。根据这两个模型，即可分别独立地预测未来台风路径的经度与度。题中已知数据编号月日时中心位置DD经度/度纬度/度1885125.923-0.9983-12886125.623-0.9971-13887125.323-0.9965-14888125.123-0.9952-15889124.823-0.9946-1.00468810124.623.1-0.994-1.007478811124.423.2-0.994-1.006488812124.323.2-0.994-1.005698813124.223.2-0.9931-1.0051088141242

9、3.2-0.9922-1.0045118815123.823.2-0.9913-1.0041128816123.623.2-0.9905-1.0038138817123.423.2-0.9896-1.0035148818123.223.2-0.9888-1.003315881912323.2-0.9872-0.9989168820122.723.1-0.9872-0.9948178821122.623-0.9863-0.9908188822122.422.9-0.9863-0.9871198823122.322.8-0.9854-0.987820890122.122.8-0.9853-0.98

10、842189112222.8-0.9844-0.988922892121.822.8-0.9819-0.989423893121.422.8-0.9803-0.989824894121.122.8-0.9787-0.990225895120.822.8-0.9779-0.990626896120.622.8-0.9764-0.99127897120.322.8-0.9748-0.99132889812022.8-0.9741-0.991629899119.822.8-0.9741-0.9876308910119.722.7-0.9725-0.988318911119.422.7-0.971-0

11、.9884328912119.122.7-0.971-0.988733891311922.7-0.9695-0.9891348914118.722.7-0.9679-0.9894358915118.422.7-0.9656-0.989736891611822.7-0.9632-0.99378917117.622.7-0.9634-0.9902388918117.522.7-0.9619-0.9905398919117.222.7-0.9612-0.990740892011722.7-0.9613-0.9866418921116.922.6-0.9622-0.974428922116.922.3

12、-0.963-0.9616438923116.922 表一建立连续变维分形模型：为简便计，在特征线度r和物体数目N 之中，最好能有一个按自然数序列1，2，3 取值。经过对比，可以将r按自然数序列取值。令特征线度r代表时间序列，t代表实际时间 表示表1中编号为i的中心位置所对应的时间。由于表1中各中心位置所对应的时间间隔均为l小时，所以可令 (5)式中亦即r=i ( i=1，2，3 ) (6)在两种模型中r值均按(8)式取值。(1)建立代表已知路径的纬度与时间序列关系的连续变维分形摸型。令表示表1中编号为i的中心位置所对应的纬度，将 值取为纬度的累计和，亦即有： (7)与r之问的连续变维分形关系

13、接(1)式处理，其中D值按(4)式计算。D结果列写在表一中。N值见附录。(2)建立代表已知路径的经度与时间序列关系的分段变维分形模型。令 表示表1中编号为i的中心位置所对应的经度，将值取为经度的累计和，亦即有： (8)与 之间的连续变维分形关系按(1)式处理，其中D值按(5)式计算。D结果列写在表一中。N值见附录。 现在讨论如何预测编号为43的下一时刻的分维数之值。在这种情况下，就需要应用已知的台风路径给出的信息，亦即表一中给出的之值。令的总增量为，则平均增量为。 根据平均增量一值，应用最简单的线性外插值公式，可近似的确定第42段及以后的分维数之值： （i42） (9) （10） 又根据 （1

14、1） （12）根据上述的公式计算得预测的经纬度列写在表中。下面给出用matlab绘制出来的图形。原始数据的图形经过插值运算后的图形MATLAB程序：jd= .; %为了节省篇幅省去了纬度的值wd=.;%为了节省篇幅省去了纬度的值j1=116.9116.9116.9;%非线性增加的部分w1=22.622.322;%非线性增加的部分j=linspace(116.9,130.5,500);w=interp1(jd,wd,j,spline); plot(jd,wd,*,j,w,r,j1,w1,r) 预测值的图形总体途径图二、台风风速的预测1灰色预测模型的问题设有原始序列对进行一次累加, 得到生成序列x

15、 (1) ( t)。根据x (1) ( t) , 建立GM (1, 1) 模型, 其白化形式的微分方程为式中a、u 为待定参数, 可根据最小二乘法来确定, 其参数向量这里由此可得时间响应函数 （1）式（1）即为GM (1, 1) 模型。对（1）求导即得还原模型 （2） 或者 （3）2、求解 利用灰色预测原理，我们假设风速构成矩阵y。预测得到的新的风速。据灰色模型的算法关系，可编写matlab程序来预测未来二十四小时的风速。程序如下：function =Greymodel(y)% 本程序主要用来计算根据灰色理论建立的模型的预测值。% 应用的数学模型是 GM(1,1)。% 原始数据的处理方法是一次

16、累加法。y=18181818202323232323232323232320202018181818232323232323202020202020181818181816161616;n=length(y);yy=ones(n,1);yy(1)=y(1);for i=2:n yy(i)=yy(i-1)+y(i);endB=ones(n-1,2);for i=1:(n-1) B(i,1)=-(yy(i)+yy(i+1)/2; B(i,2)=1;endBT=B;for j=1:n-1 YN(j)=y(j+1);endYN=YN;A=inv(BT*B)*BT*YN;a=A(1);u=A(2);t=

17、u/a;i=1:n+24;yys(i+1)=(y(1)-t).*exp(-a.*i)+t;yys(1)=y(1);for j=n+24:-1:2 ys(j)=yys(j)-yys(j-1);endx=1:n;xs=2:n+2;yn=ys(2:n+2);plot(x,y,r,xs,yn,*-b);det=0;for i=2:n det=det+abs(yn(i)-y(i);enddet=det/(n-1);disp(百分绝对误差为：,num2str(det),%);for i=1:24disp(再下个拟合值为,num2str(ys(n+i); end程序运行后的结果为图形如下图所示按预测模型进行

18、计算，并进行累减生成，然后计算原始序列与的绝对误差序列和相对误差序列： (18) (19)1)关联度检验关联系数： （20）关联度： （21）以关联度计算方法得到与原始序列的关联系数，再计算出其关联度，据经验，当=0.5时，关联度r大于0.6就满意了。2)后验差检验原始序列标准差： (22)绝对误差标准差： (23)方差为绝对误差标准差与原始序列标准差之比： (24)小误差概率：评判标准表精度等级好合格勉强合格不合格要求P0.950.800.700.70C0.350.600.650.65若此模型经残差检验、关联度检验、后验差检验都能通过，则可以用所建模型进行预测，否则，需要进行残差修正。由图形

19、观测得，数据的误差较大，这是因为风速给定的比较离散，无法准确的预测的很准，残差较大。风速具体与很多因素有关，所以根据给定数据无法较准确的预测结果。若用拟合的思想进行预测则，借助与Mathematica软件对风速和压强进行拟合，按照对数递减的规律拟合。数据进行列表后如下编程：程序如下：fxx_:=1/x;fyy_:=Logybiao=18.,998.,18.,998.,18.,998.,20.,995.,23.,990.,25.,985.,25.,985.,25.,985.,25.,985.,23.,985.,18.,995.,16.,998.,18.,996.,18.,996.,18.,996

20、.,18.,996.,20.,995.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,20.,990.,20.,990.,18.,992.,18.,992.,18.,992.,18.,992.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,23.,990.,20.,992.,20.,992.,20.,992.,20.,992.,20.,992.,20.,992.,18.,995.,18.,995.,18.,995.,1

21、8.,995.,18.,995.,16.,998.,16.,998.,16.,998.,16.,998.;nb=Tablefxbiaoi,1,fybiaoi,2,i,1,54ft3=Fitnb,1,x,x运行得到结果6.8753 +0.493939 x按照的关系进行拟合，先作变换x=1/x,y=ln y，拟合函数变形为y=lna+bx。按照上述方法得到：a=exp(6.8753) b=0.493939.从而得到风速与压强的一个预测关系：。 风速的大小影响因素很多，本文只预测了与压强的大致关系，其他的有待考证。 七、模型的结果1、 经纬度的预测：编号月日时纬度经度44810021.978116.

22、7945810121.936116.5846810221.894116.3747810321.852116.1748810421.811115.9649810521.769115.7550810621.728115.5551810721.686115.3452810821.645115.1353810921.604114.93548101021.562114.72558101121.521114.52568101221.48114.31578101321.439114.11588101421.398113.9598101521.357113.7608101621.316113.56181017

23、21.275113.29628101821.234113.09638101921.193112.89648102021.153112.69658102121.112112.49668102221.071112.28668102321.031112.082、风速的预测：时间/h01234567风速/m/s17.510817.387417.26517.143417.022616.902716.783716.6655时间/h89101112131415风速/m/s16.548116.431616.315816.200916.086815.973515.86115.7493时间/h1617181920

24、212223风速/m/s15.638415.528215.418815.310215.202415.095314.98914.8835 八、模型评价及改进方向 8.1 模型的优缺点：1. 本文的模型编程简单,可操作性强,易于推广.2. 本文在分析经纬度的时候采用了变维分形分析方法。在预测台风路径的时候获得的结果比较理想。3. 在预测台风风速的时候，使用了灰色预测的方法，这样在影响因素较多的情况下，预测出的结果较为满意。4. 但是在分析台风风速的时候，掌握的资料不多，所以分析起来较为困难，用拟合的方法作为检验手段，分析风速与压强存在一定的关系，成一定的相同趋势。 九、附录nn125.923251

25、.546376.869501.992626.7115751.3138.1875.7161.31000184.51124.2207.71248.2230.91372254.11495.6277.31619300.51742.2323.71865.2346.91987.93702110.53932232.9415.92355.2438.72477.3461.52599.3484.32721.1507.12842.5529.92963.6552.73084.4575.53205598.33325.3621.13445.3643.93565.1666.73684.8689.43804.2712.1392