matlab实现周期矩形波形信号的单边频谱

下面是MATLAB实现周期矩形波形信号的单边频谱的代码：

% 生成周期矩形波形信号
T = 1; % 周期
t = linspace(0, T, 1000); % 时间序列
f = 1/T; % 基波频率
x = square(2*pi*f*t); % 周期矩形波形信号

% 计算单边频谱
N = length(x); % 采样点数
X = fft(x)/N; % 快速傅里叶变换
X = X(1:N/2+1); % 取单边频谱
X(2:end-1) = 2*X(2:end-1); % 奇异点修正
f = (0:N/2)*(1/T)/N; % 频率序列

% 绘制结果
subplot(2,1,1);
plot(t, x);
xlabel('时间');
ylabel('幅值');
title('周期矩形波形信号');

subplot(2,1,2);
plot(f, abs(X));
xlabel('频率');
ylabel('幅值');
title('单边频谱');

运行代码后，会得到一个绘制了周期矩形波形信号和单边频谱的图像。其中，周期矩形波形信号在上半部分绘制，单边频谱在下半部分绘制。可以看到，单边频谱只包含了正频率部分，并且经过了奇异点修正，得到了正确的幅值大小。