基于LOF-SSA-PSO-LSSVM模型的旅游客流量预测.pdf

D O I:1 0.3 9 6 9/j.i s s n.1 0 0 1-5 3 3 7.2 0 2 3.3.0 3 1*收稿日期:2 0 2 3-0 3-0 7基金项目:2 0 2 2年呼伦贝尔学院科研资助项目(2 0 2 2 F D Y B 0 3);内蒙古自治区高等学校科学技术研究项目(N J Z Y 1 9 2 3 1).作者简介:张娟,女,1 9 8 6-,硕士,讲师;研究方向:经济分析与预测;E-m a i l:h l b e z h j 1 6 3.c o m.基于L O F-S S A-P S O-L S S VM模型的旅游客流量预测*张 娟(呼伦贝尔学院数学与统计学院,0 2 1 0 0 8,内蒙古自治区呼伦贝尔市)摘要:针对月度旅游客流量影响因素复杂,时间序列数据中可能存在异常值和噪声的特点,提出了一个基于局部异常因子(L O F)和奇异谱分析(S S A)的L O F-S S A-P S O-L S S VM预测模型.数据预处理阶段,对原始序列运用X 1 2季节加法模型消除季节性的影响,采用L O F方法检测季节调整后序列的异常值,并使用线性插值和最小二乘支持向量机(L S S VM)方法来修正异常值,新的待预测序列由修正后的季节调整序列与季节因子序列加和得到.预测阶段,首先采用奇异谱分析(S S A)方法来剔除新的待预测序列中的噪声,然后采用P S O-L S S VM方法对其进行预测.以上海2 0 0 4年1月至2 0 1 8年1 2月旅游客流量序列的预测为例,通过与A R I MA、L S S VM等模型比较,表明L O F-S S A-P S O-L S S VM模型能有效提高旅游客流量的预测精度和稳定性.关键词:旅游客流量预测;局部异常因子;最小二乘支持向量机;粒子群寻优;奇异谱分析中图分类号:F 5 9 2;T P 3 9 1 文献标识码:A 文章编号:1 0 0 1-5 3 3 7(2 0 2 3)0 3-0 0 3 1-0 60 引 言旅游业具有“无烟产业”和“永远的朝阳产业”的美称,是非贸易外汇收入的主要来源之一.根据国家旅游局统计数据,2 0 1 8年我国旅游总收入5.9 7万亿元,同比增长1 0.5%,旅游业对全国G D P的综合贡献为9.9 4万亿元,占全国G D P总量的1 1.0 4%.旅游客流量是旅游业发展水平的重要指标,科学合理地预测旅游客流量,对于旅游资源高效利用及地方经济的发展规划具有重要指导意义和实践价值.在旅游客流量的预测方法上,国内外学者已经进行了大量研究.涉及到的预测模型主要有时间序列模型1,2、计量经济模型3和人工智能模型4,5等.由于支持向量机模型(S VM)6在小样本、非线性条件下的明显优势,它已成为非线性时间序列预测的一种重要方法.最小二乘支持向量机(L S S VM)是基于支持向量机的改进模型,被广泛应用于多种预测领域7,8.L S S VM模型中参数的选择直接影响到其预测精度,因此,粒子群算法(P S O)9、布谷鸟算法(C S)1 0等方法被广泛应用于其参数优化.多位学者的实证结果表明,相比较L S S VM方法,基于P S O和C S方法参数寻优的L S S VM模型预测精度更高.在实际预测问题中,由于时间序列的固有复杂性,针对序列的不同特征采用相应的预处理方法便显得尤为重要.对于月度旅游客流量的数据,主要包含3个数据特征:1)由于受到旅游资源和旅游活动的影响,旅游业具有明显的季节性.对于具有季节效应的时间序列,常用的预处理方法为X 1 2季节分解1 1,即剔除季节效应影响后,提取出季节调整后序列做进一步分析.2)对于原始数据而言,由于容易受到大的经济因素(比如金融危机、世博会等事件)影响,使得时间序列数据中包含有一定的异常值,而异常值的出现会增大时间序列的预测误差,因此,异常值检测和修正成为序列预处理的一个重要步骤.常用的方法是局部异常因子算 法1 2(l o c a lo u t l i e r f a c t o r,L O F).3)数据本身可能包含噪声,如果不加处理而直接进行建模预测,模型的预测精度 第4 9卷 第3期2 0 2 3年7月 曲阜师范大学学报J o u r n a l o f Q u f u N o r m a l U n i v e r s i t y V o l.4 9 N o.3J u l y2 0 2 3 也会大大降低.常用的剔除噪声的方法有小波分解1 3,1 4(WT)、经 验 模 态 分 解1 5(e m p i r i c a lm o d ed e c o m p o s i t i o n,EMD)和 奇 异 谱 分 析1 6(s i n g u l a rs p e c t r u ma n a l y s i s,S S A)等.梁小珍等1 7通过对我国航空客流量的预测研究,发现在序列去噪方面,S S A方法的去噪效果明显优于EMD方法.以上研究的主要问题是,对于旅游客流量序列,缺乏综合季节性、异常值、噪声项3种因素的全面研究.基于上述分析,结合月度旅游客流量的数据特征,本文构建了一种以X 1 2季节分解、L O F异常值检测、S S A分解重构为序列预处理方法的P S O参数寻优下的L S S VM预测模型.针对原始数据具有季节效应的特点,首先采用X 1 2季节分解,通过L O F异常值检测方法识别出异常值并进行修正,利用S S A方法剔除了待预测序列中的噪声.最后为了保证模型精度,使用P S O方法对L S S VM中的参数进行寻优,并利用参数结果进行L S S VM预测.同时,将原始单一模型、仅S S A处理的模型、仅L O F修正的模型与L O F-S S A-L S S VM模型的预测结果做对比分析,以期找出更为精确的预测模型.本文结构安排如下:第1节介绍本文使用的主要方法和模型框架;第2节以上海月度旅游客流量为例,进行了实证研究;最后一节总结全文研究结论.1 模型设计本文基于L O F异常值检测的L O F-S S A-P S O-L S S VM模型对上海月度旅游客流量进行预测分析,模型框架如图1,具体步骤如下.(1)对原始序列进行X 1 2季节调整,得到季节因子(S F)和季节调整后序列(S A).(2)采用L O F异常值检测方法针对季节调整后序列(S A)进行异常值检测,得到异常值点所在的位置.(3)使用插值修正和L S S VM预测修正方法对序列中异常值点进行修正,得到修正后的季节调整序列.(4)修正后的季节调整序列与季节因子相加,得到新的待预测序列.(5)通过对4)中的月度数据进行S S A分解和重构,得到去噪后的新的待预测序列.(6)采用P S O-L S S VM方法对去噪后的新的待预测序列进行预测,从而得到预测值.图1 L O F-S S A-P S O-L S S VM预测模型框架2 模型预测效果检验2.1 数据来源及评价准则本文采用2 0 0 4年1月至2 0 1 8年1 2月上海市月度旅游客流量(即接待过夜游客人数)为研究数据,共1 8 0个观测值,数据来自于 上海统计年鉴.由图2可以看出,上海月度旅游客流量呈整体上升趋势,且具有明显的季节变动规律.某些月份的客流量数值明显异于周边数据,可能为异常值所在位置.为了保证实证结果的可靠性,根据模型要求,本文进行了3组滚动预测.第一组以2 0 0 4年1月2 0 1 5年1 2月数据(1 4 4个观测值)为训练样本对模型进行训练,2 0 1 6年1月1 2月的1 2个数据作为测试样本来检验模型预测精度.第二组以2 0 0 4年1月2 0 1 6年1 2月数据(1 5 6个观测值)预测2 0 1 7年1月1 2月的数据.第三组以2 0 0 4年1月2 0 1 7年1 2月数据(1 6 8个观测值)预测2 0 1 8年1月1 2月的数据.图2 2 0 0 4年1月2 0 1 8年1 2月上海旅游客流量序列图23 曲阜师范大学学报(自然科学版)2 0 2 3年为了评价不同模型之间的预测准确度,本文采用平均 绝 对 误 差(MA E)和 平 均 绝 对 百 分 误 差(MA P E)作为评价指标.2.2 预测效果检验2.2.1 L O F检测针对上海旅游客流量具有明显季节性和趋势性的特点,为了避免周期性的季节变化对模型预测准确度的影响,首先采用X 1 2季节加法模型对原序列进行分解,得到了季节调整后序列(S A)和季节因子(S F).其次依据季节调整后序列(S A)考察数据的异常变动,采用L O F异常值检测的方法,取k=1 0,计算每月客流量数据的L O F值,并以L O F临界值为1.5判断数据的异常性.图3为1 8 0个训练数据的L O F值,根据约定的阈值共检验出1 6个异常值,异常数据对应的具体位置如表1所示.图3 训练样本的L O F1 0值序列 表1 L O F异常值检测结果月份L O F1 0月份L O F1 02 0 0 4-1 21.5 2 32 0 1 0-1 04.4 5 72 0 0 5-0 21.7 1 12 0 1 0-1 12.3 3 92 0 0 5-0 61.5 1 92 0 1 1-0 12.6 2 02 0 0 6-0 12.3 3 42 0 1 3-0 52.6 0 62 0 0 9-0 82.1 3 92 0 1 4-0 12.0 5 92 0 1 0-0 21.7 1 32 0 1 5-0 22.1 1 72 0 1 0-0 42.4 1 72 0 1 6-0 21.9 8 82 0 1 0-0 95.8 8 92 0 1 8-0 22.5 6 5 为了考察1 6个异常值出现的原因,结合特殊事件对上海旅游客流量的影响,经过分析可得,2 0 0 9-0 8的异常数据受到了当年金融危机的影响;2 0 1 0-0 4、2 0 1 0-0 9、2 0 1 0-1 0、2 0 1 0-1 1的异常数据受到了上海世博会的影响;2 0 0 5-0 2、2 0 0 6-0 1、2 0 1 0-0 2、2 0 1 1-0 1、2 0 1 4-0 1、2 0 1 5-0 2、2 0 1 6-0 2、2 0 1 8-0 2这几年2月份前后的异常数据,主要是受到旅游周期的影响;2 0 0 4-1 2、2 0 0 5-0 6、2 0 1 3-0 5的异常数据可能是随机干扰所造成的.异常数据的存在会影响预测模型中的参数估计,从而增大预测误差,因此对其进行修正是非常必要的.2.2.2 线性插值及L S S VM修正季节调整序列根据表1提取出的异常数据位置,对季节调整序列进行修正,具体方法为:对于单一的异常数据点,比如2 0 0 5-0 6、2 0 0 9-0 8、2 0 1 3-0 5等,采用临近点线性插值法来修正异常数据;对于成片的异常数据点,比如2 0 1 0-0 9、2 0 1 0-1 0、2 0 1 0-1 1等,以当前异常值之前的数据序列为基础,使用L S S VM方法进行预测来修正异常数据.修正前后的数据对比如图4,可以看出,相比较原始序列,异常值修正后的季节调整序列具有更好的平稳性和规律性.图4 L O F修正前后的季节调整序列2.2.3 S S A对新的待预测序列进行降噪将异常值修正后的季节调整序列与之前得到的季节因子序列(S F)相加,即得到了新的待预测序列.在剔除了原始序列的异常值之后,对于新的待预测序列中的噪声项,使用S S A分解和重构方法来完成去噪.S S A方法的一个重要步骤是嵌入维数L的选择,由于L一般不超过序列长度的1/3,且最好为周期的整数倍,以训练集2 0 0 4年1月2 0 1 5年1 2月(共1 4 4个观测值)为例,选择嵌入维数L=4 8.图5是前2 0个特征值对应的特征值曲线图,从特征值变化趋势中可以看出,从第1 5个特征值开始,曲线下降速度缓慢,且经过计算前1 5个分量的方差累积贡献率达到9 6.1 5%.所以,提取前1 5个分量对原始序列进行重构,以得到新的重构后序列.剩余的重构分量构成原始序列的噪声序列.图6显示了原始旅游客流量序列和去噪后的新的待预测序列的对比结果,可以看出,经过异常值修正和降噪处理后的数据,序列更加光滑.33第3期 张娟:基于L O F-S S A-P S O-L S S VM模型的旅游客流量预测 图5 S S A重构分量对应的特征值 图6 原始旅游客流量序列和去噪后的新的待预测序列2.2.4 模型预测效果对比为了验证L O F-S S A-P S O-L S S VM方法预测效果的有效性,针对2 0 0 42 0 1 5年、2 0 0 42 0 1 6年、2 0 0 42 0 1 7年3组训练集,本文对比了A R I MA、L S S VM、P S O-L S S VM这3个单一模型、仅去噪的单一模型、仅异常值修正的单一模型、既异常值修正又去噪的单一模型的预测结果,各模型的MA E、MA P E汇总结果如表2所示,可以得出以下结论.表2 L O F-S S A-P S O-L S S VM与其他预测模型的预测结果对比测试集模型2 0 0 42 0 1 5年2 0 0 42 0 1 6年2 0 0 42 0 1 7年MA EMA P E/%MA EMA P E/%MA EMA P E/%A R I MA7.5 9 81 0.9 7 88.9 8 61 4.0 0 68.8 1 01 2.9 2 7L S S VM7.5 8 01 0.9 7 68.7 9 71 1.1 7 48.1 8 71 2.0 0 3P S O-L S S VM7.3 1 41 0.7 2 22.8 1 24.4 5 27.4 3 01 1.4 8 3S S A-A R I MA7.8 0 51 0.9 1 04.2 4 85.5 0 64.5 1 37.1 7 3S S A-L S S VM7.1 5 71 0.7 7 43.3 4 04.9 1 34.4 3 56.8 2 9S S A-P S O-L S S VM6.9 4 31 0.3 7 23.1 2 64.0 3 93.9 1 46.2 2 1L O F-A R I MA7.6 0 31 0.5 2 27.9 6 51 1.1 8 37.9 2 91 0.9 2 0L O F-L S S VM6.9 2 49.8 4 95.4 0 27.6 9 34.5 5 05.9 9 6L O F-P S O-L S S VM6.6 0 79.3 7 64.6 6 36.4 8 24.3 3 65.6 7 2L O F-S S A-A R I MA6.1 8 68.6 9 86.0 1 28.3 6 74.4 8 45.9 6 4L O F-S S A-L S S VM3.8 2 95.3 9 92.7 9 04.4 3 13.7 8 75.0 9 9L O F-S S A-P S O-L S S VM3.1 2 04.5 4 92.6 5 83.6 1 03.1 9 54.2 0 1 注:L S S VM和P S O-L S S VM的嵌入维数均为1 2,即用前1 2期的数据来预测下一期的数据值,S V R核函数为R B F函数.A R I MA模型阶数以A I C最小化原则来选取,最终使用的A R I MA模型为A R I MA(3,1,3).P S O算法中粒子群个数为2 0,粒子速度范围为-5,5,学习因子c1和c2均为2.0,惯性权重=1,最大迭代次数为1 0 0.用P S O算法来优化L S S VM模型中的参数和,参数的取值范围为0.0 1,1 0 0.(1)3个单一模 型(A R I MA、L S S VM、P S O-L S S VM)对比可以发现,传统的时间序列模型A R I-MA相比较人工智能模型L S S VM和P S O-L S S-VM,预测误差更大精度更低,这是因为A R I MA模型主要用来刻画序列的线性特征,对于非线性的上海旅游客流量序列的预测精度是不足的.(2)比较L S S VM和P S O-L S S VM模型,使用了P S O参数优化方法的P S O-L S S VM模型预测结果更好,说明P S O优化算法有利于寻找全局最优,能够为L S S VM找出更合适的参数并达到优化参数的目的.(3)3个基于S S A分解的集成模型(S S A-A R I-MA、S S A-L S S VM、S S A-P S O-L S S VM)与对应的单一模型A R I MA、L S S VM、P S O-L S S VM对比发现,MA E和MA P E评价指标均有了一定程度的降低,说明序列分解去 噪处理能有 效 提 升 预 测 模 型 的性能.43 曲阜师范大学学报(自然科学版)2 0 2 3年(4)3个 基 于L O F异 常 值 检 测 的 集 成 模 型(L O F-A R I MA、L O F-L S S VM、L O F-P S O-L S S VM)与对应的单一模型A R I MA、L S S VM、P S O-L S S VM对比可以发现,除了2 0 0 42 0 1 6年P S O-L S S VM和L O F-P S O-L S S VM对比组外,其他组的预测精度均得到了不同程度的提升,说明对原始序列的异常值修正能有效提高预测精度.(5)对比考虑异常值修正和S S A分解去噪的3个 模 型L O F-S S A-A R I MA、L O F-S S A-L S S VM、L O F-S S A-P S O-L S S VM与相应的单一模型、仅去噪处理的单一模型、仅异常值修正处理的单一模型,发现同时进行异常值修正和去噪处理后的模型,预测结果MA E、MA P E均是最小的,说明构建预测模型之前对时间序列进行异常值修正和去噪处理对提高模型预测精度是有效的.(6)对比L O F-S S A-A R I MA、L O F-S S A-L S S-VM、L O F-S S A-P S O-L S S VM这3个模型发现,在3组训练集上,均满足L O F-S S A-P S O-L S S VM模型预测精度最高.说明对于上海月度旅游客流量数据,L O F-S S A-P S O-L S S VM作为预测模型是最优的.3 结 论旅游客流量序列往往受到经济发展等多种因素的影响而产生异常值.本文以我国重点旅游城市上海市2 0 0 4年1月至2 0 1 8年1 2月旅游客流量月度数据为实例,采用X 1 2季节分解方法消除原始序列的季节性影响,借助基于密度的异常值检测算法L O F对序列中可能存在的异常值进行检测,同时运用插值及L S S VM进行异常值修正,此外,引入S S A分解对修正后的时间序列进行降噪处理.最后使用L O F-S S A-P S O-L S S VM模型对旅游客流量进行了滚动预测,结论如下.(1)由于受到特殊事件的影响,时间序列数据中可能包含异常值和噪声项,会影响到预测模型的参数估计,导致预测模型的预测精度下降.实证研究的预测误差表明,时间序列中异常值的检测和修正以及使用S S A方法对时间序列进行降噪处理能有效降低预测误差,改进模型的预测性能.(2)L O F-S S A-P S O-L S S VM模型综合了异常值修正、噪声提取、P S O参数寻优和L S S VM非线性处理能力的优势,相比线性模型A R I MA和未参数优化的模型L S S VM,能有效提升模型预测精度,更适合预测短期的旅游客流量.本文综合考虑了旅游客流量序列中季节性、异常值、噪声项对预测精度的干扰,并建立P S O参数寻优条件下的L S S VM人工智能模型,这为复杂旅游客流量序列的预测提供了一个新的思路和方法.根据上海旅游客流量的季节性特点,可针对不同季节制定相应的资源管理和调度措施,为游客提供优质的服务保障.模型后续的改进方向,可从互联网数据的角度拓展旅游客流量的数据来源,结合已有的旅游客流量预测模型,进一步提高模型预测精度,从而为旅游市场提供更好的指导意义.参考文献:1T S U IW,B A L L IH O,G I L B E YA,e t a l.F o r e c a s t i n go fH o n gK o n ga i r p o r t sp a s s e n g e r t h r o u g h p u tJ.T o u r i s mM a n a g e m e n t,2 0 1 4,4 2:6 2-7 6.2F I L D E SR,WE IY,I S MA I LS.E v a l u a t i n gt h ef o r e c a s-t i n gp e r f o r m a n c eo f e c o n o m e t r i cm o d e l so f a i rp a s s e n g e rt r a f f i c f l o w su s i n gm u l t i p l ee r r o rm e a s u r e sJ.I n t e r n a-t i o n a l J o u r n a l o fF o r e c a s t i n g,2 0 1 1,2 7(3):9 0 2-9 2 2.3周力.构建旅游消费价格指数对城乡居民旅游消费预测的实证分析 基于VA R模型的检验J.价格月刊,2 0 1 9(6):1-5.4梁昌勇,马银超,陈荣,等.基于S V R-A RMA组合模型的日旅游需求预测J.管理工程学报,2 0 1 5,2 9(1):1 2 2-1 2 7.5YAN G H H,CHE NSH,HUN GJY,e ta l.U t i l i z a t i o no fg e n e t i cp r o g r a mm i n gt oe s t a b l i s hd e m a n df o r e c a s t i nT a i w a n i n t e r n a t i o n a l f l i g h t sC/2 n d I n t e r n a t i o n a l C o n-f e r e n c eo nI n f o r m a t i o nE n g i n e e r i n ga n dC o m p u t e rS c i-e n c e.I E E E,2 0 1 0:1-4.6陈荣,梁昌勇.基于季节S V R的节假日旅游客流量预测J.统计与决策,2 0 1 3(9):8 2-8 4.7王晓光,张弢,周慧.基于L S-S VM的管道腐蚀速率灰色组合预测模型J.数学的实践与认识,2 0 1 4,4 4(7):8 2-8 7.8刘文颖,门德月,梁纪峰,等.基于灰色关联度与L S S VM组合的月度负荷预测J.电网技术,2 0 1 2,3 6(8):2 2 8-2 3 2.9S UN W,L I U X Y.A p p l i c a t i o nr e s e a r c ho fP S O-L S S-VMi nc a r b o ne m i s s i o np r e d i c t i o ni n H e b e iP r o v i n c eC/P r o c e e d i n g so f2 0 1 94 t hI n t e r n a t i o n a lC o n f e r e n c eo n A d v a n c e si n E n e r g y a n d E n v i r o n m e n t R e s e a r c h(I C A E E R2 0 1 9).香港环球科研协会,2 0 1 9:1 1 3 6-1 1 3 8.1 0赵坤,覃锡忠,贾振红,等.基于C S算法改进E LM的时间序列预测J.计算机工程与设计,2 0 1 8,3 9(8):53第3期 张娟:基于L O F-S S A-P S O-L S S VM模型的旅游客流量预测 2 6 4 9-2 6 5 3.1 1L I ANY,WE IH,WANGJS,e t a l.C h r o m o s o m e-l e v e lr e f e r e n c eg e n o m eo fX 1 2,ah i g h l yv i r u l e n tr a c eo ft h es o y b e a nc y s tn e m a t o d eh e t e r o d e r ag l y c i n e sJ.M o l e c-u l a rE c o l o g yR e s o u r c e s,2 0 1 9,1 9(6):1 6 3 7-1 6 4 6.1 2刘德彪,李夕兵,李响,等.基于L O F的K-m e a n s聚类方法及其在微震监测中的应用J.中国安全生产科学技术,2 0 1 9,1 5(6):8 1-8 7.1 3S UNYX,L E NGB,GUAN W.An o v e lw a v e l e t-S VMs h o r t-t i m ep a s s e n g e r f l o wp r e d i c t i o n i nB e i j i n gs u b w a ys y s t e mJ.N e u r o c o m p u t i n g,2 0 1 5,1 6 6:1 0 9-1 2 1.1 4I GHO RSDS,C A R L O SR,HA L LB,e t a l.M u l t i-r e s-o l u t i o nw a v e l e t a n a l y s i sf o rn o i s er e d u c t i o ni ni m p u l s ep u n c t u r e v o l t a g e m e a s u r e m e n t sJ.M e a s u r e m e n t,2 0 2 0,1 5 3:1 0 7 4 1 6.1 5Z HAN G X,YU L,WAN GSY,e ta l.E s t i m a t i n gt h ei m p a c to f e x t r e m ee v e n t so nc r u d eo i lp r i c e:A nEMD-b a s e de v e n ta n a l y s i sm e t h o dJ.E n e r g yE c o n o m i c s,2 0 0 9,3 1(5):7 6 8-7 7 8.1 6郭战坤,梁小珍,杨明歌,等.基于异常值检测的港口集装箱吞吐量预测模型J.数学的实践与认识,2 0 1 9,4 9(1 7):2 6-3 4.1 7梁小珍,郭战坤,张倩文,等.基于奇异谱分析的航空客运需求分析与分解集成预测模型J.系统工程理论与实践,2 0 2 0,4 0(7):1 8 4 4-1 8 5 5.T o u r i s t f l o wp r e d i c t i o nb a s e do nL O F-S S A-P S O-L S S VM m o d e lZHANGJ u a n(C o l l e g eo fM a t h e m a t i c sa n dS t a t i s t i c s,H u l u n b u i rC o l l e g e,0 2 1 0 0 8,H u l u n b u i r,I n n e rM o n g o l i a,P R C)A b s t r a c t:I nv i e wo f t h e c o m p l e x i t yo f t h e i n f l u e n c e f a c t o r s o fm o n t h l y t o u r i s t t r a f f i c a n d t h e c h a r a c t e r-i s t i c so fp o s s i b l eo u t l i e r sa n dn o i s ei nt h et i m es e r i e sd a t a,aL O F-S S A-P S O-L S S VM p r e d i c t i o n m o d e lb a s e do n l o c a l a n o m a l y f a c t o r(L O F)a n ds i n g u l a r s p e c t r a l a n a l y s i s(S S A)i sp r o p o s e d.I nt h ed a t ap r e p r o-c e s s i n gs t a g e,t h eo r i g i n a ls e q u e n c ei sd e c o m p o s e db yX 1 2s e a s o na d d i t i o n m o d e lt oe l i m i n a t es e a s o n a le f f e c t s,t h eo u t l i e r o f t h e s e a s o n a l l ya d j u s t e ds e q u e n c e i sd e t e c t e db yL O Fm e t h o d,a n d t h eo u t l i e r i s c o r r e c-t e db y i n t e r p o l a t i o na n dl e a s ts q u a r es u p p o r tv e c t o rm a c h i n e(L S S VM)m e t h o d.T h en e wp r e d i c t a b l es e-q u e n c e i sa d d e db yt h er e v i s e ds e a s o n a l a d j u s t m e n t s e q u e n c ew i t ht h es e a s o n a l f a c t o r s e q u e n c e.I nt h ep r e-d i c t i o np h a s e,t h es i n g u l a rs p e c t r a l a n a l y s i s(S S A)m e t h o d i su s e d t oe l i m i n a t e t h en o i s e i n t h en e wp r e d i c-t i o ns e q u e n c e,a n dt h e nt h eP S O-L S S VM m e t h o di su s e dt op r e d i c t i t.T a k i n gt h ef o r e c a s to ft h et o u r i s tt r a f f i c s e q u e n c e f r o mJ a n u a r y2 0 0 4t oD e c e m b e r 2 0 1 8 i nS h a n g h a i a s a ne x a m p l e,i t i s s h o w n t h a t t h eL O F-S S A-P S O-L S S VM m o d e l c a ne f f e c t i v e l y i m p r o v e t h ep r e d i c t i o na c c u r a c ya n ds t a b i l i t yo f t o u r i s t t r a f f i cb yc o m p a r i s o nw i t hA R I MAa n dL S S VM m o d e l s.K e yw o r d s:t o u r i s tt r a f f i cf o r e c a s t;l o c a la n o m a l yf a c t o r;m i n i m u mt w o-w a ys u p p o r tv e c t o rm a c h i n e;p a r t i c l eg r o u pe x c e l l e n c e;s i n g u l a rs p e c t r a l a n a l y s i s63 曲阜师范大学学报(自然科学版)2 0 2 3年