using for loop to read multiple column vectors

Question

shah nawaz 2020-6-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/546894-using-for-loop-to-read-multiple-column-vectors

编辑： Matt J 2020-6-16

Rat.m

I want to use a for loop to read different column vectors which are mentioned with different fn data in attached file and then I want to plot several figures from the output of for loop. Main problem is I don't know how to import fn vector data multiple times in loop and replace it with next fn vector when code complete for previous loop.

I have already imported data from txt file by using import tool and named the columns with variable names. fn is a varible which has name fn 10, fn22, fn32.... these are 5 column vectors.

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

shah nawaz 2020-6-12

importfile7.m

and this is the function file to import data from txt

请先登录，再进行评论。

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

Answer 1

Matt J 2020-6-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/546894-using-for-loop-to-read-multiple-column-vectors#answer_449937

This should help,

https://www.mathworks.com/matlabcentral/answers/304528-tutorial-why-variables-should-not-be-named-dynamically-eval

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

shah nawaz 2020-6-14

在 MATLAB Online 中打开

Is there anyway to use cellfun(@mean for multiple vectors of array. code above gives me mean of 1 vector in my array only. Now I have modified code to :

%using importfile1.m function file to import data by cell arrays
file199 = importfile1('DMconfaf199',9805,29531); 
file171 = importfile1('DMconfaf171',5045,14081);
file132 = importfile1('DMconfaf132',1028,3097);
file122 = importfile1('DMconfaf122',488,1443);
file110 = importfile1('DMconfaf110',104,307);
%%combinig arrays into 1 array 
fn =  {file199(:,3) file171(:,3) file132(:,3) file122(:,3) file110(:,3)}; 
Particle1 = {file199(:,1) file171(:,1) file132(:,1) file122(:,1) file110(:,1)};
Particle2 = {file199(:,2) file171(:,2) file132(:,2) file122(:,2) file110(:,2)};
mfn = cellfun(@mean,fn{1,1});
% have to compute mean and standard deviation and plot figures by using hold on as well as given below
%     sfn = std2(fn); % standard deviation 
% figure(3)
% hold on
% y2 = normpdf(mfn,mmfn,smfn);
% logy2 = log10(y2);
% plot(lgmfn,logy2, '.r', 'LineWidth' ,4);
% title('Log-log plots of the probability distributions of\newline normalized normal forces for nq=71, μ=0')
% xlabel('Log_1_0(fn/mfn)')
% ylabel('Log_1_0(P_n fn/mfn)')
% grid on
% hold on
% y2 = normpdf(mfn,mmfn,smfn);
% logy2 = log10(y2);
% plot(lgmfn,logy2, '.g', 'LineWidth' ,2);
%  hold on

shah nawaz 2020-6-15

moving between matrices and cell arrays i am stuck at how to use fn values again and again in loop to find fn/it's mean value each time and plot. fn column divided by it's mean is shown by mfn here. The code you told me above worked for me now for calculating mean of each column of fn arry. Can you please have a look. Thanks..

[a,b,fn99,~,~]= importfile('DMconfaf199', 9805, 29531); % a = particle 1 for nq = 99 & b = particle 2 for nq = 99 sample

[g,h,fn22,~,~]= importfile('DMconfaf122', 488, 1443); % ~ operator ignores the function outputs.

[e,f,fn32,~,~]= importfile('DMconfaf132', 1028, 3097);

[c,d,fn71,~,~]= importfile('DMconfaf171', 5045, 14081);

[k,l,fn10,~,~]= importfile('DMconfaf110', 104, 307);

Particles = {a b c d e f g h k l};

fn = {fn10 fn22 fn32 fn71 fn99};

hold on

for i = {fn10 fn22 fn32 fn71 fn99} %for loop for different fn values in the column vectors above

sfn = cellfun(@(c)std2(c), fn); % standard deviation

mean_fn = cellfun(@(c)mean2(c), fn); % mean of whole column matrix

mean_fn = num2cell(mean_fn, 1);

sfn = num2cell(sfn, 1);

for i = 1:length(mean_fn)

mfn = cellfun(@divide, fn{1,1},mean_fn{1,1}, 'uni', 0); % fn/its mean

smfn = std2(mfn);

lgmfn = log10(mfn);

mmfn = mean2(mfn); %in order to fill normpdf function parametres we take mean again smfn = std2(mfn);

log_fn = log10(fn);

lmfn = log10(mfn);

hold on

figure(3)

box on

y2 = normpdf(mfn,mmfn,smfn);

logy2 = log10(y2);

plot(lgmfn,logy2, '.r', 'LineWidth' ,4);

title('Log-log plots of the probability distributions of\newline normalized normal forces for nq=71, μ=0')

xlabel('Log_1_0(fn/mfn)')

ylabel('Log_1_0(P_n fn/mfn)')

grid on

hold on

end

shah nawaz 2020-6-15

Rat.m

i have attached the code file as well.

Matt J 2020-6-16

编辑：Matt J 2020-6-16

You are not using the same variable 'i' for both the inner and outer loops loop counters. Additionally, you don't seem to be using the loop counters anywhere inside the loops.

请先登录，再进行评论。