资源描述：
《matlab的语音信号频谱分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验一[y,fs,bit]=wavread('Ido')%读取音乐片段，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);%设
2、置角频率plot(w1,abs(Y1));%画频谱图title('信号频谱');xlabel('数字角频率');ylabel('幅度');gridon;sound(y,fs);实验二[y,fs,bit]=wavread('Ido')y1=y(:,1);n1=length(y1);D=2;%设置抽样间隔y2=y1(1:D:n1);%减抽样n2=length(y2);%减抽样后信号长度t2=(0:n2-1)/fs;%设置横坐标figuresubplot(2,2,1);plot(t2,y2);%绘制减抽样
3、信号波形图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('幅度');gridon;sound(y2,fs/D);D=5;%设置抽样间隔y3=y1(1:D:n1);%减抽样n3=length(y
4、3);%减抽样后信号长度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(w3,abs(Y3));%绘制减抽样信号频谱图title('5:1减抽样信号频谱');xlabel('数字角频率w');ylabel('幅
5、度');gridon;sound(y3,fs/D)实验三[y,fs,bit]=wavread('Ido')y1=y(:,1);n1=length(y1);n3=0:n1-1;b1=cos(0.75*pi*n3);%设置调制信号b2=cos(0.5*pi*n3);c1=b1'.*y1;%对原信号进行调制c2=b2'.*y1;lc1=length(c1);t=(0:lc1-1)/fs;figure%用载波对信号进行调制，并对其做fft变换subplot(2,2,1)%获取频谱，从图中可以观察到，调制后的
6、plot(t,c1);%信号频谱发生搬移xlabel('时间(s)');ylabel('幅度');title('调制后信号');subplot(2,2,2)%获取频谱，从图中可以观察到，调制后的plot(t,c2);%信号频谱发生搬移xlabel('时间(s)');ylabel('幅度');title('调制后信号');w1=2/lc1*[0:lc1-1];%设置角频率WC1=fft(c1);C2=fft(c2);subplot(2,2,3)plot(w1,abs(C1));xlabel('数字角频
7、率w');ylabel('幅度');title('调制后信号的频谱(高频率调制)');gridon;subplot(2,2,4)plot(w1,abs(C2));xlabel('数字角频率w');ylabel('幅度');title('调制后信号的频谱(低频率调制)');gridon;sound(c1,fs);实验四clearall;cla;closeall[a,fs,bit]=wavread('Ido');y1=a(:,1);%去单列数据进行分析f1=fft(y1);n=length(f1);tn
8、=(0:n-1)/fs;w=2/n*[0:n-1];%sound(y1,fs);figure(1)subplot(2,3,1);plot(tn,y1);gridon;title('音频');xlabel('时间');ylabel('幅度');subplot(2,3,4);plot(w,abs(f1));gridon;title('频谱');xlabel('频率/pi');ylabel('幅度');t=[0:n-1];y2=cos(pi*1/2*t);%载波函数y3=y
