Вы берете скалярное произведение каждого ряда phis с w. Вы можете сделать это легко, используя линейную алгебру.
out = phis * w';
Это матричное умножение сохраняет ваши звонки на sum, ones и size, что должно сделать ваш код намного быстрее Кроме того, операции линейной алгебры часто бывают очень быстрыми в Matlab, поскольку именно для этого программа исторически оптимизирована.