blob: f20c460f89e701e17e8ca2e357d89282d919fbf8 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
// Ex 65 Page 406
clc;clear;close;
// Given
C=6*10**-6;//F
L=2.5;//H
R=300;//ohm
a=R/2/L
omega = sqrt(1/L/C-R^2/4/L^2);//rad/s
//i=Im*%e**(-a*t)*sin(omega*t+fi)
//at t=0 sec
i0=0;//A
vc=10;//V
fi=asin(i0);//degree
//L*di/dt=vc at t=0
Im=poly([0],'Im')
function i=current(t)
i=Im*expm(-a*t)*sin(omega*t+fi)
endfunction
//i=Im*expm(-a*t)*sin(omega*t+fi)
LdiBYdt=L*numderivative(current,0)
temp = coeff(LdiBYdt)
Im=vc/temp(2)
Rn=2*sqrt(L/C);//ohm
Rad=Rn-R;//ohm
printf("Additional resistance required = %d ohm",Rad)
|
