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
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
#! /usr/bin/env perl
#
# Demo x04 for PLplot PDL bindings
#
# Log plot demo
# (inspired from t/x04.t of module Graphics::PLplot, by Tim Jenness)
#
# Copyright (C) 2004 Rafael Laboissiere
#
# This file is part of PLplot.
#
# PLplot is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Library Public License as published
# by the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# PLplot is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Library General Public License for more details.
#
# You should have received a copy of the GNU Library General Public License
# along with PLplot; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
use PDL;
use PDL::Graphics::PLplot;
use Math::Trig qw [pi];
# Parse and process command line arguments
plParseOpts (\@ARGV, PL_PARSE_SKIP | PL_PARSE_NOPROGRAM);
# Initialize plplot
plinit ();
plfont (2);
&plot1 (0);
&plot1 (1);
plend ();
sub plot1 {
my $type = shift;
pladv (0);
# Set up data for log plot
my $f0 = 1.0;
my ($freql, $ampl, $phase);
my $i = sequence (101);
$freql = -2.0 + $i / 20.0;
my $freq = 10 ** $freql;
$ampl = 20 * log10 (1 / sqrt (1 + ($freq / $f0) ** 2 ));
$phase = -(180/pi) * atan2 ($freq, $f0);
plvpor (0.15, 0.85, 0.1, 0.9);
plwind (-2, 3, -80, 0);
# Try different axis and labelling styles
plcol0 (1);
if ($type == 0) {
plbox (0, 0, 0, 0, "bclnst", "bnstv");
} else {
plbox (0, 0, 0, 0, "bcfghlnst", "bcghnstv");
}
# Plot ampl vs freq
plcol0 (2);
plline ($freql, $ampl);
plcol0 (1);
plptex (1.6, -30, 1, -20, 0.5, "-20 dB/decade");
# Put labels on
plcol0 (1);
plmtex (3.2, 0.5, 0.5, "b", "Frequency");
plmtex (2, 0.5, 0.5, "t", "Single Pole Low-Pass Filter");
plcol0 (2);
plmtex (5, 0.5, 0.5, "l", "Amplitude (dB)");
# For the gridless case, put phase vs freq on same plot
if ($type == 0) {
plcol0 (1);
plwind (-2, 3, -100, 0);
plbox (0, 0, 30, 3, "", "cmstv");
plcol0 (3);
plline ($freql, $phase);
plcol0 (3);
plmtex (5, 0.5, 0.5, "r", "Phase shift (degrees)");
}
}
|
