音频信号的幅度与相位

[Y，fs，bits]= wavread(Blip', N)用于读取音频，对音频信
号进行采集，采样值放在参数Y中，fs表示每秒采样
点数，即采样频率，bits表示每个采样点在编码时所
占位数.N表示采样点总数.参数’Blip’为音频所在
地址，如:'C:\yinpinl'.
    调用函数fft可对己采集音频信号进行时频转
换.通过函数abs()和angle ( )可分别得到信号频谱
的幅频图和相频图
    分别对放在C盘根目录下的女性音频信号
"yinpinl.wav'，和男性音频信号“yinpin2. wav”进行
采样和频谱分析.具体代码如下:
    %实现对音频信号的采样与频谱分析
    clear all;
    [Y，Fs0，bits]=wavread('C:\yinpinl');
    a=round(length(y0));
    y01=fft(y0);y011=abs     t0=linspace(0，a/ Fs0，a);
    figure(1);subplot<2，1，1);plot(t0，YO);
    title('(a)','fonts',10.5 ,'position ,[25，-1.5];
    xlabel(‘时间t/s','fonts',10.5 ,'position ,[45，-1.4]);
    ylabel(‘幅度n', 'fonts', 10. 5,'position',[-3,0.7]
    f0=linspace(0，Fs0，a);

    figure(2);subplot<2，1，1);plot(f0，y011);
    title(‘(a)‘，‘fonts', 10. 5，‘position',巨8000，一
}80」);
    xlabel(‘频率f/ha', 'fonts', 10. 5, ' position',
巨14000，一soot);
    ylabel(‘幅度ri ,'fonts',10. 5 ,'position ,巨一1500,
2000)
    figure(3);subplot<2，1，1);
    plot(f0(1:round(a/16))，y011(1:round(a/
16)));
    title(‘(a)‘，‘fonts‘，10. 5，‘position‘，巨500，一
}80」);
    xlabel(‘频率f/hz ,'fonts',10.5 ,'position ,巨850,
一soo」);
    ylabel(‘幅度ri ,'fonts', 10. 5,'position',巨-90,
2000一)
    巨y1，Fsl，nbitsl}=wavread('C:\yinpin2');
    b1=round(length(yl));
    yll=fft(yl);ylll=abs     t1=linspace<0，b1/ Fsl，b1);
    figure(1);subplot<2，1，2);plot(tl，y1);
    title('(b)','fonts',10.5,'positiori,[15，-1.5]);
    xlabel(‘时间t/s','fonts',10.5 ,'position ,[27，-1.4]);
    ylabel(‘幅度ri ,'fonts',10. 5 ,'position ,[-1.7,0.7]);
    fl=linspace<0，Fsl，b1);
    figure(2);subplot<2，1，2);plot(fl，ylll);
    title(‘(b)‘，'fonts‘，10. 5，‘position',[8000，-280]);
    xlabel(‘频率f/ha', 'fonts', 10. 5, ' position',[4000，-200]);
    yabel(‘幅度ri ,'fonts',10. 5 ,'position ,[-1500,660]);
    b=round(a*(Fs0/a))/(Fsl/bl);
    figure(3);subplot<2，1，2);
    plot(fl(1:round(b/16))，ylll(1:round(b/16)));
    title(‘(b)‘，‘fonts‘，10. 5，‘position‘，[500，-260]);
    xlabel(‘频率f/hz ,'fonts',10.5 ,'position ,[850,-200]）;
    ylabel(‘幅度ri ,'fonts', 10. 5,'position',[-90,-670]）；