blob: 40b43f4a8e2e0d6486976efd61dca623dd8ba937 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
clc
//Example 16.2
disp('Given')
disp('R=40Kohm L=1H C=1/64 uF w=8.2krad/s')
R=40*10^3; L=1; C=1/64 *10^-6; w=8.2*10^3;
//The value of Q0 must be at least 5
Q0=5;
w0=1/sqrt(L*C)
printf("w0= %3.1f krad/s \n",w0*10^-3);
f0=w0/(2*%pi)
B=w0/Q0
printf("Bandwidth= %3.1f krad/s \n",B*10^-3);
//Number of half bandwidths be N
N=2*(w-w0)/B
disp(N)
//Admittance Y(s)=(1+i*N)/R
//Finding the magnitude and angle
magY=sqrt(1+N^2)/R
angY=atan(N)*(180/%pi)
disp(angY,'angY=')
printf("admittance value=%3.2f uS",magY*10^6)
|
