exm040835_1.m

matlab mfiles实例分析

%完整序列卷积
a=ones(1,10);n1=3;n2=12;			%完整A(n)序列的非平凡值和区间端点
b=ones(1,8);n3=2;n4=9;				%完整B(n)序列的非平凡值和区间端点
c=conv(a,b);nc1=n1+n3;nc2=n2+n4;	%计算卷积和确定卷积非平凡区间端点
kc=nc1:nc2;							%构成非平凡区间的序号自变量
%截尾序列卷积
aa=a(1:6);nn1=3;nn2=8;				%截尾A(n)序列的非平凡值和区间端点
cc=conv(aa,b);ncc1=nn1+n3;
nx=nn2+n4;							%“非平凡”区间右端点
ncc2=min(nn1+n4,nn2+n3);			%截尾序列卷积被正确计算区间的右端点
kx=(ncc2+1):nx;kcc=ncc1:ncc2;N=length(kcc);
stem(kcc,cc(1:N),'r','filled')
axis([nc1-2,nc2+2,0,10]),grid,hold on
stem(kc,c,'b'),stem(kx,cc(N+1:end),'g'),hold off