chengxu.zip_MFC_chengxu_convolution_卷积 实现

R6=[ones(1,6)]; %产生长度为6的矩形序列 R8=[ones(1,8)]; %产生长度为6的矩形序列 M1=length(R6); M2=length(R8); L=M1+M2-1; r1=[R6 zeros(1,L-M1)]; %对R6补L-M1个零 r2=[R8 zeros(1,L-M2)]; %对R8补L-M2个零 x1=fft(r1,L); %对r1进行L点快速傅里叶变换 x2=fft(r2,L); %对r2进行L点快速傅里叶变换 x=x1.*x2; %将x1与x2相乘 y=ifft(x,L); %对x进行L点快速傅里叶逆变换 figure; %下一个绘图结果，将在一个新的绘图窗口显示 subplot(4,1,1); stem(0:L-1,y,'.');xlabel('n');ylabel('y(n)'); %绘制线性卷积卷积结果图 title('R6与R8线性卷积结果'); %对R6与R8两矩形序列分别进行8点循环卷积运算 R6=[ones(1,6)]; %产生长度为6的矩形序列 R8=[ones(1,8)]; %产生长度为6的矩形序列 x1=fft(R6,8); %对R6进行8点快速傅里叶变换 x2=fft(R8,8); %对R8进行8点快速傅里叶变换 x=x1.*x2; %将x1与x2相乘 y=ifft(x,8); %对x进行8点快速傅里叶逆变换 text(-2,12,'\color{red}14级通信2班 梁玮 29号') subplot(4,1,2); stem(0:length(y)-1,y,'.');xlabel('n');ylabel('y(n)'); %绘制8点循环卷积的结果图 title('R6与R8进行8点循环卷积结果'); %对R6与R8两矩形序列分别进行14点循环卷积运算 R6=[ones(1,6)]; %产生长度为6的矩形序列 R8=[ones(1,8)]; %产生长度为6的矩形序列 x1=fft(R6,14); %对R6进行14点快速傅里叶变换 x2=fft(R8,14); %对R8进行14点快速傅里叶变换 x=x1.*x2; %将x1与x2相乘 y=ifft(x,14); %对x进行14点快速傅里叶逆变换 subplot(4,1,3); stem(0:length(y)-1,y,'.');xlabel('n');ylabel('y(n)'); %绘制14点循环卷积的结果图 title('R6与R8进行14点循环卷积结果'); %对R6与R8两矩形序列分别进行16点循环卷积运算 R6=[ones(1,6)]; %产生长度为6的矩形序列 R8=[ones(1,8)]; %产生长度为6的矩形序列 x1=fft(R6,16); %对R6进行16点快速傅里叶变换 x2=fft(R8,16); %对R8进行16点快速傅里叶变换 x=x1.*x2; %将x1与x2相乘 y=ifft(x,16); %对x进行16点快速傅里叶逆变换 subplot(4,1,4); stem(0:length(y)-1,y,'.');xlabel('n');ylabel('y(n)'); %绘制16点循环卷积的结果图 title('R6与R8进行16点循环卷积结果');