wDSB调制与解调的MATLAB实现及调制性能分析.docVIP

%DSB调制的Matlab演示源程序。被调信号%当0tt0/3时为+1，当t0/3t2t0/3时为-2，其它情况下为零。echo ont0=.15;? ?%信号的持续时间ts=.001;??%抽样间隔fc=250;? ?%载波频率snr=20;? ?%用dB（对数的）表示的信噪比fs=1/ts;??%抽样频率df=0.3;? ?%所需的频率分辨率t=[0:ts:t0];? ?? ? %时间向量snr_lin=10^(snr/10);? ?%线性信噪比%message signalm=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),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/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(The message signal)pause??%按任意键可看到载波的曲线subplot(2,2,2)plot(t,c(1:length(t)))xlabel(Time)title(The carrier)pause 按任意键可看到调制信号的曲线subplot(2,2,3)plot(t,u(1:length(t)))xlabel(Time)title(The modulated signal)pause??%按任意键可看到消息的幅度曲线和在频域范围内的调制信号subplot(2,1,1)plot(f,abs(fftshift(M)))xlabel(frequency)title(Spectrum of the message signal)subplot(2,1,2)plot(f,abs(fftshift(U)))xlabel(frequency)title(Spectrum of the modulated signal)pause??%按任意键可看到一个噪声抽样subplot(2,1,1)plot(t,noise(1:length(t)))xlabel(Time)title(noise sample)pause??%按任意键可看到调制信号和噪声subplot(2,1,2)plot(t,r(1:length(t)))xlabel(Time)title(Signal and noise)pause %按任意键可看到频域的调制信号和噪声subplot(2,1,1)plot(f,abs(fftshift(U)))title(Signal spectrum)xlabel(frequency)subplot(2,1,2)plot(f,abs(fftshift(R)))title(Signal and noise spectrum)xlabel(frequency)fftseq.mfunction [M,m,df]=fftseq(m,ts,df)%[M,m,df]=fftseq(m,ts,df) %[M,m,df]=fftseq(m,ts) fs=1/ts; if nargin == 2 n1=0; else n1=fs/df; end n2=length(m); n=2^(max(nextpow2(n1),nextpow2(n2)));??M=fft(m,n); m=[m,zeros(1,n-n2)]; df=fs/n;===-- DSB信号的表达式、频谱及带宽????在幅度调制的一般模型中，若假设滤波器为全通网络（＝1），，