目錄

1?主要內(nèi)容

2?部分程序

3?下載鏈接

---

1?主要內(nèi)容

這個(gè)程序《考慮時(shí)空相關(guān)性的風(fēng)電功率預(yù)測誤差建模與分析》畫的圖片非常漂亮，和原文獻(xiàn)基本一致，但是實(shí)際上內(nèi)容并未實(shí)現(xiàn)出來，主要就是利用現(xiàn)有的風(fēng)電預(yù)測的數(shù)據(jù)和結(jié)果做了相關(guān)的圖，大家不要白花錢，今天把這個(gè)程序免費(fèi)分享給大家，大家可以學(xué)習(xí)一下畫圖技巧以及數(shù)據(jù)分析方面的知識。

ps：原程序的95%置信區(qū)間是采用均值加減標(biāo)準(zhǔn)差是不對的，已經(jīng)更正過來，應(yīng)該增加1.96的系數(shù)，程序包附上了置信區(qū)間和均值標(biāo)準(zhǔn)差的關(guān)系說明。

2?部分程序

    %% 清空環(huán)境變量
clc
clear all;
%% 提取數(shù)據(jù) 
data=xlsread('實(shí)驗(yàn)數(shù)據(jù).xlsx',1);
?
%% 提取對應(yīng)各段中點(diǎn)位置處的誤差值
error_fenbu_1=[];
for i=1:size(data,1)if data(i,3)>=220 && data(i,3)<=240error_fenbu_1(i)=data(i,8);elseerror_fenbu_1(i)=0;end
end
error_1=error_fenbu_1(find(error_fenbu_1~=0));
error_fenbu_2=[];
for i=1:size(data,1)if data(i,3)>=670&&data(i,3)<=690;error_fenbu_2(i)=data(i,8);else error_fenbu_2(i)=0;end
end  
error_2=error_fenbu_2(find(error_fenbu_2~=0));
?
error_fenbu_3=[];
for i=1:size(data,1)if data(i,3)>=1128 && data(i,3)<=1148;error_fenbu_3(i)=data(i,8);else error_fenbu_3(i)=0;end
end
error_3=error_fenbu_3(find(error_fenbu_3~=0));
error_fenbu_4=[];
for i=1:size(data,1)if data(i,3)>=1585&&data(i,3)<=1605;error_fenbu_4(i)=data(i,8);else error_fenbu_4(i)=0;end
end  
error_4=error_fenbu_4(find(error_fenbu_4~=0));
error_fenbu_5=[];
for i=1:size(data,1) if data(i,3)>=2040&&data(i,3)<=2060;error_fenbu_5(i)=data(i,8);else   error_fenbu_5(i)=0;end
end
error_5=error_fenbu_5(find(error_fenbu_5~=0));
?
error_fenbu_6=[];
for i=1:size(data,1) if data(i,3)>=2495 && data(i,3)<=2515;error_fenbu_6(i)=data(i,8);else   error_fenbu_6(i)=0;end
end
error_6=error_fenbu_6(find(error_fenbu_6~=0));
error_fenbu_7=[];
for i=1:size(data,1)  if data(i,3)>=2950&&data(i,3)<=2970;error_fenbu_7(i)=data(i,8);else   error_fenbu_7(i)=0;end
end 
error_7=error_fenbu_7(find(error_fenbu_7~=0));
error_fenbu_8=[];
for i=1:size(data,1)  if data(i,3)>=3406 && data(i,3)<=3426;error_fenbu_8(i)=data(i,8);   else   error_fenbu_8(i)=0;end
end 
error_8=error_fenbu_8(find(error_fenbu_8~=0));
error_fenbu_9=[];
for i=1:size(data,1)  if data(i,3)>=3860&&data(i,3)<=3880;error_fenbu_9(i)=data(i,8); else   error_fenbu_9(i)=0;end
end 
error_9=error_fenbu_9(find(error_fenbu_9~=0));
error_fenbu_10=[];
for i=1:size(data,1)  if data(i,3)>=4317&&data(i,3)<=4337;error_fenbu_10(i)=data(i,8); else   error_fenbu_10(i)=0;end
end 
error_10=error_fenbu_10(find(error_fenbu_10~=0));
?
%% 擬合分布—求取t分布參數(shù)進(jìn)行擬合
error_values=-3000:0.5:3000;
pd_1= fitdist(error_1','tLocationScale');
desity_1= pdf(pd_1,error_values);
pd_2= fitdist(error_2','tLocationScale');
desity_2= pdf(pd_2,error_values);
pd_3= fitdist(error_3','tLocationScale');
desity_3= pdf(pd_3,error_values);
pd_4= fitdist(error_4','tLocationScale');
desity_4= pdf(pd_4,error_values);
pd_5= fitdist(error_5','tLocationScale');

3?下載鏈接