%% 相关性分析 figure data2_limit_steady(:,end+1) = (data2_limit_steady(:,2) - data2_limit_steady(:,3))./data2_limit_steady(:,2)*100; n_data = length(data2_limit_steady(:,1)); n_cluster = 12; Q=zeros(n_cluster,1); step = floor(n_data/n_cluster); for i = 1:n_cluster subplot(n_cluster/2,2,i); index = find((data2_limit_steady((i-1)step+1:istep,10)==1) & (data2_limit_steady((i-1)step+1:istep,1)>480)); row = (i-1)step + index; plot(data2_limit_steady(row,6),data2_limit_steady(row,end),'s'); axis([0,500,60,100]); xlabel('喷氨量'); ylabel('效率'); title(num2str(i)); Q(i)=sum(data2_limit_steady(row,1).(data2_limit_steady(row,2)-data2_limit_steady(row,3)))./sum(data2_limit_steady(row,6));修改代码,使负荷大于480小于510
时间: 2024-02-13 22:07:33 浏览: 17
以下是修改后的代码:
```
figure
data2_limit_steady(:,end+1) = (data2_limit_steady(:,2) - data2_limit_steady(:,3))./data2_limit_steady(:,2)*100;
n_data = length(data2_limit_steady(:,1));
n_cluster = 12;
Q=zeros(n_cluster,1);
step = floor(n_data/n_cluster);
for i = 1:n_cluster
subplot(n_cluster/2,2,i);
index = find((data2_limit_steady((i-1)*step+1:i*step,10)==1) & (data2_limit_steady((i-1)*step+1:i*step,1)>480) & (data2_limit_steady((i-1)*step+1:i*step,1)<510));
row = (i-1)*step + index;
plot(data2_limit_steady(row,6),data2_limit_steady(row,end),'s');
axis([0,500,60,100]);
xlabel('喷氨量');
ylabel('效率');
title(num2str(i));
Q(i)=sum(data2_limit_steady(row,1).*(data2_limit_steady(row,2)-data2_limit_steady(row,3)))./sum(data2_limit_steady(row,6));
end
```
主要的修改有:
1. 在第9行,将 `istep` 修改为 `i*step`。
2. 在第10行,增加一个条件 `& (data2_limit_steady((i-1)*step+1:i*step,1)<510)`,限制喷氨量小于510。
3. 在第11行,将 `(i-1)step` 修改为 `(i-1)*step`。
4. 在第13行,将 `data2_limit_steady(row,1).(data2_limit_steady(row,2)-data2_limit_steady(row,3))` 修改为 `data2_limit_steady(row,1).*(data2_limit_steady(row,2)-data2_limit_steady(row,3))`。