nikcleju@14
|
1 # -*- coding: utf-8 -*-
|
nikcleju@14
|
2 """
|
nikcleju@17
|
3 Defines standard parameters for exact reconstruction simulation
|
nikcleju@17
|
4 Author: Nicolae Cleju
|
nikcleju@17
|
5 """
|
nikcleju@17
|
6 __author__ = "Nicolae Cleju"
|
nikcleju@17
|
7 __license__ = "GPL"
|
nikcleju@17
|
8 __email__ = "nikcleju@gmail.com"
|
nikcleju@14
|
9
|
nikcleju@14
|
10
|
nikcleju@14
|
11 import numpy
|
nikcleju@17
|
12
|
nikcleju@17
|
13 # Solver algorithms to run
|
nikcleju@14
|
14 from algos import *
|
nikcleju@14
|
15
|
nikcleju@17
|
16
|
nikcleju@17
|
17 # Test parameters
|
nikcleju@15
|
18 paramstest = dict()
|
nikcleju@15
|
19 paramstest['algos'] = exact_gap,exact_sl0,exact_bp,exact_ompeps,exact_tst # tuple of algorithms
|
nikcleju@15
|
20 #paramstest['algos'] = exact_bp_cvxopt, # tuple of algorithms
|
nikcleju@18
|
21 paramstest['d'] = 200.0
|
nikcleju@15
|
22 paramstest['sigma'] = 1.2
|
nikcleju@15
|
23 paramstest['deltas'] = numpy.array([0.05, 0.45, 0.95])
|
nikcleju@15
|
24 paramstest['rhos'] = numpy.array([0.05, 0.45, 0.95])
|
nikcleju@15
|
25 #deltas = numpy.array([0.6])
|
nikcleju@15
|
26 #deltas = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
27 #rhos = numpy.array([0.05])
|
nikcleju@15
|
28 paramstest['numvects'] = 10; # Number of vectors to generate
|
nikcleju@15
|
29 paramstest['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
|
nikcleju@15
|
30 paramstest['savedataname'] = 'exact_pt_stdtest.mat'
|
nikcleju@15
|
31 paramstest['saveplotbase'] = 'exact_pt_stdtest_'
|
nikcleju@15
|
32 paramstest['saveplotexts'] = ('png','pdf','eps')
|
nikcleju@14
|
33
|
nikcleju@14
|
34 # Standard parameters 1
|
nikcleju@14
|
35 # All algorithms, 100 vectors
|
nikcleju@17
|
36 # d = 200, sigma = 1.2, delta and rho full resolution (0.05 step)
|
nikcleju@17
|
37 # Virtually no noise (100db)
|
nikcleju@15
|
38 params1 = dict()
|
nikcleju@15
|
39 params1['algos'] = exact_gap,exact_sl0,exact_bp_cvxopt,exact_ompeps,exact_tst # tuple of algorithms
|
nikcleju@17
|
40 params1['d'] = 200.0;
|
nikcleju@15
|
41 params1['sigma'] = 1.2
|
nikcleju@15
|
42 params1['deltas'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
43 params1['rhos'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
44 params1['numvects'] = 100; # Number of vectors to generate
|
nikcleju@15
|
45 params1['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
|
nikcleju@15
|
46 params1['savedataname'] = 'exact_pt_std1.mat'
|
nikcleju@15
|
47 params1['saveplotbase'] = 'exact_pt_std1_'
|
nikcleju@15
|
48 params1['saveplotexts'] = ('png','pdf','eps')
|
nikcleju@14
|
49
|
nikcleju@15
|
50
|
nikcleju@14
|
51 # Standard parameters 2
|
nikcleju@14
|
52 # All algorithms, 100 vectors
|
nikcleju@17
|
53 # d = 20, sigma = 10, delta and rho full resolution (0.05 step)
|
nikcleju@17
|
54 # Virtually no noise (100db)
|
nikcleju@15
|
55 params2 = dict()
|
nikcleju@15
|
56 params2['algos'] = exact_gap,exact_sl0,exact_bp_cvxopt,exact_ompeps,exact_tst # tuple of algorithms
|
nikcleju@15
|
57 params2['d'] = 20.0
|
nikcleju@15
|
58 params2['sigma'] = 10.0
|
nikcleju@15
|
59 params2['deltas'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
60 params2['rhos'] = numpy.arange(0.05,1.,0.05)
|
nikcleju@15
|
61 params2['numvects'] = 100; # Number of vectors to generate
|
nikcleju@15
|
62 params2['SNRdb'] = 100.; # This is norm(signal)/norm(noise), so power, not energy
|
nikcleju@15
|
63 params2['savedataname'] = 'exact_pt_std2.mat'
|
nikcleju@15
|
64 params2['saveplotbase'] = 'exact_pt_std2_'
|
nikcleju@15
|
65 params2['saveplotexts'] = ('png','pdf','eps') |