基于频谱搬移的变调算法

[z fs bits] =wavread（‘F:\apple.wav‘）;
x=z（1:131072）;
x=x‘;
N=131072;%待处理序列长度
X=fft（x）;
% X=X‘;
f=0:N-1;


a=1000;%范围+5到-6之间，取整数，正为升，负为降

qian1=X（1:N/2）;
qian2=qian1‘;
qian=[zeros（1a）qian2（1:N/2-a）];


hou1=X（N/2+1:N）;
hou2=hou1‘;
hou=[hou2（a+1:N/2）zeros（1a）];
Y=[qianhou];

Y=Y‘;
y=real（ifft（Y））;

% 
% subplot（121）;
% plot（abs（X））;%画频谱图
% title（‘原始信号频谱‘）;
% subplot（122）;
% plot（abs（Y））;%画频谱图
%