summaryrefslogtreecommitdiff
path: root/testsuite/vests/vhdl-ams/ashenden/compliant/util/src_sine.vhd
blob: f5b82c839ca467b7bf1bb921d8b9324b8b44cdf4 (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

-- Copyright (C) 2002 Morgan Kaufmann Publishers, Inc

-- This file is part of VESTs (Vhdl tESTs).

-- VESTs is free software; you can redistribute it and/or modify it
-- under the terms of the GNU General Public License as published by the
-- Free Software Foundation; either version 2 of the License, or (at
-- your option) any later version. 

-- VESTs 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 General Public License
-- for more details. 

-- You should have received a copy of the GNU General Public License
-- along with VESTs; if not, write to the Free Software Foundation,
-- Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA 

library IEEE;  use IEEE.MATH_REAL.all;
library IEEE_proposed;  use IEEE_proposed.ELECTRICAL_SYSTEMS.all;

entity src_sine is

  generic ( freq : real;              -- frequency [Hertz]
            amplitude : voltage;      -- amplitude [Volts]
            phase : real := 0.0;      -- initial phase [Degrees]
            offset : voltage := 0.0;  -- DC value [Volts]
            df : real := 0.0;         -- damping factor [1/second]
            ac_mag : voltage := 1.0;  -- AC magnitude [Volts]
            ac_phase : real := 0.0);  -- AC phase [Degrees]

  port ( quantity output : out real );

end entity src_sine;


architecture ideal of src_sine is
  
  -- Declare quantity for phase in radians (calculated below)
  quantity phase_rad : real;
  -- Declare quantity in frequency domain for AC analysis
  quantity ac_spec : real spectrum ac_mag, math_2_pi * ac_phase / 360.0;

begin
  
  -- Convert phase to radians
  phase_rad == math_2_pi *(freq * now + phase / 360.0);

  if domain = quiescent_domain or domain = time_domain use
    output == offset + amplitude * sin(phase_rad) * exp(-now * df);
  else
    output == ac_spec;  -- used for Frequency (AC) analysis
  end use;

end architecture ideal;