or i=1:length(flist) fname=flist(i).name; t=readtable(fname,'VariableNamingRule', 'preserve'); index1=find(t.CI_WindSpeed1>8 | t.CI_WindSpeed1<2); index2=find(t.CI_YawError1>0.17452 | t.CI_YawError1<-0.17452); index3=find(t.CI_PcsMeasuredGeneratorSpeed>183 | t.CI_PcsMeasuredGeneratorSpeed<111); index4=find(t.CI_IprRealPower>3600000 | t.CI_IprRealPower<1000000); index5=find(t.CI_PitchPositionA1>0); index=[index1;index2;index3;index4;index5]; index=unique(index); t(index,:)=[]; temp=t{:,:}; data=[data;temp]; end t(1:size(data,1),:)=num2cell(data); t(size(data,1)+1:end,:)=[]; ws=t.CI_WindSpeed1; Power=t.CI_IprRealPower; wsi=2:1:ceil(max(ws)); for i=1:length(wsi)-1 index=find(ws>=wsi(i) & ws<wsi(i+1)); str{i}=[num2str(wsi(i)) '~' num2str(wsi(i+1))]; WS(i,1)=length(index); WS(i,2)=mean(Power(index))*3600; if isnan(WS(i,2)) WS(i,2)=0; end disp(['风速区间 ' str{i} ' 小时发电量为:' num2str(WS(i,2))]) end figure bar(WS(:,2),0.5); set(gca,'xticklabel',str) xlabel('WindSpeed') ylabel('RealPower'),想把风速变成0.5一个区间统计
时间: 2024-03-08 16:48:25 浏览: 20
将wsi=0.5:0.5:ceil(max(ws)),即将风速区间的步长从1改为了0.5。同时,在统计每个风速区间的发电量时,也需要将相应的代码进行修改,例如:
```matlab
wsi=0.5:0.5:ceil(max(ws));
for i=1:length(wsi)-1
index=find(ws>=wsi(i) & ws<wsi(i+1));
str{i}=[num2str(wsi(i)) '~' num2str(wsi(i+1))];
WS(i,1)=length(index);
WS(i,2)=mean(Power(index))*3600;
if isnan(WS(i,2))
WS(i,2)=0;
end
disp(['风速区间 ' str{i} ' 小时发电量为:' num2str(WS(i,2))])
end
```
最后,将xticklabel的值也需要进行修改,例如:
```matlab
set(gca,'xticklabel',str(1:2:end))
```
这样可以将显示的风速区间从1个变为2个,从而适应步长变为0.5的情况。
相关问题
flist.append(lambda x, i=i: x * i)这句代码什么意思
这段代码实现了一个 lambda 函数,将 x 乘以一个变量 i,然后将这个 lambda 函数添加到 flist 列表中。在添加时,使用了 i=i 的 trick,将循环中的变量 i 作为默认参数传递给 lambda 函数,以避免在循环中出现变量绑定问题。
for i=1:length(flist) fname=flist(i).name; t=readtable(fname,'VariableNamingRule', 'preserve'); index1=find(t.CI_WindSpeed1>8 | t.CI_WindSpeed1<2); index2=find(t.CI_YawError1>0.17452 | t.CI_YawError1<-0.17452); index3=find(t.CI_PcsMeasuredGeneratorSpeed>183 | t.CI_PcsMeasuredGeneratorSpeed<111); index4=find(t.CI_IprRealPower>3600000 | t.CI_IprRealPower<1000000); index5=find(t.CI_PitchPositionA1>0); index=[index1;index2;index3;index4;index5]; index=unique(index); t(index,:)=[]; temp=t{:,:}; data=[data;temp]; end t(1:size(data,1),:)=num2cell(data); t(size(data,1)+1:end,:)=[]; ws=t.CI_WindSpeed1; Power=t.CI_IprRealPower; wsi=2:1:ceil(max(ws)); for i=1:length(wsi)-1 index=find(ws>=wsi(i) & ws<wsi(i+1)); str{i}=[num2str(wsi(i)) '~' num2str(wsi(i+1))]; WS(i,1)=length(index); WS(i,2)=mean(Power(index))*3600; if isnan(WS(i,2)) WS(i,2)=0; end disp(['风速区间 ' str{i} ' 小时发电量为:' num2str(WS(i,2))]) end figure bar(WS(:,2),0.5); set(gca,'xticklabel',str) xlabel('WindSpeed') ylabel('RealPower'),分析上列代码
这段代码的作用是读取多个表格文件并进行数据处理和可视化。具体来说,它通过循环读取文件列表(flist)中的每个文件,然后读取文件中的数据表格(t),并在数据表格中筛选出符合特定条件的数据行(index1,index2,index3,index4,index5),并将这些行从数据表格中删除,最终将处理后的数据存入data中。接着,代码将data转换为cell格式的数据表格(t),并根据风速对实际功率进行分组统计,得到不同风速范围内的发电量(WS(:,2)),最后以条形图的形式将发电量可视化出来。