实验四 离散傅里叶变换的性质
离散傅里叶变换性质
X [m] X [ N m]
实序列 实部周期偶对称，虚部 周期奇对称
当x[k]是实序列周期偶对称时：X [m] X [ N m]
实序列，周期偶对称 实序列，周期偶对称
当x[k]是实序列周期奇对称时： X [m] X [ N m]
实序列，周期奇对称 纯虚序列，周期奇对称
第2章 离散傅里叶变换(DFT)
问题的提出
有限长序列的傅里叶分析
离散傅里叶变换的性质 利用DFT计算线性卷积 利用DFT分析信号的频谱
离散傅里叶变换的性质
1. 线性
2. 循环位移
3. 对称性 4. 序列的循环卷积 5. 序列DFT与z变换的关系
离散傅里叶变换的性质
1. 线性
X1[m] DFTx1[k ]
k
-5 -4 -3 -2 -1 0 1 2 3 4 5 6 7
x[(k 2) 5 ]R N [k ]
k
0 1 2 3 4
DFT时域循环位移特性
mn ~ ~ DFS {x[k n]} WN X [m]
DFTx[( k n) N ]RN [k ]
mn WN X [m]
时域的循环位移对应频域的相移
DFT{ x [k ]} X [(m)N ]RN [m] X [ N m]
时域共轭 频域周期共轭 ~ DFT{x [(k ) N ]RN [k ]} X [m] 时域周期共轭 频域共轭
DFT性质
离散傅里叶变换的性质
3. 对称性 (symmetry)
当x[k]是实序列时：
DFT性质
离散傅里叶变换的性质
序列的循环卷积步骤： （1）补零
（2）周期延拓
离散傅里叶变换的性质
DFT [ x1 (n)] X 1 (k )
DFT [ x2 (n)] X 2 (k )
若 则
Y (k ) X 1 (k ) X 2 (k )
y (n) IDFT [Y (k )] [ x1 (m) x2 ((n m)) N ]RN ( n)
m 0 N 1
圆周卷积和
——电子信息工程
纯虚序列的共轭对称性 序列
Re[ x ( n )] 0 j Im[ x ( n )]
DFT
X ep ( k ) 0 X op ( k ) X ( k )
X ( k ) X (( N k )) N R N ( k )
xep ( n ) xop ( n ) Re[ X ( k )] j Im[ X ( k )]
*
1/ 2[ X (( k )) N X (( N k )) N ]RN ( k )
*
X op (k ) X op (( N k )) N RN (k )
*
1/ 2[ X (( k )) N X (( N k )) N ]RN ( k )
*
——电子信息工程
共轭对称性
~ x (n 2)R N (n)
x2 (n)
——电子信息工程
DFT 的性质: 圆周移位(续)
N=15
x(n)
~ x (n)
~ x ( n 13 )
~ x ( n 13 ) R 15 ( n )
——电子信息工程
X m (k ) DFT [ xm (n)] DFT [ x((n m)) N RN (n)]
1 N 1
N 1
2
(n m) x2 (( n m )) N (( n m )) N
§4.11 离散傅里叶变换及其性质
f N(k m) e
i m
j
2 kn N
N m 1
[
j
f N(i) e
2 2 in j mn N ]e N
j
2 in N
都是以N为周期的函数，因此
2 in N
j
f N(i) e
N 1 i 0
f N(i) e
j
2 in N
解 F (n) = DFT[f(k)] = 仅当n=0时，
j 2n N
e
1
F (0) =N
当n=1，2，…，N-1时，
j
2n j 2n N )N j 1 (e , (e N 1) 2n j 1 (e N ) 2n j N, (e N 1)
若
则 证明
f(k)←→ F(n)
f ((k –m))NGN(k) ←→ WmnF(n)
▲ ■ 第 7页
DFT时移特性证明
DFT[ f ((k –m))NGN(k)]=DFT[ fN (k –m)GN(k)]
N 1 k 0
令i=k-m，有 DFT[ f ((k –m))NGN(k)]= 由于fN (k )和 e
■
第 1页
一．离散傅里叶变换(DFT)
借助周期序列DFS的概念导出有限长序列的DFT。 将有限长序列f(k)延拓成周期为N的周期序列fN(k)
f (k)
f N (k )
l
f (k lN )
o
(a)
fN(k)
N-1
k
o
(b)
N-1 N
2N-1 k
F (n) DFT[ f (k )]
实验离散傅里叶变换的性质及应用
实验一离散傅里叶变换地性质及应用一、实验目地1．了解DFT 地性质及其应用2．熟悉MATLAB 编程特点二、实验仪器及材料计算机，MATLAB 软件三、实验内容及要求1．用三种不同地DFT 程序计算)()(5n R n x =地256点离散傅里叶变换()X k ，并比较三种程序计算机运行时间.（1）编制用for loop 语句地M 函数文件dft1.m ，用循环变量逐点计算()X k ；（2）编写用MATLAB 矩阵运算地M 函数文件dft2.m ，完成下列矩阵运算：00000121012(1)(1)(1) (0)(0) (1)(1) (1)(1) N N N N N N N N N N N N N N N N N X x W W W W X x W W W W x N X N W W W W -----⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ （3）调用fft 库函数，直接计算()X k ；（4）分别调用上述三种不同方式编写地DFT 程序计算序列()x n 地离散傅里叶变换()X k ，并画出相应地幅频和相频特性，再比较各个程序地计算机运行时间.2．研究实序列地DFT 特点及性质.已知序列)()/2cos(5)(n R N rn n x N π=，若128=N 计算下列三种情况下序列地DFT 地幅度、实部及虚部，并用图形表示相应地[][])(Im ,)(Re ,)(),(k X k X k X n x .（1）)/2cos(5)(N rn n x π=，5=r（2）)/2cos(5)(N rn n x π=，4.4=r（3）)/2cos(5)(N rn n x π=，0=r3．利用DFT 实现两序列地卷积运算，并研究DFT 点数与混叠地关系.（1）已知两序列)()(3n R n x =)()1()(5n R n n h +=（2）用直接法（即用线性卷积地定义计算，见下式）计算线性卷积y (n )=x (n )*h (n )地结果，并以图形方式表示结果；20),()()(10-+≤≤-⨯=∑-=M N n m n h m x n y N m 其中：序列)1N n 0(),n (x -≤≤和序列)1M n 0(),n (h -≤≤（3）用MATLAB 编制利用DFT 计算线性卷积y (n )=x (n )*h (n )地程序；分别令圆周卷积地点数为L=5，6，7，8，以图形方式表示结果.（4）对比直接法和圆周卷积法所得地结果.clc;clear all;close all;N=256;x=[ones(1,8),zeros(1,N-8)];n=[0: (length(x)-1)]t=cputime;[Am1,Pha1]=dft1(x);t1=cputime-t,t=cputime;[Am2,Pha2]=dft2(x);t2=cputime-t,t=cputime;[Am3,Pha3]=dft3(x);t3=cputime-t,sizex=4;sizeh=6;x=[1,2,3,4]h=[0,1,0,0,0,-2];y=conv(x,h);figure(1)subplot(2,3,1)le1=0:1:sizex-1;stem(le1,x,'r');title('x sequence')ylabel('x(n)')subplot(2,3,2)le2=0:1:sizeh-1;stem(le2,h);title('h sequence')ylabel('h(n)')subplot(2,3,3)le3=0:1:sizex+sizeh-2;stem(le3,y,'g')title('y sequence by direct metod')ylabel('y(n)')LL=[6 9 12];for n=1:3L=LL(n);X=fft(x,L);H=fft(h,L);Y=X.*HyFFT=ifft(Y,L);subplot(2,3,n+3)le3=0:1:L-1;stem(le3,yFFT);if (n==1)title('y sequence by FFT L=6')elseif (n==2)title('y sequemce by FFT L=9')elseif (n==3)title('y sequence by FFT L=12')endendr=4.4N=128;n=0:N-1;x=5*cos(2*pi*r*n/N);f1=1000*r/N,f2=1000*0/N,f3=1000*4.4/NX=fft(x,N);figure(2)subplot(221);stem(n,x);subplot(222);stem(n,abs(X));subplot(223);stem(n,real(X));subplot(224);stem(n,imag(X));n =Columns 1 through 220 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21Columns 23 through 4422 23 24 25 26 27 28 29 30 31 32 3334 35 36 37 38 39 40 41 42 43Columns 45 through 6644 45 46 47 48 49 50 51 52 53 54 5556 57 58 59 60 61 62 63 64 65Columns 67 through 8866 67 68 69 70 71 72 73 74 75 76 7778 79 80 81 82 83 84 85 86 87Columns 89 through 11088 89 90 91 92 93 94 95 96 97 98 99100 101 102 103 104 105 106 107 108 109Columns 111 through 132110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131Columns 133 through 154132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153Columns 155 through 176154 155 156 157 158 159 160 161 162 163 164 165166 167 168 169 170 171 172 173 174 175Columns 177 through 198176 177 178 179 180 181 182 183 184 185 186 187188 189 190 191 192 193 194 195 196 197Columns 199 through 220198 199 200 201 202 203 204 205 206 207 208 209210 211 212 213 214 215 216 217 218 219Columns 221 through 242220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241Columns 243 through 256242 243 244 245 246 247 248 249 250 251 252 253254 255t1 =0.7644t2 =0.6396t3 =x =1 2 3 4Y =-10.0000 -9.5000 +11.2583i 3.5000 - 6.0622i -2.0000 3.5000 + 6.0622i -9.5000 -11.2583iY =Columns 1 through 7-10.0000 -7.4363 -21.7778i 4.5753 - 1.6805i 3.5000 -6.0622i 4.3610 + 1.9108i 4.3610 - 1.9108i 3.5000 + 6.0622iColumns 8 through 94.5753 + 1.6805i -7.4363 +21.7778iY =Columns 1 through 7-10.0000 14.7942 -17.6244i -9.5000 +11.2583i -2.0000 -2.0000i3.5000 - 6.0622i -0.7942 + 6.6244i -2.0000Columns 8 through 12-0.7942 - 6.6244i 3.5000 + 6.0622i -2.0000 + 2.0000i -9.5000-11.2583i 14.7942 +17.6244ir =4.4000f1 =34.3750f2 =f3 =34.3750版权申明本文部分内容，包括文字、图片、以及设计等在网上搜集整理.版权为个人所有This article includes some parts, including text, pictures, and design. Copyright is personal ownership.sQsAE。
实验四 离散傅里叶变换的性质
实验四离散傅里叶变换的性质一、实验目的1. 熟悉matlab软件中离散傅里叶变换的实现方法及FFT函数的使用方法；2. 通过软件仿真，加深对离散傅里叶变换性质的理解。
二、实验内容1. 验证离散傅里叶变换的线性性质；2. 掌握用matlab实现圆周移位的方法；3. 验证圆周卷积与线性卷积的关系。
三、实验步骤1. 验证线性性质设两个有限长序列分别为xn1=[3,1,-2,2,3,4]，xn2=[1,1,1,1]，做4DFT[xn1]+2DFT[xn2],及DFT[4xn1+2xn2]的运算，比较它们的结果。
代码如下：clear,N=20;n=[0:1:N-1];xn1=[3,1,-2,2,3,4];n1=0:length(xn1)-1; %定义序列xn1xn2=[1,1,1,1];n2=0:length(xn2)-1; %定义序列xn2yn1=4*xn1;yn2=2*xn2;[yn,ny]=seqadd(yn1,n1,yn2,n2); %计算4xn1+2xn2xk1=fft(xn1,N);xk2=fft(xn2,N); %分别求DFT[xn1] 和DFT[xn2]yk0=4*xk1+2*xk2; %计算4DFT[xn1]+2DFT[xn2]yk=fft(yn,N); %计算DFT[4xn1+2xn2]subplot(2,1,1);stem(n,yk0);title('傅里叶变换之和') %显示4DFT[xn1]+2DFT[xn2]subplot(2,1,2);stem(n,yk);title('序列和之傅里叶变换') %显示DFT[4xn1+2xn2]运行结果如图1所示，从图中可知，用两种方法计算的DFT完全相等，所以离散傅里叶变换的线性性质得到验证。
图1 离散傅里叶变换的线性性质2. 圆周移位设x=[7,6,5,4,3,2]，位于主值区间，现要把x分别圆周右移两位，圆周左移1位，成为新主值区间的向量，以此观察圆周移位的特点。
离散傅里叶变换及性质c
思考题解答
（1）实验1中序列1和序列2在N=8时的幅频特性是否相同，为什么？N=16呢？
由以上图像可知序列1和序列2在N=8时的幅频特性相同，因为周期延拓，它们的图像有一部分相同，当抽样段不同他们的幅频特性就可以相同。当N=16时是周期的整数倍，此时它们的幅频特性就不相同，如1.1.1图和1.1.2图所示。
N=length(xn);
subplot(211);stem([0:N-1],xn);
grid;
X=zeros(1,N);
for k=0:N-1
for n=0:N-1
X(k+1)=X(k+1)+xn(n+1)*exp(-j*2*pi*n*k/N);%按定义式计算序列的DFT
end;
end;
subplot(212);stem([0:N-1],abs(X));
2、抽样频率应为信号频率的整数倍，且倍数应不小于3；
3、由于在用快速傅里叶变换实现信号的频谱分析时往往要求数据的长度是2的整次幂，所以抽样频率应取信号频率的4倍，即每个周期抽四个点；
4、对抽样后的数据截短时，应使截短后的数据包含整周期；
5、在截短后的数据后不宜补零以避免频谱的泄露。
分析补零对频谱的影响：1、能使序列的频谱变得细致,但不能提高序列的频率分辨率,只有采集更多的有效数据,才能得到序列的高分辨率频谱；2、补零使谱线插值式密化对提高谱分析精度也有一定好处；3、插零使谱线周期性密化，可方便地实现基带带宽压缩和多采样；
3、fs和N确定后，即可确定所需相应模拟信号x(t)的长度T=N/fs.分辨率变化量f反比于T，而不是N，在给定的T情况下，靠减小Ts来增加N时不能提高分辨率的，因为T=NTs为常数。
离散傅里叶变换及其特性验证
实验名：离散傅里叶变换及其特性验证一、实验目的1、掌握离散时间傅立叶变换（DTFT ）的计算方法和编程技术。
2、掌握离散傅立叶变换（DFT ）的计算方法和编程技术。
3、理解离散傅立叶变换（DFT ）的性质并用MA TLAB 进行验证。
二、实验原理与计算方法1、离散时间傅立叶变换如果序列x (n )满足绝对可和的条件，即∞<∑∞-∞=n n x |)(|，则其离散时间傅立叶变换定义为： ∑∞-∞=-==n nj j en x n x F e X ωω)()]([)( （1）假设序列x (n )在N n n n ≤≤1(即不一定在[0, N -1])有N 个样本，要估计下列各点上的X (e j ω)：M k k Mk ...,2,1,0==， πω它们是[0,π]之间的(M +1)个等间隔频点，则（1）式可写成： M k n x ee X Nl l kn Mjj l...,2,1,0)()(1==∑=-， πω（2）将{x (n l )}和{X (e j ωk)}分别排列成向量x 和X ，则有：X=Wx （3） 其中W 是一个(M +1)×N 维矩阵：⎭⎬⎫⎩⎨⎧=≤≤=-M k n n n e N kn M j ...,2,1,0;1， πW将{k }和{n }排成列向量，则⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=n k W T M j πexp 在MA TLAB 中，把序列和下标排成行向量，对（3）式取转置得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=k n x X T T T M j πexp其中n T k 是一个N ×(M +1)维矩阵。
用MATLAB 实现如下：k=[0:M]; n=[n1:n2];X=x*(exp(-j*pi/M)).^(n ’*k); 2、离散傅立叶变换一个有限长序列的离散傅立叶变换对定义为：10,)()(10-≤≤=∑-=N k W n x k X N n nk N （4）10,)(1)(1-≤≤=∑-=-N n Wk X Nn x N k kn N （5）以列向量x 和X 形式排列x (n )和X (k )，则式（4）、（5）可写成： X =W N x可由下面的MA TLAB 函数dft 和idft 实现离散傅立叶变换运算。
离散傅里叶变换及其性质
离散傅⾥叶变换及其性质1 ⼀维与⼆维离散傅⾥叶变换以周期对函数 f(t) 采样可表⽰为，对采样函数进⾏傅⾥叶变换得，整理得。
由于对函数 f(t) 的采样周期为，采样函数的傅⾥叶变换的⼀个完整周期为，同样的，也是采样函数的傅⾥叶变换的⼀个完整周期，只是这个周期不是以原点对称的。
在区间中取 M 个点，则第 m 个点的频率为，带⼊公式得，其中，为连续函数 f(t) 对应的 M 个离散值，为取样函数的傅⾥叶变换对应的 M 个离散值，整理公式得（由于函数仅在 [0,M-1] 上有⾮零值，故真实求和区间为 [0,M-1]）。
因此，⼀维离散傅⾥叶变换对为，。
类似的，⼆维离散傅⾥叶变换对为，。
2 傅⾥叶变换的性质1）傅⾥叶变换平移特性，⽤指数项乘以 f(t) 使得傅⾥叶变换后原点移动到处，使⽤负指数乘以使得反傅⾥叶变换后原点移动到处，证明如下：，使⽤替换得，因此有，类似推导可得。
将平移特性扩展到⼆维离散变量上有。
2）离散傅⾥叶变换⼀定具有周期特性，因为离散傅⾥叶变换的频率取值在区间内，有限频率导致必然具有周期性，连续傅⾥叶变换频率取值为⽆穷⼤，所以连续傅⾥叶变换⼀般不具有周期性（但也有所有频率都⼀样的函数）。
离散傅⾥叶变换周期性可表⽰为。
观察公式 或，发现频率取值在之间，⽽⼀个完整的频率应该在之间，如下图：如果直接应⽤公式进⾏傅⾥叶变换，得到的频率为 [0,M-1]区间，这是两个半周期组成的⼀个周期。
在图像中则表现为低频信号分布在4个⾓落，这显然不便于观察频率信息。
结合傅⾥叶变换的平移特性，可以将原函数乘以⼀个正指数项，使得平移后傅⾥叶变换再 [0,M-1]区间正好是⼀个完整的周期。
将原函数平移 M/2 可以实现该⽬标，具体分析如下： 原函数平移 M/2 得 ，由于 x 为⾮负整数，,最终得到。
对于⼆维离散变量有相似结论 。
3）原函数（⼆维及以上）旋转⼀定⾓度，其傅⾥叶变换也旋转对应⾓度。
令 为原函数变量的列向量， 为傅⾥叶变换函数变量的列向量，对的傅⾥叶变换可表⽰为，对 旋转⼀定⾓度可表⽰为，其中 R 为旋转矩阵，对 的傅⾥叶变换可表⽰为 ，由 得 ，并将其带⼊上式得，由于，因此 ，使得傅⾥叶变换旋转相应⾓度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验四离散傅里叶变换的性质
一、实验目的
1. 熟悉matlab软件中离散傅里叶变换的实现方法及FFT函数的使用方法；
2. 通过软件仿真，加深对离散傅里叶变换性质的理解。
二、实验内容
1. 验证离散傅里叶变换的线性性质；
2. 掌握用matlab实现圆周移位的方法；
3. 验证圆周卷积与线性卷积的关系。
三、实验步骤
1. 验证线性性质
设两个有限长序列分别为xn1=[3,1,-2,2,3,4]，xn2=[1,1,1,1]，做4DFT[xn1]+2DFT[xn2],及DFT[4xn1+2xn2]的运算，比较它们的结果。
代码如下：
clear,N=20;n=[0:1:N-1];
xn1=[3,1,-2,2,3,4];n1=0:length(xn1)-1; %定义序列xn1
xn2=[1,1,1,1];n2=0:length(xn2)-1; %定义序列xn2
yn1=4*xn1;yn2=2*xn2;[yn,ny]=seqadd(yn1,n1,yn2,n2); %计算4xn1+2xn2
xk1=fft(xn1,N);xk2=fft(xn2,N); %分别求DFT[xn1] 和DFT[xn2]
yk0=4*xk1+2*xk2; %计算4DFT[xn1]+2DFT[xn2]
yk=fft(yn,N); %计算DFT[4xn1+2xn2]
subplot(2,1,1);stem(n,yk0);title('傅里叶变换之和') %显示4DFT[xn1]+2DFT[xn2]
subplot(2,1,2);stem(n,yk);title('序列和之傅里叶变换') %显示DFT[4xn1+2xn2]
运行结果如图1所示，从图中可知，用两种方法计算的DFT完全相等，所以离散傅里叶变换的线性性质得到验证。
图1 离散傅里叶变换的线性性质
2. 圆周移位
设x=[7,6,5,4,3,2]，位于主值区间，现要把x分别圆周右移两位，圆周左移1位，成为新主值区间的向量，以此观察圆周移位的特点。
代码如下：
x=[7,6,5,4,3,2];Nx=length(x);nx=0:Nx-1; % 给出x序列
nx1=-Nx:2*Nx-1;x1=x(mod(nx1,Nx)+1); % 延拓为周期向量x1,注意mod用法
[y1,ny1]=seqshift(x1,nx1,2); % 将x1右移两位，得到y1
[y2,ny2]=seqshift(x1,nx1,-1); % 将x1左移三位，得到y2
RN=(nx1>=0)&(nx1<Nx); % 在x1的位置向量nx1上设置主值窗
RN1=(ny1>=0)&(ny1<Nx); % 在y1的位置向量ny1上设置主值窗
RN2=(ny2>=0)&(ny2<Nx); % 在y2的位置向量ny2上设置主值窗
subplot(6,1,1),stem(nx1,RN.*x1) % 在子图上画出x1的主值
title('主值序列x')
axis([-6,15,0,10])
subplot(6,1,2),stem(nx1,x1) % 画出x1
title('周期序列x1')
axis([-6,15,0,10])
subplot(6,1,3),stem(ny1,y1) % 画出y1
title('移位周期序列y1')
axis([-6,15,0,10])
subplot(6,1,4),stem(ny2,y2) % 画出y2
title('移位周期序列y2')
axis([-6,15,0,10])
subplot(6,1,5),stem(ny1,RN1.*y1) % 画出y1的主值
title('y1的主值序列y')
axis([-6,15,0,10])
subplot(6,1,6),stem(ny2,RN2.*y2) % 画出y2的主值
title('y2的主值序列y')
axis([-6,15,0,10])
运行结果如图2所示。
图2 序列的圆周移位
3. 圆周卷积和与线性卷积和
设x1=[1,2,3,0]；x2=[5,4,-3,-2]；分别求它们的N点圆周卷积和与线性卷积和，并做比较。
代码如下：
圆周卷积和函数
function y = circonv(x,h,N)
if length(x) > N
error('N 必须>= x的长度')
end
% 检查h的长度
if length(h) > N
error('N 必须>= h的长度')
end
x=[x,zeros(1,N-length(x))]; % 将x的长度扩展至N
h=[h,zeros(1,N-length(h))]; % 将h的长度扩展至N
m = [0:N-1];
hm = h(mod(-m,N)+1); % 将h循环折叠
H = toeplitz(hm,[0,h(2:N)]); % 用toeplitz函数产生循环卷积矩阵
y = x*H; % 用向量－矩阵乘法求卷积
x1与x2的N（N分别取4、7、10）点圆周卷积和及线性卷积和
x1=[1,2,3,0];x2=[5,4,-3,-2];
y1=circonv(x1,x2,4)
y2=circonv(x1,x2,7)
y3=circonv(x1,x2,10)
y4=conv(x1,x2)
运行结果如下：
y1 =
-8 8 20 4
y2 =
5 14 20 4 -13 -
6 0
y3 =
5 14 20 4 -13 -
6 0 0 0 0
y4 =
5 14 20 4 -13 -
6 0
从计算结果中，可得，当N≥N1+N2-1时，圆周卷积和与线性卷积和相等。