基本模擬信號調(diào)制仿真1基本數(shù)字信號調(diào)制仿真2課堂設(shè)計教學(xué)內(nèi)容

掌握基本的模擬信號、數(shù)字信號的調(diào)制仿真

理解調(diào)制過程用MATLAB語言模擬信號的調(diào)制過程重點難點課堂設(shè)計教學(xué)目標(biāo)實驗?zāi)康膶嶒瀮?nèi)容（1）將不同調(diào)制方式映射為IQ正交調(diào)制；（2）領(lǐng)會模擬調(diào)制與數(shù)字調(diào)制的聯(lián)系與差異。（1）實現(xiàn)基本模擬信號的調(diào)制，如AM、FM；（2）實現(xiàn)基本數(shù)字信號的調(diào)制，如QAM。調(diào)制的基本原理調(diào)制的基本原理

極坐標(biāo)調(diào)制方式很難精確控制相位極坐標(biāo)調(diào)制調(diào)制的基本原理調(diào)制的基本原理調(diào)制的基本原理IQ正交調(diào)制模擬信號調(diào)制AM調(diào)制令m(t)為待調(diào)制基帶信號，則AM調(diào)制信號為：clearall;closeall;%%%%AM調(diào)制%%N=4000;%采樣點fs=1000;%采樣頻率t=(0:N-1)/fs;%定義時間區(qū)間(N=4000,fs=1000)f0=1;f1=2;%設(shè)置待調(diào)信號頻率f2=10;%設(shè)置載波信號頻率fai=pi/4;m=cos(2*pi*f0*t)+cos(2*pi*f1*t);%待調(diào)制信號z1=cos(2*pi*f2*t);%設(shè)置載波z1z2=sin(2*pi*f2*t);%設(shè)置載波z2z=cos(2*pi*f2*t+fai);A=2;%%映射m1=(A+m).*cos(fai);%I路基帶信號m2=(A+m).*sin(fai);%Q路基帶信號%%正交調(diào)制s=m1.*z1-m2.*z2;%已調(diào)制信號%%頻域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('待調(diào)制信號');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('已調(diào)制信號');xlabel('t');ylabel('幅度');subplot(3,2,6);plot(f,fftshift(affts));title('AM已調(diào)信號頻譜');xlabel('f');ylabel('幅頻特性');FM調(diào)制令m(t)為待調(diào)制基帶信號，則FM調(diào)制信號為：模擬信號調(diào)制參數(shù)準(zhǔn)備同AM調(diào)制%求積分m0(1)=0;dt=1/fs;forii=1:length(t)-1m0(ii+1)=m0(ii)+m(ii)*dt;endz1=cos(2*pi*f2*t);%設(shè)置載波z1z2=sin(2*pi*f2*t);%設(shè)置載波z2A=2;%設(shè)置振幅kf=10;%設(shè)置頻偏%%映射m1=A*cos(kf*m0);%I路基帶信號m2=A*sin(kf*m0);%Q路基帶信號s=m1.*z1-m2.*z2;%已調(diào)制信號%%頻域df=fs/length(t);ff=fftshift(fft(s));ww=linspace(-fs/2,fs/2,N);figure(1)subplot(3,2,1);plot(m);title('待調(diào)制信號');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('已調(diào)制信號');xlabel('t');ylabel('幅度');subplot(326);plot(ww,0.5*log10(abs(ff)));axis([0,100,-7.5,5]);title('頻譜');xlabel('頻率');ylabel('幅度');QAM調(diào)制當(dāng)Ak和φk有多種值時，根據(jù)信號星座點的數(shù)目M，可稱為MQAM數(shù)字信號調(diào)制clearall;closeall;%%%QAM調(diào)制%%N=4000;%采樣點fs=1000;%采樣頻率t=(0:N-1)/fs;%定義時間區(qū)間f2=10;%設(shè)置載頻m0=[73912151062115];%映射關(guān)系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);%設(shè)置載波z1z2=sin(2*pi*f2*t);%設(shè)置載波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路信號%已調(diào)制信號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('待調(diào)制信號');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('幅