序列的FT变换.pptVIP

几种常用序列 function [x,n] = impseq(n0,n1,n2) % 产生 x(n) = delta(n-n0); n1 = n,n0 = n2 % [x,n] = impseq(n0,n1,n2) if ((n0 n1) | (n0 n2) | (n1 n2)) error(参数必须满足 n1 = n0 = n2) end n = [n1:n2]; %x = [zeros(1,(n0-n1)), 1, zeros(1,(n2-n0))]; x = [(n-n0) == 0]; 几种常用序列 (2)单位阶跃函数 单位阶跃序列的产生函数如下： function [x,n] = stepseq(n0,n1,n2) % 产生 x(n) = u(n-n0); n1 = n,n0 = n2 % [x,n] = stepseq(n0,n1,n2) if ((n0 n1) | (n0 n2) | (n1 n2)) error(参数必须满足 n1 = n0 = n2) end n = [n1:n2]; %x = [zeros(1,(n0-n1)), ones(1,(n2-n0+1))]; x = [(n-n0) = 0];  几种常用序列 例2.1 用MATLAB产生各种离散序列。 解 MATLAB程序如下： n=[-5:5]; x1=impseq(0,-5,5); subplot(2,2,1);stem(n,x1);title(单位脉冲序列) xlabel(n);ylabel(x(n)); n=[0:10]; x2=stepseq(0,0,10); subplot(2,2,2);stem(n,x2);title(单位阶跃序列); xlabel(n);ylabel(x(n));  几种常用序列 n=[0:10]; x3=stepseq(0,0,10)-stepseq(5,0,10); subplot(2,2,3);stem(n,x3);title(矩形序列); xlabel(n);ylabel(x(n)); n=[0:20]; x4=sin(0.3*n); subplot(2,2,4);stem(n,x4);title(正弦序列); xlabel(n); ylabel(x(n));  几种常用序列  序列的运算 8. 用MATLAB实现序列的运算 (1)两个序列相加减 function [y,n] = sigadd(x1,n1,x2,n2) % 实现 y(n) = x1(n)+x2(n) % [y,n] = sigadd(x1,n1,x2,n2) % y = 在包含n1 和 n2 的n点上求序列和 % x1 = 在 n1上的第一序列 % x2 = 在 n2上的第二序列(n2可与 n1不等)  序列的运算 n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度 y1 = zeros(1,length(n)); y2 = y1; % 初始化 y1(find((n=min(n1))(n=max(n1))==1))=x1; % 具有y的长度的 x1 y2(find((n=min(n2))(n=max(n2))==1))=x2; % 具有y的长度的x2 y = y1+y2; % 序列相加.  序列的运算