MEMBUAT MENU PILIHAN DENGAN MATLAB

clear all
clc
disp ('***************************************************************')
disp ('               PENYELESAIAN SOAL TRANSMISI                     ')
disp ('***************************************************************')
disp ('NAMA              = MASHUDAN KURNIAWAN                         ')
disp ('NIM               = 061030310851                               ')
disp ('KELAS             = 5 ELB                                      ')
disp ('MATA KULIAH       = PEMROGRAMAN KOMPUTER                       ')
disp ('DOSEN PEMBIMBING  = ANTON FIRMANSYAH,S.T.,M.T.                 ')
disp ('HARI/TANGGAL      = SABTU,20 NOVEMBER 2012                     ')
disp ('---------------------------------------------------------------')
disp ('silahkan pilih penyelesaian soal di bawah ini                  ')
disp ('---------------------------------------------------------------')
disp ('PILIHAN')
disp ('---------------------------------------------------------------')
disp ('1.  soal 1')
disp ('2.  soal 2')
disp ('3.  soal 3')
disp ('4.  selesai')
disp ('---------------------------------------------------------------')
pilih=input('masukkan pilihan Anda =  ');
disp ('---------------------------------------------------------------')
while pilih ~=3
switch pilih
case 1
disp('SOAL 1')
disp('Suatu saluran transmisi fasa tiga,50 km,70 KV,mempunyai konstanta:')
disp('R=0,20 ohm per km;X=0,608 ohm per km; Y=j4,0 X 10^-6 mho per km.')
disp('Saluran transmisi tersebut mensuplai beban 30 MW dan faktor daya pf=0,9 terbelakang')
disp('Tegangan pada ujung beban 70 KV.Hitunglah:')
disp('a.tegangan pada ujung kirim ')
disp('b.daya pada ujung kirim')
disp('c.efisiensi transmisi')
disp('d.pengaturan tegangan')
disp('Penyelesaian')
disp(' ')
disp('Is = IR = I ')
disp('VS = VR + IZ')
Pr=30;
V=70;
VR=(V/sqrt(3))*(1000);
disp(['PR = ',num2str(Pr),(' MW'),' ;pf = 0,9 terbelakang'])
disp('                       __')
disp('VR = 70 KV(L-L)= 70 / /3 KV (L-N) ')
disp(['               = ',num2str(VR),' Volt(L-N)'])
S=50;
R=0.2;
X=0.608;
Y=(j*4)*1e-6;
PR=30000;
pf=0.9;
%jawab:
%Soal a
disp(' ')
disp('a. Karena panjang saluran 50 km,jadi ini adalah saluran pendek')
%
IR=(Pr)*(1000)/(sqrt(3)*(V)*(pf));
sudutIR=acosd(pf);
%
z=complex(R,X)*50;
Z=sqrt(real(z)^2+imag(z)^2);
sudutZ=atand(imag(z)/real(z));
%
IZ=(IR*Z);
sudutIZ=(sudutZ-sudutIR);
z1=complex(IZ*cosd(sudutIZ),IZ*sind(sudutIZ));
VS=VR+real(z1);
image_z1=imag(z1);
z2=complex(VS,image_z1);
Vs=(sqrt(real(z2)^2+imag(z2)^2))/1000;
sudut_VS=atand(imag(z2)/real(z2));
disp (['   Vs=',num2str(Vs),(' /__ '),num2str(sudut_VS),(' KV(L-N)')]);
disp ('   atau')
VS=((sqrt(real(z2)^2+imag(z2)^2))/1000)*(sqrt(3));
disp (['   |Vs|=',num2str(VS),(' KV(L-L)')]);
disp (['   Jadi Tegangan pada ujung kirim(Vs)=',num2str(VS),(' KV(L-L)')]);
disp (' ')
%Soal B
teta_s=(sudut_VS)+sudutIR;
Ps=((sqrt(3))*(VS)*(IR)*(cosd(teta_s)))/(1000);
disp (['(B)Daya pada ujung kirim(Ps)=',num2str(Ps),(' MW')]);
disp (' ')
%Soal C
Eff=(Pr/Ps)*100;
disp (['(C)Efisiensi Transmisi (Eff)=',num2str(Eff),('%')]);
disp (' ')
%Soal D
VRF=((VS-V)/(V))*(100);
disp (['(D)Pengaturan Tegangan VR(%)=',num2str(VRF),('%')]);
pause;
return
%
case 2
disp('Soal No 2')
disp('Suatu saluran transmisi fasa tiga, 200 km, 220 KV. Konstanta')
disp('Kawat:Z=0,64 sudut71,8 ohm/km; Y=4e-6 sudut90 mho.km.Saluran')
disp('Transmisi di atas menyalurkan daya sebesar 100 MW dan pf=1,0')
disp('Pada ujung beban.Tegangan pada ujung beban 220KV.dengan meng-')
disp('Gunakan representasi nominal PI')
disp('Tentukanlah:')
disp('a.Tegangan dan arus pada ujung kirim')
disp('b.Efisiensi Transmisi')
disp('c.Pengaturan Tegangan')
disp('Penyelesaian')
%dik:
Pr=100;
pf=1;
y=4e-6;
sudut_y=90;
V=220;
s=200;
z=0.64;
sudut_z=71.8;
%Jawab
IR=(Pr)*(1000)/(sqrt(3)*V);
sudutIR=atand(pf);
VR=V/sqrt(3);
Z=z*s;
Y=s*y;
zy_2=(Z)*(Y/2);
sudutzy=sudut_y+sudut_z;
a=complex(zy_2*cosd(sudutzy),zy_2*sind(sudutzy));
b=1+real(a);
zy2=sqrt((b)^2+imag(a)^2);
sudut_zy2=atand(imag(a)/real(b));
zy_4=(Z)*(Y/4);
c=complex(zy_4*cosd(sudutzy),zy_4*sind(sudutzy));
d=1+real(c);
zy4=sqrt((d)^2+imag(c)^2);
sudut_zy4=atand(imag(c)/real(d));
%Tegangan pada ujung kirim
VS1=(zy2)*(VR*1000);
VS2=(Z)*(IR);
z1=complex(VS2*cosd(sudut_z),VS2*sind(sudut_z))+complex(VS1*cosd(sudut_zy2),VS1*sind(sudut_zy2));
VS=((sqrt(real(z1)^2+imag(z1)^2))*sqrt(3))/1000;
sudut_VS=atand(imag(z1)/real(z1));
disp (['(A)Tegangan pada ujung kirim(Vs)=',num2str(VS),(' sudut'),num2str(sudut_VS),(' KV(L-L)')]);
%Arus pada ujung kirim
IS1=zy4*Y*(VR*1000);
IS2=zy2*IR;
sudut_IS1=sudut_y+sudut_zy4;
z2=complex(IS1*cosd(sudut_IS1),IS1*sind(sudut_IS1))+complex(IS2*cosd(sudut_zy2),IS2*sind(sudut_zy2));
IS=sqrt(real(z2)^2+imag(z2)^2);
sudut_IS=atand(imag(z2)/real(z2));
disp (['   Daya pada ujung kirim    (Ps)=',num2str(IS),(' sudut'),num2str(sudut_IS),(' Ampere')]);
teta_S=sudut_VS-sudut_IS;
Ps=(3*(sqrt(real(z1)^2+imag(z1)^2))/1000*IS*cosd(teta_S))/1000;
disp (' ')
%Efisiensi transmisi
Eff=(Pr/Ps)*100;
disp (['(B)Efisiensi Transmisi (Eff)=',num2str(Eff),('%')]);
disp (' ')
%Pengaturan Tegangan
VR_NL=((sqrt(real(z1)^2+imag(z1)^2))/1000)/zy2;
VRF=((VR_NL-VR)/(VR))*100;
disp (['(C)Pengaturan Tegangan VR(%)=',num2str(VRF),('%')]);
pause;
return
%
case 3
 disp('Soal 3')
disp('Penentuan besar beban pada rugi-rugi tertentu. Tentukanlah jumlah')
disp('Daya pada ujung beban,  bila  faktor daya beban cos teta_R , yang')
disp('Dapat ditransmisikan melalui saluran pendek dengan tegangan ujung')
disp('Kirim Vs dan rugi-rugi transmisi yang diijinkan L?')
disp('Jawab:')
%SOAL A
%dik:
l=20;
z=complex(0.263,0.765);
cos_teta=0.906;
teta_R=-25;
teta=acosd(cos_teta);
Vs_LL=30;
Vs_LN=19;
%jawab:
Zr=(100*real(z)*l)/(10*cos_teta);
a=19000+j*0;
b=0.808+j*0;
c=complex(20*cosd(71),20*sind(71));
d=complex(58.1*cosd(25),58.1*sind(25));
e=1+j*0;
z=a/((b*c/d)+1);
f=complex(real(z),imag(z));
Vr=(sqrt(real(f)^2+imag(f)^2))/1000;
sudut_Vr=atand(imag(f)/real(f));
Ir=(Vr*1000)/Zr;
%jadi
Pr_fasa=((Vr*Ir)*cosd(-25));
Pr_3fasa=((Vr*Ir)*cosd(-25))*3;
teta_S=teta_R+sudut_Vr;
Ps_fasa=(Vs_LN*Ir)*cosd(teta_S);
Ps_3fasa=(Vs_LN*Ir)*cosd(teta_S)*3;
eff=Pr_fasa/Ps_fasa;
rugi_rugi=Ps_fasa-Pr_fasa;
Vr_persen=((Vs_LN-Vr)/Vr)*100;
disp(['(A)Tegangan Pada Ujung Beban(Vr)=',num2str(Vr),(' KV(L-L)')]);
disp(['   Daya Pada Ujung Beban(Pr)    =',num2str(Pr_3fasa),(' KW/3fasa')]);
disp(['   Daya Pada Ujung kirim(Ps)    =',num2str(Ps_3fasa),(' KW/3fasa')]);
disp(['   Efisiensi                    =',num2str(eff),]);
disp(['   Rugi-rugi                    =',num2str(rugi_rugi),(' KW/fasa')]);
disp(['   Pengaturan Tegangan VR(%)    =',num2str(Vr_persen),(' %')]);
%SOAL B
Y=((5.58)*(10^4));
thet=90;
disp('Y=5.58*10^4/__90 ')
Vs=66;
disp(['VS = ',num2str(Vs),'KV (L-L)'])
Zr=(100*0.263*100)/(15*cos_teta);
sudutZr=25;
Zr1=200;
disp(['ZR = ',num2str(Zr),' /__25 ohm'])
Zs=40.4;
sudutZs=71;
disp(['ZS = ',num2str(Zs),' /__71 ohm'])
Vs_a=(66/sqrt(3));
a1=1;
r1=40.4;
theta1=71;
z1_a=cosd(theta1)*r1;
z1_b=sind(theta1)*r1;
z1=complex(z1_a,z1_b)
%
z2_a=cosd(thet)*Y;
z2_b=sind(thet)*Y;
z2=complex(z2_a,z2_b)
brew=((a1)+(z1)+(z2));
bre=(brew/4);
br=brew+bre;
%
r3=80.8;
tet=71
z3_a=cosd(tet)*r3;
z3_b=sind(tet)*r3;
z3=complex(z3_a,z3_b)
%

z4_a=cosd(sudutZr)*Zr1;
z4_b=sind(sudutZr)*Zr1;
z4=complex(z4_a,z4_b)
%
bro=(z3/z4);
%
Has=(Vs_a)/(br);
hasil=(Has*bro);
%
IR=sqrt((real(hasil))^2+(imag(hasil)^2))
sudut_Z1=atand(imag(hasil)/real(hasil))


pause;
return
end
end