[基礎(chǔ)科學(xué)]風(fēng)功率預(yù)測論文-資料下載頁
2025-01-18 13:13本頁面
【正文】 +y(i)。endB=ones(n1,2)。for i=1:(n1)B(i,1)=(yy(i)+yy(i+1))/2。B(i,2)=1。endBT=B39。for j=1:n1YN(j)=y(j+1)。endYN=YN39。A=inv(BT*B)*BT*YN。a=A(1)。u=A(2)。t=u/a。t_test=input(39。請輸入需要預(yù)測個數(shù)：39。)。i=1:t_test+n。yys(i+1)=(y(1)t).*exp(a.*i)+t。yys(1)=y(1)。for j=n+t_test:1:2ys(j)=yys(j)yys(j1)。endx=1:n。xs=2:n+t_test。yn=ys(2:n+t_test)。plot(x+1960,y,39。*r39。,xs+1960,yn,39。.b39。)。det=0。for i=2:ndet=det+abs(yn(i)y(i))。enddet=det/(n1)。%disp([39。百分絕對誤差為：39。,num2str(det),39。%39。])。disp([39。預(yù)測值為： 39。,num2str(ys(n+1:n+t_test))])。：clcclearload data%% 數(shù)據(jù)累加作為網(wǎng)絡(luò)輸入[n,m]=size(X)。for i=1:n for j=1:m y(i,j)=sum(X(1:i,j))。 endend%% 網(wǎng)絡(luò)參數(shù)初始化a=+rand(1)/4。for k=1:(m1) b(k)=+rand(1)/4。end%% 學(xué)習(xí)速率初始化for k=1:(m1) u(k)=。end%% 權(quán)值閥值初始化t=1。w11=a。w2(1,1)=y(1,1)。for k=2:mw2(k)=2*b(k1)/a。endfor k=1:mw3(k)=1+exp(a*t)。endtheta1=0。for k=1:(m1)theta1=theta1+b(1,k)*y(1,(k+1))/a。endtheta2=(1+exp(a*t))*(theta1y(1,1))。kk=1。%% 循環(huán)迭代for j=1:10%循環(huán)迭代E(j)=0。for i=1:(n6) %% 網(wǎng)絡(luò)輸出計算 t=i。 %LB層輸出 LB_b=1/(1+exp(w11*t))。 %LC層輸出 LC_c(1,1)=LB_b*w2(1,1)。 for k=2:m LC_c(k)=y(i,k)*LB_b*w2(1,k)。 end %LD層輸出 LD_d=0。 for k=1:m LD_d=LD_d+w3(1,k)*LC_c(1,k)。 end %閥值 theta3=0。 for k=2:m theta3=theta3+w2(1,k)*y(i,k)/2。 end theta4=(1+exp(w11*t))*(theta3y(1,1))。 %網(wǎng)絡(luò)輸出值 ym=LD_dtheta4。 yc(i)=ym。 %% 權(quán)值修正 error=ymy(i,1)。 %計算誤差 E(j)=E(j)+abs(error)。 %誤差求和 %計算誤差 for k=1:m errorr(k)=error*(1+exp(w11*t))。 end thetaE0=0。 for k=1:m thetaE0=thetaE0+w2(1,k)*errorr(1,k)。 end thetaE=(1/(1+exp(w11*t)))*(11/(1+exp(w11*t)))*thetaE0。 %修改權(quán)值 for k=2:m w2(1,k)=w2(1,k)u(k1)*errorr(1,k)*LB_b。 end w11=w11+a*t*thetaE。endend %畫誤差隨進(jìn)化次數(shù)變化趨勢figure(1)plot(E)title(39。訓(xùn)練誤差39。,39。fontsize39。,12)。xlabel(39。進(jìn)化次數(shù)39。,39。fontsize39。,12)。ylabel(39。誤差39。,39。fontsize39。,12)。%print dtiff r600 283%根據(jù)訓(xùn)出的神經(jīng)網(wǎng)絡(luò)進(jìn)行預(yù)測for i=(n5):n t=i。 LB_b=1/(1+exp(w11*t))。 %LB層輸出 LC_c(1,1)=LB_b*w2(1,1)。 %LC層輸出 for k=2:m LC_c(1,k)=y(i,k)*LB_b*w2(1,k)。 %LC層輸出 end LD_d=0。 for k=1:m LD_d=LD_d+w3(1,k)*LC_c(1,k)。 %LD層輸出 end theta5=0。 for k=2:m theta5=theta5+w2(1,k)*y(i,k)/2。 end theta6=(1+exp(w11*t))*(theta5y(1,1))。 %閥值 ym=LD_dtheta6。 %網(wǎng)絡(luò)輸出值 yc(i)=ym。endyc=yc*10。y(:,1)=y(:,1)。%計算預(yù)測for j=n:1:2 ys(j)=(yc(j)yc(j1))/10。endfigure(2)plot(ys((n5):n),39。*39。)。hold onplot(X((n5):n,1),39。r:o39。)。legend(39。預(yù)測值39。,39。實際值39。)title(39?；疑窠?jīng)網(wǎng)絡(luò)預(yù)測39。,39。fontsize39。,12)xlabel(39。x39。,39。fontsize39。,12)ylabel(39。y39。,39。fontsize39。,12) 23
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1