dan@0: function wmp = win_mod_poly(N, t, poly, w)
dan@0: R = size(poly,1);
dan@0: ply_val = zeros(N,R);
dan@0: w_mat = repmat (w, 1, R);
dan@0: for k = 1:R
dan@0:   ply_val(:,k) = polyval(poly(k,:),t).';
dan@0: end
dan@0: wmp = sum (w_mat.*exp(ply_val),1);
dan@0: end