jiandanhecheng.m

选取windows系统自带的ding.wav信号作为分析对象

[x,fs,bits]=wavread('ding.wav');       % 读入声音文件（*.wav）      

y=x(:,1);                        % 左声道信号波形数据
subplot(311)
plot(y);
title(['左声道声音波形']);
xlabel(['采样点数']);
%sound(y,fs,bits);                % 数据通过声卡转换为声音


[k]=fft(y,32768);
subplot(312)
plot(abs(k));
title(['左声道声音频域的幅值']);
xlabel(['采样点数']);

[m1,i1]=max(abs(k));              % 找出频域最大值            
F1=i1/32768*fs;
F1t=[num2str(F1) ' Hz'];
text(1500,450,F1t);   

[m2,i2]=max(abs(k(4000:5000)));   % 找出频域次大值
F2=(4000+i2)/32768*fs;
F2t=[num2str(F2) ' Hz'];
text(3000,40,F2t);
[m3,i3]=max(abs(k(10000:16390))); % 找出频域次次大值
F3=(10000+i3)/32768*fs;
F3t=[num2str(F3) ' Hz'];
text(8000,20,F3t);

t=(0:0.0001:1);
y1=(sin(2*pi*788*t)+sin(2*pi*3174*t)*(22.65/490)+sin(2*pi*6936*t)*(2.996/490))*0.18;


subplot(313)
plot(t,y1);