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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
|
function [increment, pm_candidates_imag, pm_candidates_real]= equations_gen(Lh)
###########################################################################
# Copyright (C) 2008 by Piotr Krysik #
# pkrysik@stud.elka.pw.edu.pl #
# #
# This program 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 3 of the License, or #
# (at your option) any later version. #
# #
# This program 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 this program; if not, write to the #
# Free Software Foundation, Inc., #
# 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. #
###########################################################################
INCREMENT_NUM=2**(Lh-1);
SYMBOLS = make_symbols(Lh);
PREVIOUS = generate_previous(Lh);
load tests/data/rhh.dat
[INCREMENT NUMBERS] = generate_increment(SYMBOLS,PREVIOUS);
function_real="real";
function_imag="imag";
sign="";
fction="";
equation_num="";
increment = { };
equation="";
for equation_num=1:INCREMENT_NUM,
for part_num=1:Lh,
coefficient=INCREMENT(equation_num,part_num);
if(coefficient==1 || coefficient==-1),
fction=function_real;
else
fction=function_imag;
coefficient=coefficient*j;
end
if(coefficient==1 && part_num != 1),
sign="+";
elseif(coefficient==-1),
sign="-";
else
sign="";
end
equation = [equation " " sign "rhh[" int2str(part_num) "]." fction "()" ];
end
increment(equation_num) = ["increment[" int2str(equation_num-1) "] =" equation];
equation="";
end
BRANCH_NUM=2**(Lh+1);
%make path metrics
pm_candidates_real={};
pm_candidates_imag={};
for symbol_num=1:BRANCH_NUM,
inc1_num=NUMBERS(PREVIOUS(symbol_num,1),symbol_num);
inc2_num=NUMBERS(PREVIOUS(symbol_num,2),symbol_num);
symbol=SYMBOLS(symbol_num,1);
if(inc1_num<0)
inc1_sign="+";
inc1_num=-inc1_num;
else
inc1_sign="-";
end
if(inc2_num<0)
inc2_sign="+";
inc2_num=-inc2_num;
else
inc2_sign="-";
end
inc1_num=inc1_num-1;
inc2_num=inc2_num-1;
if(symbol==1 || symbol==-1),
fction=function_real;
else
fction=function_imag;
symbol=symbol*(-j);
end
if(symbol==1),
symbol_sign="+";
else
symbol_sign="-";
end
branch_metric1=[" " symbol_sign " input_symbol_" fction " " inc1_sign " increment[" int2str(inc1_num) "]" ];
branch_metric2=[" " symbol_sign " input_symbol_" fction " " inc2_sign " increment[" int2str(inc2_num) "]" ];
num = ceil((symbol_num)/2);
prev1_num = ceil((PREVIOUS(symbol_num,1))/2)-1;
prev2_num = ceil((PREVIOUS(symbol_num,2))/2)-1;
if(mod(symbol_num,2)==1),
pm_candidates_imag(num, 1)=[ "pm_candidate1 = old_path_metrics[" int2str(prev1_num) "]" branch_metric1 ];
pm_candidates_imag(num, 2)=[ "pm_candidate2 = old_path_metrics[" int2str(prev2_num) "]" branch_metric2 ];
else
pm_candidates_real(num, 1)=[ "pm_candidate1 = old_path_metrics[" int2str(prev1_num) "]" branch_metric1 ];
pm_candidates_real(num, 2)=[ "pm_candidate2 = old_path_metrics[" int2str(prev2_num) "]" branch_metric2 ];
end
end
|
