Can my code be written like this?

Question

yang-En Hsiao 2019-1-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/442434-can-my-code-be-written-like-this

编辑： Jan 2019-1-31

在 MATLAB Online 中打开

i want to write this formula as code in matlab,but i am not very sure that is my code right.

The range of j and k are both 1~K

for k= 1:K
     for j=1:K
        if k==j
            a=mu(k).*( abs(h_(k)'*f_(k)^2 ) / hat_r(k);
        else
           b=mu(k).*( abs(h_(k)'*f_(j)^2 );
        end      
    end
    total=a+b;
end

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

Jan 2019-1-31

编辑：Jan 2019-1-31

在 MATLAB Online 中打开

Neither the elementwise .* does not matter, when you work with scalars. Using conj looks less confusing than using the transposition.

In abs(h_(k)'*f_(k)^2) you are squaring f only. Use this instead:

abs(conj(h_(k)) * f_(k))^2

yang-En Hsiao 2019-1-31

K_k means that the range of j is the same as k,that is 1 ~K.but j can't not be equal to k

请先登录，再进行评论。

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

Answer 1

Jan 2019-1-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/442434-can-my-code-be-written-like-this#answer_358847

编辑：Jan 2019-1-31

在 MATLAB Online 中打开

There is no "hat_r(k)" in your formula. A trailing parenthesis is missing also.

total = 0;
hc    = conj(h_);
for k = 1:K
    c1   = ;
    sumj = 0;
    for j = 1:K
        if j ~= k   % [EDITED]
          sumj = sumj + abs(hc(k) * f_(j)) ^ 2;
        end
    end
    total = total + mu(k) * (abs(hc(k) * f_(k)) ^ 2 - sumj) / r_hat(k);
end

2 个评论
显示无隐藏无

yang-En Hsiao 2019-1-31

it should be bar_r(k).there is a bar_r(k) ,1.,in the picture.

and the 2 means that the range of j is the same as k,that is 1 ~K.but j can't not be equal to k

Jan 2019-1-31

编辑：Jan 2019-1-31

If "it" should be "bar_r(k)", than simply insert this in my code. You do not have to wait until I do it. By the way, the character is a greek "gamma", not an "r". But names do not matter.

There is still a missing closing parenthesis. Please post a valid formula. A smaller zoom-level would be nice.

请先登录，再进行评论。