一十种概率密度函数.DOC

1、一十种概率密度函数function zhifangtu(x,m)%画数据的直方图,x表示要画的随机数,m表示所要画的条数%a=min(x);b=max(x);l=length(x);h=(b-a)/m;%量化xx=x/h;x=ceil(x);w=zeros(1,m);for i=1:lfor j=1:mif (x(i)=j)%x(i)落在j的区间上，则w(j)加1 w(j)=w(j)+1;else continueendendendw=w/(h*l);z=a:h:(b-h);bar(z,w);title(直方图)function y=junyun(n)%0-1的均匀分布，n代表数据量，一般要大

2、于1024%y=ones(1,n);x=ones(1,n);m=100000;x0=mod(ceil(m*rand(1,1),m);x0=floor(x0/2);x0=2*x0+1;u=11;x(1)=x0;for i=1:n-1x(i+1)=u*x(i)+0;x(i+1)=mod(x(i+1),m);x(i)=x(i)/m;end%x(n)单位化x(n)=x(n)/m;y=x;function y=zhishu(m,n)%指数分布,m表示指数分布的参数,m不能为0.n表示数据量,n一般要大于1024%x=junyun(n);for i=1;nif (x(i)=0)x(i)=0.0001;el

3、se continue;endend u=log(x);y=-(1/m)*u;function y=ruili(m,n)%瑞利分布,m是瑞利分布的参数,n代表数据量,n一般要大于1024%x=junyun(n);for i=1:nif (x(i)=0)x(i)=0.0001;else continue;endend u=(-2)*log(x);y=m*sqrt(u);function y=weibuer(a,b,n)%韦布尔分布,a,b表示参数,b不能为0.n表示数据量,一般要大于1024%a=1时,是指数分布%a=2时,是瑞利分布%x=junyun(n);for i=1:nif (x(i)=

4、0)x(i)=0.0001;elsecontinue;endendu=-log(x);y=b*u.(1/a);function y=swerling(n)%swelingII分布%r=ones(1,n);u=junyun(n);v=junyun(n);for i=1:nif (u(i)=0)u(i)=0.0001;elsecontinueendendfor i=1:nif (u(i)=v(i)u(i)=u(i)+0.0001else continueendend t=-2*log(u);h=2*pi*v;x=sqrt(t).*cos(h);z=sqrt(t).*sin(h);y=(r/2).*

5、(x.2+z.2);function y=bernoulli(p,n)%产生数据量为n的贝努利分布,其中p属于(0-1)之间。%-%u=junyun(n);y=zeros(1,n);for i=1:nif(u(i)1)p=-log(exp(1)+a)/exp(1)-y2)/a);if(x1=1)while(k=log(z)y(k)=x;k=k+1;elsecontinue;endendendy=b*y;function y=beitafenbu(a1,a2,n)%产生贝他分布的随机数，其中a1、a2是贝他分布的参数，n代表数据量%x1=dajiama(a1,1,n);x2=dajiama(a2

6、,1,n);y=x1./(x1+x2);function y1,y2=gaussian(n)%产生数据量为n的两个相互独立高斯分布y1、y2%-%k=1;y1=zeros(1,n);y2=zeros(1,n);while(k=1)continue;elseif(s=0)k=k+1;elsey1(k)=v1*sqrt(-2*log(s)/s);y2(k)=v2*sqrt(-2*log(s)/s);k=k+1;endendfunction y=canshu(x);y=ones(1,2);n=length(x);y(1)=sum(x)/n;z=x-y(1);z=z.2;y(2)=sum(z)/(n-

7、1);function y=correlation(x)%计算x的自相关函数%n=length(x);for i=1:nx1(i)=x(n+1-i);endy=conv(x,x1);二三种相关杂波function y=gaussianpu(x)%由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的高斯随机向量%N=0:20;f=20;T=1/256;c=2*f*T*sqrt(pi)*exp(-4*f2*pi2*T2*N.2);n=length(x);y=zeros(1,n);for k=1:nfor i=20:-1:0if (k-i)=0)continue;elsey(k)=y(k)+c(2

8、1-i)*x(k-i);endendfor i=20:40if (k-i)=0)continue;elsey(k)=y(k)+c(i-19)*x(k-i);endendendy=0.5*y;function y=weibuerpu(a,b,n)%由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的韦布尔分布的随机向量% z1,z2=gaussian(n);z1=5*z1;z2=5*z2;y1=sqrt(ba/2)*z1;y2=sqrt(ba/2)*z2;x1=gaussianpu(y1);x2=gaussianpu(y2);x1=sqrt(ba/2)*x1;x2=sqrt(ba/2)*x2;

9、y=x1.2+x2.2;b=canshu(y);y=y-b(1);function y=duishuzhengtaipu(a,b,n)%由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的对数正态随机向量%a表示标准方差,b表示均值%z1=gaussian(n);x=gaussianpu(z1);y=a*x;y=exp(y);y=b*y;b=canshu(y);y=y-b(1); %去掉直流分量function y=swerling2pu(n)%由数据量为n的高斯白噪声产生向量为n，功率谱为高斯型的斯维凌II型随机向量%r=6;z1,z2=gaussian(n);x1=gaussianpu(

10、z1);x2=gaussianpu(z2);y=x1.2+x2.2;y=r*y;b=canshu(y);y=y-b(1); %去掉直流分量function y=kexipu(m,n)%由数据量为n的高斯白噪声产生向量为n，功率谱为柯西谱的高斯随机向量wc=2*pi*256;T0=1/(256*m);x=gaussian(n);y=zeros(1,n);y(1)=wc*T0*x(1);for i=2:ny(i)=wc*T0*x(i)+exp(-wc*T0)*y(i-1);endb=canshu(y);%y=y-b(1); %去掉直流分量 y=conv(y,y);y=fft(y);y=abs(y)

11、;i=1:2*n-1;plot(i,y)function plotpu(x)%绘出随机数的功率谱密度函数频域的图形。%w=fft(x);w=abs(w);v=2*pi/length(w);i=0:v:(2*pi-v);plot(i,w);三.雷达系统仿真function t,s,g,f0,fs,f1=huibo%产生目标回波信号x,系统噪声y，地物杂波z以及回波p%f0=3*107; %发射信号频率w=0; %发射信号初始相位c=3*108; %光速l=c/f0; %雷达信号波长(载波波长)R=40000; %目标范围Vd=200; %雷达与目标之间的径向速度fd=2*Vd/l; %多普勒频率

12、Tr=600/f0; %脉冲重复周期N=10; %雷达脉冲串长度f1=f0/5; %调频带宽是发射信号频率的1/5k=10*f1/Tr;fs=3*f0; %仿真采样频率Ts=1/fs;%Tt=2*R/c;Btar=4*pi*R/l; %M=floor(Tr*fs); %一个脉冲重复周期内的采样点数M=600mt1=floor(2*Tr*fs/5);%mt1=720mt2=floor(3*Tr*fs/10);%mt2=540mt3=floor(7*Tr*fs/10);%mt3=1260mt4=floor(3*Tr*fs/5);%mt4=1080mt5=mt1-mt2;%mt5=180Vgain=

13、6;s=zeros(M,N); %回波幅度起伏 for m=1:Mfor n=1:Nv(m,n)=(u(mt1-m)-u(mt2-m)*Vgain*cos(2*pi*f0*(m-mt2)*Ts+2*pi*k*Ts*(m-mt2)/2*(m-mt2)*Ts+2*pi*fd*n*Tr);g(m,n)=(u(mt3-m)-u(mt4-m)*Vgain*cos(2*pi*(f0*(m-mt4)*Ts+k*Ts*(m-mt4)/2*(m-mt4)*Ts);%u(t)是发射信号包络endendfor i=1:M/10t(i)=Vgain*cos(2*pi*(k*Ts*(M/10-i)/2*(M/10-i)

14、*Ts);endfor n=1:Ny(1:M,n)=gaussian(M); %系统噪声服从高斯分布endfor n=1:Nz(1:M,n)=swerling2pu(M); %地物杂波服从swerling2型分布ends=v+g+y+z;%s=v+g;m,n=size(s);x=zeros(1,m*n);q=zeros(1,m*n);x=s(:);q=g(:);p=x;i=0:length(p)-1;subplot(2,1,1);plot(i,x),title(目标回波信号); %目标回波信号xsubplot(2,1,2); l=canshu(x);b=x-l(1);plotpu(b),tit

15、le(目标回波频谱); function y=gaofang(s,f0,fs,f1)%高频放大器%M,N=size(s);Vgain=3; %高放增益w1=(f0-4*f1/3)*2*pi/fs; %w1=26pi/45w3=(f0-f1/3)*2*pi/fs; %w3=29pi/45w2=(f0+7*f1/3)*2*pi/fs; %w2=37pi/45w4=(f0+4*f1/3)*2*pi/fs; %w4=34pi/45th=min(w3-w1),(w2-w4);%w3-w1=pi/15,w2-w4=pi/15M1=ceil(6.6*pi/th)+1; %M1=99%n=0:M1-1;w5=(w1+w3)/2;w6=(w4+w2)/2;%w5=11pi/18,w6=71pi/90h=wide(w6,M1)-wide(w5,M1);w=(hamming(M1);