MATLAB实现移位寄存器产生M序列，要求采用四级寄存器，初始值任选且不全为0，产生序列长度为60.给出.m注释的代码

在MATLAB中，我们可以使用`bitshift`和`xor`操作符来模拟四级移位寄存器（LFSR）并生成M序列。以下是一个简化的.M文件示例，它接受一个初始状态向量作为输入，并生成长度为60的M序列：

% 导入所需功能
function m_seq = generate_M_sequence(initial_state)
    % 定义移位寄存器的反馈多项式系数
    feedback_poly = [15, 17, 23, 31]; % 四级寄存器系数

    % 确保初始状态不是全零向量，且长度为4
    if all(initial_state == 0)
        initial_state = ones(1, 4); % 随机初始化一个非全零状态
    end
    
    % 初始化移位寄存器
    register_states = initial_state;
    seq_length = 60;

    % 生成M序列
    m_seq = zeros(seq_length, 1);
    for i = 1:seq_length
        next_register = bitshift(register_states(1), -feedback_poly(1)) xor register_states(2);
        next_register = bitshift(next_register, -feedback_poly(2)) xor register_states(3);
        next_register = bitshift(next_register, -feedback_poly(3)) xor register_states(4);
        
        % 更新寄存器
        register_states = [next_register, register_states(1:end-1)];
        
        % 保留序列的最后一项
        m_seq(i) = register_states(1);
    end
end

要运行此函数，首先保存代码为`generate_M_sequence.m`，然后在MATLAB命令窗口中调用它：

initial_state = randi([0 1], 1, 4); % 示例初始状态，可以任意设置
m_sequence = generate_M_sequence(initial_state);
disp(m_sequence(1:60)); % 显示前60个序列元素

这个代码会生成一个长度为60的M序列，初始状态由用户指定，但不能全部为0。如果需要其他帮助或详细解释，请告诉我。