ch7随机信号的产生与处理选编.ppt

取N＝3000点的情况执行该程序，所得结果如图7-12所示。3000点中，接受2301点，舍弃699点，效率达到76.70%，接近理论值（N-∞时）78.54%。■ % File: c7_rejex1.m R = 7; % default value of R M = 4/pi; % value of M N = input(Input number of points N ); % set N fx = zeros(1,N); % array of output samples u1 = rand(1,N); u2 = rand(1,N); % generate u1 and u2 v1 = R*u1; % generate v1 v2 = (M/R)*rand(1,N); % generate v2 (g(x)) kpts = 0; % initialize counter for k=1:N if v2(k)(M/(R*R))*sqrt(R*R-v1(k)*v1(k)); kpts=kpts+1; % increment counter fx(kpts)=v1(k); % save output sample end end fx = fx(1:kpts); [N_samp,x] = hist(f