【正文】 ntrol(39。Style39。,39。Edit39。,39。Units39。,39。Normal39。,39。Position39。,... [.89 .5 .1 .04],39。Back39。,[1 1 1],39。String39。,39。38e+0639。,39。Call39。,39。am fc。39。)。 H6=uicontrol(39。Style39。,39。Push39。,39。Units39。,39。Normal39。,39。Position39。,... [.84 .4 .15 .05],39。String39。,39。RESTORE39。,39。Call39。,39。am restore。39。)。 H7=uicontrol(39。Style39。,39。Push39。,39。Units39。,39。Normal39。,39。Position39。,... [.84 .3 .15 .05],39。String39。,39。DEMOD39。,39。Call39。,39。am demod。39。)。 32 幅度調(diào)制程序（續(xù) 2） f=+03。fc=38e+06。m=1。Kw=1。Kmod=1。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1) elseif strcmp(action,39。method39。)。 Kmod=get(gco,39。Value39。)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 if Kmod1。set(H3,39。value39。,1)。end。 elseif strcmp(action,39。wave39。)。 Kw=get(gco,39。Value39。)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 elseif strcmp(action,39。mr39。)。 k=get(gco,39。Value39。)。mt=[1 .9 .875 .7 .6 .5]。m=mt(k)。 if Kmod==1。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 else。 set(H3,39。value39。,1)。 end。 elseif strcmp(action,39。f39。)。 f=get(gco,39。string39。)。f=str2num(f)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 33 幅度調(diào)制程序（續(xù) 3） elseif strcmp(action,39。fc39。)。 fc=get(gco,39。string39。)。fc=str2num(fc)。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 elseif strcmp(action,39。restore39。)。 f=+03。fc=38e+06。m=1。Kw=1。Kmod=1。 [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,1)。 set(H42,39。string39。,39。+0339。)。set(H52,39。string39。,39。38e+0639。)。 set(H1,39。value39。,1)。set(H2,39。value39。,1)。set(H3,39。value39。,1)。 elseif strcmp(action,39。demod39。)。 x1=demod(y,fc,fs,39。am39。)。x1=x1mean(x1)。 if strcmp(v(1),39。439。)。 subplot(1,1,1)。set(gca,39。vis39。,39。off39。)。 subplot(39。position39。,[.065 .5811 .7 .3439])。 plot(t,y,39。b39。,t,x,39。r39。)。zoom xon。xlabel([39。t(39。 tstr 39。)39。])。 end。 subplot(39。position39。,[.065 .11 .7 .3439])。 i=1:length(t)8。plot(t(i),x1(i),39。b39。)。 xlabel([39。t(39。 tstr 39。)39。])。title(39。Demodulated Signal39。)。 end。 %======================================================================================== 34 幅度調(diào)制程序（續(xù) 4） [y,x,t,fc,fs,tstr]=am1(f,fc,m,Kmod,Kw,Ks)。 % % THEE 2021/10 %======================================================================================== if nargin6。Ks=0。 end。 if nargin5。Kw=1。 end。 if nargin4。Kmod=1。 end。 if nargin3。m=1。 end。 if nargin2。fc=38e+06。 end。 if nargin1。f=+03。end。 fs=4*fc。dt=1/fs。T=2/f。t=0:dt:Tdt。f1=1/T。 if Kw==1。 x=sin(2*pi*f*tpi/2)。 elseif Kw==2。 x=square(2*pi*f*t,10)。m=m(m==1)*eps。 elseif Kw==3。 x=sawtooth(2*pi*f*t,.98)。 elseif Kw==4。 x=diric(2*pi*f*t,32)。 end。 35 幅度調(diào)制程序（續(xù) 5） if Kmod==1。 y=(1+m*x).*cos(2*pi*fc*t)。astr=39。AM,39。 elseif Kmod==2。 y=x.*cos(2*pi*fc*t)。astr=39。AMDSB,39。 elseif Kmod==3。 y=x.*cos(2*pi*fc*t)+imag((hilbert(x))).*sin(2*pi*fc*t)。astr=39。AMSSB lower band,39。 I=find(y)。y(I)=ones(size(I))。I=find(y)。 y(I)=ones(size(I))。 elseif Kmod==4。 y=x.*cos(2*pi*fc*t)imag((hilbert(x))).*sin(2*pi*fc*t)。astr=39。AMSSB upper band,39。 I=find(y)。y(I)=ones(size(I))。I=find(y)。 y(I)=ones(size(I))。 end。 if fc=1e+06。 t=t*1e+06。f1=f1/1e+06。tstr=39。us39。fstr=39。MHz39。Kf=1e06。 elseif fc=1000。 t=t*1000。f1=f1/1000。tstr=39。ms39。fstr=39。kHz39。Kf=1e03。 else。 tstr=39。s39。fstr=39。Hz39。Kf=1。 end。 36 幅度調(diào)制程序（續(xù) 6） if Ks==1。subplot(39。position39。,[.065 .5811 .7 .3439])。 else。subplot(211)。end。 plot(t,y,39。b39。,t,x,39。r39。)。zoom xon。xlabel([39。t(39。 tstr 39。)39。])。 n=length(y)。yy=abs(fft(y)/n)。yy=yy(1:fix(n/2))。ymax=max(yy)。yy=yy/ymax。yy=20*log10(yy)。 I=find(yy40)。fn=(I1)*f1。yn=yy(I)。m=length(I)。 if Ks==1。subplot(39。position39。,[.065 .11 .7 .3439])。 else。subplot(212)。end。 for i=1:m。 plot([0 0] +fn(i),[40 yn(i)],39。b39。)。hold on。 end。 xlabel([39。f(39。 fstr 39。)39。])。ylabel(39。dB39。)。v=axis。zoom xon。 dv=(v(2)v(1))*.05。axis([v(1)dv v(2)+dv v(3:4)])。hold off。 set(gcf,39。num39。,39。off39。,39。Name39。,[astr blanks(10) 39。THEE 2021/839。])。 %========================================================================================