具体内容:利用快速傅里叶变换函数(fft)计算循环卷积,找出时域混叠的条件;
两原信号:1.中心点在1,宽度为3 2.中心点在4,宽度为2;
时域混叠条件:时间长度T<宽度1+宽度2时,会产生混叠,代码中T采用4(<3+2);
%1.6循环卷积
clc;
clear;
close all;
Ts=0.01;
T=4;
T1S=-1;
T2S=2;
w1=3;
w2=2;
c1=1;
c2=4;
T1E=T1S+T;
T2E=T2S+T;
t1=T1S:Ts:T1E;
t2=T2S:Ts:T2E;
x1=rectpuls(t1-c1,w1);
x2=rectpuls(t2-c2,w2);
T3S=T1S+T2S;
T3E=T1E+T2S;
t3=T3S:Ts:T3E;
F1=fft(x1)*Ts;
F2=fft(x2)*Ts;
y=ifft(F1.*F2)/Ts;
subplot(3,1,1);
plot(t1,x1);
axis([T1S-1,T1E+1,-0.5,1.5]);
subplot(3,1,2);
plot(t2,x2);
axis([T2S-1,T2E+1,-0.5,1.5]);
subplot(3,1,3);
plot(t3,y);
axis([T3S-1,T3E+1,-0.5,4]);