【正文】 句 。對(duì)于這類(lèi)錯(cuò)誤，通常能在運(yùn)行時(shí)發(fā)現(xiàn)，終止執(zhí)行，并給出相應(yīng)的錯(cuò)誤原因以及所在行號(hào)。 ? 利用 echo 指令，使運(yùn)行時(shí)在屏幕上逐行顯示文件內(nèi)容。 Single Step： 單步執(zhí)行函數(shù)。 12 2． 5 程序設(shè)計(jì)實(shí)例 2． 5． 1 音頻信號(hào)的分析與合成 （韓利竹 ， ） ? 采用 MATLAB分析 WAV文件。 ? WAV文件的輸入方式采用 uigetfile 函數(shù)，該函數(shù)打開(kāi)一個(gè)標(biāo)準(zhǔn)的對(duì)話(huà)框，然后可以選擇聲音文件的路徑和文件的名稱(chēng)。 strcmp(action,39。 ? 假定調(diào)制信號(hào)為 m(t)， 其幅度為 1v ， 載波信號(hào)為 ，各種調(diào)幅方式的已調(diào)波的表達(dá)式如下 : t?sin? SAM: ttmmtsA ?s i n)](1[)( ????是調(diào)制度。 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))。張志涌 （ P299） 程序文本。 （ 張志涌 P432） 20 音頻信號(hào)的分析與合成程序 （韓利竹 ， ） 返回 y=sigspec(action)。 [fname,pname]=uigetfile(39。 % 數(shù)據(jù)通過(guò)聲卡轉(zhuǎn)換為聲音 %==========選擇聲道 ========================================== %y=x(:)。)。,10)。按任意鍵聽(tīng)左聲道信號(hào)三次??！ 39。 [k]=fft(y,32768)。39。FontSize39。 [m1,i1]=max(abs(k))。 F2t=[num2str(F2) 39。]。39。FontSize39。 disp(39。按任意鍵顯示左聲道簡(jiǎn)單合成聲音的頻域幅值?。?39。 pause plot(abs(k1))。 title([39。 F1_rt=[num2str(F1_r) 39。]。 if strcmp(action,39。Open Wave File39。 % nbits是每次采樣比特?cái)?shù) file=[file 39。file=[file 39。439。 end。 H_osc=figure(39。,39。,39。39。Callback39。)。,39。 if strcmp(v(1),39。c2=39。b39。 % dt=1/fs， fs是采樣速率（ Hz） if T1e03。 end。,10)。])。 % 返回當(dāng)前軸坐標(biāo)的設(shè)置 del=t(2)t(1)。sound(xx,fs)。Spectrum39。n2=length(x)。 H2=figure(39。,39。,39。39。 N2=fix(N/2)。fn=fn(I)。,39。,39。ylabel(ystr)。 plot(fn,y,c2)。YLim39。 xlabel([39。 %=================================================================================== 30 幅度調(diào)制程序 （ 車(chē)晴 P203 ） 返回 function []=am(action) % AM Modulate THEE 2021/10 % []=am(action)。)。,39。,[1 1 1])。,39。w39。Units39。 H1=uicontrol(39。,39。,... 39。,39。,39。,str2mat(39。),... 39。 H3=uicontrol(39。,39。,39。),... 39。 H41=uicontrol(39。,39。,[.8 .8 .8])。Normal39。+0339。Style39。Position39。 H52=uicontrol(39。,39。,39。,39。,... [.84 .4 .15 .05],39。39。,39。,39。m=1。)。wave39。)。 else。 f=get(gco,39。 fc=get(gco,39。 f=+03。,39。 set(H1,39。,1)。 if strcmp(v(1),39。)。)。 subplot(39。t(39。 %======================================================================================== 34 幅度調(diào)制程序（續(xù) 4） [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,Ks)。Kmod=1。 if nargin1。 if Kw==1。 x=diric(2*pi*f*t,32)。astr=39。I=find(y)。I=find(y)。fstr=39。fstr=39。Kf=1。end。t(39。yy=yy/ymax。position39。)。ylabel(39。 set(gcf,39。])。,39。zoom xon。f(39。end。yn=yy(I)。 n=length(y)。r39。position39。 tstr=39。 t=t*1000。 t=t*1e+06。astr=39。astr=39。astr=39。m=m(m==1)*eps。dt=1/fs。 end。 if nargin5。title(39。plot(t(i),x1(i),39。 tstr 39。 plot(t,y,39。set(gca,39。 x1=demod(y,fc,fs,39。value39。string39。Kmod=1。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 end。mt=[1 .9 .875 .7 .6 .5]。)。value39。method39。39。,... [.84 .3 .15 .05],39。,39。,39。,39。39。,[1 1 1],39。,39。fc(Hz):39。Units39。am f。Back39。Edit39。,39。,39。,39。,39。,str2mat(39。,39。,39。,39。,39。 31 幅度調(diào)制程序（續(xù) 1） H2=uicontrol(39。),... 39。,str2mat(39。,39。Position39。Style39。,39。 figure(39。,39。 figure(39。initialized39。])。Arial39。Ygrid39。)。frequency (39。,[60 0],... % gca是返回當(dāng)前軸的句柄并設(shè)置 39。(dB)39。y(I)=1e3*ones(size(I))。 % N是離散傅氏變換用的點(diǎn)數(shù) % N應(yīng)該選擇 2的乘方的數(shù)（如： 16， 128...） f1=1/(dt*N)。39。Label39。,... 39。)。del=del*.5。 % 放大 zoom(H_osc,39。 if isempty(n2)。 27 音頻頻譜分析儀程序（續(xù) 2） elseif strcmp(action,39。xlabel([39。,39。 elseif T1。t=(0:N1)*dt。)。k39。39。amp。Play39。,39。39。Label39。 blanks(6) file],... 39。Play39。fs=39。 % 若是雙聲道，將數(shù)據(jù)合并為一列 %x=x(:,1)。 int2str(nbits) 39。539。 % 清除工作區(qū) v=version。 % AUDIO SPECTRUM ANALYSER THEE 2021/10 %=================================================================================== global x t dt fs c1 c2 H_osc H2 fn y ystr Ktype if nargin1。 [m2_r,i2_r]=max(abs(k1(10000:12021)))。采樣點(diǎn)數(shù) 39。Arial39。)。 % 數(shù)據(jù)通過(guò)聲卡轉(zhuǎn)換為聲音 pause sound(y1,fs,bits)。])。FontName39。 disp(39。 [m3,i3]=max(ab