0001 function Ar = M2Ar(M) 0002 nArtic = DIVA('VocalTract','nArtic'); 0003 nMotor = DIVA('MotorCortex','nMotor'); 0004 0005 W = zeros(nArtic,nMotor); 0006 pb = diag(DIVA('VocalTract','MaedaSpread')*ones(nArtic,1)); 0007 0008 W(:,1:2:end) = pb; 0009 W(:,2:2:end) = -pb; 0010 0011 Ar = W*M; 0012 Ar(8,:) = Ar(8,:)+DIVA('VocalTract','JawLip')*Ar(1,:); 0013 Ar = max(min(Ar,DIVA('VocalTract','MaedaSpread')),-DIVA('VocalTract','MaedaSpread'));