summaryrefslogtreecommitdiff
path: root/3802/CH1/EX1.10/Ex1_10.sce
blob: 49d7070a234ef0976e3d7255fec6fd7ce4dd0e76 (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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
//Book Name:Fundamentals of Electrical Engineering
//Author:Rajendra Prasad
//Publisher: PHI Learning Private Limited
//Edition:Third ,2014

//Ex1_10.sce.

clc;
clear;
subplot(2,2,1)
t=[0:0.001:8];
x=length(t);
v=ones(1,x);
for n=1:x;
    L=5;
    if t(n)<=2
       v(n)=6.25;
    else if t(n)>=6 & t(n)<8
           v(n)=-6.25;
        else
            v(n)=0;
            end
    end 
    end
xlabel("Time in seconds")
ylabel("voltage in volts")
title("voltage waveform")
plot(t,v)
subplot(2,2,2)
t=[0:0.001:8];
x=length(t);
p=ones(1,x);
for n=1:x;
    if t(n)<=2
       v(n)=6.25;
       i(n)=1.25;
       p(n)=v(n)*t(n)*i(n);
    else if t(n)>=6 & t(n)<8
           v(n)=-6.25;
           i(n)=10;
           p(n)=(i(n)-(1.25*t(n)))*v(n);
        else
            v(n)=0;
            i(n)=2.5;
            p(n)=v(n)*t(n)*i(n);
            end
    end 
    end
xlabel("Time in seconds")
ylabel("power in watts")
title("power waveform")
plot(t,p)
subplot(2,2,3)
t=[0:0.001:8];
x=length(t);
e=ones(1,x);
L=5;
for n=1:x;
    if t(n)<=2
       i(n)=1.25;
       e(n)=(1/2)*L*(t(n)*i(n))^2;
    else if t(n)>=6 & t(n)<8
           i(n)=10;
           e(n)=(1/2)*L*(i(n)-(1.25*t(n)))^2;
        else
            i(n)=2.5;
            e(n)=(1/2)*L*(i(n))^2;
            end
    end 
    end
xlabel("Time in seconds")
ylabel("Energy in joules")
title("Energy waveform")
plot(t,e)