【正文】 % v(1)和 v(2)是 X軸坐標(biāo)的設(shè)置 n1=find(abs(tv(1))del)。 if isempty(n2)。end。sound(xx,fs)。Zoom39。 % 放大 zoom(H_osc,39。)。Spectrum39。 % 頻譜分析儀 v=axis。del=del*.5。n2=find(abs(tv(2))del)。n2=length(x)。 28 音頻頻譜分析儀程序（續(xù) 3） % 增加菜單項(xiàng)（ Log）， 生成信號(hào)頻譜分析圖形，回應(yīng)（ Callback） 調(diào)用函數(shù) audspec(39。)。delete(H2)。 H2=figure(39。,39。,... 39。,39。,39。,[155 32 640 480])。Label39。amp。,39。,39。39。39。39。 xx=x(n1:n2)。 % N是離散傅氏變換用的點(diǎn)數(shù) % N應(yīng)該選擇 2的乘方的數(shù)（如： 16， 128...） f1=1/(dt*N)。 變換點(diǎn)的實(shí)際頻率是：（變換點(diǎn)的位置序號(hào) /N） *fs y=fft(xx/N,32768)。 N2=fix(N/2)。y=y/max(y)。y(I)=1e3*ones(size(I))。 29 音頻頻譜分析儀程序（續(xù) 4） fn=f1*(0:N21)。fn=fn(I)。 plot(fn,y,c2)。(dB)39。 set(gca,39。,39。,39。,[60 0],... % gca是返回當(dāng)前軸的句柄并設(shè)置 39。,39。,39。,10)。frequency (39。Hz)39。ylabel(ystr)。 elseif strcmp(action,39。)。 semilogx(fn,y,c2)。 plot(fn,y,c2)。 Ktype=Ktype。Ygrid39。on39。YLim39。FontName39。Arial39。FontSize39。 xlabel([39。 39。])。zoom xon。 %=================================================================================== 30 幅度調(diào)制程序 （ 車晴 P203 ） 返回 function []=am(action) % AM Modulate THEE 2021/10 % []=am(action)。if nargin1。initialized39。 if strcmp(action,39。)。539。 figure(39。,39。,39。,39。,39。,... [5 29 792 530],39。,[1 1 1])。439。 figure(39。,39。,39。,39。,39。,[5 29 792 530])。w39。 end。Style39。Frame39。Units39。Normal39。Position39。Back39。 H1=uicontrol(39。,39。,39。,39。,39。,... [.84 .9 .15 .05],39。,str2mat(39。,39。,... 39。,39。),... 39。,[1 1 1],39。,39。39。 31 幅度調(diào)制程序（續(xù) 1） H2=uicontrol(39。,39。,39。,39。,39。,... [.84 .8 .15 .05],39。,str2mat(39。,... 39。,39。,39。),... 39。,[1 1 1],39。,39。39。 H3=uicontrol(39。,39。,39。,39。,39。,... [.84 .7 .15 .05],39。,str2mat(39。,... 39。,39。,39。,39。,39。),... 39。,[1 1 1],39。,39。39。 H41=uicontrol(39。,39。,39。,39。,39。,... [.83 .595 .05 .04],39。,39。,39。,[.8 .8 .8])。Style39。Edit39。Units39。Normal39。Position39。Back39。String39。+0339。Call39。am f。)。Style39。Text39。Units39。Normal39。Position39。String39。fc(Hz):39。Back39。 H52=uicontrol(39。,39。,39。,39。,39。,... [.89 .5 .1 .04],39。,[1 1 1],39。,39。,39。,39。39。 H6=uicontrol(39。,39。,39。,39。,39。,... [.84 .4 .15 .05],39。,39。,39。,39。39。 H7=uicontrol(39。,39。,39。,39。,39。,... [.84 .3 .15 .05],39。,39。,39。,39。39。 32 幅度調(diào)制程序（續(xù) 2） f=+03。m=1。Kmod=1。method39。 Kmod=get(gco,39。)。 if Kmod1。value39。end。wave39。 Kw=get(gco,39。)。 elseif strcmp(action,39。)。Value39。mt=[1 .9 .875 .7 .6 .5]。 if Kmod==1。 else。value39。 end。f39。 f=get(gco,39。)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。fc39。 fc=get(gco,39。)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。restore39。 f=+03。m=1。Kmod=1。 set(H42,39。,39。)。string39。38e+0639。 set(H1,39。,1)。value39。set(H3,39。,1)。demod39。 x1=demod(y,fc,fs,39。)。 if strcmp(v(1),39。)。set(gca,39。,39。)。position39。 plot(t,y,39。,t,x,39。)。xlabel([39。 tstr 39。])。 subplot(39。,[.065 .11 .7 .3439])。plot(t(i),x1(i),39。)。t(39。)39。title(39。)。 %======================================================================================== 34 幅度調(diào)制程序（續(xù) 4） [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,Ks)。Ks=0。 if nargin5。 end。Kmod=1。 if nargin3。 end。fc=38e+06。 if nargin1。end。dt=1/fs。t=0:dt:Tdt。 if Kw==1。 elseif Kw==2。m=m(m==1)*eps。 x=sawtooth(2*pi*f*t,.98)。 x=diric(2*pi*f*t,32)。 35 幅度調(diào)制程序（續(xù) 5） if Kmod==1。astr=39。 elseif Kmod==2。astr=39。 elseif Kmod==3。astr=39。 I=find(y)。I=find(y)。 elseif Kmod==4。astr=39。 I=find(y)。I=find(y)。 end。 t=t*1e+06。tstr=39。fstr=39。Kf=1e06。 t=t*1000。tstr=39。fstr=39。Kf=1e03。 tstr=39。fstr=39。Kf=1。 36 幅度調(diào)制程序（續(xù) 6） if Ks==1。position39。 else。end。b39。r39。zoom xon。t(39。)39。 n=length(y)。yy=yy(1:fix(n/2))。yy=yy/ymax。 I=find(yy40)。yn=yy(I)。 if Ks==1。position39。 else。end。 plot([0 0] +fn(i),[40 yn(i)],39。)。 end。f(39。)39。ylabel(39。)。zoom xon。axis([v(1)dv v(2)+dv v(3:4)])。 set(gcf,39。,39。,39。,[astr blanks(10) 39