【文章內容簡介】 (F2) 39。 Hz39。]。 text(3000,40,F2t)。 [m3,i3]=max(abs(k(10000:12021)))。 % 找出頻域次次大值 F3=(10000+i3)/32768*fs。 F3t=[num2str(F3) 39。 Hz39。]。 text(8000,20,F3t)。 %=========================================================== 23 音頻信號的分析與合成程序（續(xù) 3） %=========================================================== t=(0::1)。 y1=(sin(2*pi*788*t)+sin(2*pi*3174*t)*()+sin(2*pi*6936*t)*())*。 disp(39。按任意鍵顯示左聲道簡單合成聲音波形！！ 39。)。 disp(39。39。)。 pause plot(t,y1)。 set(gca,39。FontName39。,39。Arial39。,39。FontSize39。,10)。 title([39。左聲道簡單合成聲音波形 39。])。 disp(39。按任意鍵聽簡單合成左聲道信號三次?。?39。)。 disp(39。39。)。 pause sound(y1,fs,bits)。 % 數(shù)據(jù)通過聲卡轉換為聲音 pause sound(y1,fs,bits)。 pause sound(y1,fs,bits)。 [k1]=fft(y1,32768)。 disp(39。按任意鍵顯示左聲道簡單合成聲音的頻域幅值??！ 39。)。 disp(39。39。)。 24 音頻信號的分析與合成程序（續(xù) 4） disp(39。39。)。 pause plot(abs(k1))。 set(gca,39。FontName39。,39。Arial39。,39。FontSize39。,10)。 title([39。左聲道簡單合成聲音的頻域幅值 39。])。 xlabel([39。采樣點數(shù) 39。])。 [m1_r,i1_r]=max(abs(k1))。 % 找出簡單合成聲音頻域最大值 F1_r=i1_r/32768*10000。 F1_rt=[num2str(F1_r) 39。 Hz39。]。 text(3000,850,F1_rt)。 [m2_r,i2_r]=max(abs(k1(10000:12021)))。 F2_r=(10000+i2_r)/32768*10000。 % 找出簡單合成聲音頻域次大值 F2_rt=[num2str(F2_r) 39。 Hz39。]。 text(9000,60,F2_rt)。 end。 %=========================================================== 25 音頻頻譜分析儀程序 （ 車晴 ） 返回 y=audspec(action)。 % AUDIO SPECTRUM ANALYSER THEE 2021/10 %=================================================================================== global x t dt fs c1 c2 H_osc H2 fn y ystr Ktype if nargin1。action=39。initialized39。end。 if strcmp(action,39。initialized39。)。 % 示波器 clc。 % 清除工作區(qū) v=version。 % MATLAB版本 [fname,pname]=uigetfile(39。*.wav39。,39。Open Wave File39。)。% 打開對話框輸入文件 file=[pname,fname]。 % 文件名稱和路徑 if strcmp(v(1),39。539。)。 % 讀入聲音文件（ *.wav） [x,fs,nbits]=wavread(file)。 % x是采樣數(shù)據(jù)， fs是采樣速率（ Hz） [m,n]=size(x)。 % nbits是每次采樣比特數(shù) file=[file 39。 fs=39。 int2str(fs) 39。Hz,39。 int2str(nbits) 39。bits39。]。 if n==2。file=[file 39。,stero39。]。 x=x(:)。 % 若是雙聲道，將數(shù)據(jù)合并為一列 %x=x(:,1)。 % 取左聲道數(shù)據(jù) %x=x(:,2)。 % 取右聲道數(shù)據(jù) end。 elseif strcmp(v(1),39。439。)。 [x,fs]=wavread(file)。 file=[file 39。fs=39。 int2str(fs) 39。Hz,39。]。 end。 26 音頻頻譜分析儀程序（續(xù) 1） % 增加菜單項（ Zoom， Play， Spectrum）， 生成信號波形，回應（ Callback） 分別調用 % 函數(shù) audspec(39。Zoom39。)， audspec(39。Play39。)和 audspec(39。Spectrum39。)。 H_osc=figure(39。Name39。,[39。OSCILLOSCOPE THEE 2021/839。 blanks(6) file],... 39。NumberTitle39。,39。off39。,39。Position39。,[10 78 640 480])。 uimenu(39。Label39。,39。amp。Zoom39。,39。Callback39。,39。audspec(39。39。Zoom39。39。)。39。)。 uimenu(39。Label39。,39。amp。Play39。,39。Callback39。,39。audspec(39。39。Play39。39。)。39。)。 uimenu(39。Label39。,39。amp。Spectrum39。,39。Callback39。,39。audspec(39。39。Spectrum39。39。)。39。)。 if strcmp(v(1),39。439。)。 whitebg(39。k39。)。c1=39。y39。c2=39。g39。 elseif strcmp(v(1),39。539。)。 c1=39。r39。c2=39。b39。 end。 dt=1/fs。N=length(x)。t=(0:N1)*dt。T=N*dt。ss=39。s39。 % dt=1/fs， fs是采樣速率（ Hz） if T1e03。t=t*1e+06。ss=39。us39。 elseif T1。t=t*1000。ss=39。ms39。 end。 plot(t,x,c1)。 % 畫左右聲道信號波形 set(gca,39。FontName39。,39。Arial39。,39。FontSize39。,10)。 % gca是返回當前軸的句柄并設置 title([39。畫左右聲道聲音波形 39。])。xlabel([39。T (39。 ss 39。)39。])。 n1=1。n2=N。H2=1。 27 音頻頻譜分析儀程序（續(xù) 2） elseif strcmp(action,39。Play39。)。 % 播放聲音文件 v=axis。 % 返回當前軸坐標的設置 del=t(2)t(1)。del=del*.5。 % v(1)和 v(2)是 X軸坐標的設置 n1=find(abs(tv(1))del)。n2=find(abs(tv(2))del)。 if isempty(n2)。n2=length(x)。end。 xx=x(n1:n2)。sound(xx,fs)。 % 以采樣速率 fs輸出聲音文件數(shù)據(jù) elseif strcmp(action,39。Zoom39。)。 % 放大 zoom(H_osc,39。xon39。)。