From 7ebd609b4f4139f4dac8627ac00678de89ef4575 Mon Sep 17 00:00:00 2001
From: Andreas Bogk <andreas@pt141.(none)>
Date: Thu, 15 Jan 2009 17:17:29 +0100
Subject: Viterbi generator by Piotr Krysik.

---
 .../utils/lower_utils/equations_gen.m              | 124 +++++++++++++++++++++
 1 file changed, 124 insertions(+)
 create mode 100644 viterbi_generator/utils/lower_utils/equations_gen.m

(limited to 'viterbi_generator/utils/lower_utils/equations_gen.m')

diff --git a/viterbi_generator/utils/lower_utils/equations_gen.m b/viterbi_generator/utils/lower_utils/equations_gen.m
new file mode 100644
index 0000000..af7d37e
--- /dev/null
+++ b/viterbi_generator/utils/lower_utils/equations_gen.m
@@ -0,0 +1,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
+
-- 
cgit v1.2.3