matlab怎么實(shí)現(xiàn)風(fēng)電功率

本文小編為大家詳細(xì)介紹“matlab怎么實(shí)現(xiàn)風(fēng)電功率”，內(nèi)容詳細(xì)，步驟清晰，細(xì)節(jié)處理妥當(dāng)，希望這篇“matlab怎么實(shí)現(xiàn)風(fēng)電功率”文章能幫助大家解決疑惑，下面跟著小編的思路慢慢深入，一起來學(xué)習(xí)新知識(shí)吧。

部分代碼

%% ARMA 預(yù)測(cè)clc,clear,close allload data%導(dǎo)入數(shù)據(jù)T=30;buchang=size(unnamed,1)-T;%預(yù)測(cè)步長(zhǎng)y = unnamed(1:T);[m,n]=size(y);%% %% 3.確定ARMA模型階數(shù)% ACF和PACF法，確定階數(shù)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即殘差% 模型驗(yàn)證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)% 預(yù)測(cè)[yF,yMSE] = forecast(EstMdl,buchang,'Y0',y);UB = yF + 1.96*sqrt(yMSE); %95置信區(qū)間下限LB = yF - 1.96*sqrt(yMSE); %95置信區(qū)間下限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('實(shí)際幅值','預(yù)測(cè)幅值');xlabel('時(shí)間序列')ylabel('幅值')title('arma預(yù)測(cè)圖')bp_mse = mean((yF-unnamed).^2);%msedisp(['ARMA預(yù)測(cè)的mse=',num2str(bp_mse)])bp_mae = mean(abs(yF-unnamed));%maedisp(['ARMA預(yù)測(cè)的mae=',num2str(bp_mae)])bp_rmse = sqrt(mean((yF-unnamed).^2));%均方差disp(['ARMA預(yù)測(cè)的rmse=',num2str(bp_rmse)])

仿真結(jié)果

編輯

讀到這里，這篇“matlab怎么實(shí)現(xiàn)風(fēng)電功率”文章已經(jīng)介紹完畢，想要掌握這篇文章的知識(shí)點(diǎn)還需要大家自己動(dòng)手實(shí)踐使用過才能領(lǐng)會(huì)，如果想了解更多相關(guān)內(nèi)容的文章，歡迎關(guān)注億速云行業(yè)資訊頻道。