【正文】 ? 步驟五：根據(jù)該聲音的頻譜，反演時(shí)域圖形 。（有 失真 ） ? 常用函數(shù)： ? [x， fs， bits]=waveread(?filename?) ? [d]= FFT (w， l) ? sound(w， fs， bits) ? 步驟六：進(jìn)行 付立葉逆變換 IFFT并畫頻域圖形 。 13 2． 5． 2 音頻頻譜分析儀 （ 車晴 ） ? M函數(shù)文件 用于對(duì)音頻頻譜分析儀進(jìn)行仿真，它通過鼠標(biāo)來選擇一個(gè) WAV文件，它可以播放該聲音文件、顯示音頻輸入信號(hào)的波形和頻譜。 ? WAV文件的輸入方式采用 uigetfile 函數(shù)，該函數(shù)打開一個(gè)標(biāo)準(zhǔn)的對(duì)話框，然后可以選擇聲音文件的路徑和文件的名稱。 ? 信號(hào)波形生成使用繪圖函數(shù)就可以完成，而信號(hào)的頻譜是通過快速 付立葉變換 而得到的。 ? 常用函數(shù)： strcmp(action,39。initialized39。)。 strcmp(action,39。Play39。)。 strcmp(action,39。Zoom39。)。 strcmp(action,39。Spectrum39。)。 14 2． 5． 3 幅度調(diào)制 （ 車晴 P203 ） ? 幅度調(diào)制， 即載波幅度隨調(diào)制信號(hào)變化的調(diào)制（ AM）。幅度調(diào)制有多種實(shí)現(xiàn)方式，包括：標(biāo)準(zhǔn)幅度調(diào)制（ SAM），雙邊帶幅度調(diào)制（ DSBAM）， 單邊帶幅度調(diào)制（ SSB），殘留邊帶調(diào)制（ VSB） 和平衡正交幅度調(diào)制（ QAM） 等。 ? 假定調(diào)制信號(hào)為 m(t)， 其幅度為 1v ， 載波信號(hào)為 ，各種調(diào)幅方式的已調(diào)波的表達(dá)式如下 : t?sin? SAM: ttmmtsA ?s i n)](1[)( ????是調(diào)制度。 Am? DSBAM： ttmts ?s i n)()( ??? SSB： 是調(diào)制信號(hà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腳本文件入門 返回 通過 M腳本文件，畫出下列分段函數(shù)所表示的曲面。 ?????????????????????????),(21216212112122212212122xxexxexxexxpxxxxxxxx （ ） 16 M腳本文件入門（續(xù)） 返回 % This is my first example. a=2。b=2。 clf。 x=a::a。y=b::b。 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))。 elseif x(j)+y(i)=1 z(i,j)=*exp(*y(i)^*x(j)^2+*x(j))。 else z(i,j)=*exp(y(i)^26.*x(j)^2)。 end end end axis([a,a,b,b,min(min(z)),max(max(z))])。 colormap(flipud(winter))。surf(x,y,z)。 17 例 2_2_1： M函數(shù)文件示例。 返回 sa = circle(r,s) % CIRCLE plot a circle of radii r in the line specified by s. % r 指定半徑的數(shù)值 % s 指定線色的字符串 % sa 圓面積 % circle(r) 利用藍(lán)實(shí)線畫半徑為 r 的圓周線 . % circle(r,s) 利用串 s 指定的線色畫半徑為 r 的圓周線 . % sa=circle(r) 計(jì)算圓面積，并畫半徑為 r 的藍(lán)色圓面 . % sa=circle(r,s) 計(jì)算圓面積，并畫半徑為 r 的 s 色圓面 . % 編寫于 1999年 4月 7日，修改于 1999年 8月 27日。張志涌 （ P299） 程序文本。 18 控件的 Callback 屬性 返回 （張志涌 P426） ? 對(duì)于傳遞函數(shù)為 的歸一化二階系統(tǒng)，制作一 個(gè)能繪制該系統(tǒng)單位階躍響應(yīng)的圖形用戶界面。 1212 ??? ssG ?? 在編輯框輸入 ， 可以看到 一根曲線 。 ? 在編輯框輸入 ::1， 可以看到十根曲線 。 ? 本例設(shè)置斷點(diǎn)演示：（ A） 圖形界面的大致生成 過程； （ B） 靜態(tài)文本和編輯框的生成； （ C） 坐標(biāo)方格控制鍵的形成； （ D） 如何使用該界面。 ? 運(yùn)行程序 ， 驗(yàn)證 Callback 屬性定義的子程序。 19 菜單的 Callback 屬性 返回 本例的目標(biāo)是：在圖形窗上自制一個(gè)名為【 Test】 的“頂層菜單項(xiàng)”；當(dāng)用鼠標(biāo)點(diǎn)動(dòng)該菜單項(xiàng)時(shí)，將產(chǎn)生一個(gè)帶分格的封閉坐標(biāo)軸。通過本例說明： （ A） 回調(diào)屬性的運(yùn)作機(jī)理； （ B） 用戶頂層菜單項(xiàng)的制作； （ C） uimenu屬性的設(shè)置方法； （ D） 復(fù)雜字符串的構(gòu)成方法和注意事項(xiàng)。 （ 張志涌 P432） 20 音頻信號(hào)的分析與合成程序 （韓利竹 ， ） 返回 y=sigspec(action)。 % SIGNAL SPECTRUM ANALYSER THEE 2021/10 % y=sigspec(action)。 % Usage:y=sigspec。 %=========================================================== clear。 if nargin1。action=39。initialized39。end。 [fname,pname]=uigetfile(39。*.wav39。,39。Open Wave File39。)。 file=[pname,fname]。 [x,fs,bits]=wavread(file)。 % 讀入聲音文件（ *.wav） sound(x,fs,bits)。 % 數(shù)據(jù)通過聲卡轉(zhuǎn)換為聲音 %==========選擇聲道 ========================================== %y=x(:)。 % 雙聲道信號(hào)波形數(shù)據(jù) y=x(:,1)。 % 左聲道信號(hào)波形數(shù)據(jù) %y=x(:,2)。 % 右聲道信號(hào)波形數(shù)據(jù) %=========================================================== 21 音頻信號(hào)的分析與合成程序（續(xù) 1） %=========================================================== disp(39。按任意鍵顯示左聲道聲音波形！！ 39。)。 disp(39。39。)。 pause plot(y)。 set(gca,39。FontName39。,39。Arial39。,39。FontSize39。,10)。 title([39。左聲道聲音波形 39。])。 xlabel([39。采樣點(diǎn)數(shù) 39。])。 disp(39。按任意鍵聽左聲道信號(hào)三次??！ 39。)。 disp(39。39。)。 pause sound(y,fs,bits)。 % 數(shù)據(jù)通過聲卡轉(zhuǎn)換為聲音 pause sound(y,fs,bits)。 pause sound(y,fs,bits)。 [k]=fft(y,32768)。 disp(39。按任意鍵顯示左聲道聲音頻域的幅值?。?39。)。 disp(39。39。)。 22 音頻信號(hào)的分析與合成程序（續(xù) 2） disp(39。39。)。 pause plot(abs(k))。 set(gca,39。FontName39。,39。Arial39。,39。FontSize39。,10)。 title([39。左聲道聲音頻域的幅值 39。])。 xlabel([39。采樣點(diǎn)數(shù) 39。])。 [m1,i1]=max(abs(k))。 % 找出頻域最大值 F1=i1/32768*fs。 F1t=[num2str(F1) 39。 Hz39。]。 text(1500,450,F1t)。 [m2,i2]=max(abs(k(4000:5000)))。 % 找出頻域次大值 F2=(4000+i2)/32768*fs。 F2t=[num2str