initial ci of NLMS object without feedback subtraction

svn path=/trunk/externals/iem/iem_adaptfilt/; revision=14738

1 files changed, 214 insertions, 0 deletions

diff --git a/help/NLMSerr_in~-help.pd b/help/NLMSerr_in~-help.pd
new file mode 100644
index 0000000..328a888
--- /dev/null
+++ b/help/NLMSerr_in~-help.pd
@@ -0,0 +1,214 @@
+#N canvas 378 66 925 637 10;
+#N canvas 0 22 450 300 (subpatch) 0;
+#X array W 40 float 0;
+#X coords 0 1 39 -1 200 140 1;
+#X restore 545 340 graph;
+#X msg 25 245 update \$1;
+#X msg 102 255 beta \$1;
+#X obj 272 24 tgl 15 0 empty empty empty 0 -6 0 8 -225280 -1 -1 0 1
+;
+#X obj 272 45 dsp;
+#X floatatom 272 91 5 0 0 0 - - -;
+#X floatatom 285 70 5 0 0 0 - - -;
+#X obj 276 563 bng 15 150 20 0 empty empty empty 0 -6 0 8 -262144 -1
+-1;
+#X text 453 230 2.arg: <float> learn-rate = beta;
+#X obj 53 275 noise~;
+#X obj 25 58 vradio 15 1 0 8 empty empty empty 0 -6 0 8 -225280 -1
+-1 0;
+#N canvas 165 115 464 314 (subpatch) 0;
+#X obj 76 61 inlet;
+#X msg 32 163 0;
+#X msg 63 159 1;
+#X msg 97 158 2;
+#X msg 132 156 4;
+#X msg 159 157 8;
+#X msg 191 162 16;
+#X msg 219 164 32;
+#X msg 248 165 64;
+#X obj 76 84 sel 0 1 2 3 4 5 6 7;
+#X obj 32 217 outlet;
+#X connect 0 0 9 0;
+#X connect 1 0 10 0;
+#X connect 2 0 10 0;
+#X connect 3 0 10 0;
+#X connect 4 0 10 0;
+#X connect 5 0 10 0;
+#X connect 6 0 10 0;
+#X connect 7 0 10 0;
+#X connect 8 0 10 0;
+#X connect 9 0 1 0;
+#X connect 9 1 2 0;
+#X connect 9 2 3 0;
+#X connect 9 3 4 0;
+#X connect 9 4 5 0;
+#X connect 9 5 6 0;
+#X connect 9 6 7 0;
+#X connect 9 7 8 0;
+#X restore 25 205 pd;
+#X text 34 37 internal downsampling of update;
+#X msg 96 337 gamma \$1;
+#X text 75 363 input signal;
+#X text 190 362 desired signal;
+#N canvas 0 22 450 300 (subpatch) 0;
+#X array IR 40 float 0;
+#X coords 0 1 39 -1 200 140 1;
+#X restore 367 11 graph;
+#X obj 54 540 unsig~;
+#X floatatom 54 566 9 0 0 0 - - -;
+#X obj 174 539 unsig~;
+#X floatatom 174 565 9 0 0 0 - - -;
+#X text 26 363 x(n);
+#X text 243 377 d(n);
+#X text 29 508 y(n) = W * x(n);
+#X obj 276 291 FIR~ IR 32;
+#X obj 618 73 loadbang;
+#X text 503 220 (array-sizes have to be >= then FIR_size);
+#X text 453 210 1.arg: <float> number of order of FIR-filter;
+#X text 42 71 update every sample;
+#X text 42 56 stop \, no update;
+#X text 42 86 update every 2nd sample;
+#X text 42 101 update every 4th sample;
+#X text 42 116 update every 8th sample;
+#X text 42 131 update every 16th sample;
+#X text 42 146 update every 32nd sample;
+#X text 42 161 update every 64th sample;
+#N canvas 0 22 474 324 (subpatch) 0;
+#X obj 48 46 inlet;
+#X obj 205 47 inlet;
+#X msg 205 71 \; IR const 0;
+#X msg 48 120 \; IR 0 0 0 0 0.314287 0.8 0.75 0 0 0 0 -0.7 -0.65 0
+0 0 0.157143 0.128572 0 0 -0.128572 -0.1 0 0 0 0 0 0 0 0 0 0 0 0 0
+0 0 0 0 0 0;
+#X obj 92 84 loadbang;
+#X connect 0 0 3 0;
+#X connect 1 0 2 0;
+#X connect 4 0 3 0;
+#X restore 619 48 pd;
+#X obj 632 30 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1
+-1;
+#X text 649 27 clear;
+#X obj 619 11 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1
+-1;
+#X text 636 8 load;
+#X text 644 50 IR;
+#X text 4 495 filtered output signal;
+#X text 178 495 error signal;
+#X text 152 507 e(n) = d(n) - W * x(n);
+#X text 293 557 constrain;
+#X text 293 567 coefficients;
+#X msg 276 522 0;
+#X obj 276 542 speedlim 100;
+#X obj 83 512 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0;
+#X obj 247 511 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0;
+#N canvas 0 22 470 320 (subpatch) 0;
+#X obj 137 73 tgl 15 0 empty empty empty 0 -6 0 8 -262144 -1 -1 1 1
+;
+#X msg 137 115 39;
+#X msg 137 156 \$1 39;
+#X msg 136 53 1;
+#X obj 137 92 metro 200;
+#X obj 136 24 inlet;
+#X obj 137 135 tabread W;
+#X obj 137 180 tabwrite W;
+#X connect 0 0 4 0;
+#X connect 1 0 6 0;
+#X connect 2 0 7 0;
+#X connect 3 0 0 0;
+#X connect 4 0 1 0;
+#X connect 5 0 3 0;
+#X connect 6 0 2 0;
+#X restore 618 94 pd;
+#X text 643 100 update of W;
+#X text 642 89 graphical;
+#X text 453 250 4.arg: <symbol> table-name of W;
+#X text 453 240 3.arg: <float> minimum input value gamma;
+#X text 132 317 minimum input value;
+#X msg 102 218 0.1;
+#X msg 96 299 1e-05;
+#X text 155 256 beta [0 .. 2];
+#X text 152 336 gamma [0 .. 1];
+#X obj 455 187 cnv 15 68 17 empty empty empty 20 12 0 14 -225280 -66577
+0;
+#N canvas 0 22 499 295 FORMULAS 0;
+#X obj 167 52 cnv 15 150 40 empty empty empty 20 12 0 14 -225280 -66577
+0;
+#X obj 167 123 cnv 15 270 50 empty empty empty 20 12 0 14 -225280 -66577
+0;
+#X obj 167 205 cnv 15 260 30 empty empty empty 20 12 0 14 -225280 -66577
+0;
+#X text 280 129 beta;
+#X text 231 146 >;
+#X text 232 144 _;
+#X text 232 138 _;
+#X text 283 141 2;
+#X text 346 141 2;
+#X text 311 146 gamma * blocksize);
+#X text 300 146 +;
+#X text 243 146 x[n-i];
+#X text 174 135 my(n) =;
+#X text 223 130 _________________________________;
+#X text 16 135 normalized learn rate:;
+#X text 173 56 y(n) = W * x(n);
+#X obj 237 54 cnv 11 1 1 empty empty * 0 9 0 14 -225280 -1 0;
+#X text 173 72 e(n) = d(n) - W * x(n);
+#X obj 286 71 cnv 11 1 1 empty empty * 0 9 0 14 -225280 -1 0;
+#X text 119 73 error:;
+#X text 116 55 output:;
+#X text 26 210 coefficient iteration:;
+#X text 188 211 W(k+1 \, i) = W(k \, i) + my(n)* e(n)* x(n);
+#X restore 455 187 pd FORMULAS;
+#X obj 102 238 cnv 14 31 14 empty empty empty 20 12 0 14 -225280 -66577
+0;
+#X floatatom 102 238 5 0 2 0 - - -;
+#X text 138 237 learn-rate;
+#X obj 96 319 cnv 14 31 14 empty empty empty 20 12 0 14 -225280 -66577
+0;
+#X floatatom 96 319 5 0 2 0 - - -;
+#X floatatom 25 226 5 0 0 0 - - -;
+#X text 455 169 Normalized Least Mean Square (linear adaptive FIR-filter)
+;
+#X text 451 268 (C) 2005 \, m.noisternig & t.musil \, IEM \, Graz \,
+Austria;
+#X text 479 281 [noisternig \, musil]_AT_iem.at;
+#X obj 53 450 NLMSerr_in~ 32 0.1 1e-05 W;
+#X obj 20 390 delwrite~ causal_del 2;
+#X msg 756 387 64 44.1;
+#X obj 760 415 /;
+#X obj 760 442 print;
+#X obj 20 411 delread~ causal_del 1.45126;
+#X obj 226 475 delwrite~ feedback_del 2;
+#X obj 351 358 delread~ feedback_del 0;
+#X obj 248 415 -~;
+#X text 249 437 err;
+#X text 351 376 filt out;
+#X connect 1 0 72 0;
+#X connect 2 0 72 0;
+#X connect 3 0 4 0;
+#X connect 4 0 5 0;
+#X connect 4 1 6 0;
+#X connect 9 0 24 0;
+#X connect 9 0 73 0;
+#X connect 10 0 11 0;
+#X connect 11 0 68 0;
+#X connect 13 0 72 0;
+#X connect 17 0 18 0;
+#X connect 19 0 20 0;
+#X connect 24 0 80 0;
+#X connect 25 0 51 0;
+#X connect 37 0 36 1;
+#X connect 39 0 36 0;
+#X connect 47 0 48 0;
+#X connect 48 0 7 0;
+#X connect 57 0 64 0;
+#X connect 58 0 67 0;
+#X connect 64 0 2 0;
+#X connect 67 0 13 0;
+#X connect 68 0 1 0;
+#X connect 72 0 17 0;
+#X connect 72 0 78 0;
+#X connect 74 0 75 0;
+#X connect 75 0 76 0;
+#X connect 77 0 72 0;
+#X connect 79 0 80 1;
+#X connect 80 0 72 1;