【正文】 ll figure windows n = 5:5。選項(xiàng)，因此繪制的圖形中每根棒條線的頂端是一個(gè)實(shí)心點(diǎn)。) 由于在程序的stem(n,x,39。) % Open a figure window and draw the plot of x[n] grid on, title (39。 程序Program1_3用來(lái)仿真下面形式的離散時(shí)間信號(hào)： x[n]={...., , , , 0, , ….} ↑n=0 % Program1_3 % This program is used to generate a discretetime sequence % and draw its plot clear, % Clear all variables close all, % Close all figure windows n = 5:5。Sinusoidal signal x[n]39。 離散時(shí)間信號(hào)的仿真 程序Program1_2用來(lái)產(chǎn)生離散時(shí)間信號(hào)x[n]=sin()。產(chǎn)生單位沖激信號(hào)的擴(kuò)展函數(shù)為： function y = delta(t) dt = 。) xlabel(39。 % Specify the step of time variable t = 2:dt:2。 注意，此時(shí)imp的點(diǎn)數(shù)應(yīng)與g或G相同。卷積后s的長(zhǎng)度是length(imp)+length(g)1?，F(xiàn)規(guī)定M、N都是2的整冪，于是L自然也是2的整冪。2) 產(chǎn)生沖擊序列信號(hào)imp(t)imp=zeros(1,N)。PAM信號(hào)是指所有形如的信號(hào)。語(yǔ)句sign(rand(1,M))產(chǎn)生M個(gè)取值177。設(shè)仿真系統(tǒng)的取樣點(diǎn)數(shù)是N，系統(tǒng)帶寬為Bs，矢量的總功率為，最高頻率分量為Bs，并且各樣點(diǎn)的值互不相關(guān)，故它代表雙邊功率譜密度為（W/Hz）的白噪聲。這兩個(gè)函數(shù)分別如附錄所示。對(duì)應(yīng)的其他參數(shù)為Bs＝1MHz。要提高仿真的精度，就必須降低時(shí)域取樣間隔Δt及頻域取樣間隔Δf，也就是要加大總?cè)狱c(diǎn)數(shù)N。考慮到系統(tǒng)帶寬為Bs，便把頻譜的截短區(qū)間設(shè)計(jì)為[Bs, Bs]然后再按間隔Δf均勻取樣，得到取樣點(diǎn)數(shù)為： (33)將式（32）代入式（33）得 (34)同樣，信號(hào)在頻域被離散后，對(duì)應(yīng)到時(shí)域也是一個(gè)周期信號(hào)，其周期為。如果信號(hào)的最高頻率為fH，那么必須有fH≤才能保證不發(fā)生頻域混疊失真。 信號(hào)及系統(tǒng)在計(jì)算機(jī)中的表示 時(shí)域取樣及頻域取樣一般來(lái)說(shuō)，任意信號(hào)s(t)是定義在時(shí)間區(qū)間（∞，+∞）上的連續(xù)函數(shù)，但所有計(jì)算機(jī)的CPU都只能按指令周期離散運(yùn)行，同時(shí)計(jì)算機(jī)也不能處理（∞，+∞）這樣一個(gè)時(shí)間段。而從數(shù)學(xué)的角度來(lái)看，信息從一地傳送到另一地的整個(gè)過(guò)程或者各個(gè)環(huán)節(jié)不外乎是一些碼或信號(hào)的交換過(guò)程。通信原理實(shí)驗(yàn)報(bào)告 實(shí)驗(yàn)名稱(chēng)： 實(shí)驗(yàn)一 數(shù)字基帶傳輸系統(tǒng)的MATLAB仿真 實(shí)驗(yàn)二 模擬信號(hào)幅度調(diào)制仿真實(shí)驗(yàn) 班 級(jí)： 10通信工程三班 學(xué) 號(hào)： 2010550920 姓 名： 彭龍龍 指導(dǎo)老師： 王仕果 實(shí)驗(yàn)一 數(shù)字基帶傳輸系統(tǒng)的MATLAB仿真一、實(shí)驗(yàn)?zāi)康氖煜ず驼莆粘Ｓ玫挠糜谕ㄐ旁頃r(shí)域仿真分析的MATLAB函數(shù)； 掌握連續(xù)時(shí)間和離散時(shí)間信號(hào)的MATLAB產(chǎn)生； 牢固掌握沖激函數(shù)和階躍函數(shù)等函數(shù)的概念，掌握卷積表達(dá)式及其物理意義，掌握卷積的計(jì)算方法、卷積的基本性質(zhì)； 掌握利用MATLAB計(jì)算卷積的編程方法，并利用所編寫(xiě)的MATLAB程序驗(yàn)證卷積的常用基本性質(zhì)； 掌握MATLAB描述通信系統(tǒng)中不同波形的常用方法及有關(guān)函數(shù)，并學(xué)會(huì)利用MATLAB求解系統(tǒng)功率譜，繪制相應(yīng)曲線。例如信源壓縮編碼、糾錯(cuò)編碼、AMI編碼、擾碼等屬于碼層次上的變換，而基帶成形、濾波、調(diào)制等則是信號(hào)層坎上的處理。為此將把s(t)按區(qū)間截短為sT(t)，再對(duì)sT(t)按時(shí)間間隔Δt均勻取樣，得到取樣點(diǎn)數(shù)為： (31)仿真時(shí)用這個(gè)樣值集合來(lái)表示信號(hào)s(t)。設(shè) (32)則稱(chēng)Bs為仿真系統(tǒng)的系統(tǒng)帶寬。如果時(shí)域截短時(shí)間為T(mén)，那么必須T≤1/df才能保證不發(fā)生頻域混疊失真。這說(shuō)明仿真的精度與仿真系統(tǒng)的運(yùn)算量直接有關(guān)。 頻域分析限于篇幅，在此不介紹MATLAB中關(guān)于傅立葉變換的有關(guān)函數(shù)。我們關(guān)心的另一個(gè)指標(biāo)是信號(hào)的功率譜密度，任意信號(hào)s(t)的功率譜的定義是，其中是s(t)截短后所得信號(hào)的傅立葉變化，是的能量譜，是在截短時(shí)間T內(nèi)的功率譜。 隨機(jī)碼序列的產(chǎn)生語(yǔ)句round(rand(1,M))產(chǎn)生M個(gè)取值0等概的隨機(jī)碼。1等概的隨機(jī)碼。它可以用如下的等效模型來(lái)表示：故對(duì)所有不同的數(shù)字信號(hào)，都可以用相同的方法來(lái)產(chǎn)生。imp(1:L:M)=a/dt。3) 產(chǎn)生PAM數(shù)字信號(hào)s(t)代表脈沖波形，例如，升余弦滾降要求的波形有：設(shè)的傅立葉變換為，的傅立葉變換為，則有： %t2f為傅立葉變換函數(shù) %f2t為傅立葉反變換函數(shù) 設(shè)矢量s代表數(shù)字信號(hào)，矢量g代表脈沖波形，矢量G代表其頻譜。扣除延遲時(shí)間及拖尾時(shí)間后，數(shù)字信號(hào)為： Ii=find(g=max(g))。若g的寬度小于imp，則應(yīng)用零補(bǔ)齊。 % Specify the interval of time x = sin(2*pi*t)。Time t (sec)39。 y = (u(t)u(tdt))/dt。 % Program1_2 % This program is used to generate a discretetime sinusoidal signal and draw its plot clear, % Clear all variables close all, % Close all figure windows n = 10:10。) xlabel (39。 % Specify the interval of time, the number of points of n is 11. x = [0, 0, 0, 0, , , , 0, , 0, 0]。A discretetime sequence x[n]39。.39。 如果需要在序列的前后補(bǔ)較多的零的話，可以利用函數(shù)zeros()，其語(yǔ)法為zeros(1, N)：圓括號(hào)中的1和N表示該函數(shù)將產(chǎn)生一個(gè)一行N列的矩陣，矩陣中的所有元素均為零。 % Specify the interval of time x = [zeros(1,4), , , , 0, , zeros(1,2)]。r39。（2）修改program1_1,，存盤(pán)程序名為Q1_2，生成實(shí)指數(shù)信號(hào)x(t)=e2t。 % Specify the step of time variable t = 2:dt:2。) xlabel(39。 　y = (u(t)u(tdt))/dt。 % Specify the interval of time axis{[2,5,]}x = [zeros(1,4), , , , 0, , zeros(1,2)]。r39。 grid ontitle(39。輸入?yún)?shù)t表示輸入數(shù)組的時(shí)間變量，