用MATLAB编写rect函数表示结构

clear all
close all
clc;
TransPoint=[0.0509586,0.1013709,0.4542960];        
TransPoint_position=TransPoint;
a=(TransPoint_position(3)+TransPoint_position(1))/2;
b=TransPoint_position(3)-TransPoint_position(1);
x=-0.4033374:0.0001:0.4033374;
t=rectpuls((x-a)/b);
y=fft(t); %傅里叶变换函数
plot(abs(y));

想生成周期型的矩形函数，高度设置为pi

使用repmat函数重复矩阵，生成周期方波即可。代码：


clear all
close all
clc;
TransPoint=[0.0509586,0.1013709,0.4542960];        
TransPoint_position=TransPoint;
a=(TransPoint_position(3)+TransPoint_position(1))/2;
b=TransPoint_position(3)-TransPoint_position(1);
x=-0.4033374:0.0001:0.4033374;
t=rectpuls((x-a)/b);
% Repeat
X0 = x(end) - x(1);
X = [x X0+x 2*X0+x  3*X0+x 4*X0+x];
T = pi*repmat(t,1,5);
plot(X,T)
ylim([-0.5 4])

结果：

img