m个流atlab交通事故预测.docVIP

交通事故预测问题 问题 采用1982~1998年全国交通事故统计数据作为预测样本，1999~2002年的事故数据作为检验样本。 要用不同的方法得出不同的平面图线，找到最接近1999到2002数据的图形。 年份 事故次数 年份 事故次数 1982 103777 1993 242343 1983 107758 1994 253537 1984 118886 1995 271843 1985 202394 1996 287685 1986 295136 1997 304217 1987 298147 1998 346129 1988 276071 1999 412860 1989 258030 2000 616974 1990 250297 2001 754919 1991 264817 2002 773137 1992 228278 二、模型的建立与求解 1.直接作图。 MATLAB 程序： 得到图形如下： 可以看出，这个图形几乎是没办法得到预测答案的。 2.带常数项的一次线性回归模型。 设 ，（x为表示年份的编号，1982年时x=1） 由此，编写MATLAB程序如下： 得到b1=1.4794e+005,b2=1.0420e+004,则整个模型为： 代入x2=[18,19,20,21]，即1999~2002年，得到事故数为： fy2 =[335510 345930 356350 366770] 该预测结果与检验样本相比明显偏低。 模型图像为： 3.在2的基础上增加一个指数项。 设 ，（x为表示年份的编号，1982年时x=1） 由此，编写MATLAB程序如下： 得到b1= 1.5116e+005,b2= 9.8618e+003,b3= 8.0385e-004,则整个模型为： 代入x2=[18,19,20,21]，即1999~2002年，得到事故数为： fy2 =[ 381500,482000,738400,1418400] 可以看到，新的模型拟合度有所提高，但2001年之后结果明显过高。 模型图像为： 4.综合以上结果重新修正的模型。 设 ，（x为表示年份的编号，1982年时x=1） 由此，编写MATLAB程序如下： 得到b1= 7.9172e+004,b2= -1.8220e+004,b3= 3.3280e-004,b4=1.5484e+005,则整个模型为： 代入x2=[18,19,20,21]，即1999~2002年，得到事故数为： fy2 =[ 402500,486390,604610,767360] 与检验样本对比： 年份 事故次数 预测值 实际值 1999 402500 412860 2000 486390 616974 2001 604610 754919 2002 767360 773137 可以看到，新的模型与检验样本有一定的拟合度。 模型图像为： 三、总结 一般的线性回归是建立在一些基础影响值之上的，在本例中，只有单一的结果数据。单纯根据结果进行线性回归具有较大的不准确性，通过不断地修正模型公式，得到了具有一定拟合度的结果，但如果进行D-W检验（本例中未进行），结果很可能无法让人满意。 综上所述，在解答此类问题时还需要换用不同于线性回归的其他方法。 下面是“十个小故事大道理”不需要的朋友可以下载后编辑删除！！！谢谢！！！ 小故事1、《扁鹊的医术》 　　魏文王问名医扁鹊说：“你们家兄弟三人，都精于医术，到底哪一位最好呢? 　　扁鹊答：“长兄最好，中兄次之，我最差。 　　文王再问：“那么为什么你最出名呢? 　　扁鹊答：“长兄治病，是治病于病情发作之前。由于一般人不知道他事先能铲除病因，所以他的名气无法传出去;中兄治病，是治病于病情初起时。一般人以为他只 能治轻微的小病，所以他的名气只及本乡里。而我是治病于病情严重之时。一般人都看到我在经脉上穿针管放血、在皮肤上敷药等大手术，所以以为我的医术高明， 名气因此响遍全国。 　　大道理：事后控制不如事中控制，事中控制不如事前控制。 　　小故事2、危险的森林里 　　一个人在森林中漫游时，突然遇见了一只饥饿的老虎，老虎大吼一声就扑了上来。他立刻用最快的速度逃开，但是老虎紧追不舍，他一直跑一直跑，最后被老虎逼到了断崖边。 　　站在悬崖边上，他想：“与其被老虎捉到，活活被咬死，还不如跳入悬崖，说不定还有一线生机。” 　　他纵身跳入悬崖，非常幸运地卡在一棵树上。那是长在断崖边的梅树，树上结满了梅子。 　　正在庆幸之时，他听到断崖深处传来巨大的吼声，往崖底望去，原来有一只凶猛的狮子正抬头看着他，狮子的声音使他心颤，但转念一想：“狮子与老虎是相同的猛兽，被什么吃掉，都是一样的。”