【正文】 tep： 單步執(zhí)行函數(shù)。 Quit Debugging： 退出調(diào)試狀態(tài)。 Clear All Breakpoints： 清除程序中的所有斷點(diǎn) 。 Stop if Warning： 運(yùn)行至警告消息或結(jié)束。 12 2． 5 程序設(shè)計(jì)實(shí)例 2． 5． 1 音頻信號(hào)的分析與合成 （韓利竹 ， ） ? 采用 MATLAB分析 WAV文件。（ ） ? 步驟二：讀 WAV文件數(shù)據(jù)并畫時(shí)域圖形 。 ? 步驟四：進(jìn)行該聲波主要頻譜的分析 。（有 失真 ） ? 常用函數(shù)： ? [x， fs， bits]=waveread(?filename?) ? [d]= FFT (w， l) ? sound(w， fs， bits) ? 步驟六：進(jìn)行 付立葉逆變換 IFFT并畫頻域圖形 。 ? WAV文件的輸入方式采用 uigetfile 函數(shù)，該函數(shù)打開(kāi)一個(gè)標(biāo)準(zhǔn)的對(duì)話框，然后可以選擇聲音文件的路徑和文件的名稱。 ? 常用函數(shù)： strcmp(action,39。)。Play39。 strcmp(action,39。)。Spectrum39。 14 2． 5． 3 幅度調(diào)制 （ 車晴 P203 ） ? 幅度調(diào)制， 即載波幅度隨調(diào)制信號(hào)變化的調(diào)制（ AM）。 ? 假定調(diào)制信號(hào)為 m(t)， 其幅度為 1v ， 載波信號(hào)為 ，各種調(diào)幅方式的已調(diào)波的表達(dá)式如下 : t?sin? SAM: ttmmtsA ?s i n)](1[)( ????是調(diào)制度。 ? QAM： ttmttmts ?? c o s)(2s i n)(1)( ????ttmttmts ?? s i n)( o s)()( ?????? ?ttmttmts ?? s i n)( o s)()( ?????? ?)(tm? )(tm15 M腳本文件入門 返回 通過(guò) M腳本文件，畫出下列分段函數(shù)所表示的曲面。b=2。 x=a::a。 for i=1:length(y) for j=1:length(x) if x(j)+y(i)1 z(i,j)=*exp(*y(i)^*x(j)^*x(j))。 else z(i,j)=*exp(y(i)^26.*x(j)^2)。 colormap(flipud(winter))。 17 例 2_2_1： M函數(shù)文件示例。張志涌 （ P299） 程序文本。 1212 ??? ssG ?? 在編輯框輸入 ， 可以看到 一根曲線 。 ? 本例設(shè)置斷點(diǎn)演示：（ A） 圖形界面的大致生成 過(guò)程； （ B） 靜態(tài)文本和編輯框的生成； （ C） 坐標(biāo)方格控制鍵的形成； （ D） 如何使用該界面。 19 菜單的 Callback 屬性 返回 本例的目標(biāo)是：在圖形窗上自制一個(gè)名為【 Test】 的“頂層菜單項(xiàng)”；當(dāng)用鼠標(biāo)點(diǎn)動(dòng)該菜單項(xiàng)時(shí)，將產(chǎn)生一個(gè)帶分格的封閉坐標(biāo)軸。 （ 張志涌 P432） 20 音頻信號(hào)的分析與合成程序 （韓利竹 ， ） 返回 y=sigspec(action)。 % Usage:y=sigspec。 if nargin1。initialized39。 [fname,pname]=uigetfile(39。,39。)。 [x,fs,bits]=wavread(file)。 % 數(shù)據(jù)通過(guò)聲卡轉(zhuǎn)換為聲音 %==========選擇聲道 ========================================== %y=x(:)。 % 左聲道信號(hào)波形數(shù)據(jù) %y=x(:,2)。按任意鍵顯示左聲道聲音波形??！ 39。 disp(39。)。 set(gca,39。,39。,39。,10)。左聲道聲音波形 39。 xlabel([39。])。按任意鍵聽(tīng)左聲道信號(hào)三次??！ 39。 disp(39。)。 % 數(shù)據(jù)通過(guò)聲卡轉(zhuǎn)換為聲音 pause sound(y,fs,bits)。 [k]=fft(y,32768)。按任意鍵顯示左聲道聲音頻域的幅值??！ 39。 disp(39。)。39。 pause plot(abs(k))。FontName39。Arial39。FontSize39。 title([39。])。采樣點(diǎn)數(shù) 39。 [m1,i1]=max(abs(k))。 F1t=[num2str(F1) 39。]。 [m2,i2]=max(abs(k(4000:5000)))。 F2t=[num2str(F2) 39。]。 [m3,i3]=max(abs(k(10000:12021)))。 F3t=[num2str(F3) 39。]。 %=========================================================== 23 音頻信號(hào)的分析與合成程序（續(xù) 3） %=========================================================== t=(0::1)。 disp(39。)。39。 pause plot(t,y1)。FontName39。Arial39。FontSize39。 title([39。])。按任意鍵聽(tīng)簡(jiǎn)單合成左聲道信號(hào)三次??！ 39。 disp(39。)。 % 數(shù)據(jù)通過(guò)聲卡轉(zhuǎn)換為聲音 pause sound(y1,fs,bits)。 [k1]=fft(y1,32768)。按任意鍵顯示左聲道簡(jiǎn)單合成聲音的頻域幅值！！ 39。 disp(39。)。39。 pause plot(abs(k1))。FontName39。Arial39。FontSize39。 title([39。])。采樣點(diǎn)數(shù) 39。 [m1_r,i1_r]=max(abs(k1))。 F1_rt=[num2str(F1_r) 39。]。 [m2_r,i2_r]=max(abs(k1(10000:12021)))。 % 找出簡(jiǎn)單合成聲音頻域次大值 F2_rt=[num2str(F2_r) 39。]。 end。 % AUDIO SPECTRUM ANALYSER THEE 2021/10 %=================================================================================== global x t dt fs c1 c2 H_osc H2 fn y ystr Ktype if nargin1。initialized39。 if strcmp(action,39。)。 % 清除工作區(qū) v=version。*.wav39。Open Wave File39。% 打開(kāi)對(duì)話框輸入文件 file=[pname,fname]。539。 % 讀入聲音文件（ *.wav） [x,fs,nbits]=wavread(file)。 % nbits是每次采樣比特?cái)?shù) file=[file 39。 int2str(fs) 39。 int2str(nbits) 39。]。file=[file 39。]。 % 若是雙聲道，將數(shù)據(jù)合并為一列 %x=x(:,1)。 % 取右聲道數(shù)據(jù) end。439。 [x,fs]=wavread(file)。fs=39。Hz,39。 end。Zoom39。Play39。Spectrum39。 H_osc=figure(39。,[39。 blanks(6) file],... 39。,39。,39。,[10 78 640 480])。Label39。amp。,39。,39。39。39。39。 uimenu(39。,39。Play39。Callback39。audspec(39。Play39。)。)。Label39。amp。,39。,39。39。39。39。 if strcmp(v(1),39。)。k39。c1=39。c2=39。 elseif strcmp(v(1),39。)。r39。b39。 dt=1/fs。t=(0:N1)*dt。ss=39。 % dt=1/fs， fs是采樣速率（ Hz） if T1e03。ss=39。 elseif T1。ss=39。 end。 % 畫左右聲道信號(hào)波形 set(gca,39。,39。,39。,10)。畫左右聲道聲音波形 39。xlabel([39。 ss 39。])。n2=N。 27 音頻頻譜分析儀程序（續(xù) 2） elseif strcmp(action,39。)。 % 返回當(dāng)前軸坐標(biāo)的設(shè)置 del=t(2)t(1)。