Program Penyelesaian persamaan kuadrat (ABC)

disp('---------------------');
disp('Program Penyelesaian persamaan kuadrat (ABC)')
disp('dibuat Oleh : Mashudan Kurniawan')
disp('Hari/tanggal: sabtu/10 Okt 2012')
%-------------------------
%persamaan ax^2+bx+c, dimana:
a=1     %harga a
b=-5     %harga b
c=1     %harga c
%Deskripsi dan Perhitungan
%------------------------------
d=b^2-4*a*c       %D=diskriminan
x1=(-b+sqrt(d))/(2*a) %Akar x1
x2=(-b-sqrt(d))/(2*a) %Akar x2

MENU INPUT

disp('---------------------');
disp('Program Penyelesaian persamaan kuadrat (ABC)')
disp('dibuat Oleh : Mashudan Kurniawan')
disp('Hari/tanggal: sabtu/10 Okt 2012')
disp('---------------------------------')
%-------------------------
%persamaan ax^2+bx+c, dimana:
a=input('masukan nilai a =.....?')     %harga a
b=input('masukan nilai b =.....?')     %harga b
c=input('masukan nilai c =.....?')     %harga c
%Deskripsi dan Perhitungan
%------------------------------
d=b^2-4*a*c       %D=diskriminan
x1=(-b+sqrt(d))/(2*a) %Akar x1
x2=(-b-sqrt(d))/(2*a) %Akar x2

Read more >>

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

Read more >>

VEKTOR DAN MATRIK DENGAN MATLAB

VEKTOR DAN MATRIK DENGAN MATLAB

1.     FORMAT DATA

Dalam matlab dikenal dua jenis format data yaitu vektor dan matrik. Vektor atau matrik di definisikan dengan tanda kurung siku “ [ ] “untuk memisahkan kolom digunakan tanda koma ” , ” atau spasi tanda “ ; ” digunakan untuk memisahkan baris.

a.      Vektor

Vektor adalah sekelompok bilangan yang tersusun atas satu dimensi, dimensi ini bisa berupa baris atau kolom.

Contoh vektor baris:

>>a = [1,2,3]

a =

     1     2     3

>> b = [1 2 3]

b =

     1     2     3

Contoh vektor kolom:

>> a = [1;2;3]

a =

     1

     2

     3

Untuk mendefinisikan suatu vektor dengan jumlah elemen yang banyak tetapi memiliki interval yang teratur digunakan operator ” : “

>> a = 1:10

a =

     1     2     3     4     5     6     7     8     9    10

Untuk membuat suatu vektor baris dengan interval 2 sebagai contoh berikut:

>> a = 1:2:10

a =

     1     3     5     7     9

2.      MATRIK

            Matrik adalah sekelompok bilangan yang tersusun dalam segi empat dua dimensi,matrik bisa matrik bujur sangkar dan matrik segiempat.

a.      Menginputkan matrik

Ada empat cara untuk menginputkan matrik

Cara 1:

>> a=[1 2 3;4 5 6;7 8 9]

Cara 2:

>> a=[1 2 3

5 5 6

7 8 9]

Cara 3:

>> a1=[1 2 3];

>> a2=[4 5 6];

>> a3=[7 8 9];

>> a=[a1;a2;a3];

>> a

Cara 4:

>> a=input('masukan matrik= ');

masukan matrik= [1 2 3;4 5 6;7 8 9]

>> disp(a)

3.      Perkalian Matrik

1.      Perkalian Skalar

Untuk mengalikan satu matrik dengan bilangan tunggal (skalar) masing-masing elemen matrik harus dikalikan dengan faktor tersebut.

Contoh:

>> a = 4

a =

     4

>> b = [3 2 5;6 1 7]

b =

     3     2     5

     6     1     7

>> c = a*b

c =

    12     8    20

    24     4    28

2.      Perkalian Dua Matrik

Dua matrik dapat dikalikan satu sama lain apabila jumlah kolom dalam matrik pertama sama dengan jumlah baris matrik kedua.

Contoh:

>> a=[1 2 3;4 5 6;7 8 0]

a =

     1     2     3

     4     5     6

     7     8     0

>> b=[366;804;351]

b =

   366

   804

   351

>> det(a)

ans =

    27

>> inv(a)

ans =

   -1.7778    0.8889   -0.1111

    1.5556   -0.7778    0.2222

   -0.1111    0.2222   -0.1111

Penyelasaian dengan statement inv

Cara 1:

>> x=inv(a)*b

x =

   25.0000

   22.0000

   99.0000

Cara 2:

>> x=a\b

x =

   25.0000

   22.0000

   99.0000

Cara 3:

>> [L,U]=lu(a)

L =

    0.1429    1.0000         0

    0.5714    0.5000    1.0000

    1.0000         0         0

U =

    7.0000    8.0000         0

         0    0.8571    3.0000

         0         0    4.5000

>> x=inv(U)*inv(L)*b

x =

   25.0000

   22.0000

   99.0000

Cara lain:

>>  ABC = [6 -3 4;12 5 -7;-6 2 6]

ABC =

     6    -3     4

    12     5    -7

    -6     2     6

>> K=[41;-26;14]

K =

    41

   -26

    14

>> xyz=inv(ABC)*K

xyz =

    2.0035

   -2.6842

    5.2316

MODIFIKASI MATRIK

Kita dapat mengganti elemen dari matrik yang telah didefinisikan hal ini dapat berguna apabila kita bekerja dengan matrik yang sangat besar. Kemudian kita ingin mengganti eleman matrik  pada baris dua kolom tiga dengan angka 10.

>> d=[1 2 3;4 5 6]

d =

     1     2     3

     4     5     6

>> d(2,3)=10

d =

     1     2     3

     4     5    10

TRANSPOSE MATRIK

Transpose matrik merupakan salah satu operasi penting dalam matrik, dengan mengetik tanda kutip tunggal “ ‘ “ bisa digunakan untuk mengubah bentuk elemen baris menjadi kolom.

>> d

d =

     0     5     3

     4     5    10

>> d'

ans =

     0     4

     5     5

     3    10

BILANGAN KOMPLEK

Bilangan komplek adalah suatu bilangan yang terdiri dari atas dua suku yang terpisah sering ditulis sebagai:

z=a+jb

dimana:

a=bilangan real dari z

b=bilangan imajiner dari z

bilangan imajiner sering ditulis juga Me^jᶿ

dimana:

M adalah magnitude, M=

0 adalah sudut, =tan^-1(b/a)

BENTUK POLAR SUATU BIL. KOMPLEKS

Contoh:

>> z

z =

   3.0000 + 4.0000i

>> abs(z)

ans =

     5

>> angle(z)

ans =

    0.9273

>> real(z)

ans =

     3

>> imag(z)

ans =

     4

>> conj(z)

ans =

   3.0000 - 4.0000i

Read more >>

Teknik Pengoperasian kubikel 20 KV

Penamaan Kubikel 20 KV

Disebut kubikel karena peralatan-peralatan tersebut di kemas plat blok berbentuk almari dengan pintu di bagian yang bisa dibuka dan ditutup menurut standar operasi yang diminta.

            Dikatakan kubikel 20 KV atau switchgear 20 KV, karena peralatan tersebut bekerja pada tegangan nominal fasa-fasa 20 KV, yang termasuk kategori teganga menengah (> 1 KV hingga 35 KV) dan kubikel 20 KV banyak digunakan di lapangan.

Teknik Pengoperasian kubikel 20 KV

Penyaluran daya

            Sumber tegangan menengah yang mensuplai kubikel 20 KV di laboratorium mekanik teknik listrik diperoleh dari transformator step-up (380 V/20 KV). Daya listrik yang disalurkan  dari panel hubung bagi sentral politeknik sebesar 380 V masuk ke box panel PMS (pemisah). Dari panel PMS ini keluarannya masuk ke ...(380 V) dari transformator step-up, yang sekundernya (tegangan 20 KV) dimasukkan ke incoming pada kubikel.

            Keluaran kubikel (outgoing) dihubungkan ke transformator step-down (20 KV/380 V), yang selanjutnya akan digunakan ke beban-beban melalui keluaran tegangan rendah transformator step-down atau tegangan menengah langsung dari keluaran kubikel.

Bagian-bagian utama kubikel 20 KV

            Bagian-bagian utama alat ini terbagi 3 yaitu:

Incoming, Metering dan Outgoing. Masing-masing bagian memiliki kontak DS (Disconecting switch) da bagian kontak pentanahan (grounding).

Incoming

Incoming adalah tempat masuknya tegangan dari sumber sebesar 20 KV yang mana pada incoming trsebut terdapat 3 buah lampu indikator neon (R,S,T) yang akan menyala apabila tegangan dari luar masuk ke incoming.

Pada incoming juga terdapat 3 buah trafo arus yangkan digunakan oleh KWH meter pada incoming dengan kode P7 untuk mengukur daya yang dikeluarkan . Pada beberapa kubikel terdapat jenis pemasangan transformator arus (CT) yaitu:

1.      Current Transformator Double Secondary

2.      Current Transformator Single Secondary

Selain itu pada incoming juga terdapat heater (pemanas) yang berfungsi untuk memanaskan ruang atau kompartemen pada kubikel sehingga tidak lembab dan dapat menghindari terjadinya percikan bunga api.

Metering

      Metering adalah suatu bagian dari kubikel yang memiliki kemampuan pengukuran besarn-besarn listrik. Pada metering terdapat alat-alat ukur pengukuran seperti voltmeter,amperemeter,dan KV meter. Pada metering terdapat fuse-fuse pengaman sebanyak 3 buah buatan merlin gerin dengan data fuse HV 20 KV;6,3 A; W/O Striker. Tepat di bawah fuse-fuse ini terdapat 3 buah transformator tegangan yang menurunkan tegangan dari 20 KV menjadi 100 V AC yang akan digunakan oleh alat-alat ukur seperti kv meter dan amperemeter serta voltmeter. Selain fuse-fuse juaga terdapat MCB 3 fasa yang akan mengamankan hubungan ke material-material pengukuran.

Outgoing

      Outgoing berfungsi sebagai tempat keluarnya tegangan menengah setelah melalui incoming dan metering. Ougoing memiliki kompartemen yang paling besar dimana pada kompartemen bagian atas terdapat kontak grounding dan juga terdapat kontak disconecting switch DS. Pada bagian atas juga terdapat 2 buah kunci dan satu kunci lagi pada bagian bawah kompartemen.

Kunci-kunci bekerja secara interlock dangan tipe kunci 50. Keterangan mengenai ketiga kunci itu antara lain:

1.      Missing key

2.      Free key

3.      Described key

Pada bagian bawah terdapat tuas yang akan menekan pegas yang ada di dalamnya. Selain menggunakan pegas juga bisa dilakukan otomatis dengan motorized.

Bagian-bagian kubikel

Read more >>