ARIMA股价预测.docx
这个人逃单，查了地址，应该是上海大学宝山区的，结果最后白嫖，程序就给大家免费分享T<长君 amy(4)NHASMJ400元，已付款，可以做了<长君a my (4)好的前后你把数据发出来吧长名老师你是用matlab还是eviews做啊6月23日晚上20：23matlab6月21日上午11 44长昌智飞生物(1).xls528 0 KBy城信电髯版长君D好的哈好喻长追老师这些数据可以吗，还需要什么跟我说我再去找ARIMA股价预测20220623.rar159.0 KBq微信电脑版数据分析代做ttMfita 同分析««s囹0分析 用关分析生存分析“同序列«*«及度学习 交间计时间分析经网第机届学习i俱仿Bttft9香发洛舞色优化电子电气SPSS数据分析python服务mat06月23日11:35前发货*100.oo×4wnx,卖了换钱申请售后实付款v4OO.oo收货信息：刘宸君，86-15221033615,上海上海市宝山区杨行镇江杨北路1568弄54号901室订单编号:1596151202394461699 ：复制这是用ns2去对AIC和BIC进行排序，选择最优p, q参数，差分用的是一阶完整程序如下，其中需要用到老版本的econ函数包，可以去下载个matlab2011版，在toolbox文件夹内clearclc%读取数据，选取2021年1月4 -2022年6月2 口数据放在X.txt中，需要检验的2022年6月6 0-2022年6月10日数据放在Y.txt中load X.txt %历史收盘价数据，一维load Y.txt %用于预测检测的数据，一维figureplot(X)xlim(l,length(X)xticks(lJength(X)xticklabels(,2022-l-42022-6-2,)tit®收盘价-时间序列%ARIMA预测nn=5;%预测步数rll=autocorr(X); %计算自相关函数rl2=parcorr(X); %计算偏相关函数da=diff(X); %计算1阶差分r21=autocorr(da); %计算自相关函数r22=parcorr(da); %计算偏相关函数figuresubplot(211)zautocorr(X);subplot(212),parcorr(X);figuresubplot(211),autocorr(da);subplot(212),parcorr(da);n=length(da); %计算差分后的数据个数AICBIC=;for i=0:8for j=0:8spec= garchsetCR'J/M'J/DisplayVoff'); % 指定模型的结构coeffX,errorsX,LLFX = garchfit(spec,da); %拟合参数num=garchcount(coeffX); %计算拟合参数的个数%compute Akaike and Bayesian Information Criteriaaic,bic=aicbic(LLFX/num,n);fprintf('R=%d,M=%d,AIC=%f/BIC=%fn，,i,j,aic/bic); %显示计算结果AICBIC=AICBIC;ij/aic,bic;endend%非支配排序选出最佳R、M系数AICBIC(:/3：4)/AICBIC(:/l:2)=ns2(AICBIC(:,l:2)/AICBIC(:/3),AICBIC(:,4)；r=AICBIC(l,l)jm=AICBIC(lz2spec2= garchset('R'zr,'MmDisplayoff)j % 指定模型的结构coeffX,errorsX,LLFX = garchfit(spec2,da); %拟合参数sigmaForecast,w_Forecast = garchpred(coeffX,da,nn); %计算 5 步预报值x_pred=X(end-nn+l:end)+cumsum(w_Forecast)% 预测值%平均百分比误差err=round(mean(abs(x_pred-Y)./Y)*100,2);figurehold onplot(x-pred,'r-o,)plot(X'b-*')xlabel('时间')ylabel(收盘价)title('ARIMA时间序列预测，平均误差百分比为：num2str(ej,%)legend。预测值？检验值,'Location','SouthEast)ylim(0,100)xticks(l:5)xticklabels(,2022-6-6'z,2022-6-7,z'2022-6-872022-6-92022-6-10')function TT,chrom=ns2(NN,FlzF2)%快速非支配排序a = 0;T1 =；T2 = ;chrom=NN;chroml =;chrom2 =;while a = 0 %根据被支配数进行分级和排序M = ;for i = llength(Fl)M(izl) = length(find(Fl<Fl(i,l)+length(find(F2<F2(i,l)依 目标函数最小化这里为 <,最大化改成>endbl =；b2 =;bl,b2 = sort(M); %bl返回从小到大排序，b2返回原始序号if length(chrom)>0 && bl(l) = 0%无被支配数进入一级用T1矩阵保存Tl = Tl;Fl(b2(l)zF2(b2(l);chroml = chroml;chrom(b2(l),:);Fl(b2(l)=;F2(b2(l)=;chrom(b2(l) =;else %有被支配数进入二级用T2矩阵保存a = 1;T2 = F1,F2;chrom2 = chrom;endendT2=T2(b2/：);chrom2 = chrom2(b2,:);ifsize(Tl,l)>2 %T1矩阵不用进行拥挤度调整排序，直接对T2进行排序调整即可y = yongji(Tl)% 拥挤度for i = 2:size(Tl,l)ify(i-l) >y(i)Tl(i-l:l:i,:) = Tl(i:-l:i-l,:); %根据拥挤度调整排序,如果后者优于前者则反转顺序chroml(i-l:l:i,:) = chroml(i:-l:i-l,:);endendendiflength(T2)>0 %T1矩阵不用进行拥挤度调整排序，直接对T2进行排序调整即可y = yongji(T2)% 拥挤度for i = 2:size(T2,l)if bl(i)=bl(i-l)ify(i-l) >y(i)T2(i-l:l:i,:) = T2(i:-l:i-lz:); %根据拥挤度调整排序，如果后者优于前者则反转顺序chrom2(i-l:l:i,:) = chrom2(i:-l:i-l,:);endendendend%排序重组 = T1;T2;chrom = chroml;chrom2;function y=yongji(H)%计算拥挤度yi=H(：/i);y2=H(:,2);yyl,al=sort(yl);yy2,a2=sort(y2);L=；puf(3z)uns=Ze)T(T'e)l=1l =puaJ(3X)u!iu-(3)xe)(z-!)3-(z+!)3X),()u!UJ-()xe)(z-!)X-(z+!)f=1-(AA)8u3N=!oj,l =
