blob: 8cfbb9abde278a5a59b680d8f73a71f8f4c1ab12 (
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
29
30
31
32
33
34
35
36
37
38
39
40
41
|
clc
disp("Example 2.2")
printf("\n")
printf("Given")
disp("Inductance used is 30mH")
disp("Current flow is i=10*sin(50*t)")
L=30*10^-3;iamp=10;
t=0:0.01:0.06283;
i=10*sin(50*t)
//v=L*d/dt(i)
//d/dt(sin 50t)=50*cos t
vamp=L*iamp*50;
v=vamp*cos(50*t)
//sinA*cosB=(sin(A+B)+sin(A-B))/2
pamp=vamp*iamp/2;
p=pamp*sin(100*t)
//On integrating 'p' w.r.t t
wL=0.75*(1-cos(100*t));
subplot(221)
plot(t,i)
xtitle ('i vs wt','wt','i');
subplot(222)
plot(t,v)
xtitle ('v vs wt','wt','v ');
subplot(223)
plot(t,p)
xtitle ('p vs wt','wt','p ');
subplot(224)
plot(t,wL)
xtitle ('wL vs wt','wt','wL ');
|
