第一章实验报告

1、数字信号处理实验实验报告第一章：离散时间信号的时域分析学 院： 专 业： 班 级： 学 号： 学生姓名： 2015年 月 日1.1运行程序P1.1，以产生单位样本序列u(n)并显示它。% 程序 P1_1% 单位采样序列的产生 clf;% 产生向量n=-10:20n = -10:20;% 单位采样序列的产生u = zeros(1,10) 1 zeros(1,20);% Plot the unit sample sequencestem(n,u);xlabel(时间指标 n);ylabel(振幅);title(单位采样序列);axis(-10 20 0 1.2);图形显示如下：1.3修改程序P1.

2、1，以产生带有延时11个样本的延迟单位样本序列udn。运行修改的程序并显示产生的序列。 % Program P1_3% 产生单位采样序列 clf;% 产生一个向量n=-10:20n = -10:20;% 产生单位采样序列u = zeros(1,21) 1 zeros(1,9); %改变序列位置% 绘制单位采样序列stem(n,u);xlabel(时间指标n);ylabel(振幅);title(单位采样序列);axis(-10 20 0 1.2);图形显示如下：1.5修改程序P1.1，以产生带有超前7个样本的延时单位阶跃序列sdn。运行修改后的程序并显示产生的序列。% Program P1_5%

3、 产生一个延时的单位阶跃序列clf;% 产生一个向量n=-10:20n = -10:20;% 产生一个超前的单位序列u = zeros(1,17) ones(1,14); %将0序列延长至n=6时即可% 绘制单位阶跃序列stem(n,u);xlabel(时间指标 n);ylabel(振幅);title(单位阶跃延迟序列);axis(-10 20 0 1.2);图形显示如下：1.7哪个参数控制该序列的增长或衰减率？哪个参数控制该序列的振幅？答：在1.6的程序中，是real(c)控制该序列的增长或者衰减率。K控制该序列的振幅。1.9运算符real和imag的作用是什么？答：运算符real的作用是确

4、定一个复数或者矩阵的实部。运算符imag的作用是确定一个复数或者矩阵的虚部。1.11运行程序P1.3，以产生实数值得指数序列。% Program P1_3% 产生一个实的指数序列clf;n = 0:35; a = 1.2; K = 0.2;x = K*a.n;stem(n,x);xlabel(时间指标 n);ylabel(振幅);图形显示如下：1.13算术运算符和.之间的区别是什么？答：算术运算符表示矩阵乘法，也就是，必须满足前一个矩阵的列数等于后一个矩阵的行数。算术运算符.表示元素相乘，也就是向量或者矩阵中对应元素相乘。1.15该序列的长度是多少？怎样才能改变他？答：运行1.14程序后，得到

5、该序列的长度是36；通过改变“n=0:35”中的“0”或者“35”这个数字可以改变序列的长度。1.17运行程序P1.4，以产生正弦序列并显示它。% Program P1_4% 产生一个正弦序列n = 0:40; f = 0.1; phase = 0; A = 1.5; arg = 2*pi*f*n - phase; x = A*cos(arg);clf; % Clear old graphstem(n,x); % Plot the generated sequenceaxis(0 40 -2 2);grid; title(Sinusoidal Sequence);xlabel(Time ind

6、ex n);ylabel(Amplitude);axis;图形显示如下：1.19该序列的长度是多少？怎样可以改变它？答：该序列长度为41；可以通过改变n序列的取值来改变它，例如n=0:50,则序列长度为51，同时应注意axis的x域的取值范围。1.21 axis和grid命令的作用是什么？答：axis的命令的作用是手工设置图中坐标的尺度。grid命令是在当前图形上增加或者减少网格线。1.23修改上述程序，以产生一个长度为50、频率为0.08、振幅为2.5、相移为90度的一个正弦序列并显示它。该序列的周期是多少？n = 0:49; f =0.08; phase = pi/2; A = 2.5;

7、arg = 2*pi*f*n - phase; x = A*cos(arg);clf; % Clear old graphstem(n,x); % Plot the generated sequenceaxis(0 49 -3 3);grid; title(正弦序列 Q1.23 f=0.08);xlabel(时间指标 n);ylabel(振幅);axis;图形显示如下：1.25在程序P1.4中用stairs命令代替stem命令，运行新程序。新图形与习题 Q1.17和习题Q1.24中产生的图形有什么区别？n = 0:49; f =0.08; phase = pi/2; A = 2.5; arg

8、= 2*pi*f*n - phase; x = A*cos(arg);clf; % Clear old graphstairs(n,x); % Plot the generated sequenceaxis(0 49 -3 3);grid; title(正弦序列 Q1.23 f=0.08);xlabel(时间指标 n);ylabel(振幅);axis;图形显示如下：这是一个杆状图，stairs是将离散的点通过直线连接起来，stem是直接画离散图，而plot是曲线联合在一起，一般画离散图采用stem进行。1.27编写一个MATLAB程序，以产生并显示一个长度为75的高斯随机信号，该信号正态分布且

9、均值为0，方差为3.clf;N=75;x=0+3*randn(1,N);n=1:N;stem(n-1,x(n);title(随机信号);xlabel(时间指标 n);ylabel(x);axis;图形显示为：1.29运行程序P1.5，以产生所有相关的信号。% Program P1_5%通过平均的信号平滑clf;R = 51;d = 0.8*(rand(R,1) - 0.5); % 产生随机噪声m = 0:R-1;s = 2*m.*(0.9.m); % 产生为污染的信号x = s + d; % 产生被噪音污染的信号subplot(2,1,1);plot(m,d,r-,m,s,g-,m,x,b-.

10、);xlabel(时间序列 n);ylabel(振幅);legend(dn ,sn ,xn );x1 = 0 0 x;x2 = 0 x 0;x3 = x 0 0;y = (x1 + x2 + x3)/3;subplot(2,1,2);plot(m,y(2:R+1),r-,m,s,g-);legend( yn ,sn );xlabel(时间序列 n);ylabel(振幅);图形显示为：1.31使用语句x=s+d能产生被噪声污染的信号吗？若不能，为什么？答：不能，因为d产生的序列是长度为51的列向量，而s是长度为51的行向量，必须转置后才能进行运算，可以将d的表达式变为d = 0.8*(rand(

11、1,R) - 0.5);1.33 legend命令的作用是什么？答：用一个给定字符串作为标记在当前图上插上说明。1.35算术运算符*和.*之间的区别是什么？答：算术运算符*表示矩阵相乘。算术运算符.*表示元素相乘。1.37该信号的最大频率和最小频率是多少？答：假设自变量的形式是，则角频率是；最大频率为pi；最小频率为0。1.39在命令窗口中键入who。此时在命令窗口会显示什么信息？答：who：列出内存中的当前变量。1.41编写matlab程序，以产生指定方波和锯齿波序列，并将序列绘制出来。t=0:30;x=square(2*pi*0.1.*t,60); %占空比 60%stem(t,x);axis(0,30,-1.5,1.5)