基于matlab的数字变声器和滤波器

学习的时候做的代码,可实现男女童声变换,还有一些滤波器,程序简单已调试无误

clear all;
[xfs] = audioread（‘66.wav‘）;
B1=resample（x11）;  
%sound（B1fs）;
X1=fft（B1）;
aX1=abs（X1）;
figure
subplot（321）
plot（B1）
title（‘时域波形图fs=48000‘）;
xlabel（‘时间‘）
yo1=（0:length（X1）-1）‘*fs/length（X1）;
subplot（322）
plot（yo1（1:length（yo1）/2）aX1（1:length（yo1）/2））;
title（‘原始信号频谱fs=48000‘）;
axis（[0 1500 0 500]）;
xlabel（‘频率/HZ‘）
ylabel（‘幅度‘）
B2=resample（x14）; 
%sound（B2fs）;
X2=fft（B2）;
aX2=abs（X2）;
subplot（323）
plot（B2）
title（‘时域波形图fs=12000‘）;
ylabel（‘幅度‘）
yo2=（0:length（X2）-1）‘*12000/length（X2）;
subplot（324）
plot（yo2（1:length（yo2）/2）aX2（1:length（yo2）/2））;
title（‘原始信号频谱fs=24000‘）;
axis（[0 1500 0 300]）;
xlabel（‘频率/HZ‘）
ylabel（‘幅度‘）
B3=resample（x21）; 
%sound（B3fs）;
X3=fft（B3）;
aX3=abs（X3）;
subplot（325）
plot（B3）
title（‘时域波形图fs=96000‘）;
ylabel（‘幅度‘）
yo3=（0:length（X3）-1）‘*96000/length（X3）;
subplot（326）
plot（yo3（1:length（yo3）/2）aX3（1:length（yo3）/2））;
title（‘原始信号频谱fs=96000‘）;
axis（[0 1500 0 1000]）;
xlabel（‘频率/HZ‘）
ylabel（‘幅度‘）
%m2=mean（B）;
%v2=var（B）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1077  2018-05-18 21:01  三级项目\a.m

     文件        697  2018-05-18 21:44  三级项目\b.m

     文件        986  2018-05-18 22:09  三级项目\daitong.m

     文件        563  2010-10-23 20:04  三级项目\findpitch.m

     文件        821  2018-05-14 17:22  三级项目\FIR_H.m

     文件        956  2018-05-14 17:26  三级项目\IIR_L.m

     文件       5185  2018-05-13 09:21  三级项目\nantonv.m

     文件       3255  2018-05-14 13:04  三级项目\nantotong.m

     文件       1067  2018-05-14 12:55  三级项目\nvtonan.m

     文件        359  2018-05-12 09:52  三级项目\pinpuyasuo.m

     文件      16044  2009-05-04 10:18  三级项目\sunday_2.wav

     目录          0  2018-07-06 19:29  三级项目

----------- ---------  ---------- -----  ----

                31010                    12