十种概率密度函数Word下载.docx

1、x（n）=x（n）/m;y=x;function y=zhishu（m,n）%旨数分布，m表示指数分布的参数，m不能为O.n表示数据量,n 一般要大于1024 % x=junyun（n）;for i=1;nif （x（i）=0）x（i）=0.0001;continue ;u=log（x）;y=-（1/m）*u;function y=ruili（m,n）%瑞利分布，m是瑞利分布的参数,n代表数据量,n 一般要大于1024 % x=junyun（n）;u=（-2）*log（x）; y=m*sqrt（u）;function y=weibuer（a,b,n）%韦布尔分布 ,a,b 表示参数 ,b 不能

2、为 0.n 表示数据量 , 一般要大于 1024%a=1时,是指数分布%a=2时,是瑞利分布 % x=junyun（n）;end u=-log（x）;y=b*u.A（1/a）;function y=swerling（n）%swelingII 分布% r=ones（1,n）;u=junyun（n）;v=junyun（n）;if （u（i）=0）u（i）=0.0001;if （u（i）=v（i）u（i）=u（i）+0.0001else continuet=-2*log（u）;h=2*pi*v;x=sqrt（t）.*cos（h）;z=sqrt（t）.*sin（h）;y=（r/2）.*（x.A2+z.

3、A2）;function y=bernoulli（p,n）沪生数据量为n的贝努利分布，其中p属于（0-1）之间% %y=zeros（1,n）;n if （u（i）=p） y（i）=1;y（i）=0;function y=duishuzhengtai（a,b,n）%产生对数正态分布， a,b 为随机分布的参数， n 为数据量% x=gaussian（n）;u=sqrt（b）*x+a;y=exp（u）;function y=kaifeng（m,n）沪生开丰分布，其中 m弋表开丰分布的自由度， n表示产生的点数量% y=zeros（1,n）;if （floor（m/2）=m/2）m/2x1,x2=g

4、aussian（n）;y（j）=x1（j）A2+x2（j）A2+y（j）;floor（m/2）x=gaussian（n）;y（j）=y（j）+x（j）A2;end function y=dajiama（a,b,n）沪生伽马随机分布的数据， a、b为随机分布的参数，数据量为 n%k=1;if （a1）while （kp=-log（exp（1）+a）/exp（1）-y2）/a）;if （x1=1） v=（2*a-1）A（-0.5）*log（x1/（1-x2）; x=a*exp（v）;z=x1A2*x2; w=a-log（4）+（a+sqrt（2*a-1）*v-x;if （w=log（z） y（k）

5、=x;end end y=b*y;function y=beitafenbu（a1,a2,n）沪生贝他分布的随机数，其中 al、a2是贝他分布的参数，n代表数据量% x1=dajiama（a1,1,n）;x2=dajiama（a2,1,n）;y=x1./（x1+x2）;function y1,y2=gaussian（n）沪生数据量为n的两个相互独立高斯分布 y1、y2 y1=zeros（1,n）; y2=zeros（1,n）;u1=junyun（1）;u2=junyun（1）;v1=2*u1-1;v2=2*u2-1;s=v1A2+v2A2;if （s=1） continue ;elseif （

6、s=0）k=k+1;y1（k）=v1*sqrt（-2*log（s）/s）;y2（k）=v2*sqrt（-2*log（s）/s）;end function y=canshu（x）;y=ones（1,2）;n=length（x）;y（1）=sum（x）/n;z=x-y（1）;z=z.A2;y（2）=sum（z）/（n-1）;function y=correlation（x）%计算x的自相关函数%x1（i）=x（n+1-i）;y=conv（x,x1）;二三种相关杂波function y=gaussianpu（x）%由数据量为n的高斯白噪声产生向量为 n,功率谱为高斯型的高斯随机向量%N=0:20;f

7、=20;T=1/256;c=2*f*T*sqrt（pi）*exp（-4*fA2*p22*TA2*N.A2）;for k=1:for i=20:-1:if （k-i）=0）y（k）=y（k）+c（21-i）*x（k-i）;40y（k）=y（k）+c（i-19）*x（k-i）;y=0.5*y;function y=weibuerpu（a,b,n）%由数据量为n的高斯白噪声产生向量为 n,功率谱为高斯型的韦布尔分布的随机向量%z1,z2=gaussian（n）;z1=5*z1;z2=5*z2;y1=sqrt（bAa/2）*z1;y2=sqrt（bAa/2）*z2;x1=gaussianpu（y1）;

8、x2=gaussianpu（y2）;x1=sqrt（bAa/2）*x1;x2=sqrt（bAa/2）*x2;y=x1.A2+x2.A2;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; %去掉直流分量function y=swerling2pu（n）%由数据量为n的高斯白噪声产生向量为 n,功率谱为高斯型的斯维凌 II型随机

9、向量r=6;x1=gaussianpu（z1）;x2=gaussianpu（z2）;y=r*y;function y=kexipu（m,n）%由数据量为n的高斯白噪声产生向量为 n,功率谱为柯西谱的高斯随机向量wc=2*pi*256;T0=1/（256*m）;y（1）=wc*T0*x（1）;for i=2:y（i）=wc*T0*x（i）+exp（-wc*T0）*y（i-1）;y=conv（y,y）y=fft（y）;y=abs（y）; i=1:2*n-1;plot（i,y）function plotpu（x）%绘出随机数的功率谱密度函数频域的图形。% w=fft（x）;w=abs（w）;v=2*

10、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*10A7;w=0;c=3*10A8;l=c/f0;R=40000;Vd=200;fd=2*Vd/l;Tr=600/f0;%发射信号频率%发射信号初始相位%光速%雷达信号波长 （ 载波波长 ）%目标范围%雷达与目标之间的径向速度%多普勒频率%脉冲重复周期N=10; %雷达脉冲串长度f1=f0/5; %调频带宽是发射信号频率的 1/5k=10*f1/Tr;fs=3*f0;Ts=1/f

11、s;%Tt=2*R/c;Btar=4*pi*R/l; M=floor（Tr*fs）;%仿真采样频率%一个脉冲重复周期内的采样点数 M=600mt1=floor（2*Tr*fs/5）;mt2=floor（3*Tr*fs/10）;mt3=floor（7*Tr*fs/10）;mt4=floor（3*Tr*fs/5）;%mt1=720%mt2=540%mt3=1260%mt4=1080mt5=mt1-mt2; %mt5=180Vgain=6;s=zeros（M,N）; %回波幅度起伏for m=1:Mfor n=1:Nv（m,n）=（u（mt1-m）-u（mt2-m）*Vgain*cos（2*pi*f

12、0*（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） 是发射信号包络M/10t（i）=Vgain*cos（2*pi*（k*Ts*（M/10-i）/2*（M/10-i）*Ts）;y（1:M,n）=gaussian（M）; %系统噪声服从高斯分布z（1:M,n）=swerling2pu（M） %地物杂波服从 swerling2 型分布s=v+g+y+z;%

13、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）,title（ 目标回波频谱 function y=gaofang（s,f0,fs,f1） %高%频%放%大%器% % M,N=size（s）;Vgain=3;%高放增益w1=（f0-4*f1/3）*2*pi/fs;w3=（f0-f1/3）*2*pi

14、/fs;w2=（f0+7*f1/3）*2*pi/fs;w4=（f0+4*f1/3）*2*pi/fs; th=min（w3-w1）,（w2-w4）;M1=ceil（6.6*pi/th）+1;%n=0:M1-1;%w1=26pi/45%w3=29pi/45%w2=37pi/45%w4=34pi/45%w3-w1=pi/15,w2-w4=pi/15%M1=99w5=（w1+w3）/2;w6=（w4+w2）/2; %w5=11pi/18,w6=71pi/90h=wide（w6,M1）-wide（w5,M1）;w=（hamming（M1）h=h.*w;L=length（h）;s=Vgain*s;z=ze

15、ros（M+L-1,N）;y=zeros（M,N）;K=ceil（L/2）;10z（1+（j-1）*M/10:j*M/10+L-1,i）=conv（h,s（1+（j-1）*M/10:j*M/10,i）y（1+（j-1）*M/10:j*M/10,i）=z（1+（j-1）*M/10+K:j*M/10+K,i）;w=y（:subplot（3,1,1）;length（w）-1;plot（i,w）,title（ 高放 subplot（3,1,2）;plotpu（w）,title（ 频谱 subplot（3,1,3）;plotpu（h）,title（ 幅频特性function y=hunpin（s,f0,

16、fs,f1）%进%行%混%频%，%输%出%为%中频信号%M,N=size（s）;M-1;z1=cos（2*pi*2*f0*i/3/fs）;y1=zeros（M,N）;y1（1:M,i）=（z1.*s（1:M,i）w1=（f0-5*f1）*2*pi/（3*fs）;w3=（f0-2*f1）*2*pi/（3*fs）; w2=（f0+8*f1）*2*pi/（3*fs）;%w1=pi/9%w3=1.6pi/9%w2=0.4piw4=（f0+5*f1）*2*pi/（3*fs）; M1=ceil（6.6*pi/th）+1;%w4=pi/3%w3-w1=0.2pi/3,w2-w4=0.2pi/3%w5=1.3

17、pi/9,w6=1.1pi/3wth=（hamming（M1）h=h.*wth; z=zeros（M+L-1,N）; y=zeros（M,N）;10 z（1+（j-1）*M/10: y（1+（j-1）*M/10:v=y1（:subplot（4,1,1）; plot（i,w）,title（ 混频 subplot（4,1,2）;plotpu（v）,title（ subplot（4,1,3）;plotpu（w）;subplot（4,1,4）; plotpu（h）,title（ 幅频特性 function y=zhongfang（s,f0,fs,f1）%中%频%放%大%器% M,N=size（s）;%

18、中放增益Vgain=5;w2=（f0+8*f1）*2*pi/（3*fs）;th=min（w3-w1）,（w2-w4）; %n=0: w5=（w1+w3）/2; h=wide（w6,M1）-wide（w5,M1）; w=（hamming（M1） s=Vgain*s; for i=1:中放 频谱幅频特性 end end w=y（: subplot（3,1,1）; plot（i,w）,title（ subplot（3,1,2）; plotpu（w）,title（ subplot（3,1,3）; plotpu（h）,title（ function I,Q=xiangganjianbo（s,fs,f0,f1） %相%位%相%干%检%波% % M,N=size（s）;z1=cos（2*pi*f0*i/（3*fs）; z2=sin（2*pi*f0*i/（3*fs）;y2=zeros（M,N）;y2（1:M,i）=（z2.*s（1:wt1=2*pi*（4*f1）/（3*fs）;%wt1=4pi/45wt2=2*pi*（7*f1）/（3*fs