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
|
clc;
//page no 43
//problem 2.5
//All voltages are in V
//All power in watts
R=5; //ohm
C0=5; //dc value
C1=8;
C2=6;
C3=3; //volts
Vrms=sqrt(C0^2+0.5*(C1^2+C2^2+C3^2)); //rms voltage
disp(Vrms,'(a) The rms value of voltage is');
P=Vrms^2/R; //watts
disp('W',P,'(b) The average power dissipated in resistor is')
P0=C0^2/R;
disp(P0,'(c) The dc power is ');
P1=C1^2/(2*R);
disp(P1,'The power in fundamental is');
P2=C2^2/(2*R);
P3=C3^2/(2*R);
disp(P3,P2,'The second and third harmonics are');
//for plot
fHz=[0:1600];
f1=500; //fundamental freq.
f2=1000; f3=1500;
Pn=[P0 zeros(1:f1-1) P1 zeros(f1+1:f2-1) P2 zeros(f2+1:f3-1) P3 zeros(f3+1:1600)]
clf
plot2d(fHz,Pn,[3],rect=[0,0,1600,8])
xtitle('Power spectrum','f,Hz','Pn(W)')
xgrid
|
