请给我一段matlab代码，输入为一组采样自一个周期三角波的离散采样点数据，输出原三角波的基波与谐波

在MATLAB中，你可以使用傅里叶变换（FFT）来分析离散时间信号，包括提取基波（直流分量）和高阶谐波。以下是一个简单的示例，假设你有一个名为`sample_data`的一维数组，代表周期三角波的离散采样数据：

% 假设 sample_data 是你的离散样本数据
Fs = 100; % 采样频率
T = 1/Fs; % 采样间隔

% 对数据进行离散傅立叶变换 (DFT)
N = length(sample_data);
Y = fft(sample_data); % DFT结果
Y_mag = abs(Y/N); % 计算幅度谱，除以长度纠正幅度

% 提取基波 (DC成分) 和第一、二阶谐波
frequencies = (0:N-1)*(Fs/N); % 频率轴
fundamental_freq = Fs / 2; % 三角波基频

% 选取基波及其前两个谐波
[~, index] = min(abs(frequencies - fundamental_freq)); % 找到基频索引
harmonics_index = [index, 2*index+1, 2*index+2]; % 第一、二、三阶谐波索引

% 提取基波和谐波部分
base_wave = Y_mag(harmonics_index);
harmonic_components = Y_mag([0, harmonics_index]);

% 可能还需要对结果进行归一化或其他处理，具体取决于需求
% ...

% 显示或保存结果
figure;
plot(frequencies, Y_mag, 'b', frequencies(harmonics_index), harmonic_components, 'ro');
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Spectrum of the sampled triangle wave');
grid on;

% 或者将结果存储为变量供后续分析
freq_axis = frequencies(harmonics_index);
output = struct('BaseWave', base_wave, 'HarmonicComponents', harmonic_components, 'FreqAxis', freq_axis);