北邮数字信号处理实验报告

1、2011级数字信号处理实验报告实验名称： 实验一 数字信号的产生和基本运算1实验要求因为现实世界里存在的是模拟信号，因此数字信号处理的第一个问题是将信号离散化，得到一个数字信号，然后再进行数字处理。(1) 常用数字信号序列的产生：熟悉 Matlab 产生数字信号的基本命令，加深对数字信号概念的理解，并能够用 Matlab 产生和绘制出一些常用离散信号序列。请用 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) 数字信号的基本运算：a)

2、 加、减、尺度（乘除）和移位是数字信号处理中最基本的算术运算，将上述基本序列进行这些基本运算，得到多个序列构成的组合序列。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) 的时域波形；³ 提交实验报告时，请把声音文件转换为 mp3 格式，图像转换为 JPEG 格式，以节省存储空间。&#

3、179; 通过本次实验，掌握 Matlab 中这些基本运算命令，对数字信号处理有一个基本概念，为后面的数字信号分析和滤波打下基础。2. 程序分析第一部分：1、 (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、 阶跃函数u（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, 窗函数R(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、运行结果第二部分：1、各个函数求和>>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) ');运行结果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,&