2、,zeros(1,t0/(3*ts)+1)];c=cos(2*pi*fc.*t); %载波信号u=m.*c; %调制信号[M,m,df1]=fftseq(m,ts,df); %傅立叶变换M=M/fs; %缩放[U,u,df1]=fftseq(u,ts,df);U=U/fs;[C,c,df1]=fftseq(c,ts,df);f=[0:df1:df1*(length(m)-1)]-fs/2;%频率向量signal_power=power(u(1:length(t)));%调制信号的功率noise_power=signal_power/
3、snr_lin;%计算噪声功率noise_std=sqrt(noise_power);%计算噪声标准偏差noise=noise_std*randn(1,length(u));%产生噪声r=u+noise; %给调制信号加入噪声[R,r,df1]=fftseq(r,ts,df);%信号加噪声的频谱R=R/fs;pause %按任意键可看到调制信号功率signal_powerpause %按任意键可看到消息的曲线clfsubplot(2,2,1)plot(t,m(1:length(t)))xlabel('Time')title('Themessa
4、gesignal')pause %按任意键可看到载波的曲线subplot(2,2,2)plot(t,c(1:length(t)))xlabel('Time')title('Thecarrier')pause按任意键可看到调制信号的曲线subplot(2,2,3)plot(t,u(1:length(t)))xlabel('Time')title('Themodulatedsignal')pause %按任意键可看到消息的幅度曲线和在频域范围内的调制信号subplot(2,1,1)plot(f,abs(fftshift(M)))xlabel('fr
5、equency')title('Spectrumofthemessagesignal')subplot(2,1,2)plot(f,abs(fftshift(U)))xlabel('frequency')title('Spectrumofthemodulatedsignal')pause %按任意键可看到一个噪声抽样subplot(2,1,1)plot(t,noise(1:length(t)))xlabel('Time')title('noisesample')pause %按任意键可看到调制信号和噪声subplot(2,1,2)plot(t,r
6、(1:length(t)))xlabel('Time')title('Signalandnoise')pause%按任意键可看到频域的调制信号和噪声subplot(2,1,1)plot(f,abs(fftshift(U)))title('Signalspectrum')xlabel('frequency')subplot(2,1,2)plot(f,abs(fftshift(R)))title('Signalandnoisespectrum')xlabel('frequency')fftseq.mfunction[M,m,df]=fftseq(m,t
7、s,df)%[M,m,df]=fftseq(m,ts,df)%[M,m,df]=fftseq(m,ts)fs=1/ts;ifnargin==2n1=0;elsen1=fs/df;endn2=length(m);n=2^(max(nextpow2(n1),nextpow2(n2))); M=fft(m,n);m=[m,zeros(1,n-n2)];df=fs/n;===--DSB信号的表达式、频谱及带宽 在幅度调制的一般模型中，若假设滤波器为全通网络（＝1），，调制信号m(t)中无直流分量，则输出的已调信号就是无载波分量的双边带调制信号，或称
8、抑制载波双边带（DSB-SC）调制信号，简称双边带（DSB）信号。 DSB调制器模型如图3-7所示。可见DSB信号实
