MATLAB的语音信号频谱分析

0123456-1-0.500.51幅 值s时 间 （ ）信 号 波 形00.20.40.60.811.21.41.61.820200040006000信 号 频 谱数 字 角 频 率幅 度0123-1-0.500.51幅度s时 间 （ ）2:1减 抽 样 信 号 波 形00.511.5201000200030002:1减 抽 样 信 号 频 谱w数 字 角 频 率幅度00.511.5-1-0.500.51幅度s时 间 （ ）5:1减 抽 样 信 号 波 形00.511.52050010005:1减 抽 样 信 号 频 谱w数 字 角 频 率幅度实验一［y,fs,bit］=wavread('I do’)％读取音乐片段,fs是采样率size(y）%求矩阵的行数和列数y1=y(:，1);%对信号进行分列处理n1=length(y1)；％取y的长度t1=(0：n1—1)/fs;％设置波形图横坐标figuresubplot（2,1，1)；plot（t1，y1）; ％画出时域波形图ylabel（'幅值’)；xlabel(’时间（s）’)；title（’信号波形');subplot（2，1，2)；Y1=fft(y1）；w1=2/n1＊（0：n1—1)；%设置角频率plot(w1，abs（Y1))；％画频谱图title（'信号频谱’);xlabel(’数字角频率');ylabel（'幅度'）;grid on;sound（y,fs); 实验二[y,fs，bit]=wavread（'I do'）y1=y（：,1）;n1=length(y1）;D=2；%设置抽样间隔 y2=y1（1：D:n1)；％减抽样n2=length（y2)；％减抽样后信号长度t2=(0:n2—1）/fs;％设置横坐标figure subplot（2，2，1）；plot(t2，y2)； %绘制减抽样信号波形图ylabel（'幅度');xlabel(’时间(s）')；title（'2：1减抽样信号波形’）；Y2=fft（y2）； ％对y2进行n2点fft谱分析w2=2/n2*[0:n2—1];subplot(2，2，3)；plot(w2,abs（Y2））;％绘制减抽样信号频谱图title(’2:1减抽样信号频谱’)；xlabel（'数字角频率w');ylabel(’幅度’）；grid on;sound(y2,fs/D）；D=5;％设置抽样间隔 y3=y1（1:D：n1);%减抽样n3=length(y3)；％减抽样后信号长度0246-1-0.500.51(s)时 间幅度调 制 后 信 号0246-1-0.500.51(s)时 间幅度调 制 后 信 号00.511.520100020003000w数 字 角 频 率幅度()调 制 后 信 号 的 频 谱 高 频 率 调 制00.511.520100020003000w数 字 角 频 率幅度()调 制 后 信 号 的 频 谱 低 频 率 调 制t3=（0:n3-1)/fs；％设置横坐标subplot（2，2，2)；plot(t3，y3）; ％绘制减抽样信号波形图ylabel('幅度’）；xlabel('时间（s）’);title（’5：1减抽样信号波形');Y3=fft（y3）; ％对y2进行n2点fft谱分析w3=2/n3*[0：n3-1］；subplot(2，2,4)；plot(...