《數(shù)字信號處理實(shí)驗(yàn)一.ppt》由會員分享，可在線閱讀，更多相關(guān)《數(shù)字信號處理實(shí)驗(yàn)一.ppt（26頁珍藏版）》請?jiān)谘b配圖網(wǎng)上搜索。
1、下載地址：,ftp://192.168.1.253 用戶名：student 密碼：student 目錄：數(shù)字信號處理-第一次實(shí)驗(yàn),數(shù)字信號處理實(shí)驗(yàn),測控技術(shù)與儀器研究室 指導(dǎo)教師：林曉鋼 ,實(shí)驗(yàn)一 離散時間信號分析,2011年10月,實(shí)驗(yàn)?zāi)康?1掌握各種常用的序列，理解其數(shù)學(xué)表達(dá)式和波形表示。 2掌握在計(jì)算機(jī)中生成及繪制數(shù)字信號波形的方法。 3掌握序列的相加、相乘、移位、反轉(zhuǎn)等基本運(yùn)算及計(jì)算機(jī)實(shí)現(xiàn)與作用。 4掌握線性卷積軟件實(shí)現(xiàn)的方法。 5掌握計(jì)算機(jī)的使用方法和相關(guān)應(yīng)用軟件的使用。 6通過編程，上機(jī)調(diào)試程序，進(jìn)一步增強(qiáng)使用計(jì)算機(jī)解決問題的能力。,實(shí)驗(yàn)原理,1序列的基本概念,離散時間信號在數(shù)學(xué)上
2、可用時間序列x(n)來表示，其中x(n)代表序列的第n個數(shù)字，n代表時間的序列。注意： x(n)只 在 n 為 整 數(shù) 時 才 有 意 義, n 不 是 整 數(shù) 時 無 定 義, 但 不 能 認(rèn) 為 是 0。 離散時間信號可以是由模擬信號通過采樣得到，例如對模擬信號xa(t)進(jìn)行等間隔采樣，采樣間隔為T，得到一個有序的數(shù)字序列xa(nT)就是離散時間信號，簡稱序列。,實(shí)驗(yàn)原理,2常用序列,單位脈沖序列（單位抽樣）、單位階躍序列、矩形序列、實(shí)指數(shù)序列、復(fù)指數(shù)序列、正弦型序列等。,方法1 n1 = -5:5; x1=(n1-0)==0; %顯然,用邏輯式是比較高明的方法 stem(n1,x1)
3、; %作用是繪制以n1為橫坐標(biāo)，x1為縱坐標(biāo)的離散桿圖 title(單位脈沖序列); %繪圖標(biāo)題 axis(-5,5,0,1); %橫坐標(biāo)從-5到5，縱坐標(biāo)從0到1 方法2 x=1 zeros(1,N-1); % 產(chǎn)生長度為N的單位脈沖序列,實(shí)驗(yàn)原理,2常用序列,單位脈沖序列（單位抽樣）、單位階躍序列、矩形序列、實(shí)指數(shù)序列、復(fù)指數(shù)序列、正弦型序列等。,clear; n1 = -5:5; x1=(n1-0)=0; % 用邏輯式是比較高明的方法 stem(n1,x1); title(階躍序列); axis(-5,5,0,1.1);,實(shí)驗(yàn)原理,2常用序列,clear; n1 =
4、 -5:5; x1=(n1-0)=0; x2=(n1-3=0); x3=x1-x2; stem(n1,x3);title(矩形序列) axis(-5,5,0,1.1),實(shí)驗(yàn)原理,2常用序列,，式中0為數(shù)字頻率,clear; n1 = 0:40; x1=exp((0.1+(pi/6)*i)*n1); subplot(2,1,1); stem(n1,real(x1));title(復(fù)指數(shù)序列); ylabel(實(shí)部); subplot(2,1,2); stem(n1,imag(x1)); ylabel(虛部);,實(shí)驗(yàn)原理,2常用序列,clear; n = 0:40; x1=3*cos(0.1*pi
5、*n+pi/3); x2=2*sin(0.2*pi*n); subplot(2,1,1); stem(n,x1);title(余弦序列); subplot(2,1,2); stem(n,x2);title(正弦序列);,實(shí)驗(yàn)原理,2序列的基本運(yùn)算,2.1 序列的加法 ：z(n)=x (n)+y(n) 加法用于序列的合成。它把兩個序列中位置序號相同的樣本相加，形成新的樣本序列。,function y,n = seqadd(x1,n1,x2,n2) % y,n = seqadd(x1,n1,x2,n2) % ----------------------------- % 實(shí)現(xiàn)y(n) = x1(
6、n)+x2(n) % y = 在包含n1和n2的n點(diǎn)上求序列和, % x1 = 在位置向量n1上的第一序列 % x2 = 在位置向量n2上的第二序列(n2可與 n1不同) n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的長度 y1 = zeros(1,length(n)); y2 = y1; % 初始化 y1(find((n=min(n1)) % 序列相加,實(shí)驗(yàn)原理,2序列的基本運(yùn)算,2.2 序列的乘法 ：z(n)=x(n)y(n) 序列的乘法是一種非線性運(yùn)算，它用于信號的調(diào)制。它把兩個序
7、列中位置序號相同的樣本相乘，形成新的樣本序列。例如：,在MATLAB中序列加法可用算術(shù)運(yùn)算符“+”實(shí)現(xiàn)，序列乘法要用點(diǎn)乘運(yùn)算符“.*”實(shí)現(xiàn)。然后首先必須使x1和x2的位置向量的起點(diǎn)、終點(diǎn)相同，長度相等。 Matrix dimensions must agree.,,實(shí)驗(yàn)原理,2序列的基本運(yùn)算 （乘法）,function y,n = seqmult(x1,n1,x2,n2) % y,ny = seqmult(x1,n1,x2,n2) % ----------------------------- % 實(shí)現(xiàn) y(n) = x1(n)*x2(n) % y = 在n區(qū)間上的乘積序列,ny 包含 n1
8、和 n2 % x1 = 在 n1上的第一序列 % x2 = 在 n2上的第二序列(n2可與 n1不同) n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的位置向量 y1 = zeros(1,length(n)); y2 = y1; % y1,y2初始置零 y1(find((n=min(n1)) % 序列相乘,例：求序列x1和x2的和與乘積。,x1=0,1,2,3,4,3,2,1,0;ns1=-2;% 給定x1及ns1 x2=2,2,0,0,0,-2,-2; ns2=2;% 給定x2及ns
9、2 nf1=ns1+length(x1)-1; nf2=ns2+length(x2)-1; ny= min(ns1,ns2):max(nf1,nf2); % y(n)的位置向量 y1 = zeros(1,length(ny)); y2 = y1; % 延拓序列初始化 y1(find((ny=ns1) % 序列相乘,畫圖,subplot(4,2,1), stem(ns1:nf1,x1,.)% 繪圖 xlabel(nx1),ylabel(x1),axis(-5,10,0,4) subplot(4,2,3), stem(ns2:nf2,x2,.),axis(-5,10,-2,2)
10、 xlabel(nx2),ylabel(x2) subplot(4,2,2), stem(ny,y1,.)% 繪圖 xlabel(ny),ylabel(y1) subplot(4,2,4), stem(ny,y2,.) xlabel(ny),ylabel(y2) line(ny(1),ny(end),0,0)% 畫x軸 subplot(4,2,6), stem(ny,ya,.) xlabel(ny),ylabel(ya) line(ny(1),ny(end),0,0)% 畫x軸 subplot(4,2,8), stem(ny,yp,.) xlabel(ny),ylabel(yp) line(n
11、y(1),ny(end),0,0)% 畫x軸 set(gcf,color,w) % 置圖形背景色為白,,實(shí)驗(yàn)原理,2序列的基本運(yùn)算,2.3 序列的移位 ：y(n)=x(n-k) 在這個運(yùn)算中，x(n)的每一個樣本都向右移動k個采樣周期。實(shí)際上， x(n-k)取的是序列過去的值，具有物理可實(shí)現(xiàn)性；而x(n+k) 要知道序列未來的值，物理上無法直接實(shí)現(xiàn)。所以數(shù)字信號處理中通常都用序列右移。,function y,ny = seqshift(x,nx,n0) % y,ny = seqshift(x,nx,n0) % ------------------------- % 實(shí)現(xiàn) y(n) = x
12、(n-n0) % n0為平移樣本數(shù) ny = nx + n0; % 位置向量移位 y = x; % 序列的值不變,實(shí)驗(yàn)原理,2序列的基本運(yùn)算,2.4 序列的反轉(zhuǎn) ：y(n)=x(-n) 在這個運(yùn)算中，x(n)的每一個樣本都對縱坐標(biāo)軸n=0做對稱翻轉(zhuǎn)(flip-flop)得到反轉(zhuǎn)后的序列。在MATLAB中，這一運(yùn)算要兩次調(diào)用fliplr函數(shù)來實(shí)現(xiàn)。y=fliplr(x)的作用是把行向量x中元素排列的次序左右翻轉(zhuǎn)。ny=-fliplr(nx)既將其位置向量左右翻轉(zhuǎn)，又變了正負(fù)號。 程序： y = fliplr(x); % 將序列數(shù)值左右翻轉(zhuǎn) ny = -fliplr(nx); %
13、將序列位置對零位置左右翻轉(zhuǎn)，故同時改變正負(fù)號,實(shí)驗(yàn)原理,3序列的卷積,上式的運(yùn)算關(guān)系稱為卷積運(yùn)算，式中代表兩個序列卷積運(yùn)算。兩個序列的卷積是一個序列與另一個序列反褶后逐次移位乘積之和，故稱為離散卷積，也稱兩序列的線性卷積。 四個步驟： （1）反轉(zhuǎn)：先將x(n)和h(n)的變量n換成m，變成x(m)和h(m) ，再將h(m)以縱軸為對稱軸反轉(zhuǎn)成h(-m) 。 （2）移位：將h(-m)移位n，得h(n-m) 。當(dāng)n為正數(shù)時，右移n位；當(dāng)n為負(fù)數(shù)時，左移n位。 （3）相乘：將h(n-m)和x(m)的對應(yīng)點(diǎn)值相乘。 （4）求和：將以上所有對應(yīng)點(diǎn)的乘積累加起來，即得y(n) 。,function y,
14、ny = seqfold(x,nx) y = fliplr(x); ny = -fliplr(nx);,function y,ny = seqshift(x,nx,n0) ny = nx + n0; y = x;,function y,ny = seqmult(x,nx,h2,nh2),y = sum(y1),函數(shù)調(diào)用,x1=0,1,2,3,4,3,2,1,0; n1=-2:6;% 給定x1及ns1 x2=2,2,0,0,0,-2,-2; n2=2:8; y,n = seqadd(x1,n1,x2,n2); stem(n,y);,實(shí)驗(yàn)內(nèi)容,1、用matlab編制信號產(chǎn)生子程序，產(chǎn)生以下典型信號
15、序列。,,實(shí)驗(yàn)內(nèi)容,2、利用MATLAB編程完成序列(5)和(6)的移位、反轉(zhuǎn)、加法、乘法等運(yùn)算，并繪制運(yùn)算后序列的波形。,,3、利用MATLAB編制一個計(jì)算兩個序列線性卷積的通用程序，計(jì)算上述兩序列，并繪制卷積后序列的波形。,實(shí)驗(yàn)報(bào)告要求,1簡述實(shí)驗(yàn)原理及目的。 2給出上述序列的實(shí)驗(yàn)結(jié)果。 3列出計(jì)算卷積的公式，畫出程序框圖，并列出實(shí)驗(yàn)程序清單（包括必要的程序說明）。 4記錄調(diào)試運(yùn)行情況及所遇問題的解決方法。 5給出實(shí)驗(yàn)結(jié)果，并對結(jié)果做出分析。 6簡要回答思考題。 如何產(chǎn)生方波信號序列和鋸齒波信號序列？ 實(shí)驗(yàn)中所產(chǎn)生的正弦序列的頻率是多少？是否是周期序列？,實(shí)驗(yàn)報(bào)告封面格式,實(shí)驗(yàn)課程名稱：數(shù)字信號處理 開 課 實(shí)驗(yàn) 室：測控技術(shù)與儀器研究室 學(xué) 院：光電工程學(xué)院 年級 2009 專業(yè)班 XXX 學(xué) 生 姓 名： XXX 學(xué)號 2009XXXX 開 課 時 間： 2011 至 2012學(xué)年 第 1 學(xué)期,實(shí)驗(yàn)報(bào)告封面格式,開課學(xué)院實(shí)驗(yàn)室：光電工程學(xué)院、測控技術(shù)與儀器研究室 實(shí) 驗(yàn) 時 間：2011年10月25日 學(xué)生 ： XXX 專業(yè)班、年級 XXX 2009,測控技術(shù)與儀器研究室,下載地址：,ftp://192.168.1.253 用戶名：student 密碼：student 目錄：數(shù)字信號處理-第一次實(shí)驗(yàn),
數(shù)字信號處理實(shí)驗(yàn)一.ppt
