matlab中用FFT实现线性卷积循环卷积
3星 · 超过75%的资源 需积分: 43 178 浏览量
2014-06-12
09:37:58
上传
评论 8
收藏 44KB DOC 举报 
用FFT实现序列 和 的线性卷积和
N1=-5:11;
x1=(0.8.^N1).*[zeros(1,5),ones(1,12)];
x2=[zeros(1,5),ones(1,6),zeros(1,6)];
lenx=length(x1);
lenh=length(x2);
N=lenx+lenh-1;
x1=[x1 zeros(1,N-lenx)];
subplot(221);
stem(x1,'.')
x2=[x2 zeros(1,N-lenh)];
subplot(222);
stem(x2,'.')
N=2^(ceil(log2(N)));
X2=fft(x2);
X1=fft(x1);
Y=X2.*X1;
y=real(ifft(Y));
subplot(223);
stem(y,'.')
z=conv(x1,x2);
subplot(224);
stem(z,'.')
axis([-5,35,0,4]);
0 10 20 30 40
0
0.2
0.4
0.6
0.8
1
0 10 20 30 40
0
0.2
0.4
0.6
0.8
1
0 10 20 30 40
0
1
2
3
4
0 10 20 30
0
1
2
3
4
X1(n)={1,2,2,1} x2(n)={1,-1,-1,1}
1 求它们的线性卷积 2 计算循环卷积
x1=[1,2,2,1];
x2=[1,-1,-1,1];
lenx=length(x1);
lenh=length(x2);
N=lenx+lenh-1;
x1=[x1 zeros(1,N-lenx)];
subplot(221);
12
( ) 0.8 ( )
n
x n R n
6
( ) ( )h n R n
资源评论
Sawakott2014-11-23要怎么用唉？？
