文档介绍：数字信号处理实验报告(一) 实验一、<1> A=[,3,,-,],n=[-2,-1,0,1,2];B= m,m=0~10; 请编程完成下面操作并画图表示结果 A×B,A-B ,A.*B ,A(n-2 ) 解答: 1、A×B 这是矩阵的叉乘, 所以要保证矩阵乘法的规则。即A 的列数等于 B 的行数。故A 要补零 A2=[,3,,-,]; --------- A3=[,3,,-,,0,0,0,0,0,0]; 同时要转置成列矩阵。 A1=A3'; C=A1*B; 2、A-B D=A3-B; A和B 要保持长度一致,所以用补零的 A3 来和 B 作差。 3、A.*B E=A3.*B; 点乘要求长度一致,计算时对应的元素相乘。 4、A(n-2 ) A右移两位 F=[0,0,A2(1:end)]; 程序: %1 、 A=[,3,,-,],n=[-2,-1,0,1,2];B=,m=0~10; % 请编程完成下面操作并画图表示结果%A ×B, A-B , A.*B ,A( n-2 ) m=0:10; % A1=[;3;;-;;0;0;0;0;0;0]; A2=[,3,,-,]; A3=[,3,,-,,0,0,0,0,0,0];A1=A3'; n1=[-2,-1,0,1,2]; B=.^m; C=A1*B; D=A3-B; subplot(311); stem(m,D); title( 'A-B' ); E=A3.*B; subplot(312); stem(m,E); F=[0,0,A2(1:end)]; subplot(313); u=[-2,-1,0,1,2,3,4]; stem(u,F); <2> p138 , . 解答: 1、求频率响应的幅度相应和相位响应可以直接用 freqz 来求, 也可以是用公式计算出来 clear all ; n=0:1:1000; h=.^n; w=-pi:2*pi/N:pi; H=freqz(h,1,w); % 绘制定制区域, 则直接给定范围, 该范围为给定自变量向量, 例如 freqz(b,a,[-pi:2*pi/n:pi]) plot(w,abs(H)); 注意:用 freqz 做 dtft 时需要求特定范围的 dtft ,这时可以在 freqz 中直接写上给定的范围,比如我们这次要求-pi 到+pi 范围内的 dtft ,所以在 freqz 函数中直接写上范围就可以了。这是一种方法。 x=exp(i*w)./(exp(i*w)-*ones(1,1000)); y=exp(i*w)./(exp(i*w)-*ones(1,1000))+exp(i*w)./(exp(i*w)-*ones( 1,1000)); 也可以像上面那样直接算出来。 2、当输入为 x(n)时,可以用 conv 和 filter 两种方法。 Conv 在计算卷积时,得到的序列的长度是 length MAX([LENGTH(A)+LENGTH(B)-1,LENGTH(A),LENGTH(B)]). 而 filter 卷积时得到的序列长度是 MAX(LENGTH(A),LENGTH(B))-1 ); 可以说 filter 做卷积是 conv 卷积的截断。程序: w=linspace(-pi,pi,1000); % LINSPACE(X1, X2, N) generates N points between X1 and X2. x=exp(i*w)./(exp(i*w)-*ones(1,1000)); % 这是公式 ONES(M,N) or ONES([M,N]) is an M-by-N matrix of ones. % A./B denotes element-by-element division. A and B subplot(4,2,1); plot(w/pi,abs(x)); xlabel( ' 频率 f' ); title( ' 幅度特性' ); ylabel( ' 幅值') grid; %GRID ON adds major grid lines to the current axes. %GRID OFF removes major and minor grid lines from the current axes. subplot(4,2,2); plot(w/pi,angle(x)); grid; xlabel( ' 频率 f' ); title( ' 相位特性' ); ylabe