通信系统计算机模拟第十讲.ppt

通信系统的计算机模拟 第十讲 第九讲 回顾 需要随机信号的产生与处理的原因 平稳与遍历性过程 均匀随机数发生器 线性同余 混合同余算法 具有素模数的乘性算法 具有非素模数的乘性算法 随机数发生器的测试 散点图 Durbin-Watson测试 最低标准 Lewis, Goodman和Miller最低标准 Wichmann-Hill算法 MATLAB实现 7.3将均匀分布随机变量映射成任意pdf 1．已知目标随机变量的累积分布函数具有闭合形式?逆变换法。 2．已知目标随机变量的pdf具有闭合形式，但CDF不具有闭合形式。高斯随机变量，专用方法，可以用舍弃法。 3．已知pdf和CDF都不具有闭合形式。在寻求一个随机数发生器，使其pdf与实验数据的pdf相一致时，经常会出现这种情况。 7.3.1 逆变换法 将一个不相关的均匀分布的随机序列U变换为一个具有概率分布函数Fx(x)的不相关的（独立的样本）序列。 无记忆的非线性器 由于非线性器是无记忆的，在输入序列不相关时，能保证输出序列也是不相关的。根据维纳——辛钦定理，不相关随机数序列的功率谱密度是常数（白噪声）。逆变换方法简单地设定： 根据维纳——辛钦定理，不相关随机数序列的功率谱密度是常数（白噪声）。逆变换方法简单地设定： 逆变换方法能产生具有所需分布函数的随机变量，分布函数是自变量的一个非减函数，如图7-5所示。 根据分布函数的定义 例7－7 将一个均匀分布的随机变量变换成具有单边指数分布的随机变量 其中ux是单位阶跃函数 Solution 我们首先是要求出CDF，也就是 实现均匀分布－》指数分布变换的MATLAB代码 clear all % be safe n = input(Enter number of points ); b = 3; % set pdf parameter u = rand(1,n); % generate U y_exp = -log(u)/b; % transformation [N_samp,x] = hist(y_exp,20) % get histogram parameters subplot(2,1,1) bar(x,N_samp,1) % plot histogram ylabel(Number of Samples) xlabel(Independent Variable - x) subplot(2,1,2) Resultβ=3,N＝100 Result N＝2000 例7-8 作为第二个例子，考虑瑞利随机变量，其pdf表示为 Solution U与1－U等价 N=10000次的结果 clear all % be safe n = input(Enter number of points ); varR = 3; % set pdf parameter u = rand(1,n); % generate U y_exp = sqrt(-2*varR*log(u)); % transformation [N_samp,r] = hist(y_exp,20); % get histogram parameters subplot(2,1,1) bar(r,N_samp,1) % plot histogram ylabel(Number of Samples) xlabel(Independent Variable - x) subplot(2,1,2) term1 = r.*r/2/varR; % exponent ray = (r/varR).*exp(-term1); % Rayleigh pdf del_r = r(3)-r(2); % determine bin width p_hist = N_samp/n/del_r; % probability from histogram plot(r,ray,k,r,p_hist,ok) % compare results ylabel(Probability Density) xlabel(Independent Variable - x) legend(true pdf,samples from histogram,1) 7.3.2 直方图法 假设有一组通过实验测得的数据 尽管pdf可由数据的直方图近似，但pdf和CDF都未知 ？？？——》 如何开发出一种算法，来产生一组样本点使之与实验数据具有相似的pdf。 第一步是产生实验数据的直方图，就大体知道了pdf和CDF，运用逆变换方法。 这里变换方法的简单扩展