%%?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個文件信息