通信原理實(shí)驗(yàn)報(bào)告 實(shí)驗(yàn)一二

1、通信原理實(shí)驗(yàn)一實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)日期：2014年12月11日 學(xué) 院：信息工程學(xué)院 班 級(jí)：12級(jí)電子信息工程二班 學(xué) 號(hào)：2012600119 姓 名：張正潔 指導(dǎo)老師：彭思齊 實(shí)驗(yàn)一 數(shù)字基帶傳輸系統(tǒng)的MATLAB仿真一、實(shí)驗(yàn)?zāi)康?1、熟悉和掌握常用的用于通信原理時(shí)域仿真分析的MATLAB函數(shù)； 2、掌握連續(xù)時(shí)間和離散時(shí)間信號(hào)的MATLAB產(chǎn)生； 3、牢固掌握沖激函數(shù)和階躍函數(shù)等函數(shù)的概念，掌握卷積表達(dá)式及其物理意義，掌握卷積的計(jì)算方法、卷積的基本性質(zhì)； 4、掌握利用MATLAB計(jì)算卷積的編程方法，并利用所編寫(xiě)的MATLAB程序驗(yàn)證卷積的常用基本性質(zhì)； 5、掌握MATLAB描述通信系統(tǒng)中不同波形

2、的常用方法及有關(guān)函數(shù)，并學(xué)會(huì)利用MATLAB求解系統(tǒng)功率譜，繪制相應(yīng)曲線。 基本要求：掌握用MATLAB描述連續(xù)時(shí)間信號(hào)和離散時(shí)間信號(hào)的方法，能夠編寫(xiě)MATLAB程序，實(shí)現(xiàn)各種常用信號(hào)的MATLAB實(shí)現(xiàn)，并且以圖形的方式再現(xiàn)各種信號(hào)的波形。 二、實(shí)驗(yàn)內(nèi)容 1、編寫(xiě) MATLAB 程序產(chǎn)生離散隨機(jī)信號(hào) 2、編寫(xiě) MATLAB 程序生成連續(xù)時(shí)間信號(hào) 3、編寫(xiě) MATLAB 程序?qū)崿F(xiàn)常見(jiàn)特殊信號(hào)三、實(shí)驗(yàn)原理 從通信的角度來(lái)看，通信的過(guò)程就是消息的交換和傳遞的過(guò)程。而從數(shù)學(xué)的角度來(lái)看，信息從一地傳送到另一地的整個(gè)過(guò)程或者各個(gè)環(huán)節(jié)不外乎是一些碼或信號(hào)的交換過(guò)程。例如信源壓縮編碼、糾錯(cuò)編碼、AMI編碼、擾

3、碼等屬于碼層次上的變換，而基帶成形、濾波、調(diào)制等則是信號(hào)層次上的處理。碼的變換是易于用軟件來(lái)仿真的。要仿真信號(hào)的變換，必須解決信號(hào)與信號(hào)系統(tǒng)在軟件中表示的問(wèn)題。四、實(shí)驗(yàn)步驟 （1）分析程序 program1_1 每條指令的作用，運(yùn)行該程序，將結(jié)果保存，貼在下面的空白處。然后修改程序，將 dt 改為 0.2，并執(zhí)行修改后的程序，保存圖形，看看所得圖形的效果怎樣。 dt=0.01 時(shí)的信號(hào)波形 dt=0.2 時(shí)的信號(hào)波形 請(qǐng)問(wèn)：上述的兩幅圖形有什么區(qū)別，哪一副圖形看起來(lái)更接近于實(shí)際信號(hào)波形？為什么會(huì)有這種區(qū)別？ 答： 取點(diǎn)變少，間距變大，發(fā)生失真（2）修改program1_1,，存盤(pán)程序名為Q1_

4、2，生成實(shí)指數(shù)信號(hào)x(t)= exp(-2*t)。 要求在圖形中加上網(wǎng)格線，并使用函數(shù)axis()控制圖形的時(shí)間范圍在02秒之間。然后執(zhí)行該程序，保存所的圖形。 修改Program1_1后得到的程序Q1_2如下： %program1_1% This program is used to generate a sinusoidal signal and draw its plot clear, % Clear all variables close all, % Close all figure windows dt = 0.01; % Specify the step of time vari

5、able t = -2:dt:2;% Specify the interval of time x = exp(-2*t); % Generate the signal plot(t,x) grid on,axis(0,2,-0.5,1.5)% Open a figure window and draw the plot of x(t) title() xlabel(Time t (sec) 圖形結(jié)果如下：(3)將前文中所給的單位沖激信號(hào)和單位階躍信號(hào)的函數(shù)文件在 MATLAB 文件編輯器中編寫(xiě)好，并分別以文件名 delta和u 存入 work文件夾中以便于使用。 抄寫(xiě)函數(shù)文件 delta 如

6、下： % delta function function y = delta(t) dt = 0.01; y = (u(t)-u(t-dt)/dt; 抄寫(xiě)函數(shù)文件u 如下：% Unit step function function y = u(t) y = (t=0); % y = 1 for t 0, else y = 0 (4) 修改程序Program1_4，并以Q1_4為文件名存盤(pán)，利用axis()函數(shù)，將圖形窗口的橫坐標(biāo)范圍改為-2n5，縱坐標(biāo)范圍改為-1.5 x 1.5。 修改Program1_4后得到的程序Q1_4如下： % Program1_4 % This program is

7、 used to generate a discrete-time sinusoidal signal % and draw its plot clear, % Clear all variables close all, % Close all figure windows n = -5:5; % Specify the interval of time x = zeros(1,4), 0.1, 1.1, -1.2, 0, 1.3, zeros(1,2); % Generate the sequence stem(n,x,filled,r) % Open a figure window an

8、d draw the plot of xntitle() grid on,axis(-2,5,-1.5,1.5)信號(hào)的波形圖 (5)根據(jù)示例程序的編寫(xiě)方法，編寫(xiě)一個(gè)MATLAB程序，以Q1_5文件名存盤(pán)，給給定信號(hào) x = exp(-0.5*t).*u(t);求信號(hào)y(t)=x(1.5t+3)，并繪制出x(t)和y(t)的圖形。 編寫(xiě)的程序 Q1_5 如下： clear,close all, t = -5:0.01:5; x = exp(-0.5*t).*u(t); y=exp(-0.5*(1.5*t+3).*u(1.5*t+3); subplot(211) plot(t,x) grid on

9、,title (x(t)subplot (212) plot (t,y) grid on,title ( y(t)xlabel (Q1_5)x(t)和 y(t)的圖形分別為：通信原理實(shí)驗(yàn)二實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)日期：2014年12月11日 學(xué) 院：信息工程學(xué)院 班 級(jí)：12級(jí)電子信息工程二班 學(xué) 號(hào)：2012600119 姓 名：張正潔 指導(dǎo)老師：彭思齊 實(shí)驗(yàn)二 模擬信號(hào)幅度調(diào)制仿真實(shí)驗(yàn)一、實(shí)驗(yàn)?zāi)康?1. 加深對(duì)模擬線性調(diào)制的原理與實(shí)現(xiàn)方法的理解； 2. 掌握 AM、DSB、SSB 功率譜密度函數(shù)的特點(diǎn)，并進(jìn)行對(duì)比； 3. 掌握 MATLAB基本指令的使用； 4. 掌握 MATLAB中 M 文件的調(diào)試以

10、及子函數(shù)調(diào)用的方法。 二、實(shí)驗(yàn)內(nèi)容 1. 復(fù)習(xí) AM、DSB 和SSB 調(diào)制的相關(guān)原理 2. 編寫(xiě) MATLAB 程序?qū)崿F(xiàn) AM 調(diào)制； 3. 編寫(xiě) MATLAB 程序?qū)崿F(xiàn) DSB 調(diào)制； 4. 編寫(xiě) MATLAB 程序?qū)崿F(xiàn) SSB 調(diào)制； 三、實(shí)驗(yàn)原理 調(diào)制是一個(gè)將信號(hào)變換成適于在信道傳輸?shù)倪^(guò)程。由于信源的特性與信道的特性可能不匹配，直接傳輸可能嚴(yán)重影響傳輸質(zhì)量。模擬調(diào)制針對(duì)的信源為模擬信號(hào)，常用的模擬調(diào)制有調(diào)幅、調(diào)相、調(diào)頻。本次實(shí)驗(yàn)進(jìn)行的是模擬信號(hào)的幅度調(diào)制。 幅度調(diào)制是由調(diào)制信號(hào)去控制高頻載波的幅度，使之隨調(diào)制信號(hào)做線性變化的過(guò)程。由于已調(diào)信號(hào)的幅度隨基帶信號(hào)的規(guī)律呈正比地變化，這一特點(diǎn)

11、反映在頻譜結(jié)構(gòu)上，表現(xiàn)為已調(diào)信號(hào)的頻譜完全是基帶信號(hào)頻譜在頻域內(nèi)的簡(jiǎn)單搬移。所以，幅度調(diào)制通常又稱為線性調(diào)制。幅度調(diào)制包括 AM、DSB和 SSB 調(diào)制。四、實(shí)驗(yàn)步驟 0（1）按照 3.1 所提供的 AM 調(diào)制的思路，運(yùn)行提供的范例程序，存檔為 Q2_1，并將所得的結(jié)果存盤(pán)，貼在下面空格處。（2）程序 1 中定義加入的直流分量為A0=2，請(qǐng)?jiān)贏0的值分別改為 1 和 10，看得到的調(diào)制波形會(huì)有什么變化？ A0=1 的調(diào)制波形 A0=10 時(shí)的調(diào)制波形請(qǐng)問(wèn)，調(diào)制波形為什么會(huì)有這種變化，這種變化會(huì)造成什么影響？ 答：（3）按照 3.2 所提供的DSB 調(diào)制的思路，運(yùn)行提供的范例程序，存檔為 Q2_

12、2，并將所得的結(jié)果存盤(pán)，貼在下面空格處。（4） 按照3.3所提供的相移法進(jìn)行SSB調(diào)制的思路，運(yùn)行提供的范例程序， 存檔為Q2_3，并將所得的結(jié)果存盤(pán)，貼在下面空格處。 （5）按照 3.3 所提供的濾波法進(jìn)行 SSB 調(diào)制的思路，編寫(xiě)用濾波法實(shí)現(xiàn) SSB 調(diào)制的程序，存檔為 Q2_4，并將所得的結(jié)果存盤(pán)，貼在下面空格處。 （提示：使用 LPF 函數(shù)實(shí)現(xiàn)濾波） 程序?yàn)椋篶lose all;clear all; dt=0.001; fm=1;fc=10;T=5;t=0:dt:T;mt=sqrt(2)*cos(2*pi*fm*t);s_dsb=mt.*cos(2*pi*fc*t); B=2*fm;f

13、igure(1);subplot(311);plot(t,s_dsb); hold on; plot(t,mt,r-);title(DSB調(diào)制信號(hào));xlabel(t);f_dsb=fft(s_dsb); temp=f_dsb;temp(50:4953)=0; s_ssb=ifft(temp); subplot(312);plot(t,s_ssb);hold on;plot(t,mt,r-);title(SSB上邊帶);xlabel(t);temp=f_dsb;temp(1:49)=0;temp(4953:end)=0; s_ssb=ifft(temp); subplot(313);plot(

14、t,s_ssb);hold on;plot(t,mt,r-);title(SSB下邊帶);xlabel(t); 實(shí)驗(yàn)結(jié)果為：（6）按照實(shí)驗(yàn)原理中介紹的功率譜的計(jì)算公式，在同一圖形的四個(gè)子圖中，分別畫(huà)出基帶信號(hào)、AM 調(diào)制信號(hào)、DSB 調(diào)制信號(hào)和 SSB 調(diào)制信號(hào)的功率譜，要求寫(xiě)出響應(yīng)的程序，畫(huà)出圖形，并在圖中標(biāo)出相應(yīng)的標(biāo)題和坐標(biāo)軸。 程序?yàn)椋?close all; clear all;dt=0.001;fm=1;fc=10;T=5;t=0:dt:T;mt=sqrt(2)*cos(2*pi*fm*t);A=2;s_am=(A+mt).*cos(2*pi*fc*t); s_dsb=mt.*cos(

15、2*pi*fc*t); s_ssb=real(hilbert(mt).*exp(j*2*pi*fc*t); subplot(221)f,sf=T2F(t,mt); psf=(abs(sf).2)/T; plot(f,sf);axis(-2*fc 2*fc 0 max(psf);title(基帶信號(hào)功率譜);xlabel(f);subplot(222)f,sf=T2F(t,s_am);psf=(abs(sf).2)/T; plot(f,psf,r-);hold on;axis(-2*fc 2*fc 0 max(psf);title(AM信號(hào)功率譜);xlabel(f);subplot(223)f,sf=T2F(t,s_dsb);psf=(abs(sf).2)/T; plot(f,sf);a