# Printing loop(element wise multiplications) output

1 view (last 30 days)
Arif Sayyad on 25 Jan 2020
Commented: Arif Sayyad on 25 Jan 2020
clear all
clc
vr=40;
u=0;
S=12;
P=4;
m=5;
q=1;
V=700;
f=50;
Q=666666.666
Btooth=1.8;
Biron=1.6;
s=0.1;
Kw=0.955;
pf_n = 0.35:0.05:0.5;
ac=10000:10000:40000;
Bavg=0.45:0.05:0.6;
C = (1.11 .* Bavg .* Kw .* pf_n.*ac);
vs=vr/(1-s);
L_W = Q./(C.*vs);
W = sqrt(L_W./P);
L = P .* W;
fprintf('\n\n\npf_n=%f,\nac=%f,\nBavg=%f,\nL_W=%f',pf_n,ac,Bavg,L_W)
plot(L_W,pf_n,'Color','k','LineWidth',2.0);
xlabel('L_W','Fontweight','Bold');
ylabel('pf_n','Fontweight','Bold');
title('L_W vs pf_n');
hold on

KALYAN ACHARJYA on 25 Jan 2020
Edited: KALYAN ACHARJYA on 25 Jan 2020
vr=40;
u=0;
S=12;
P=4;
m=5;
q=1;
V=700;
f=50;
Q=666666.666;
Btooth=1.8;
Biron=1.6;
s=0.1;
Kw=0.955;
pf_n=0.35:0.05:0.5;
ac=10000:10000:40000;
Bavg=0.45:0.05:0.6;
vs=vr/(1-s);
for i=1:length(Bavg)
C(i)=(1.11.*Bavg(i).*Kw.*pf_n(i).*ac(i));
L_W(i)=Q/(C(i)*vs);
W(i)=sqrt(L_W(i)/P);
L(i)=P*W(i);
fprintf('\n\n\npf_n=%f,\nac=%f,\nBavg=%f,\nL_W=%f',pf_n(i),ac(i),Bavg(i),L_W(i))
L_W(i)
pf_n(i)
plot(L_W(i),pf_n(i),'*','linewidth',3);
hold on;
end
xlabel('L_W','Fontweight','Bold');
ylabel('pf_n','Fontweight','Bold');
title('L_W vs pf_n');
grid on;

Arif Sayyad on 25 Jan 2020
Thank you very much sir,
it is working.
but i didnt got actually you have done...?
and my another question is if i want to plot Bavg also..
then what shoud i do ?
KALYAN ACHARJYA on 25 Jan 2020
Add or use same plot command
plot(Bavg(i),pf_n(i),'--','linewidth',3);
Please note, generally variable in x-axis and variant in y-axis
Arif Sayyad on 25 Jan 2020
thank you sir ... :)