Welcome I have a CF equation that performs a three-dimensional summation process and gives a two-dimensional result I did this process using (For Loops) in MATLAB. I want to c

Question

hisham mohammed 2023-3-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1937714-welcome-i-have-a-cf-equation-that-performs-a-three-dimensional-summation-process-and-gives-a-two-dim

编辑： Dyuman Joshi 2023-3-30

采纳的回答： Dyuman Joshi

clc

xm(:,:,1)= [1 2 3 4;2 3 4 5;3 4 5 6;4 5 6 7 ]

xm(:,:,2)= [2 3 4 5;3 4 5 6;4 5 6 7;5 6 7 8]

tic

[cf] =coherence_factor(xm)

toc

function [cf] = coherence_factor(xm)

for i=1:size(xm,1);

for j=1:size(xm,2);

n=size(xm,3);

cf(i,j)=(abs (sum(xm(i,j,:))).^2) / (n* sum(abs(xm(i,j,:).^2)));

end

2 个评论
显示无隐藏无

Dyuman Joshi 2023-3-30

"I want to c"

You want to what? Change the way questions on MATLAB Answers are asked?

hisham mohammed 2023-3-30

I want to remove the (for loop)in this code and use matrix oprtion for this code how i can do it !

请先登录，再进行评论。

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

Answer 1

Dyuman Joshi 2023-3-30

编辑：Dyuman Joshi 2023-3-30

在 MATLAB Online 中打开

%Data
xm(:,:,1)= [1 2 3 4;2 3 4 5;3 4 5 6;4 5 6 7 ];
xm(:,:,2)= [2 3 4 5;3 4 5 6;4 5 6 7;5 6 7 8];
n=size(xm,3);
%Loop Method
for i=1:size(xm,1);
    for j=1:size(xm,2);
        cf1(i,j)=(abs (sum(xm(i,j,:))).^2) / (n* sum(abs(xm(i,j,:).^2)));
    end
end
cf1
cf1 = 4×4
    0.9000    0.9615    0.9800    0.9878
    0.9615    0.9800    0.9878    0.9918
    0.9800    0.9878    0.9918    0.9941
    0.9878    0.9918    0.9941    0.9956
%Vectorized method
cf2 = abs(sum(xm,3).^2)./(n*sum(abs(xm.^2),3))
cf2 = 4×4
    0.9000    0.9615    0.9800    0.9878
    0.9615    0.9800    0.9878    0.9918
    0.9800    0.9878    0.9918    0.9941
    0.9878    0.9918    0.9941    0.9956