I'm trying to find the RMSD of the numbers by taking 1-2,2-3,3-4 in a sequence ,how can i loop as such ??

Question

Ramesh Bala 2019-6-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/469319-i-m-trying-to-find-the-rmsd-of-the-numbers-by-taking-1-2-2-3-3-4-in-a-sequence-how-can-i-loop-as-su

评论： Guillaume 2019-6-30

I'm right now doing as such,how to form a loop by taking in the previous nos?

N = [ 55378 55344 55310 55276 55242 55208];
a = N(1,1);
b = N(1,2);
RMSDD1 = sqrt ( sum(( b-a).^2 /sum(a).^2) );
fprintf ('the RMSD value between two nos :    %f\n' ,RMSDD1);
RMSDDP = RMSDD1*100;
fprintf ('the RMSD value between two nos in % :    %f\n' ,RMSDDP);
a = N(1,2);
b = N(1,3);
RMSDD2 = sqrt ( sum(( b-a).^2 /sum(a).^2) );
fprintf ('the RMSD value between two nos :    %f\n' ,RMSDD2);
RMSDDP = RMSDD2*100;
fprintf ('the RMSD value between two nos in % :    %f\n' ,RMSDDP);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Guillaume 2019-6-28

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/469319-i-m-trying-to-find-the-rmsd-of-the-numbers-by-taking-1-2-2-3-3-4-in-a-sequence-how-can-i-loop-as-su#answer_381164

在 MATLAB Online 中打开

No idea what RMSD is. In your formula,

RMSD = sqrt(sum(b-a).^2 / sum(a).^2);

none of the sum make any sense, since both a and b are scalar. sum has only one number to sum, so it is equivalent to:

RMSD = sqrt((b-a)^2 / a^2); %what was the point of the sums?

Assuming that it is the correct formula, to apply it to each element of your vector it would be:

N = [55378, 55344, 55310, 55276, 55242, 55208];
RMSD = sqrt((N(2:end) - N(1:end-1)).^ 2 ./ N(1:end-1).^2);

5 个评论
显示 3更早的评论隐藏 3更早的评论

Ramesh Bala 2019-6-28

在 MATLAB Online 中打开

Thanks Guillaume,

I added sum as the numbers are not actually single numerals in my case ,they are matrices with 1024 values,for which one has to load the files.

That's why to make it easier I added it as single numers in the question section.But as you've said sum isn't necessary here.

Now,lemme add the small mat files here and the code

A= load('18502.mat')
B= load('22172.mat')
C= load('24677.mat')
D=load('25892.mat')
a = ((A.ans));
b = ((B.ans));
RMSDD1 = sqrt ( sum(( b-a).^2 /sum(a).^2) );
fprintf ('the RMSD value between two sensor signals 1,2 @ 5,200 kHz :    %f\n' ,RMSDD1);
RMSDDP = RMSDD1*100;
fprintf ('the RMSD value between two sensor signals DC @ 5,200 kHz :    %f\n' ,RMSDDP);
a = ((B.ans));
b = ((C.ans));
RMSDD2 = sqrt ( sum(( b-a).^2 /sum(a).^2) );
fprintf ('the RMSD value between two sensor signals 1,2 @ 5,200 kHz :    %f\n' ,RMSDD2);
RMSDDP = RMSDD2*100;
fprintf ('the RMSD value between two sensor signals DC @ 5,200 kHz :    %f\n' ,RMSDDP);
  
a = ((C.ans));
b = ((D.ans));
RMSDD2 = sqrt ( sum(( b-a).^2 /sum(a).^2) );
fprintf ('the RMSD value between two sensor signals 1,2 @ 5,200 kHz :    %f\n' ,RMSDD2);
RMSDDP = RMSDD2*100;
fprintf ('the RMSD value between two sensor signals DC @ 5,200 kHz :    %f\n' ,RMSDDP);
% so can it made as such as per your previous comment ?
A= load('18502.mat')
B= load('22172.mat')
C= load('24677.mat')
D=load('25892.mat')
for N = load( '18502.mat','22172.mat','24677.mat','25892.mat')
RMSD = sqrt ((N (2: end)-N (1: end-1)).^ 2 ./N (1: end-1).^2);
RMSD
end

Ramesh Bala 2019-6-29

Thank you for the comment.

Is this

%calculate RMSD

RMSD = sqrt(sum((data(2:end, :) - data(1:end-1, :)) .^2, 2) ./ sum(data(1:end-1, :), 2).^2)

Actually this : RMSDD = sqrt ( sum(( b-a).^2) /sum(a.^2) );

It seems the values aren't matching??

Guillaume 2019-6-30

在 MATLAB Online 中打开

"Comments from guilame gives me : 833,1341,1711 as RMSD results ? But the actual results shld be 19,11,5."

My code gives you exactly the same output as the one you've posted, except that it is a lot more generic (will work with any number of files) and simpler.

Output from the code you've given above:

the RMSD value between two sensor signals 1,2 @ 5,200 kHz :    833.775403
the RMSD value between two sensor signals DC @ 5,200 kHz :    83377.540327
the RMSD value between two sensor signals 1,2 @ 5,200 kHz :    1341.762683
the RMSD value between two sensor signals DC @ 5,200 kHz :    134176.268293
the RMSD value between two sensor signals 1,2 @ 5,200 kHz :    1717.228148
the RMSD value between two sensor signals DC @ 5,200 kHz :    171722.814834

You get 833, 1341 and 1717 as well.

Of course, if the formula in your code was wrong, then the formula in mine is as well. I've just implemented exactly the same in more efficient way.

As Walter pointed out, it can be implemented a bit simpler as well, by using diff.

filenumbers = [18502, 22172, 24677, 25892];
%load data as numel(filenumbers) x N matrix
filecontents = arrayfun(@(fn) load(sprintf('%d.mat', fn)), filenumbers);  %this will error if the files don't all have exactly the same variable names
data = vertcat(filecontents.ans);  %this will error if the files don't have an ans variable or if the ans vectors are not the same length
%calculate RMSD
RMSD = sqrt(sum(diff(data, 1) .^2, 2) ./ sum(data(1:end-1, :), 2).^2)

请先登录，再进行评论。

Answer 2

Walter Roberson 2019-6-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/469319-i-m-trying-to-find-the-rmsd-of-the-numbers-by-taking-1-2-2-3-3-4-in-a-sequence-how-can-i-loop-as-su#answer_381286

在 MATLAB Online 中打开

sum(diff(data, 1).^2,1)

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

Ramesh Bala 2019-6-29

在 MATLAB Online 中打开

Alright still I get varying results ??

Comments from guilame gives me : 833,1341,1711 as RMSD results ?

Comment from Walter gives me : 1*1024 doubles ?

But the actual results shld be 19,11,5.

N = [18502, 22172, 24677, 25892]; % data
% filecontents = arrayfun(@(fn) load(sprintf('%d.mat', fn)), N);  %this will error if the files don't all have exactly the same variable names
% data = vertcat(filecontents.ans);  %this will error if the files don't have an ans variable or if the ans vectors are not the same length
% %calculate RMSD - guillame comments
% RMSD = sqrt(sum((data(2:end, :) - data(1:end-1, :)) .^2, 2) ./ sum(data(1:end-1, :), 2).^2);
% 
% %walter's comments
% % RMSD = sum(diff(data, 1).^2,1) % gives a 1024 double results ? instead of a scalar
a = (N(1));
b = (N(2));
RMSDD1 = sqrt ( sum(( b-a).^2) /sum(a.^2) );
RMSDDP1 = RMSDD1*100;
a = (N(2));
b = (N(3));
RMSDD2 = sqrt ( sum(( b-a).^2) /sum(a.^2) );
RMSDDP2 = RMSDD2*100;
a = (N(3));
b = (N(4));
RMSDD3 = sqrt ( sum(( b-a).^2) /sum(a.^2) );
RMSDDP3 = RMSDD3*100;

请先登录，再进行评论。

I'm trying to find the RMSD of the numbers by taking 1-2,2-3,3-4 in a sequence ,how can i loop as such ??

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

5 个评论
显示 3更早的评论隐藏 3更早的评论

更多回答（1 个）

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

另请参阅

类别

标签

产品

Community Treasure Hunt

I'm trying to find the RMSD of the numbers by taking 1-2,2-3,3-4 in a sequence ,how can i loop as such ??

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

5 个评论 显示 3更早的评论隐藏 3更早的评论

更多回答（1 个）

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

另请参阅

类别

标签

产品

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论

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