loop in a column to calculate max and min in a certain interval

Question

sud 2018-5-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/402863-loop-in-a-column-to-calculate-max-and-min-in-a-certain-interval

评论： jonas 2018-5-28

i have a column that has 200 values of length for each volume. and there can be n no. of volumes. so i need to calculate the maximum and minimum of lengths for each volume and plot volume(vs)minlength and volume(vs)maxlength in the same plot. please suggest a loop and how to apply it. help is much appreciated.

2 个评论
显示无隐藏无

jonas 2018-5-28

What do you mean by volume? Can you provide an example? The problem is not clear to me.

Stephen23 2018-5-28

编辑：Stephen23 2018-5-28

ROSETTE_S2_R5_GeometricData.csv

sud's "Answer" moved here:

hi! i will attach a file. so its 1st column is volume and 16th column is length.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

KSSV 2018-5-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/402863-loop-in-a-column-to-calculate-max-and-min-in-a-certain-interval#answer_322224

在 MATLAB Online 中打开

Let V be your volume and A be other column.

[c,ia,ib] = unique(V) ;
N = length(c) ;
iwant = zeros(N,1) ;
for i = 1:length(c)
    iwant(i) = mean(A(ib==i)) ;
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

sud 2018-5-28

在 MATLAB Online 中打开

 Hi. It is working. thanks but maybe i didnt put my problem correctly. I want to apply a loop for thia program written below! (after every 200 lengths):
 IdlCatalogue_Dir='C:\Users\yadav\Desktop\IDL_routine\OUTPUT_FILE\ROSETTE_S2_r5\';
  File = 'ROSETTE_S2_r5_GeometricData.csv';  
  FF = fullfile(IdlCatalogue_Dir,File); 
  Type = strrep(File, '_GeometricData_TOT.csv', ''); 
  Type = strrep(Type,'_','\_'); 
  data = load(FF);
  deq = data(end-199:end,16);
  for i=1:length(deq) 
  minimum(i)=min(deq(1:i)); 
  % i'th of maximum would be the max upto ith interval of deq . 
  end
x=1:1:200; 
figure() 
plot(x,maximum); 
hold on; 
plot(x,minimum); set(gca,'xscale','log'); 
xlabel('number of projections') 
ylabel('Deq(µm)') 
title('Variation of Max Deq and Min Deq with number of projections')

请先登录，再进行评论。

Answer 2

jonas 2018-5-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/402863-loop-in-a-column-to-calculate-max-and-min-in-a-certain-interval#answer_322226

在 MATLAB Online 中打开

Perhaps this will help...

data=readtable('ROSETTE_S2_R5_GeometricData.csv');
[volumes]=unique(data{:,1});
out=nan(2,length(volumes));
for i=1:length(volumes);
    out(1,i)=min(data{data{:,1}==volumes(i),16});
    out(2,i)=max(data{data{:,1}==volumes(i),16});
end
figure;
plot(volumes,out(1,:),...
    volumes,out(2,:));

9 个评论
显示 7更早的评论隐藏 7更早的评论

sud 2018-5-28

在 MATLAB Online 中打开

Hi. It is working. thanks but maybe i didnt put my problem correctly. I want to apply a loop for thia program written below! (after every 200 lengths):

IdlCatalogue_Dir = 'C:\Users\yadav\Desktop\IDL_routine\OUTPUT_FILE\ROSETTE_S2_r5\'; %path of folder File = 'ROSETTE_S2_r5_GeometricData.csv'; %name of file FF = fullfile(IdlCatalogue_Dir,File); Type = strrep(File, '_GeometricData_TOT.csv', ''); Type = strrep(Type,'_','\_'); data = load(FF);

deq = data(end-199:end,16);

% n=1:1:length(deq(1:100)); % n becomes [1,10,20,...400] for i=1:length(deq) maximum(i)=max(deq(1:i)); % i'th of maximum would be the max upto ith interval of deq . end

for i=1:length(deq) minimum(i)=min(deq(1:i)); % i'th of maximum would be the max upto ith interval of deq . end

x=1:1:200; figure() plot(x,maximum); hold on; plot(x,minimum); set(gca,'xscale','log'); xlabel('number of projections') ylabel('Deq(µm)') title('Variation of Max Deq and Min Deq with number of projections')

if true
  % code
end

sud 2018-5-28

hi. so sorry for the bad format but i did try to select it and code it. i dont know why is it not taking the form of code!

jonas 2018-5-28

matlabhelp.m

I did some formatting and wrote this for you (see attached). All values are stored in maximum and minimum. You can choose which curve to plot by changing the variable id.

If you need anything else, please format the post correctly.

请先登录，再进行评论。