风电场一年内的风速实测数据

包括风电场一年内的风速实测数据，并根据实测数据对风电场风速和输出功率进行预测和建模仿真，希望对大家有帮助

%% 2006-07-19
%% 此程序主要功能是在风向固定，风速变化的情况下，求总风场输出的功率

clc
clear
%format long
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 风机参数拟合时的原始数据
x1=[4:25]‘;
%%风力机的功率系数曲线
yy1=[0 0.230 0.352 0.406 0.431 0.440 0.438 0.419 0.384 0.335 0.281 0.232 0.192 0.160 0.135 0.115...
        0.098 0.085 0.074 0.065 0.057 0.05]‘;
yy4=[0.226 0.32 0.4	0.42 0.423 0.425 0.408 0.375 0.325 0.275 0.236 0.185 0.16 0.13 0.1 0.085 0.07 0.065	0.058 0.048	0.038 0.038]‘;


x2=[4:16]‘;
%%风电机组的功率曲线
yy2=[0 55 145 266 421 613 836 1065 1267 1407 1474 1495 1500 1500 1500 1500 1500 1500 1500 1500 1500]‘;

yy3=[0 55 145 266 421 613 836 1065 1267 1407 1474 1495 1500 ]‘;

%%%进行多项式拟合
%功率系数曲线多项式系数
p_cp=polyfit（x1yy45）;

%%以下两语句可显示拟合的效果
% pop_cp=polyval（p_cpx1）;
% plot（x1yy1‘*‘x1pop_cp‘-‘）

%功率曲线系数
p_pv=polyfit（x2yy36）;
%%以下两语句可显示拟合的效果
% pop_pv=polyval（p_pvx2）
% plot（x1yy2‘*‘x2pop_pv‘-‘）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%风机总台数
n=16;

%风机坐标 

wg=[100110; 420110;740110; 1060110;
    100430;420430;740430;1060430;
    100750;420750;740750;1060750;
    1001100;4201100;7401100;10601100;];
wg（:3）=0;

wg（:4）=0;

%计算出风机与原点的正切值
for ii=1:n
   wg（ii4）=（wg（ii2）/wg（ii1））;
end
 
%假定风机的半径
wg（:5）=31.5;

% plot（wg（11）wg（12）‘*‘）
% hold on
% plot（wg（21）wg（22）‘*‘）
% plot（wg（31）wg（32）‘*‘）
% plot（wg（41）wg（42）‘*‘）
% plot（wg（51）wg（52）‘*‘）
% plot（wg（61）wg（62）‘*‘）
% plot（wg（71）wg（72）‘*‘）
% plot（wg（81）wg（82）‘*‘）
% plot（wg（91）wg（92）‘*‘）
% plot（wg（101）wg（102）‘*‘）
% plot（wg（111）wg（112）‘*‘）
% plot（wg（121）wg（122）‘*‘）
% plot（wg（131）wg（132）‘*‘）
% plot（wg（141）wg（142）‘*‘）
% plot（wg（151）wg（152）‘*‘）
% plot（wg（161）wg（162）‘*‘）

%自然风的风向  与轴的夹角单位  度
sita=78;

   
%自然风的风速

wind_time=1:100;
wind_s_low=-0.1;
wind_s_high=0.1;
wind_speed_n=8;
wind_speed=wind_speed_n+wind_s_low+（wind_s_high-wind_s_low）*rand（1100）;
wind_number=size（wind_speed2）;
wind_time=1:wind_number;
%判断出迎着风向的第一台风机
sita1=sita*pi/180;
si_f=tan（sita1）;
%plot（wg（161）wg（161）*tan（sita1）‘+‘）
for ii=1:n
    if wg（ii4）>=si_f
        wg（ii3）=cos（pi/2-sita1-atan（wg（ii1）/wg（ii2）））*（sqrt（wg（ii1）*wg（ii1）+wg（ii2）*wg（ii2）））;
   
    else
        wg（ii3）=cos（sita1-atan（wg（ii2）/wg（ii1）））*（sqrt（wg（ii1）*wg（ii1）+wg（ii2）*wg（ii2）））;
    end
end
%%对的列进行调整
wg（:6）=wg（:1）;
wg（:7）=wg（:2）;
wg（:1）=wg（:3）;
wg（:2）=0;
wg（:3）=0;
wg（:8）=0;
wg（:9）=[1:n]‘;
wg;
wg_1=sortrows（wg）;

    %第一台风机的原始坐标
    wt_first_x=wg_1（16）;
    wt_first_y=wg_1（17）;
    %计算出以第一台风机为坐标原点顺风向为x轴后各风机的新坐标
    wg_new=wg_1;
    % 先坐标平移
    for ii=1:n
        wg_new（ii6）=wg_1（ii6）-wt_first_x;
        wg_new（ii7）=wg_1（ii7）-wt_first_y;
    end
    wg_new1=wg_new;
    % 坐标旋转
    for ii=1:n
        wg_new1（ii6）=wg_new（ii6）*cos（sita1）+wg_new（ii7）*sin（sita1）;
        wg_new1（ii7）=-wg_new（ii6）*sin（sit

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-11-07 11:45  windfarm\
     文件       54784  2013-11-07 11:02  windfarm\Book1.xls
     文件       55296  2013-11-07 11:07  windfarm\Book2.xls
     文件       39936  2013-11-07 11:08  windfarm\Book3.xls
     文件       13800  2013-11-07 12:02  windfarm\fengsu.txt
     文件       10033  2013-11-07 12:02  windfarm\windfarm.m
     文件       10859  2013-11-07 11:23  windfarm\windspeed.txt
     文件      160768  2013-11-07 10:53  windfarm\windspeed.xls
     目录           0  2013-11-07 10:53  windfarm\实测风速数据\
     目录           0  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\
     文件       58888  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\11.xls
     文件       58617  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\12.xls
     文件       57934  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\4.xls
     文件       56874  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\5.xls
     文件       58021  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\6.xls
     文件       27435  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\QQ截图未命名.jpg
     文件      160768  2013-11-07 10:53  windfarm\实测风速数据\10分钟间隔\合并.xls
     目录           0  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\
     文件       27435  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\QQ截图未命名.jpg
     文件       51712  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\TB004-08：40.xls
     文件       61440  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\TB005-08：40.xls
     文件       61440  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\TB006-08：40.xls
     文件       61440  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\TB011-08：40.xls
     文件       61440  2013-11-07 10:53  windfarm\实测风速数据\1s间隔\TB012-08：40.xls
     目录           0  2013-11-07 10:53  windfarm\实测风速数据\ws\
     文件      114688  2013-11-07 10:53  windfarm\实测风速数据\ws\TB004-07：00.xls
     文件      114688  2013-11-07 10:53  windfarm\实测风速数据\ws\TB005-07：00.xls
     文件      114176  2013-11-07 10:53  windfarm\实测风速数据\ws\TB006-07：00.xls
     文件      114176  2013-11-07 10:53  windfarm\实测风速数据\ws\TB011-07：00.xls
     文件      114176  2013-11-07 10:53  windfarm\实测风速数据\ws\TB012-07：00.xls
     文件       53857  2013-11-07 10:53  windfarm\实测风速数据\ws.rar
............此处省略0个文件信息