How to fix looping?

1 view (last 30 days)
Devia Rafika Putri
Devia Rafika Putri on 15 May 2013
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 Comments
Devia Rafika Putri
Devia Rafika Putri on 15 May 2013
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..

Sign in to comment.

Accepted Answer

Andrei Bobrov
Andrei Bobrov on 15 May 2013
Edited: Andrei Bobrov on 15 May 2013
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 Comments
Devia Rafika Putri
Devia Rafika Putri on 15 May 2013
Hi Andrei, thank you so much for your explain.. i'm so grateful.
but actually I have another problem. is about 1 phase motors, can you help me?
if you can, I will send my problem to your email..
thanks so much, Andrei..
with sincerity Devia

Sign in to comment.

More Answers (1)

Yao Li
Yao Li on 15 May 2013
It seems the length of array s is 50, but you wanna call s(51) in the for loop
  4 Comments
Devia Rafika Putri
Devia Rafika Putri on 15 May 2013
lol.. thanks so much Yao li, I've sent an email to you .. if you're not busy please check it..

Sign in to comment.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by