我现在有个命名为'1.txt'的txt文本，里面有我的原始光谱数据，一共155列，第一列第一个单元格是波长命名，第一列剩余的单元格表示的是350到2500nm的波长编号，剩下的2到155列的每一列数据都表示一个样本在350到2500nm的原始光谱反射率，第一行的2到155单元格表示的是每条原始光谱数据的编号，现在我需要用matlb对'1.txt'进行分数阶微分处理，我需要输出每个样本的原始光谱反射率的0.2,0.4,0.6,0.8，1,1.2,1.4,1.6,1.8,2阶分数阶微分数据结果，请给我一段MATLAB代码处理这个问题，并且代码里面如果有自定义函数，请帮我定义好

以下是一个处理您描述的具体问题的Matlab代码：

首先，我们可以定义一个用于计算分数阶微分的函数 `fracdiff()`，该函数在上一个问题中已经提到过：

function df = fracdiff(f,alpha)
% Compute the fractional derivative of a function
syms x;
df = simplify(diff(f,x,alpha),'Steps',100);
end

然后，我们可以读取并解析 '1.txt' 文件，提取出每个样本的原始光谱反射率数据：

filename = '1.txt';
delimiter = '\t';
headerlines = 1;
data = importdata(filename, delimiter, headerlines);
wavelength = data.data(:,1);
reflectance = data.data(:,2:end);

接下来，我们可以计算每个样本的各阶分数阶微分数据，并将结果保存在一个矩阵中：

alpha = [0.2,0.4,0.6,0.8,1,1.2,1.4,1.6,1.8,2];
n_samples = size(reflectance,2);
n_alpha = length(alpha);
diff_data = zeros(length(wavelength), n_samples, n_alpha);
for i = 1:n_samples
    for j = 1:n_alpha
        f = reflectance(:,i);
        df = fracdiff(f,alpha(j));
        diff_data(:,i,j) = eval(df);
    end
end

最后，我们可以将结果保存在一个新的文本文件中：

output_filename = 'diff_data.txt';
output_data = [wavelength, diff_data(:)];
output_header = {'Wavelength', 'Sample', 'Alpha', 'Reflectance_Diff'};
output_fmt = repmat('%f\t', 1, size(output_data,2));
output_fmt = [output_fmt(1:end-1), '\n'];
fid = fopen(output_filename, 'w');
fprintf(fid, '%s\t', output_header{:});
fprintf(fid, '\n');
fprintf(fid, output_fmt, output_data');
fclose(fid);

完整代码如下：