t = -80018;g2 = zerossizet;g2abst=2 = 1;f = g2;subplot211plottftitle门信号ft=g2txlabeltylabelftF = fftflengthf;frequencies = linspace-11lengthf;subplot212plotfrequenciesabsfftshiftFtitle傅立叶变换xlabel频率ylab

t = -8:0.01:8; %生成时间序列t，从-8到8，步长为0.01 g2 = zeros(size(t)); %生成一个与t相同大小的全0向量g2 g2(abs(t)<=2) = 1; %将g2在绝对值小于等于2的位置上赋值为1，形成矩形门信号 f = g2; %将门信号赋值给f

subplot(2,1,1) %将图像分为2行1列，当前绘制第1个 plot(t,f) %绘制门信号f(t)随时间的变化图像 title('门信号f(t)=g2(t)') %设置标题 xlabel('t') %设置x轴标签 ylabel('f(t)') %设置y轴标签

F = fft(f)/length(f); %对门信号进行傅里叶变换，并除以信号长度作为归一化系数，得到频域信息F frequencies = linspace(-1,1,length(f)); %生成频率序列，从-1到1，长度与门信号长度相同 subplot(2,1,2) %将图像分为2行1列，当前绘制第2个 plot(frequencies,abs(fftshift(F))) %绘制傅里叶变换后频域信息的幅度谱随频率的变化图像 title('傅立叶变换') %设置标题 xlabel('频率') %设置x轴标签 ylabel('|F(f)|') %设置y轴标