for sig=1:length(snr)
 for iK=1:K
 [x,fs]=audioread('fileaudio.m4a');
 .......

err=(x-x_est).^2;
  MSE(iK,sig)=sum(err(:))/numel(x);
  end
end
MSE=mean(MSE);