重庆分公司,新征程启航
为企业提供网站建设、域名注册、服务器等服务
本文小编为大家详细介绍“matlab怎么实现风电功率”,内容详细,步骤清晰,细节处理妥当,希望这篇“matlab怎么实现风电功率”文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。
石狮网站建设公司创新互联建站,石狮网站设计制作,有大型网站制作公司丰富经验。已为石狮数千家提供企业网站建设服务。企业网站搭建\成都外贸网站制作要多少钱,请找那个售后服务好的石狮做网站的公司定做!
%% ARMA 预测clc,clear,close allload data%导入数据T=30;buchang=size(unnamed,1)-T;%预测步长y = unnamed(1:T);[m,n]=size(y);%% %% 3.确定ARMA模型阶数% ACF和PACF法,确定阶数figuresubplot(211),autocorr( y );subplot(212),parcorr( y );figuredy = diff( y );subplot(211),autocorr( dy );subplot(212),parcorr( dy ); %% ARIMA 模型Mdl = arima(5,1,0);EstMdl = estimate(Mdl,y);res = infer(EstMdl,y); %res即残差% 模型验证figuresubplot(2,2,1)plot(res./sqrt(EstMdl.Variance))title('Standardized Residuals')subplot(2,2,2),qqplot(res)subplot(2,2,3),autocorr(res)subplot(2,2,4),parcorr(res)% 预测[yF,yMSE] = forecast(EstMdl,buchang,'Y0',y);UB = yF + 1.96*sqrt(yMSE); %95置信区间下限LB = yF - 1.96*sqrt(yMSE); %95置信区间下限yF=[unnamed(1:T);yF];figure(4)h5 = plot(unnamed,'b');hold onh6 = plot(yF,'r','LineWidth',2);h7 = plot(m+1:m+buchang,UB,'k--','LineWidth',1.5);plot(m+1:m+buchang,LB,'k--','LineWidth',1.5); legend('实际幅值','预测幅值');xlabel('时间序列')ylabel('幅值')title('arma预测图')bp_mse = mean((yF-unnamed).^2);%msedisp(['ARMA预测的mse=',num2str(bp_mse)])bp_mae = mean(abs(yF-unnamed));%maedisp(['ARMA预测的mae=',num2str(bp_mae)])bp_rmse = sqrt(mean((yF-unnamed).^2));%均方差disp(['ARMA预测的rmse=',num2str(bp_rmse)])
编辑
读到这里,这篇“matlab怎么实现风电功率”文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注创新互联行业资讯频道。