北郵數(shù)字信號處理實驗報告

1、2011級數(shù)字信號處理實驗報告實驗名稱： 實驗一 數(shù)字信號的產(chǎn)生和基本運算1實驗要求因為現(xiàn)實世界里存在的是模擬信號，因此數(shù)字信號處理的第一個問題是將信號離散化，得到一個數(shù)字信號，然后再進行數(shù)字處理。(1) 常用數(shù)字信號序列的產(chǎn)生：熟悉 Matlab 產(chǎn)生數(shù)字信號的基本命令，加深對數(shù)字信號概念的理解，并能夠用 Matlab 產(chǎn)生和繪制出一些常用離散信號序列。請用 Matlab 畫出下列序列的波形（-10<n<10）:a) (n)b) 單位階躍序列 2 u(n-5)c) 矩形序列 R(n) d) y(n)2sin(0.3n)+ 0.5cos(0.6n) (2) 數(shù)字信號的基本運算：a)

2、 加、減、尺度（乘除）和移位是數(shù)字信號處理中最基本的算術(shù)運算，將上述基本序列進行這些基本運算，得到多個序列構(gòu)成的組合序列。b) 請用您的計算機聲卡采用一段您自己的聲音 x(n)，長度為 45秒，單聲道，取樣頻率 44.1kHz，16bit/樣值，然后與給定的一段背景音樂 y(n) 按下式疊加為一個聲音信號 z(n)：z(n) = 0.7x(n) + 0.3y(n) 要求：³ 在同一個 Figure 中,畫出采集聲音 x(n)、背景音樂 y(n)和混音 z(n) 的時域波形；³ 提交實驗報告時，請把聲音文件轉(zhuǎn)換為 mp3 格式，圖像轉(zhuǎn)換為 JPEG 格式，以節(jié)省存儲空間。&#

3、179; 通過本次實驗，掌握 Matlab 中這些基本運算命令，對數(shù)字信號處理有一個基本概念，為后面的數(shù)字信號分析和濾波打下基礎(chǔ)。2. 程序分析第一部分：1、 (n)函數(shù)實現(xiàn)% create impulse function functionx,n=impseq(n0,n1,n2)%This function should be used to create a impseq array%It return 1 at n0 ,and return 0 at other place%You should set the arguments as: n1<n0<n2if(n0<n

4、1)|(n0>n2)|(n1>n2) error('Arguments must satisfy n1<=n0<=n2')endn=n1:n2;x=(n-n0)=0;2、 階躍函數(shù)u（n）的實現(xiàn)%create u(n)function x,n=stepseq(n0,n1,n2)%This function should be used to create a u(n) array%It return 1 when n>=n0 ,and return 0 at other place%You should set the arguments as:

5、n1<n0<n2if(n0<n1)|(n0>n2)|(n1>n2)error('Arguments must satisfy n1<=n0<=n2')endn=n1:n2;x=(n-n0)>=0;3, 窗函數(shù)R(n)的實現(xiàn)%create R(n)function y,n=R(n1,left,right,n2)%This function should be used to create a R(n) array%It return 1 when left<=n<right ,and return 0 at other

6、place%You should set the arguments as: n1<=left<=right<=n2if(n1>left|left>right|right>n2) error('Arguments must satisfy ns<=left<=right<=nf')endn=n1:n2;y=n>=left-n>=right;4、程序主代碼<<cleary1=impseq(0,-10,10);y2=stepseq(5,-10,10)*2;y3=R(-10,-5,5,10);n=-10:1

7、:10;temp=-10:0.1:10;y4=2*sin(0.3*pi*n)+0.5*cos(0.6*pi*n);y5=2*sin(0.3*pi*temp)+0.5*cos(0.6*pi*temp);subplot(2,2,1),stem(n,y1),title('(n)');subplot(2,2,2),stem(n,y2),title('u(n)');subplot(2,2,3),stem(n,y3),title('R(n)');subplot(2,2,4),stem(n,y4),hold on,plot(temp,y5,':'

8、;),title('2sin(0.3pi*n)+0.5cos(0.6pi*n)');5、運行結(jié)果第二部分：1、各個函數(shù)求和>>cleary1=impseq(0,-10,10);y2=stepseq(5,-10,10)*2;y3=R(-10,-5,5,10);n=-10:1:10;y4=2*sin(0.3*pi*n)+0.5*cos(0.6*pi*n);z1=y1+y2;z2=y3+y4;z3=y1*2;subplot(3,1,1),stem(n,z1),title('(n)+2u(n-5)');subplot(3,1,2),stem(n,z2),ti

9、tle('2sin(0.3pi*n)+0.5cos(0.6pi*n)+R(n)');subplot(3,1,3),stem(n,z3),title('2*(n) ');運行結(jié)果2、 音頻處理clearrecord,fs,bits=wavread('record.wav');background,fs,bits=wavread('background.wav');x=record(fs*0+1:fs*45),:);%錄音y=background(fs*0+1:fs*45),:); z=0.1*x+0.9*y; subplot(3,1,1),plot(x),title('record');subplot(3,1,2),plot(y),title('background');subplot(3,1,3),plot(z),title('music');wavwrite(z,fs,bits,&