书写一个matlab程序，要求如下：有一个文件其中包含很多小文件需要提取每个文件的第一、二、四列并且根据文件的一系列的区别进行分类，同一类的写入一个文件。提取结束后计算平均值，要求如下：提取好的文件是每15分钟记录的一个数据共一年的数据，其中有缺失的数据在计算的时候需要去掉以后再计算，需要计算每日的平均值和每月的平均值，计算结束后的平均值需要表明是哪一天的日平均值和那个月的月平均值，计算完的根据最初的文件分类方法再分类写入文件，文件名是第一列的名字

以下是一份示例 Matlab 程序，它实现了上述要求：

% 读取文件
data = readtable('data.csv');

% 提取前 4 列，并对文件分类
col1 = data{:, 1};
col2 = data{:, 2};
col4 = data{:, 4};

% 分组，每组存在同一个文件
unique_col1 = unique(col1);
for i = 1:length(unique_col1)
    index = strcmp(col1, unique_col1(i));
    file_data = [col2(index), col4(index)];
    
    % 计算平均值
    [daily_mean, monthly_mean] = calculate_mean(file_data);
    
    % 计算完的平均值再分类写入文件
    filename = unique_col1{i};
    save_to_file(daily_mean, monthly_mean, filename);
end

% 计算每日的平均值和每月的平均值
function [daily_mean, monthly_mean] = calculate_mean(file_data)
    % 处理缺失值
    file_data(isnan(file_data)) = [];
    
    % 计算每日平均值
    daily_mean = mean(file_data, 1);
    
    % 计算每月平均值
    monthly_mean = mean(file_data, 2);
end

% 写入文件
function save_to_file(daily_mean, monthly_mean, filename)
    fid = fopen(strcat(filename, '.txt'), 'w');
    fprintf(fid, '日平均值：%f，月平均值：%f\n', daily_mean, monthly_mean);
    fclose(fid);
end

以上程序读入一个名为 "data.csv" 的文件，提取其中的第一、二、四列。然后，根据第一列的不同值对文件分类，对每一组的数据计算日平均值和月平均值，最后将计算结果