yading@12: function step3_heavy(arg)
yading@12: 
yading@12: x=wavread(arg);          
yading@12: [y,yh,ys] = hpsmodel(x,44100,hamming(2025),4096,-100,200,0,400,5,.2,4)       
yading@12: wavwrite(ys,44100,'ysbaby_voice_heavy.wav');