dsp實(shí)驗(yàn)matlab參考模板

1、實(shí)驗(yàn)一 數(shù)字信號(hào)的產(chǎn)生與基本運(yùn)算1 / 71實(shí)驗(yàn)要求(1) 常用數(shù)字信號(hào)序列的產(chǎn)生：熟悉 Matlab 產(chǎn)生數(shù)字信號(hào)的基本命令，加深對(duì)數(shù)字信號(hào)概念的理解，并能夠用 Matlab 產(chǎn)生和繪制出一些常用離散信號(hào)序列。請(qǐng)用 Matlab 畫出下列序列的波形（-10<n<10）:a) (n)b) 單位階躍序列 2 u(n-5)c) 矩形序列 R(n)d) y(n)2sin(0.3n)+ 0.5cos2(0.6n)（2）加、減、尺度（乘除）和移位是數(shù)字信號(hào)處理中最基本的算術(shù)運(yùn)算，將上述基本序列進(jìn)行這些基本運(yùn)算，得到多個(gè)序列構(gòu)成的組合序列。（3）請(qǐng)用您的計(jì)算機(jī)聲卡采用一段您自己的聲音 x(n)

2、，長(zhǎng)度為 45秒，單聲道，取樣頻率 44.1kHz，16bit/樣值，然后與給定的一段背景音樂 y(n) 按下式疊加為一個(gè)聲音信號(hào) z(n)：z(n) = 0.7x(n) + 0.3y(n)要求在同一個(gè) Figure 中,畫出采集聲音 x(n)、背景音樂 y(n)和混音 z(n) 的時(shí)域波形；2、實(shí)驗(yàn)結(jié)果（1）常用數(shù)字信號(hào)序列的產(chǎn)生單位取樣序列x=impseq(0,-10,10);n=-10:10;stem(n,x)單位階躍序列 2 u(n-5 x=2*stepseq(5,-10,10);n=-10:10;stem(n,x)矩形序列：x1=stepseq(-10,-20,20);n1=-20:

3、20;x2=stepseq(10,-20,20);n2=-20:20;y=x1-x2;n=-20:20;stem(n,y)y(n)2sin(0.3n)+ 0.5cos2(0.6n)：n=-10:0.1:10;x=2*sin(0.3*pi*n)+0.5*cos(2*(0.6*pi*n);stem(n,x)（2）數(shù)字信號(hào)的基本運(yùn)算加法：(n)+ 2 u(n-5)x1=impseq(0,-10,10);n1=-10:10;x2=2*stepseq(5,-10,10);n2=-10:10;y=x1+x2;n=-10:10;stem(n,y)減法(n)- 2 u(n-5)x1=impseq(0,-10,

4、10);n1=-10:10;x2=2*stepseq(5,-10,10);n2=-10:10;y=x1-x2;n=-10:10;stem(n,y)乘法: u(n)* u(n-5)x1=stepseq(0,-10,10);n1=-10:10;x2=2*stepseq(5,-10,10);n2=-10:10;y=x1.*x2;n=-10:10;stem(n,y)移位：聲音波形record,fs,bits=wavread('voice.wav');background,fs,bits=wavread('background.wav');x=record(fs*0+1:fs*45),:);y=background(fs*0+1:fs*45),:);z=0.7*x+0.3*y;wavwrite(z,fs,bits,'10210607.wav');subplot(3,1,1),plot(x),title('voice');subplot(3,1,2),plot(y),title('background');subplot(3,1,3),plot(z