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
|
clc;
clear all;
//page no 54
//problem 2.9
f0=0;
f1=500; //fundamental freq.
f2=1000; f3=1500; //harmonics
//Values from ex 2.4
C=[5 8 6 3]// Values in Volts
//Values from ex 2.5
P=[5 6.4 3.6 .9]; //poweer in watts
clf;
// plot two sided linear amplitude spectrum
fHz=-1510:10^-2:1510; //x-axis matrix
//Y-axis matrix
Cn=[C(1)]
for i=2:4
Cn=[zeros(-500+10^-2:10^-2:0-10^-2) Cn zeros(0+10^-2:10^-2:500-10^-2)]
Cn=[C(i)/2 Cn C(i)/2];
end
Cn=[zeros(-10+10^-2:10^-2:0) Cn zeros(0:10^-2:10-10^-2)]
subplot(211)
plot2d(fHz,Cn,[2],rect=[-2000,0,2000,6])
xtitle('Two-sided Linear amplitude spectrum','f,Hz','Vn(V)')
xgrid
// plot two power spectrum
fHz=-1510:10^-2:1510; //x-axis matrix
//Y-axis matrix
Pn=[P(1)]
for i=2:4
Pn=[zeros(-500+10^-2:10^-2:0-10^-2) Pn zeros(0+10^-2:10^-2:500-10^-2)]
Pn=[P(i)/2 Pn P(i)/2];
end
Pn=[zeros(-10+10^-2:10^-2:0) Pn zeros(0:10^-2:10-10^-2)]
subplot(212)
plot2d(fHz,Pn,[6],rect=[-2000,0,2000,6])
xtitle('Two-sided power] spectrum','f,Hz','Pn(W)')
xgrid
|
