Have implemented for loop corectly in my code?

Question

Nisar Ahmed 2022-5-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1720285-have-implemented-for-loop-corectly-in-my-code

编辑： KALYAN ACHARJYA 2022-5-16

采纳的回答： KALYAN ACHARJYA

在 MATLAB Online 中打开

Hi,

I want to put a for loop in my code: can you please confirm if I have done it correctly? Here is my code:

f_min = 1;                          % Minimum frequency
f_max = 200;                        % Maximum frequency
f_int = 1;                          % Interval
%f = f_min:f_int:f_max;              % Frequency range
 for ii = f_min:f_int:f_max
     Logf = log10(f);
     w = 2*pi*f;
     
     I1 = sqrt(1i.*w.*s1).*coth(1i.*w.*s1./2);
     I2 = sqrt(1i.*w.*s2).*coth(1i.*w.*s2./2);
     
    b = 1./(1+1./(I1.*g1 + I2.*g2)); 
    
    E = Eo.*b;
    Q = real(E)./imag(E);
    
    Qinv = 1./Q;
    
    v=sqrt(E./rho_b); 
    
    VP = (real(1./v)).^-1;     
 end

2 个评论
显示无隐藏无

Star Strider 2022-5-16

I do not understand wanting to use a loop when everything already appears to be vectorised. (To use a loop, all the vectors need to be subscripted using the ‘ii’ subscript reference.)

Also this:

VP = (real(1./v)).^-1

does not appear to be appropriate, since it inverts the inversion. Why not just:

VP = real(v);

.

Nisar Ahmed 2022-5-16

@Star Strider I want to change VP as a function of f and w = 2 pi f, So to compute as a function of w, I need to apply loop.

请先登录，再进行评论。

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

Answer 1

KALYAN ACHARJYA 2022-5-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1720285-have-implemented-for-loop-corectly-in-my-code#answer_964805

编辑：KALYAN ACHARJYA 2022-5-16

在 MATLAB Online 中打开

I don't think, is there any loop needed, it can be avoided, just define the following parameters

s1=..
s2=..
g1=..
g2=..
Eo=..
rho_b=..

More the code:

f_min = 1;                          % Minimum frequency
f_max = 200;                        % Maximum frequency                         
f = f_min:f_max;                    % Frequency range
Logf= log10(f);
w=2*pi*f;
I1 = sqrt(1i.*w.*s1).*coth(1i.*w.*s1./2);
I2 = sqrt(1i.*w.*s2).*coth(1i.*w.*s2./2);
b = 1./(1+1./(I1.*g1 + I2.*g2)); 
E = Eo.*b;
Q = real(E)./imag(E);
Qinv=1./Q;
v=sqrt(E./rho_b); 
VP=(real(1./v)).^-1; 