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