%方法一利用square函数构造周期方波利用fft实现卷积
A=2;
TS=0.01;
f0=1/4;
w0=2*pi*f0;
t=0:TS:(5/f0-TS);
y1=A*square(w0*t)+A;%4为周期的周期方波
y2=rectpuls(t-0.5,1);%单个方波
F1=fft(y1)*TS;
F2=fft(y2)*TS;
yy1=ifft(F1.*F2)/TS;
%两个图像呈现
subplot(3,3,1);
plot(t,y1);
axis([0,(5/f0-TS+1),-1,5]);
title('周期为4有值部分为2幅值为4的周期方波')
subplot(3,3,2);
plot(t,y2);
axis([0,(5/f0-TS+1),-1,5]);
title('周期为4有值部分为1的主值方波')
%图像显示
subplot(3,3,3);
plot(t,yy1);
axis([-1,(5/f0-TS+1),-1,5]);
title('周期卷积结果')
%方法一利用square函数构造周期方波利用fft实现卷积临界混叠
A=2;
TS=0.01;
f0=1/4;
w0=2*pi*f0;
t=0:TS:(5/f0-TS);
y1=A*square(w0*t)+A;%2为周期的周期方波
y2=rectpuls(t-1,2);%单个方波
F1=fft(y1)*TS;
F2=fft(y2)*TS;
yy1=ifft(F1.*F2)/TS;
%两个图像呈现
subplot(3,3,4);
plot(t,y1);
axis([0,(5/f0-TS+1),-1,5]);
title('周期为4有值部分为2幅值为4的周期方波')
subplot(3,3,5);
plot(t,y2);
axis([0,(5/f0-TS+1),-1,5]);
title('周期为4有值部分为2的主值方波')
%图像显示
subplot(3,3,6);
plot(t,yy1);
axis([-1,(5/f0-TS+1),-1,10]);
title('周期卷积临界混叠结果')
%方法一利用square函数构造周期方波利用fft实现卷积混叠
A=2;
TS=0.01;
f0=1/4;
w0=2*pi*f0;
t=0:TS:(5/f0-TS);
y1=A*square(w0*t)+A;%2为周期的周期方波
y2=rectpuls(t-1.5,3);%单个方波
F1=fft(y1)*TS;
F2=fft(y2)*TS;
yy1=ifft(F1.*F2)/TS;
%两个图像呈现
subplot(3,3,7);
plot(t,y1);
axis([0,(5/f0-TS+1),-1,5]);
title('周期为4有值部分为2幅值为4的周期方波')
subplot(3,3,8);
plot(t,y2);
axis([0,(5/f0-TS+1),-1,5]);
title('周期为4有值部分为3的主值方波')
%图像显示
subplot(3,3,9);
plot(t,yy1);
axis([-1,(5/f0-TS+1),-1,10]);
title('周期卷积混叠结果')
