《软件无线电原理与技术》课件-18 调制技术仿真

基本模拟信号调制仿真1基本数字信号调制仿真2课堂设计教学内容

掌握基本的模拟信号、数字信号的调制仿真

理解调制过程用MATLAB语言模拟信号的调制过程重点难点课堂设计教学目标实验目的实验内容（1）将不同调制方式映射为IQ正交调制；（2）领会模拟调制与数字调制的联系与差异。（1）实现基本模拟信号的调制，如AM、FM；（2）实现基本数字信号的调制，如QAM。调制的基本原理调制的基本原理

极坐标调制方式很难精确控制相位极坐标调制调制的基本原理调制的基本原理调制的基本原理IQ正交调制模拟信号调制AM调制令m(t)为待调制基带信号，则AM调制信号为：clearall;closeall;%%%%AM调制%%N=4000;%采样点fs=1000;%采样频率t=(0:N-1)/fs;%定义时间区间(N=4000,fs=1000)f0=1;f1=2;%设置待调信号频率f2=10;%设置载波信号频率fai=pi/4;m=cos(2*pi*f0*t)+cos(2*pi*f1*t);%待调制信号z1=cos(2*pi*f2*t);%设置载波z1z2=sin(2*pi*f2*t);%设置载波z2z=cos(2*pi*f2*t+fai);A=2;%%映射m1=(A+m).*cos(fai);%I路基带信号m2=(A+m).*sin(fai);%Q路基带信号%%正交调制s=m1.*z1-m2.*z2;%已调制信号%%频域df=fs/length(t);f=-fs/2:df:fs/2-df;ffts=fft(s);affts=abs(ffts);%%作图figure(1)subplot(3,2,1);plot(m);title('待调制信号');xlabel('t');ylabel('幅度');subplot(3,2,2);plot(z);title('载波信号');xlabel('t');ylabel('幅度');subplot(3,2,3);plot(m1);title('I路基带信号');xlabel('t');ylabel('幅度');subplot(3,2,4);plot(m2);title('Q路基带信号');xlabel('t');ylabel('幅度');subplot(3,2,5);plot(s);title('已调制信号');xlabel('t');ylabel('幅度');subplot(3,2,6);plot(f,fftshift(affts));title('AM已调信号频谱');xlabel('f');ylabel('幅频特性');FM调制令m(t)为待调制基带信号，则FM调制信号为：模拟信号调制参数准备同AM调制%求积分m0(1)=0;dt=1/fs;forii=1:length(t)-1m0(ii+1)=m0(ii)+m(ii)*dt;endz1=cos(2*pi*f2*t);%设置载波z1z2=sin(2*pi*f2*t);%设置载波z2A=2;%设置振幅kf=10;%设置频偏%%映射m1=A*cos(kf*m0);%I路基带信号m2=A*sin(kf*m0);%Q路基带信号s=m1.*z1-m2.*z2;%已调制信号%%频域df=fs/length(t);ff=fftshift(fft(s));ww=linspace(-fs/2,fs/2,N);figure(1)subplot(3,2,1);plot(m);title('待调制信号');xlabel('t');ylabel('幅度');subplot(3,2,2);plot(z1);title('载波信号');xlabel('t');ylabel('幅度');subplot(3,2,3);plot(m1);title('I路基带信号');xlabel('t');ylabel('幅度');subplot(3,2,4);plot(m2);title('Q路基带信号');xlabel('t');ylabel('幅度');subplot(3,2,5);plot(s);title('已调制信号');xlabel('t');ylabel('幅度');subplot(326);plot(ww,0.5*log10(abs(ff)));axis([0,100,-7.5,5]);title('频谱');xlabel('频率');ylabel('幅度');QAM调制当Ak和φk有多种值时，根据信号星座点的数目M，可称为MQAM数字信号调制clearall;closeall;%%%QAM调制%%N=4000;%采样点fs=1000;%采样频率t=(0:N-1)/fs;%定义时间区间f2=10;%设置载频m0=[73912151062115];%映射关系0[-3-3],1[-3-1],2[-3,1],3[-3,3],4[-1,-3],5[-1,-1],6[-1,1],%7[-1,3],8[1,-3],9[1,-1],10[1,1],11[1,3],12[3,-3],13[3,-1],14[3,1],15[3,3]m1=[-1-31331-1-3-33];%串并变换m2=[33-1-33111-13];z1=cos(2*pi*f2*t);%设置载波z1z2=sin(2*pi*f2*t);%设置载波z2n1=1;n2=1;%将输入信号加入到方波中whilen1<11M0(n2:n2+399)=m0(n1);M1(n2:n2+399)=m1(n1);M2(n2:n2+399)=m2(n1);n2=n2+400;n1=n1+1;endI=M1;%I路信号Q=M2;%Q路信号%已调制信号s=I.*z1-Q.*z2;%%频域ff=fftshift(fft(s));ww=linspace(-fs/2,fs/2,N);%%作图figure(1)subplot(3,2,1);plot(M0);title('待调制信号');xlabel('t');ylabel('幅度');axis([04000-0.515.5]);subplot(3,2,2);plot(z1);title('载波信号');xlabel('t');ylabel('幅度');subplot(3,2,3);plot(I);title('I路基带信号');xlabel('t');ylabel('幅度');axis([04000-3.53.5]);subplot(3,2,4);plot(Q);title('Q路基带信号');xlabel('t');ylabel('幅