单边带调制MATLAB源代码

%Am调制的Matlab演示源程序。被调信号 %当0<t<t0/3时为+1，当t0/3<t<2t0/3为-2，其它情况下为零 echo on t0= .15; %信号的持续时间 ts = 0.001; %抽样间隔 fc = 250; %载波频率 snr = 10; %用dB表示信噪比 fs = 1/ts; %抽样频率 t = [0:ts:t0]; %时间向量 df = 0.2; %所需的频率分辨率 Lt=length(t); snr_lin = 10^(snr/10); %信噪比 %被调信号 m = [ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)]; L=2*min(m); R=2*max(abs(m)); pause clf figure(1) subplot(4,2,1); plot(t,m(1:length(t))); axis([0 t0 -R/2 R/2]); xlabel('t'); ylabel('调制信号'); pause %按任意键可看到调制信号的曲线 subplot(4,2,2); [M,m,df1,f]=T2F(m,ts,df,fs); [Bw_eq]=signalband(M,df,t0); f_start=fc-Bw_eq; f_cutoff=fc+Bw_eq; plot(f,fftshift(abs(M))); xlabel('f'); ylabel('调制信号频谱'); pause subplot(4,2,3); c1=cos(2*pi*fc*t); plot(t,c1); axis([0 t0 -1.2 1.2]); xlabel('t'); ylabel('余弦载波'); pause subplot(4,2,4); c2=sin(2*pi*fc*t); plot(t,c2); axis([0 t0 -1.2 1.2]); xlabel('t'); ylabel('正弦载波'); pause subplot(4,2,5); u1=m(1:Lt).*c1(1:Lt)+m(1:Lt).*c2(1:Lt); plot(t,u1); axis([0 t0 -R R]); xlabel('t'); ylabel('单边带信号-下边带'); subplot(4,2,6); [U1,u1,df1,f]=T2F(u1,ts,df,fs); plot(f,fftshift(abs(U1))); xlabel('f'); ylabel('单边带信号-下边频谱'); pause subplot(4,2,7); u2=m(1:Lt).*c1(1:Lt)-m(1:Lt).*c2(1:Lt); plot(t,u2); axis([0 t0 -R R]); xlabel('t'); ylabel('单边带信号-上边带'); subplot(4,2,8); [U2,u2,df1,f]=T2F(u2,ts,df,fs); plot(f,fftshift(abs(U2))); xlabel('f'); ylabel('单边带信号-上边频谱'); pause