第一章实验报告

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(-

2、10 20 0 1.2);图形显示如下：1.3修改程序P1.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修改程序P

3、1.1，以产生带有超前7个样本的延时单位阶跃序列sdn。运行修改后的程序并显示产生的序列。% Program P1_5% 产生一个延时的单位阶跃序列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哪个参数控制该序列的增长或衰减率？哪个

4、参数控制该序列的振幅？答：在1.6的程序中，是real(c)控制该序列的增长或者衰减率。K控制该序列的振幅。1.9运算符real和imag的作用是什么？答：运算符real的作用是确定一个复数或者矩阵的实部。运算符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算术运算符和.之间的区别是

5、什么？答：算术运算符表示矩阵乘法，也就是，必须满足前一个矩阵的列数等于后一个矩阵的行数。算术运算符.表示元素相乘，也就是向量或者矩阵中对应元素相乘。1.15该序列的长度是多少？怎样才能改变他？答：运行1.14程序后，得到该序列的长度是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

6、 graphstem(n,x); % Plot the generated sequenceaxis(0 40 -2 2);grid; title('Sinusoidal Sequence');xlabel('Time index n');ylabel('Amplitude');axis;图形显示如下：1.19该序列的长度是多少？怎样可以改变它？答：该序列长度为41；可以通过改变n序列的取值来改变它，例如n=0:50,则序列长度为51，同时应注意axis的x域的取值范围。1.21 axis和grid命令的作用是什么？答：axis的命令的作用是手工

7、设置图中坐标的尺度。grid命令是在当前图形上增加或者减少网格线。1.23修改上述程序，以产生一个长度为50、频率为0.08、振幅为2.5、相移为90度的一个正弦序列并显示它。该序列的周期是多少？n = 0:49; f =0.08; phase = pi/2; A = 2.5; 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');xl

8、abel('时间指标 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 = 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(

9、9;正弦序列 Q1.23 f=0.08');xlabel('时间指标 n');ylabel('振幅');axis;图形显示如下：这是一个杆状图，stairs是将离散的点通过直线连接起来，stem是直接画离散图，而plot是曲线联合在一起，一般画离散图采用stem进行。1.27编写一个MATLAB程序，以产生并显示一个长度为75的高斯随机信号，该信号正态分布且均值为0，方差为3. clf;N=75;x=0+3*randn(1,N);n=1:N;stem(n-1,x(n);title('随机信号');xlabel('时间指标

10、 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-.');xlabel('时间序列 n'

11、);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能产生被噪声污染的

12、信号吗？若不能，为什么？答：不能，因为d产生的序列是长度为51的列向量，而s是长度为51的行向量，必须转置后才能进行运算，可以将d的表达式变为d = 0.8*(rand(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);title('方波 1');xlabel('时间指