junyun.m

来自「雷达swerling 目标的仿真代码」· M 代码 · 共 18 行

function y=junyun(n) 
%0-1的均匀分布，n代表数据量，
%
y=ones(1,n); %设单位序列x
x=ones(1,n); %设单位序列y
m=100000;    %取模值
x0=mod(ceil(m*rand(1,1)),m); %对随机数进行取模
x0=floor(x0/2); %除2运算
x0=2*x0+1;    
u=21;%混合同余法 长周期相关，高度不均匀性
x(1)=x0; 
for i=1:n-1 
x(i+1)=u*x(i)+3; 
x(i+1)=mod(x(i+1),m); 
x(i)=x(i)/m; %x(n)单位化 
end 
y=x;