mprbs.m

来自「系统辨识的M文件」· M 代码 · 共 17 行

M=[1,0,0,0,0,0,1,0,0,];              % 设置初始状态，9 级移位寄存器   
Np=511*2;                            % 逆重复 M 序列的周期
u=zeros(1,Np);  a=5;             
for i=1:511;                     
    m=xor(M(9),M(5));                % 反馈逻辑运算
    for j=9:-1:2                
          M(j)=M(j-1);               % 移位
    end                     
    M(1)=m;  u(i)=M(1);    u(511+i)=M(1);  % 周期循环
end                 
for i=1:2:Np          
    u(i)=xor(u(i),1);                % 与 1 模 2 加运算，相当于异或 
end                   
for i=1:Np            
    u(i)=(1-2*u(i))*a;               % 由逻辑值到物理值
end