From f35ea80659b6a49d1bb2ce1d7d002583f3f40947 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:38:01 +0530 Subject: updated the code --- 339/CH1/EX1.4/ex1_4.sce | 31 ++++++++++++++++--------------- 1 file changed, 16 insertions(+), 15 deletions(-) (limited to '339/CH1/EX1.4/ex1_4.sce') diff --git a/339/CH1/EX1.4/ex1_4.sce b/339/CH1/EX1.4/ex1_4.sce index 58784b047..8386964c0 100755 --- a/339/CH1/EX1.4/ex1_4.sce +++ b/339/CH1/EX1.4/ex1_4.sce @@ -1,16 +1,17 @@ -f=10^6:10^7:10^10; -rs=(4.8*10^-6).*sqrt(f); -re=(33.9*10^12) ./f; -c=47*10^-12; -w=2*%pi.*f; -l=2*1.25*10^-2; -a=2.032*10^-4; -temp=log(2*l/a)/log(%e); -lex=mu0*l*(temp-1)/(2*%pi); //external inductance -z=1 ./(1 ./re +w*c*%i)+rs+w.*lex*%i; // impedance of frequency dependent capacitor -zideal=1 ./(w*c*%i); //impedance of an ideal capacitor -plot2d("gll",f,abs(z)); -plot2d(f,abs(zideal)); -title("Frequency responce of a high frequency capacitor"); -xlabel('Frequency (f) in Hz'); +f=10^6:10^7:10^10; +mu0=4*%pi*10^-7; +rs=(4.8*10^-6).*sqrt(f); +re=(33.9*10^12) ./f; +c=47*10^-12; +w=2*%pi.*f; +l=2*1.25*10^-2; +a=2.032*10^-4; +temp=log(2*l/a)/log(%e); +lex=mu0*l*(temp-1)/(2*%pi); //external inductance +z=1 ./(1 ./re +w*c*%i)+rs+w.*lex*%i; // impedance of frequency dependent capacitor +zideal=1 ./(w*c*%i); //impedance of an ideal capacitor +plot2d("gll",f,abs(z)); +plot2d(f,abs(zideal)); +title("Frequency responce of a high frequency capacitor"); +xlabel('Frequency (f) in Hz'); ylabel('Absolute impedance (|Z|) in ohms'); \ No newline at end of file -- cgit