matlab仿真程序源码 音频加解扰

clearclcclose all % 清屏
N=176400; % 采样点数 一共4秒，每秒44100个点，采集176400个点
n=[0:N-1];
[x1fs1nbits]=wavread（‘D:\bin\au.wav‘）;% 采语音信号
%sound（x1fs1nbits）; % 播放歌曲  用来试听
x=x1（1:N）;
Ts=1/fs1; % 抽样时间间隔
X=fft（x）; % 傅里叶变换
f=n/N*fs1; % 把点数转换成频率
b=abs（X）;
subplot（211）;
plot（fb）;
title（‘ 输入语音的幅度谱 ‘）;
xlabel（‘ 频率为 /Hz‘）;
%翻转频谱
f1=22050;   %中点就是44100的一半
y=x.*exp（1i*2*pi*f1.*n*Ts）; % 进行翻转
Y=fft（y）;
subplot（212）;
plot（fabs（Y））; % 画出翻转后的幅度图
title（‘ 频谱倒置语音的幅度谱 ‘）;
xlabel（‘