MATLAB在通信仿真中应用.pptVIP

第八讲 MATLAB在通信仿真中的应用;一、模拟调制;在DSB-AM中，已调信号的时域表示如下式：;;t0=2; ts=0.001; fs=1/ts Fc=250; t=-t0/2:ts:t0/2; m=sin(200*t)./(200*t); m(1001)=1 subplot(2,2,1) plot(t,m) c=cos(2*pi*Fc*t); u=m.*c; subplot(2,2,2) plot(t,u,r) F_m=ts*fft(m); F_u=ts*fft(u); f=linspace(-fs/2,fs/2,length(m)) subplot(2,2,3) plot(f, abs(fftshift(F_m))) subplot(2,2,4) plot(f, abs(fftshift(F_u)));DSB-AM解调;调制与解调源代码;y=u.*c; F_y=ts*fft(y); subplot(2,3,3) plot(f, abs(F_y)) H=zeros(size(f)); f_cutoff=100; l=fs/length(m); n=floor(f_cutoff/l); H(1:n)=2*ones(1,n); l1=length(m)-n; H(length(m)-n+1:length(m))=2*ones(1,n); Z=F_y.*H; y=fs*real(ifft(Z)); subplot(2,3,6) plot(t,y) ;去掉DSB-AM的一边就得到SSB-AM。依据所保留的边带是上边，还是下边，可以分为USSA和LSSB两种不同的方式，此时信号的时域表示为：;;t0=2; ts=0.001; fs=1/ts Fc=250; t=-t0/2:ts:t0/2; m=sin(200*t)./(200*t); m(1001)=1; m=m.*m; subplot(2,3,1); plot(t,m); axis([-0.2 0.2 -1.2 1.2]); c=cos(2*pi*Fc.*t); d=sin(2*pi*Fc.*t); u=(m.*c-imag(hilbert(m)).*d)/2; v=(m.*c+imag(hilbert(m)).*d)/2; subplot(2,3,2); plot(t,u,r); axis([-0.1 0.1 -1.2 1.2]); subplot(2,3,3); plot(t,v); axis([-0.1 0.1 -1.2 1.2]); F_m=ts*fft(m); F_u=ts*fft(u); F_v=ts*fft(v); f=linspace(-fs/2,fs/2,length(m)); subplot(2,3,4); plot(f,abs(fftshift(F_m))); subplot(2,3,5); plot(f,abs(fftshift(F_u))); subplot(2,3,6); plot(f,abs(fftshift(F_v))); ;SSB-AM解调;SSB-AM源程序;这里a是调制指数， ;;AM调制与解调源程序;角度调制;其中u(t)是调制后的信号，fc是载波的频率(单位赫兹Hz)， 是初始相位， 是瞬时相位， 随着输入信号的振幅变化。;示例;;调频源程序;其中u(t)是调制后的信号，fc是载波的频率(单位赫兹Hz)， 是初始相位， 是瞬时相位，随着输入信号的振幅变化 ;;模拟调制总结;(1) 带通模拟调制函数amod( ) 格式：y=amod(x,Fc,Fs,method…) 功能：用载波为Fc(Hz)的信号来调制模拟信号x，采样频率为Fs(Hz)，Fc Fs。根据采样定理，采样频率必须大于或等于调制信号最高频率的两倍。字符串变量method指定所用的调制方式;带通模拟解调函数ademod( );Method:调制方式;1．基带模拟调制函数amodce( ) 格式：y=amodce(x,Fs,method,···) 功能：对输入信号x进行调制，输出复包络信号。输入输出的采样频率为Fs(Hz)。输出y为一个复矩阵。Fs可以是标量也可以为二维的矢量。二维矢量中第一个值为采样频率，第二个值为载波信号的初相，初相以弧度表示，默认值为0。字符串变量method指定所用的调制方式;2．基带模拟解调函数ademodce( );[例] 利用MATLAB对一信号进行基带、通带调制解调。 抽样频率：Fs=1000; 载波频率： Fc=200; 信号：x=sinc(10*t); 区间：[-1，1];PC