雷达系统建模与仿真报告

1、设计报告一十种随机数的产生一概述.概论论是在已知随机变量的情况下，研究随机变量的统计特性及其参量，而随机变量的仿真正好与此相反，是在已知随机变量的统计特性及其参数的情况下研究如何在计算机上产生服从给定统计特性和参数随机变量。下面对雷达中常用的模型进行建模：均匀分布高斯分布指数分布广义指数分布瑞利分布广义瑞利分布Swerling分布t分布对数一正态分布韦布尔分布二随机分布模型的产生思想及建立.产生随机数最常用的是在（0,1）区间内均匀分布的随机数，其他分布的随机数可利用均匀分布随机数来产生。均匀分布1（0，1）区间的均匀分布：用混合同余法产生（0,1）之间均匀分布的随机数，伪随机数通常是利用递推

2、公式产生的，所用的混和同余法的递推公式为：x_X.Pn1一n+C(Modm)其中，C是非负整数。通过适当选取参数C可以改善随机数的统计性质。一般取作小于M的任意奇数正整数，最好使其与模M互素。其他参数的选择无M的选取与计算机的字长有关。x(1)一般取为奇数。用Matlab来实现，编程语言用Matlab语言，可以用hist函数画出产生随机数的直方图(即统计理论概率分布的一个样本的概率密度函数)，直观地看出产生随机数的有效程度。其产生程序如下：c=3;lamade=4*200+1;x(1)=11;M=2八36;fori=2:1:10000;x(i)=mod(lamade*x(i-1)+c,M);e

3、nd;x=x./M;hist(x,10);mean(x)var(x)运行结果如下：均值=方差=2(a,b)区间的均匀分布：利用已产生的(0,1)均匀分布随机数的基础上采用变换法直接产生(a,b)均匀分布的随机数。其概率密度函数如下：1p(x)=彳b-aaxb、0 xb其产生程序如下：c=3;lamade=4*201+1;a=6;b=10;x(1)=11;M=2八36;fori=2:1:10000;x(i)=mod(lamade*x(i-1)+c,M);end;x=x./M;%i=2:1:10000;y(i)=(b-a)*x(i)+a;n=5:11;hist(y,n),axis(a-1b+10m

4、ax(hist(y,n)+20);mean(y)var(y)上面程序中取a=6,b=10.即(6，10)区间上的均匀分布。运行结果如下：均值=方差=高斯分布：高斯分布的概率密度函数如下;1_(x-u)2p(x)=e262兀G其产生方法是在均匀分布随机数的基础上通过函数变换法来产生。产生步骤是产生均匀分布的随机数。产生服从标准正态分布的随机数。由标准正态分布产生一般正态分布。1标准正态分布其部分程序如下：%i=1:1:10000;u(i)=sqrt(-2*log(x(i).*cos(2.*pi.*y(i);v(i)=sqrt(-2*log(x(i).*sin(2.*pi.*y(i);n1=-5:

5、5;n2=-5:5;subplot(1,2,1);hist(u,n1);subplot(1,2,2);hist(v,n2);mean(u)var(u)mean(v)var(v)运行结果如下：均值=方差=2一般正态分布其部分程序如下：a=2;b=2;i=1:1:10000;u(i)=sqrt(-2*log(x(i).*cos(2.*pi.*y(i);v=b*u+a;n=-10:10;hist(v,n);mean(v)var(v)运行结果如下：250250均值=方差=指数分布：服从正态分布的信号通过线性检波器后其包络强度(功率)服从指数分布。其概率密度函数为：p(x)二九e入x0其产生方法亦有：在

6、均匀分布随机数的基础上产生指数分布。在正态分布随机数的基础上产生该分布。产生程序分别如下：程序1(部分)lamade1=1i=1:1:10000;y(i)=-log(x(i)./lamade1;n=0:10;hist(y,n);mean(y)var(y)运行结果：均值=方差=程序2(部分)i=1:1:10000;s(i)=(u(i).*u(i)+v(i).*v(i);n=0:25;hist(s,n);mean(s)ar(s)运行结果：瑞利分布：在雷达系统中载带信号的包络服从瑞利分布。正态随机过程在其杂波载频(f0)上可以表示为:c(t)=x(t)cost-y(t)sintcc其中x(t)、y(

7、t)是服从N(PQ2)的相互独立的随机过程，检波器的包络幅度(电压):v(t)=yx(t)2+y(t)2服从瑞利分布R&)。瑞利分布的概率密度函数为：x/x2f(x)斗頁exp(-石7)0,其产生方法亦有：在均匀分布随机数的基础上产生瑞利分布。在正态分布随机数的基础上产生该分布。其产生程序如下：程序1(部分)：segma=2;i=1:1:10000;y(i)=segma*sqrt(-2*log(x(i);n=0:10;hist(y,n);mean(y)var(y)运行结果：方差均值程序2（部分）：i=1:1:10000;s(i)=sqrt(u(i).*u(i)+v(i).*v(i);n=0:1

8、0;hist(s,n);mean(s)var(s)运行结果：TOi1111lECOI14CO13002CO-1OTi._JL.112n24601012均值=方差=广义指数分布概率密度函数为：p(x)=e-(x+s)I(2xs)o式中：s-信噪比部分程序如下：s=8;i=1:1:10000;h(i)=u(i)+sqrt(2*s);z(i)二h(i).*h(i)+y(i)2；n=0:1:60;hist(z,n);运行结果：方差=均值=广义瑞利分布x2+A2p(x)XAX、x0e26I()G20G2Aa=G2-信噪比部分程序如下：a=1；i=2:1:10240;s(i)二sqrt(u(i)+a).八

9、2+v(i).八2);n=-1:15;hist(s,n);mean(s)var(s)运行结果如下：均值=方差=韦布尔分布韦布尔分布模型的性能介于瑞利分布模型与对数一正态分布模型之间.海浪杂波和地面杂波都可以用它来表示，并且在一个相当宽的条件范围内它能精确地表示实际的杂波分布。韦布尔分布的概率密度函数为：p(x)=a式中：a-形状参数;b-比例参数;xx0 xO-位置参数;该分布是在服从瑞利分布随机数的基础上用变换法产生的，其产生源程序部分）及直方图如下：a=3;b=2;m=5;i=2:1:10000;y(i)二m+b*(Tog(x(i)八(1/a);y二m+b*(Tog(x).八(1/a);hist(y,60);mean(y)var(y)均值=方差=对数正态分布对数一正态分布模型可以用来表示高分辨率雷达在观察角小于5时，观察到的海浪杂波，在低观察角时观察到的地面杂波也可用对数一正态分布模型，这类杂波通常是形状不规则的大反射体，例如远洋舰船，较大的空间飞行器，或者SAR雷达观察到的城市等等。其概率密度函数是：11n(x/u)p(x)=e2b22兀ox均值=eu+o2/2，方差=e2u+o2(eo2-1)其产生源程序及直方图如下：i=1:1:10000;u(i)=sqrt(-2*log(x(i).*cos(2.*pi.*y(i);