[求助]Matlab仿真雷达信号检测

我是新手,刚刚用MATLAB编写了一个雷达检测的小程序,确遇到了问题,想请高手们指点一下: 比如一个高斯白噪声叠加一个恒值信号我们知道经检波后得到的是广义瑞利分布的视频信号,这样可以通过产生随机数来仿真捕获概率; 但是如果不知道检波后的目标模型统计特性,是不是能先产生瑞利噪声序列,再直接加一个常数作为恒值信号,来仿真视频输出信号呢? 下面是我的程序, 发现俩种方法结果差异很大function Rddetection;% 纯噪声虚警固定门限处理，视频输出为瑞利包络分布；tic;k=1.38e-23;T0=290;Fn=2.5;Fn=10^(Fn/10); %接收机噪声参数；Bn=2e6;Pn=k*T0*Fn*Bn;Gsigma=sqrt(Pn) % 高斯均方差值。Rsigma=sqrt(0.5*Pn); % 检波后瑞利均方差值；Pf=1e-2;T=sqrt(2*log(1/Pf))*Rsigma % 虚警率，理论门限值；epsilon=Pf/10% 计算精度；alfa=0.95; % 置信度；N=Pf*(1-Pf)/((1-alfa)*epsilon^2)N=ceil(N); %计算得出所需实验次数，向无穷大取整；N=1e5; R = raylrnd(Rsigma,1,N); %N=1e8时已超出随机矩阵最大存储量！A=(R>T);n=sum(a)pf=n/N%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 得出虚警概率后，计算发现概率,定义信噪比为幅度比。% 方法一，构造广义瑞利分布目标模型；SNR=3;a=10^(SNR/10)*Gsigma;N1=1e3;Ri=Rice(a,Gsigma,N1);b=(Ri>T);m=sum(b)pd=m/N1toc % 结果基本符合！% Pf=1e-6;N=1e7;SNR=7;N1=1e3时，得pf=1.5e-6,pd=0.944;% Pf=1e-5;N=1e7;SNR=7;N1=1e3时，得pf=1.5e-6,pd=0.974;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% .. Kuzy&NI^w  
&].1[&M]  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册