【正文】 采用 uigetfile 函數(shù)，該函數(shù)打開(kāi)一個(gè)標(biāo)準(zhǔn)的對(duì)話框，然后可以選擇聲音文件的路徑和文件的名稱(chēng)。 ? 步驟四：進(jìn)行該聲波主要頻譜的分析 。 12 2． 5 程序設(shè)計(jì)實(shí)例 2． 5． 1 音頻信號(hào)的分析與合成 （韓利竹 ， ） ? 采用 MATLAB分析 WAV文件。 Clear All Breakpoints： 清除程序中的所有斷點(diǎn) 。 Single Step： 單步執(zhí)行函數(shù)。 ? 通過(guò)將原 M函數(shù)文件的函數(shù)申明行注釋掉，可使一個(gè)中間變量難于觀察的 M函數(shù)文件變?yōu)橐粋€(gè)所有變量都保留在基本工作空間中的 M腳本文件。 ? 利用 echo 指令，使運(yùn)行時(shí)在屏幕上逐行顯示文件內(nèi)容。 ? 有兩種調(diào)試方法：直接調(diào)試法和工具調(diào)試法。對(duì)于這類(lèi)錯(cuò)誤，通常能在運(yùn)行時(shí)發(fā)現(xiàn)，終止執(zhí)行，并給出相應(yīng)的錯(cuò)誤原因以及所在行號(hào)。若有，先運(yùn)行 Callback 屬性定義的子程序，再顯示子菜單。 ? 控件的 Callback 屬性 ： Callback 屬性的取值是字符串，可以是某個(gè) M文件名或一小段 MATLAB語(yǔ)句 。 ? 警示指令： error， warning 。 for—end 循環(huán)和 whileend循環(huán)。 具有接受 “任意多輸入 ” 、返回 “任意多輸出 ”的能力 。 ? 函數(shù)體內(nèi)的 nargout 給出調(diào)用該函數(shù)時(shí)的輸出參數(shù)數(shù)目。 ? 全局（ Global） 變量：通過(guò) global 指令， MATLAB也允許幾個(gè)不同的函數(shù)空間以及基本工作空間共享同一個(gè)變量，這種被共享的變量稱(chēng)為全局變量。 ? 編寫(xiě)和修改記錄：與在線幫助文本區(qū)相隔一個(gè)“空”行，也以 %開(kāi)頭，標(biāo)志編寫(xiě)及修改該 M文件的作者和日期等 。 ? 典型 M函數(shù)文件的結(jié)構(gòu)如下 ： ? 函數(shù)申明行：位于函數(shù)文件的首行，以關(guān)鍵字 functio 開(kāi)頭，函數(shù)名以及函數(shù)的輸入輸出宗量都在這一行被定義。 ?從運(yùn)行上看 ，與腳本文件運(yùn)行不同 ，每當(dāng)函數(shù)文件運(yùn)行， MATLAB就會(huì)專(zhuān)門(mén)為它開(kāi)辟一個(gè)臨時(shí)工作空間，稱(chēng)為 函數(shù)工作空間 （ Function workspace） 。 這兩種文件的擴(kuò)展名，均為“ . m” 。 ? M腳本文件的特點(diǎn)是： ? 腳本文件的構(gòu)成比較簡(jiǎn)單，只是一串按用戶(hù)意圖排列而成的（包括控制流向指令在內(nèi)的） MATLAB指令集合。1 第二講 MATLAB的 2． 1 腳本文件和函數(shù)文件 2． 2 函數(shù)調(diào)用和參數(shù)傳遞 2． 3 MATLAB的程序結(jié)構(gòu)和控制流 2． 4 M文件的調(diào)試 2． 5 程序設(shè)計(jì)實(shí)例 2． 5． 1 音頻信號(hào)的分析與合成 2． 5． 2 音頻頻譜分析儀 2． 5． 3 幅度調(diào)制 程序設(shè)計(jì) 2 2． 1 腳本文件和函數(shù)文件 2． 1． 1 M腳本文件 ? 對(duì)于一些比較簡(jiǎn)單的問(wèn)題 ，在指令窗中直接輸入指令計(jì)算 。 ? 腳本文件運(yùn)行后 ，所產(chǎn)生的所有變量都駐留在 MATLAB基本工作空間 （ Base workspace） 中 。 3 2． 1 腳本文件和函數(shù)文件（ 續(xù) 1） 2． 1． 2 M函數(shù)文件 ? 與腳本文件不同 ，函數(shù)文件猶如一個(gè)“黑箱”，把一些數(shù)據(jù)送進(jìn)并經(jīng)加工處理，再把結(jié)果送出來(lái)。當(dāng)執(zhí)行文件最后一條指令時(shí) ，就結(jié)束該函數(shù)文件的運(yùn)行，同時(shí)該臨時(shí)函數(shù)空間及其所有的中間變量就立即被清除。 ? 笫一注釋行：緊隨函數(shù)申明行之后以 %開(kāi)頭笫一注釋行。 ? 函數(shù)體：為清晰起見(jiàn)，它與前面的注釋以“空”行相隔。 2． 2． 2 函數(shù)調(diào)用 ? 在 MATLAB中，調(diào)用函數(shù)的常用形式是： [輸出參數(shù) 1,輸出參數(shù) 2,…] = 函數(shù)名 (輸入?yún)?shù) 1,輸入?yún)?shù) 2, …) ? 函數(shù)調(diào)用可以嵌套，一個(gè)函數(shù)可以調(diào)用別的函數(shù)，甚至調(diào)用它自己 （遞歸調(diào)用）。 ? 只要在函數(shù)文件中包括這兩個(gè)變量，就可以知道該函數(shù)文件調(diào)用時(shí)的輸入?yún)?shù)和輸出參數(shù)數(shù)目。 ? 跨空間變量傳遞： evalin。 2． 3． 2 程序流控制 ? 常用指令： return， echo， input， pause， keyboard， break。 8 2． 3 MATLAB的程序結(jié)構(gòu)和控制流（ 續(xù)） 2． 3． 3 圖形用戶(hù)界面（ GUI） 編程 ? 現(xiàn)代的主流應(yīng)用程序已經(jīng)從命令行的交互方式轉(zhuǎn)變?yōu)橐詧D形界面為主的交互方式，這主要是由于它給用戶(hù)帶來(lái)了操作和控制的方便與靈活性。 當(dāng)用戶(hù)激活控件對(duì)象（例如 ：在控件對(duì)象圖標(biāo)上單擊鼠標(biāo)左鍵 ）時(shí)，應(yīng)用程序就運(yùn)行 Callback 屬性定義的子程序。 9 2． 4 M文件的調(diào)試 ? 編寫(xiě) M文件時(shí)，錯(cuò)誤（ Bug） 在所難免。 ? 運(yùn)行錯(cuò)誤是算法本身引起的，發(fā)生在運(yùn)行過(guò)程中。 10 2． 4 M文件的調(diào)試（ 續(xù) 1） ? 直接調(diào)試法：可以用下面方法發(fā)現(xiàn)某些運(yùn)行錯(cuò)誤。echo on 能顯示 M腳本文件； echo FunNsme on 能顯示名為FunNsme 的 M函數(shù)文件。 11 2． 4 M文件的調(diào)試（ 續(xù) 2） ? GUI 界面調(diào)試法： MATLAB 版提供了一個(gè)基于 GUI界面的調(diào)試。 Step In： 深入下層局部工作區(qū) 。 Stop if Error： 運(yùn)行至出錯(cuò)或結(jié)束。 ? 步驟一：選擇一個(gè) WAV文件作為分析的對(duì)象 。 ? 步驟五：根據(jù)該聲音的頻譜，反演時(shí)域圖形 。 ? 信號(hào)波形生成使用繪圖函數(shù)就可以完成，而信號(hào)的頻譜是通過(guò)快速 付立葉變換 而得到的。 strcmp(action,39。Zoom39。)。 Am? DSBAM： ttmts ?s i n)()( ??? SSB： 是調(diào)制信號(hào) 的希爾伯特 變換 。 clf。 elseif x(j)+y(i)=1 z(i,j)=*exp(*y(i)^*x(j)^2+*x(j))。surf(x,y,z)。 18 控件的 Callback 屬性 返回 （張志涌 P426） ? 對(duì)于傳遞函數(shù)為 的歸一化二階系統(tǒng)，制作一 個(gè)能繪制該系統(tǒng)單位階躍響應(yīng)的圖形用戶(hù)界面。 ? 運(yùn)行程序 ， 驗(yàn)證 Callback 屬性定義的子程序。 % SIGNAL SPECTRUM ANALYSER THEE 2021/10 % y=sigspec(action)。action=39。*.wav39。 file=[pname,fname]。 % 雙聲道信號(hào)波形數(shù)據(jù) y=x(:,1)。)。 pause plot(y)。Arial39。 title([39。采樣點(diǎn)數(shù) 39。)。 pause sound(y,fs,bits)。 disp(39。39。)。,39。,10)。 xlabel([39。 % 找出頻域最大值 F1=i1/32768*fs。 text(1500,450,F1t)。 Hz39。 % 找出頻域次次大值 F3=(10000+i3)/32768*fs。 text(8000,20,F3t)。按任意鍵顯示左聲道簡(jiǎn)單合成聲音波形??！ 39。)。,39。,10)。 disp(39。39。 pause sound(y1,fs,bits)。)。 24 音頻信號(hào)的分析與合成程序（續(xù) 4） disp(39。 set(gca,39。,39。左聲道簡(jiǎn)單合成聲音的頻域幅值 39。])。 Hz39。 F2_r=(10000+i2_r)/32768*10000。 text(9000,60,F2_rt)。action=39。initialized39。 % MATLAB版本 [fname,pname]=uigetfile(39。)。)。 fs=39。bits39。,stero39。 % 取左聲道數(shù)據(jù) %x=x(:,2)。)。 int2str(fs) 39。 26 音頻頻譜分析儀程序（續(xù) 1） % 增加菜單項(xiàng)（ Zoom， Play， Spectrum）， 生成信號(hào)波形，回應(yīng)（ Callback） 分別調(diào)用 % 函數(shù) audspec(39。)和 audspec(39。Name39。NumberTitle39。Position39。,39。Callback39。Zoom39。)。amp。,39。39。 uimenu(39。Spectrum39。audspec(39。)。439。)。g39。 c1=39。 end。T=N*dt。t=t*1e+06。t=t*1000。 plot(t,x,c1)。Arial39。 % gca是返回當(dāng)前軸的句柄并設(shè)置 title([39。T (39。 n1=1。Play39。del=del*.5。n2=length(x)。 % 以采樣速率 fs輸出聲音文件數(shù)據(jù) elseif strcmp(action,39。xon39。)。 % v(1)和 v(2)是 X軸坐標(biāo)的設(shè)置 n1=find(ab