第二章 数字信号处理.docVIP

第二章 所用函数： function [ y,n ] sigadd x1,n1,x2,n2 %implements y n x1 n +x2 n % [y,n] sigadd x1,n1,x2,n2 % y sum sequence over n,which includes n1 and n2 % x1 first sequence over n1 % x2 second sequence over n2 n2 can b different from n1 n min min n1 ,min n2 : max max n1 ,max n2 ; %duration of y n y1 zeros 1,length n ; y2 y1; %initialization y1 find n min n1 n max n1 1 x1; %x1 with duration of y y2 find n min n2 n max n2 1 x2; %x2 with duration of y y y1 + y2; %sequence addition function [ y,n ] sigshift x,m,n0 %implements y n x n-n0 % [y,n] sigshift x,m,n0 % n m+n0;y x; function [ y,n ] sigmult x1,n1,x2,n2 %implements y n x1 n *x2 n % [y,n] sigmult x1,n1,x2,n2 % y product sequence over n,which includes n1 and n2 % x1 first sequence over n1 % x2 second sequence over n2 n2 can be different from n1 % n min min n1 ,min n2 : max max n1 ,max n2 ; %duration of y n y1 zeros 1,length n ; y2 y1; % y1 find n min n1 n max n1 1 x1; %x1 with duration of y y2 find n min n2 n max n2 1 x2; %x2 with duration of y y y1 .* y2; %sequence multiplication function [ y,ny ] conv_m x,nx,h,nh %modified convolution routine for signal processing % ------------------------------------------------- % [y,ny] conv_m x,nx,h,nh % [y,ny] convolution result % [x,nx] first signal % [h,nh] second signal % nyb nx 1 +nh 1 ;nye nx length x +nh length h ; ny [nyb:nye]; y conv x,h ; P2.6设x n 2,4,-3,-1,-5,4,7 ,产生并画出下列序列的样本 用stem函数 。 2、x2 n 4x 4+n +5x n+5 +2x n 程序： n [-3:3]; x [2,4,-3,-1,-5,4,7]; [x21,n21] sigshift x,n,-4 ; [x22,n22] sigshift x,n,-5 ; [x2,n2] sigadd 4*x21,n21,5*x22,n22 ; [x2,n2] sigadd x2,n2,2*x,n ; stem n2,x2,filled ; 结果： 4、x4 n 2*e^ 0.5*n *x n +cos 0.1*pi*n *x n+2 ,-10 n 10 程序： n1 [-10:10]; n2 [-3:3]; x [2,4,-3,-1,-5,4,7]; y1 exp 0.5*n1 ; y2 cos 0.1*pi*n1 ; [x41,n41] sigshift x,n2,-2 ; [x42,n42] sigmult y1,n1,x,n2 ; [x43,n43] sigmult y2,n1,x41,n41 ; [x4,n4] sigadd 2*x42,n42,x43,n43 ;