My Project
0.0.16
QUCS Mapping
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
DownLoad
QUCS-src
qucs-0.0.16
qucs-core
src
components
vpulse.cpp
Go to the documentation of this file.
1
/*
2
* vpulse.cpp - pulse voltage source class implementation
3
*
4
* Copyright (C) 2004, 2005, 2006, 2008 Stefan Jahn <stefan@lkcc.org>
5
*
6
* This is free software; you can redistribute it and/or modify
7
* it under the terms of the GNU General Public License as published by
8
* the Free Software Foundation; either version 2, or (at your option)
9
* any later version.
10
*
11
* This software is distributed in the hope that it will be useful,
12
* but WITHOUT ANY WARRANTY; without even the implied warranty of
13
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14
* GNU General Public License for more details.
15
*
16
* You should have received a copy of the GNU General Public License
17
* along with this package; see the file COPYING. If not, write to
18
* the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
19
* Boston, MA 02110-1301, USA.
20
*
21
* $Id: vpulse.cpp 1825 2011-03-11 20:42:14Z ela $
22
*
23
*/
24
25
#if HAVE_CONFIG_H
26
# include <config.h>
27
#endif
28
29
#include "
component.h
"
30
#include "
vpulse.h
"
31
32
vpulse::vpulse () :
circuit
(2) {
33
type =
CIR_VPULSE
;
34
setVSource (
true
);
35
setVoltageSources (1);
36
}
37
38
void
vpulse::initSP
(
void
) {
39
allocMatrixS
();
40
setS
(
NODE_1
,
NODE_1
, 0.0);
41
setS
(
NODE_1
,
NODE_2
, 1.0);
42
setS
(
NODE_2
,
NODE_1
, 1.0);
43
setS
(
NODE_2
,
NODE_2
, 0.0);
44
}
45
46
void
vpulse::initDC
(
void
) {
47
allocMatrixMNA
();
48
voltageSource
(
VSRC_1
,
NODE_1
,
NODE_2
);
49
setE
(
VSRC_1
,
getPropertyDouble
(
"U1"
));
50
}
51
52
void
vpulse::initAC
(
void
) {
53
initDC
();
54
setE
(
VSRC_1
, 0);
55
}
56
57
void
vpulse::initTR
(
void
) {
58
initDC
();
59
}
60
61
void
vpulse::calcTR
(nr_double_t
t
) {
62
nr_double_t u1 =
getPropertyDouble
(
"U1"
);
63
nr_double_t u2 =
getPropertyDouble
(
"U2"
);
64
nr_double_t t1 =
getPropertyDouble
(
"T1"
);
65
nr_double_t t2 =
getPropertyDouble
(
"T2"
);
66
nr_double_t tr =
getPropertyDouble
(
"Tr"
);
67
nr_double_t tf =
getPropertyDouble
(
"Tf"
);
68
nr_double_t ut = 0;
69
nr_double_t
s
=
getNet
()->
getSrcFactor
();
70
71
if
(t < t1) {
// before pulse
72
ut = u1;
73
}
74
else
if
(t < t1 + tr) {
// rising edge
75
ut = u1 + (u2 - u1) / tr * (t - t1);
76
}
77
else
if
(t < t2 - tf) {
// during full pulse
78
ut = u2;
79
}
80
else
if
(t < t2) {
// falling edge
81
ut = u2 + (u1 - u2) / tf * (t - (t2 - tf));
82
}
83
else
{
// after pulse
84
ut = u1;
85
}
86
setE
(
VSRC_1
, ut * s);
87
}
88
89
// properties
90
PROP_REQ
[] = {
91
{
"U1"
,
PROP_REAL
, { 0,
PROP_NO_STR
},
PROP_NO_RANGE
},
92
{
"U2"
,
PROP_REAL
, { 1,
PROP_NO_STR
},
PROP_NO_RANGE
},
93
{
"T1"
,
PROP_REAL
, { 0,
PROP_NO_STR
},
PROP_POS_RANGE
},
94
{
"T2"
,
PROP_REAL
, { 1
e
-3,
PROP_NO_STR
},
PROP_POS_RANGE
},
95
PROP_NO_PROP
};
96
PROP_OPT
[] = {
97
{
"Tr"
,
PROP_REAL
, { 1
e
-9,
PROP_NO_STR
},
PROP_POS_RANGE
},
98
{
"Tf"
,
PROP_REAL
, { 1
e
-9,
PROP_NO_STR
},
PROP_POS_RANGE
},
99
PROP_NO_PROP
};
100
struct
define_t
vpulse
::cirdef =
101
{
"Vpulse"
, 2,
PROP_COMPONENT
,
PROP_NO_SUBSTRATE
,
PROP_LINEAR
,
PROP_DEF
};
Generated on Tue Dec 25 2012 14:30:42 for My Project by
1.8.2