MATLAB基础与应用教程(人民邮电出版社-蔡旭辉)第七章a课件.ppt

第 7章 MATLAB在数字信号处理中的应用;7.1 时域离散信号的产生及时域处理;单位脉冲序列δ(n-n0)的生成函数impseq function [x,n]=impseq(n0,ns,nf) n=[ns:nf];x=[(n-n0)==0]; 单位阶跃序列u(n-n0)的生成函数stepseq function [x,n]=stepseq(n0,ns,nf) n=[ns:nf]; x=[(n-n0)=0]; P173表7.1;例7.1 序列的相加和相乘;%第七章例7.1q701离散信号序列的相加和相乘 x1=[0,1,2,3,4,3,2,1,0];ns1=-2; % 给定x1及ns1 x2=[2,2,0,0,0,-2,-2]; ns2=2; % 给定x2及ns2 nf1=ns1+length(x1)-1; nf2=ns2+length(x2)-1; ny= min(ns1,ns2):max(nf1,nf2) % y(n)的时间变量 xa1 = zeros(1,length(ny)); xa2 = xa1; % 延拓序列初始化 xa1(find((ny=ns1)(ny=nf1)==1))=x1 % 给xa1赋值x1 xa2(find((ny=ns2)(ny=nf2)==1))=x2 % 给xa2赋值x2 ya = xa1 + xa2 % 序列相加 yp = xa1.* xa2 % 序列相乘 subplot(4,1,1), stem(ny,xa1,.) % 绘图 subplot(4,1,2), stem(ny,xa2,.) line([ny(1),ny(end)],[0,0]) % 画x轴 subplot(4,1,3), stem(ny,ya,.) line([ny(1),ny(end)],[0,0]) % 画x轴 subplot(4,1,4), stem(ny,yp,.) line([ny(1),ny(end)],[0,0]) % 画x轴;ny = -2 -1 0 1 2 3 4 5 6 7 8xa1 = 0 1 2 3 4 3 2 1 0 0 0xa2 = 0 0 0 0 2 2 0 0 0 -2 -2ya = 0 1 2 3 6 5 2 1 0 -2 -2yp = 0 0 0 0 8 6 0 0 0 0 0;例7.2 序列的合成和截取 ;%第七章例7.2程序q702 离散信号序列的相加和相乘 clear;close all n0=input(输入序列起点:n0=); %书上取-6 N=input(输入序列长度:N=); %书上取15 n1=input(输入位移:n1=); %书上取3 n=n0:n1+N+5; % 生成自变量数组 u=[(n-n1)=0]; % 产生单位阶跃序列（u(n-n1)） x1=[(n-n1)=0]-[(n-n1-N)=0]; % 用阶跃序列差产生矩形序列 x2=[(n=n1)(n(N+n1))]; % 用逻辑式产生矩形序列 x3=exp(j*n*pi/8).*x2; % 对复正弦序列加矩形窗（元素群乘） subplot(2,2,1);stem(n,x1,.); xlabel(n);ylabel(x1(n)); % 标注 axis([n0,max(n),0,1]); %定坐标范围 subplot(2,2,3);stem(n,x2,.); xlabel(n);ylabel(x2(n)); % 标注 axis([n0,max(n),0,1]); %定坐标范围;subplot(2,2,2);stem(n,real(x3),.); xlabel(n);ylabel(x3(n)的实部); % 标注 line([n0,max(n)],[0,0]); % 画横轴 axis([n0,max(n),-1,1]); %定坐标范围 subplot(2,2,4);stem(n,imag(x3),.); xlabel(n);ylabel(x3(n)的虚部); % 标注 line([n0,max(n)],[0,0]); % 画横轴 a