#N canvas 170 22 1021 882 10; #N canvas 0 0 450 300 graph1 0; #X array 1_W 40 float 0; #X coords 0 1 39 -1 200 140 1; #X restore 64 519 graph; #X msg 40 247 update \$1; #X msg 117 257 beta \$1; #X obj 456 26 tgl 15 0 empty empty empty 0 -6 0 8 -225280 -1 -1 0 1 ; #X obj 456 47 dsp; #X floatatom 456 93 5 0 0 0 - - -; #X floatatom 469 72 5 0 0 0 - - -; #X obj 452 463 bng 15 150 20 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X obj 68 277 noise~; #X obj 40 60 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 40 207 pd; #X text 22 41 internal downsampling of update; #X msg 75 339 gamma \$1; #X text 90 365 input signal; #X text 335 364 desired signal; #N canvas 0 0 450 300 graph1 0; #X array IR 40 float 0; #X coords 0 1 39 -1 200 140 1; #X restore 240 25 graph; #X obj 69 462 unsig~; #X floatatom 69 488 9 0 0 0 - - -; #X obj 155 462 unsig~; #X floatatom 155 488 9 0 0 0 - - -; #X text 426 364 d(n); #X obj 516 270 loadbang; #X text 704 47 (array-sizes have to be >= then FIR_size); #X text 654 37 1.arg: number of order of FIR-filter; #X text 57 73 update every sample; #X text 57 58 stop \, no update; #X text 57 88 update every 2nd sample; #X text 57 103 update every 4th sample; #X text 57 118 update every 8th sample; #X text 57 133 update every 16th sample; #X text 57 148 update every 32nd sample; #X text 57 163 update every 64th sample; #N canvas 0 0 486 336 /SUBPATCH/ 0; #X obj 48 46 inlet; #X obj 205 47 inlet; #X msg 205 71 \; IR const 0; #X obj 92 84 loadbang; #X msg 48 120 \; IR 0 0 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; #X connect 0 0 4 0; #X connect 1 0 2 0; #X connect 3 0 4 0; #X restore 517 245 pd; #X obj 530 227 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X text 547 224 clear; #X obj 517 208 bng 15 250 50 0 empty empty empty 0 -6 0 8 -262144 -1 -1; #X text 534 205 load; #X text 579 214 IR; #X text 19 397 filtered output signal; #X text 339 396 error signal; #X text 469 457 constrain; #X text 469 467 coefficients; #X msg 452 422 0; #X obj 452 442 speedlim 100; #X obj 110 414 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0; #N canvas 0 0 482 332 /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 92 190 \$1 39; #X msg 136 53 1; #X obj 137 92 metro 200; #X obj 136 24 inlet; #X obj 92 169 tabread 1_W; #X obj 92 214 tabwrite 1_W; #X msg 177 192 \$1 39; #X msg 260 191 \$1 39; #X obj 260 170 tabread 3_W; #X obj 260 215 tabwrite 3_W; #X obj 177 171 tabread 2_W; #X obj 177 216 tabwrite 2_W; #X connect 0 0 4 0; #X connect 1 0 6 0; #X connect 1 0 12 0; #X connect 1 0 10 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 connect 8 0 13 0; #X connect 9 0 11 0; #X connect 10 0 9 0; #X connect 12 0 8 0; #X restore 516 291 pd; #X text 541 297 update of W; #X text 540 286 graphical; #X text 111 319 minimum input value; #X msg 117 220 0.1; #X msg 75 301 1e-05; #X text 170 258 beta [0 .. 2]; #X text 131 338 gamma [0 .. 1]; #X obj 574 38 cnv 15 68 17 empty empty empty 20 12 0 14 -225280 -66577 0; #N canvas 0 0 568 372 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 277 141 2; #X text 340 141 2; #X text 311 146 gamma * blocksize); #X text 294 146 +; #X text 243 146 x[n-i]; #X text 174 135 my(n) =; #X text 223 130 _________________________________; #X text 6 135 normalized learn rate:; #X text 173 54 y(n) = W * x(n); #X obj 225 56 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 267 74 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 7 212 coefficient iteration:; #X text 176 213 W(k+1 \, i) = W(k \, i) + my(n)* e(n)* x(n); #X restore 573 39 pd FORMULAS; #X obj 117 240 cnv 14 31 14 empty empty empty 20 12 0 14 -225280 -66577 0; #X floatatom 117 240 5 0 2 0 - - -; #X text 153 239 learn-rate; #X obj 75 321 cnv 14 31 14 empty empty empty 20 12 0 14 -225280 -66577 0; #X floatatom 75 321 5 0 1 0 - - -; #X floatatom 40 228 5 0 0 0 - - -; #X text 572 14 Normalized Least Mean Square (linear adaptive FIR-filter) ; #N canvas 0 0 450 300 graph1 0; #X array 2_W 40 float 0; #X coords 0 1 39 -1 200 140 1; #X restore 314 519 graph; #N canvas 0 0 450 300 graph1 0; #X array 3_W 40 float 0; #X coords 0 1 39 -1 200 140 1; #X restore 564 519 graph; #X obj 68 382 n_CNLMS~ 3 32 0.1 1e-05 1 W _______________________________ ; #X obj 244 463 unsig~; #X floatatom 244 489 9 0 0 0 - - -; #X obj 334 461 unsig~; #X floatatom 334 487 9 0 0 0 - - -; #X text 304 233 minimum input value; #X obj 268 235 cnv 14 31 14 empty empty empty 20 12 0 14 -225280 -66577 0; #X floatatom 268 235 5 0 10000 0 - - -; #X msg 268 253 kappa \$1; #X text 324 252 kappa [0 .. 1]; #X msg 268 215 1; #X obj 423 293 FIR~ IR 32; #X obj 186 293 z~ 2; #X obj 304 293 z~ 4; #X text 41 365 x1(n); #X text 272 365 x3(n); #X text 191 364 x2(n); #X text 39 410 y1(n) = 1_W * x1(n); #X obj 196 429 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0; #X obj 286 403 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0; #X text 125 425 y2(n) = 2_W * x2(n); #X text 215 399 y3(n) = 3_W * x3(n); #X obj 375 425 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0; #X text 341 421 - 1_W * x1(n); #X text 312 410 e(n) = d(n); #X obj 375 435 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0; #X text 341 431 - 2_W * x2(n); #X obj 375 445 cnv 8 1 1 empty empty * 0 7 0 14 -262144 -1 0; #X text 341 441 - 3_W * x3(n); #X text 572 3 Multi Channel Constraint; #X msg 463 139 \; 1_W const 0 \; 2_W const 0 \; 3_W const 0; #X text 462 117 RESET TABLES; #X text 654 150 (C) 2005 \, m.noisternig & t.musil \, IEM \, Graz \, Austria; #X text 682 163 [noisternig \, musil]_AT_iem.at; #X text 655 64 2.arg: number of order of FIR-filter; #X text 655 78 3.arg: learn-rate = beta; #X text 656 103 5.arg: constraint = kappa; #X text 655 90 4.arg: regularization parameter = gamma; #X text 656 118 6.arg: table-name of W; #X connect 1 0 64 0; #X connect 2 0 64 0; #X connect 3 0 4 0; #X connect 4 0 5 0; #X connect 4 1 6 0; #X connect 8 0 75 0; #X connect 8 0 64 0; #X connect 8 0 76 0; #X connect 8 0 77 0; #X connect 9 0 10 0; #X connect 10 0 60 0; #X connect 12 0 64 0; #X connect 16 0 17 0; #X connect 18 0 19 0; #X connect 21 0 45 0; #X connect 33 0 32 1; #X connect 35 0 32 0; #X connect 42 0 43 0; #X connect 43 0 7 0; #X connect 49 0 56 0; #X connect 50 0 59 0; #X connect 56 0 2 0; #X connect 59 0 12 0; #X connect 60 0 1 0; #X connect 64 0 16 0; #X connect 64 1 18 0; #X connect 64 2 65 0; #X connect 64 3 67 0; #X connect 64 4 42 0; #X connect 65 0 66 0; #X connect 67 0 68 0; #X connect 71 0 72 0; #X connect 72 0 64 0; #X connect 74 0 71 0; #X connect 75 0 64 3; #X connect 76 0 64 1; #X connect 77 0 64 2;