内容简介作者简介 - 科学与工程计算系

In D:\MATLAB6p5\work\double_int.m at line 8
SS =
9.5810
(4)
Ssym=vpa(int(int('x^2+y^2','x','sqrt(y)',2),'y',1,4))
Ssym =
9.5809523809523809523809523809524
4.8.3 卷积
4.8.3.1 “完整”离散序列的数值卷积
4.8.3.2 “截尾”离散序列的数值卷积
4.8.3.3 多项式乘法与离散卷积的算法同构
4.8.3.4 连续时间函数的数值卷积
4.8.3.5 卷积的 MATLAB 实现
⎧ 1 n = 3,
4,
, 12 ⎧ 1 n = 2,
3,
, 9
【例 4.8.3.5-1】有序列 A(
n)
= ⎨
和 B(
n)
= ⎨
。
⎩0
else
⎩0
else
(A)
%
a=ones(1,10);n1=3;n2=12;
b=ones(1,8);n3=2;n4=9;
c=conv(a,b);nc1=n1+n3;nc2=n2+n4;
kc=nc1:nc2;
%
aa=a(1:6);nn1=3;nn2=8;
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
10
9
8
7
6
5
4
3
2
1
0
4 6 8 10 12 14 16 18 20 22
图 4.8-2
16