振动函数的波形读取误差matlab

Question

我正试图编译我在互联网上发现的颤音效果代码。代码有一个函数调用wavread，在这个函数中matlab显示错误，我搜索了另一个代码来做到这一点，每个代码都使用这个函数来打开de wav文件，有人发生了什么?以下代码:

颤音脚本:

clear all;
close all;
clc;
infile = 'musica.wav';
% read in wav sample
[ x, Fs, N ] = wavread(infile);

%set Parameters for vibrato
% Change these to experiment with vibrato
Modfreq = 10; %10 Khz
Width = 0.0008; % 0.8 Milliseconds
% Do vibrato
yvib = vibrato(x, Fs, Modfreq, Width);
% write output wav files
wavwrite(yvib, Fs,  'out_vibrato.wav');
% plot the original and equalised waveforms
figure(1)
hold on
plot(x(1:500),'r');
plot(yvib(1:500),'b');
title('Vibrato First 500 Samples');

颤音功能:

% Vibrato 
function y=vibrato(x,SAMPLERATE,Modfreq,Width)
ya_alt=0;
Delay=Width; % basic delay of input sample in sec
DELAY=round(Delay*SAMPLERATE); % basic delay in # samples
WIDTH=round(Width*SAMPLERATE); % modulation width in # samples
if WIDTH>DELAY 
  error('delay greater than basic delay !!!');
  return;
end
MODFREQ=Modfreq/SAMPLERATE; % modulation frequency in # samples
LEN=length(x);        % # of samples in WAV-file
L=2+DELAY+WIDTH*2;    % length of the entire delay  
Delayline=zeros(L,1); % memory allocation for delay
y=zeros(size(x));     % memory allocation for output vector
for n=1:(LEN-1)
   M=MODFREQ;
   MOD=sin(M*2*pi*n);
   ZEIGER=1+DELAY+WIDTH*MOD;
   i=floor(ZEIGER);
   frac=ZEIGER-i;
   Delayline=[x(n);Delayline(1:L-1)]; 
   %---Linear Interpolation-----------------------------
   y(n,1)=Delayline(i+1)*frac+Delayline(i)*(1-frac); 
   %---Allpass Interpolation------------------------------
   %y(n,1)=(Delayline(i+1)+(1-frac)*Delayline(i)-(1-frac)*ya_alt);  
   %ya_alt=ya(n,1);
end 

出现的错误:

Undefined function or variable 'wavread' .

在这一行:

[ x, Fs, N ] = wavread(infile);

Answer 1

从Matlab R2015b开始不支持wavread函数。

此功能已被audioread取代，原型略有变化。

请用

替换有故障的行

% read in wav sample
[x, Fs] = audioread(infile);

那么wavwrite被audiowrite取代的情况也是一样的。

应该改成

% write output wav files
audiowrite('out_vibrato.wav', yvib, Fs);