基于MATLAB的抑制载波双边带调频仿真.docVIP

基于MATLAB的抑制载波双边带调频仿真 图1 抑制载波双边带调幅波形 三、仿真实现??? 未调制信号为S（t）= sinc（200t）（|t|=t0），其中t0取2s，载波频率fc = 100Hz。用抑制载波调幅来调制信号，给出调制信号M（t）波形，画出S（t）和M（t）的频谱。??? 环境为Matlab 7.6（R2008a），安装后占3.6GB硬盘空间。安装完之后，可以看到ID号，为破解版。 图3 About Matlab ??? 编写了脚本文件dsbmod.m来完成带限信号的抑制载波双边带调频的例子。 ? 图6 Matlab环境下的脚本编写 ??? 脚本关键代码： ?%模拟线性调制 - 抑制载波双边带条幅（DSB-SC） - 天水浪客clearecho ont0 = 2; %信号持续时间ts = 0.001; %采样时间间隔fc = 100; %载波频率fs = 1/ts;df = 0.3; %频率分辨力t = [-t0/2:ts:t0/2]; %定义时间序列x = sin(200*t);m = x./(200*t);m(1001) = 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; %频率压缩f = [0:df1:df1*(length(m)-1)]-fs/2;pauseclfsubplot(2,2,1)plot(t,m(1:length(t)))axis([-0.4,0.4,-0.5,1.1])xlabel(时间);title(未调信号)pause;subplot(2,2,3)plot(t,u(1:length(t)))axis([-0.1,0.1,-1.5,1.5]);xlabel(时间)title(载波);pausesubplot(2,2,2)plot(t,u(1:length(t)))axis([-0.2,0.2,-1,1.2]);xlabel(时间)title(已调信号);pausesubplot(2,1,1)plot(f,abs(fftshift(M)))xlabel(频率);title(未调信号的频谱)subplot(2,1,2)plot(f,abs(fftshift(U)))title(已调信号的频谱);xlabel(频率) ?%傅立叶变换函数 - 天水浪客function[M,m,df] = fftseq(m,tz,df)fz = 1/tz;if nargin == 2??? n1 = 0;else n1 = fz/df;endn2 = length(m);n = 2^(max(nextpow2(n1),nextpow2(n2)));M = fft(m,n);m = [m,zeros(1,n-n2)]; 图7抑制载波调幅频谱 图8 抑制载波调幅频谱（续） ??? 本来sinc函数的频谱应该是一个矩形波谱，但是由于S（t）只是sinc函数的一段，并且在计算机上使用离散的数字来存储，因为计算精度的要求，是频谱结果与矩形波谱不完全一致，但是差强人意。因此，本实验的仿真还是成功的。 五、参考资料1、《通信原理》，孙学军、王秉军编著，电子工业出版社2、《Matlab通信仿真与应用》，刘敏、魏玲编著，国防工业出版社