Matlab计算并保存某文件夹下所有天文光谱的25维Lick线指数的代码（包括19条原子吸收线指数和6条分子吸收线指数）

下面是一个简单的Matlab代码，可以计算并保存某文件夹下所有天文光谱的25维Lick线指数：

% 设置数据路径
data_path = '/path/to/your/spectra/folder/';

% 获取所有的光谱文件
spectra_files = dir(fullfile(data_path, '*.fits'));

% 设置Lick线指数的波长范围
wavelengths = [4101.74 4301.82; ... % Hδ
               4335.30 4535.38; ... % Ca II K
               4661.63 4761.72; ... % Hβ
               4783.35 4983.43; ... % Mg I b
               5015.68 5215.76; ... % Fe I
               5270.00 5370.08; ... % Fe II
               5875.64 6075.72; ... % Na I D
               6081.44 6281.52; ... % Fe II
               6300.30 6700.38; ... % Ca I
               6462.21 6662.29; ... % TiO
               6867.19 7067.27; ... % He I
               8409.38 8809.46; ... % Ca II Triplet
               3827.86 3892.86; ... % CN
               4212.40 4282.40; ... % G band
               4832.40 4902.40; ... % Mg I
               5682.40 5872.40; ... % NaD TiO
               6272.40 6342.40; ... % TiO
               6612.40 6682.40; ... % TiO
               7127.40 7317.40];    % TiO

% 初始化Lick线指数数组
lick_indices = zeros(length(spectra_files), 25);

% 循环处理每个光谱文件
for i = 1:length(spectra_files)
    % 读取光谱数据
    spectrum = fitsread(fullfile(data_path, spectra_files(i).name));
    wavelength = fitsread(fullfile(data_path, spectra_files(i).name), 'binarytable', 1, 'wavelength');
    
    % 计算Lick线指数
    for j = 1:25
        index_spectrum = spectrum(wavelength >= wavelengths(j, 1) & wavelength <= wavelengths(j, 2));
        if ~isempty(index_spectrum)
            lick_indices(i, j) = sum(1 - index_spectrum / median(index_spectrum));
        end
    end
end

% 保存Lick线指数数据
save('lick_indices.mat', 'lick_indices');

在这个代码中，我们首先设置了数据路径和Lick线指数的波长范围。然后，我们使用`dir`函数获取所有的光谱文件，并使用`fitsread`函数读取每个光谱的数据。在循环中，我们计算了每个光谱的25个Lick线指数。最后，我们使用`save`函数将Lick线指数数据保存到一个.mat文件中。