将100条DNA序列txt文件one-hot编码的MATLAB代码

假设你有一个包含100条DNA序列的txt文件，每个序列占一行。下面是一个MATLAB示例代码，用于将这些DNA序列进行one-hot编码：

% 读取DNA序列文件
dna_file = 'dna_sequences.txt';
dna_seqs = readcell(dna_file);

% 将DNA序列转换为数字序列
num_seqs = cellfun(@nt2int, dna_seqs, 'UniformOutput', false);

% 定义one-hot编码矩阵
one_hot = zeros(length(num_seqs), length(num_seqs{1}), 4);

% 对每个数字进行one-hot编码
for i = 1:length(num_seqs)
    for j = 1:length(num_seqs{1})
        if num_seqs{i}(j) == 1
            one_hot(i, j, :) = [1 0 0 0]; % A
        elseif num_seqs{i}(j) == 2
            one_hot(i, j, :) = [0 1 0 0]; % C
        elseif num_seqs{i}(j) == 3
            one_hot(i, j, :) = [0 0 1 0]; % G
        elseif num_seqs{i}(j) == 4
            one_hot(i, j, :) = [0 0 0 1]; % T
        end
    end
end

这个示例中，我们首先读取DNA序列文件，并将每个序列转换为数字序列。接着，我们定义了一个3维矩阵，用于存储one-hot编码后的结果，其中第一维表示序列的编号，第二维表示碱基的位置，第三维表示one-hot编码的4个维度。最后，我们遍历每个数字序列中的每个数字，根据数字对应的碱基，将对应的one-hot编码写入one_hot矩阵中。这样，我们就得到了100条DNA序列的one-hot编码。