How to fix looping?

Question

Devia Rafika Putri 2013-5-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75804-how-to-fix-looping

采纳的回答： Andrei Bobrov

this is my program:

v = 120;

r1 = 1.8;

x1 = 2.4;

r2 = 3.5;

x2 = 1.2;

xm = 60;

ns = 1800;

ws = 188.5;

s = (0.5:1:50)/50;

s (1) = 0.0001;

nm = (1-s)*ns;

for i = 1:51;

zf(i) =(((r2/s(i))+(j*x2))*(j*xm))/(((r2/s(i))+(j*x2))+(j*xm));

zb(i) = (((r2/(2-s(i)))+(j*x2))*(j*xm))/(((r2/(2-s(i)))+(j*x2))+(j*xm));

I(i) = v/(r1+(j*x1)+(0.5*zf(i))+(0.5*zb(i)));

PagF(i) = (abs(I(i)^2))*(0.5*real(zf(i)));

PagB(i) = (abs(I(i)^2))*(0.5*real(zb(i)));

Pag(i)=PagF(i)-PagB(i);

Tind(i) = Pag(i)/ws;

end

figure (1);

plot(nm,Tind,'Color','b','LineWidth',2.0);

grid on;

hold off;

*but after i run this program an error occurred with argument: Attempted to access s(51); index out of bounds because numel(s)=50.

Error in ==> fasbel3 at 15 zf(i) =(((r2/s(i))+(j*x2))*(j*xm))/(((r2/s(i))+(j*x2))+(j*xm));

how to fix it? please help thanks.. with sincerity :)

4 个评论
显示 2更早的评论隐藏 2更早的评论

Craig Cowled 2013-5-15

Devia, Andrei's code is much neater. No need for a loop, just use element wise divide './'. And no, I'm not an electrical engineer. I'm a structural engineer working on experimental structural dynamics problems.

Devia Rafika Putri 2013-5-15

ok craig.. i see.. but some code from andrei, i still dont understand.. but i'll find it..

ooh, I thought you were an electrical engineer

because I'm having problems with programs relating to single phase motor..

thanks so much craig..you are very helpful..

请先登录，再进行评论。

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

Answer 1

Andrei Bobrov 2013-5-15

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/75804-how-to-fix-looping#answer_85474

编辑：Andrei Bobrov 2013-5-15

在 MATLAB Online 中打开

v = 120;
r1 = 1.8;
x1 = 2.4;
r2 = 3.5;
x2 = 1.2;
xm = 60;
ns = 1800;
ws = 188.5;
s = [.0001;(.5:50).'/50];
s1 = [2-s,s];
z =(r2./s1+1i*x2)*1i*xm./( r2./s1+1i*(x2+xm) );
I = v./(r1+1i*x1+mean(z,2));
Tind = diff(bsxfun(@times,abs(I.^2),real(z)*.5),1,2)/ws;
plot((1-s)*ns,Tind,'Color','b','LineWidth',2.0);
grid on;

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

Andrei Bobrov 2013-5-15

在 MATLAB Online 中打开

Hi Devia!

Please read about functions of the MATLAB: bsxfun, diff.

Rewrite code:

s = [.0001;(.5:50).'/50];
s1 = [2-s,s];
z =(r2./s1+1i*x2)*1i*xm./( r2./s1+1i*(x2+xm) );
I = v./(r1+1i*x1+mean(z,2));
Tind = diff(bsxfun(@times,abs(I.^2),real(z)*.5),1,2)/ws;

as:

s = [.0001;(.5:50).'/50];
s1 = [2-s,s];
z =(r2./s1+1i*x2)*1i*xm./( r2./s1+1i*(x2+xm) );
I = v./(r1+1i*x1+mean(z,2));
Ia2 = abs(I.^2);
r = real(z)*.5;
P = bsxfun(@times,Ia2,r); % as in your code: P = [PagB PagF];
Pag = diff(P,1,2); % as Pag = PagF - PagB;
Tind = Pag/ws;