判決式仿真代碼及仿真圖

《判決式仿真代碼及仿真圖》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在教育資源-天天文庫。

1、判決式仿真代碼及仿真圖mx=[-2:0.01:2];a=0;%隨機變量A取一確定值sigma=1;%發(fā)射隨機信號x的標準差asigma=2;%隨機變量A的標準差gama=0.1;%虛警概率約束為10%N=20;%接收隨機信號采樣點數(shù)%H0為真的接收信號采樣均值的概率密度函數(shù)pdh0=normpdf(mx,0,sigma/sqrt(N));plot(mx,pdh0);holdon;%H1為真的接收信號采樣均值的概率密度函數(shù)pdh1=normpdf(mx,a,sqrt((sigma^2+asigma^2)/N));plot(mx

2、,pdh1,'r');holdon;axis([-2,2,0,2])xlabel('mx');ylabel('mx概率密度');legend('p(mx

3、H0)','p(mx

4、H1)');%接收信號采樣均值的判決門限值beta=sqrt(sigma^2/N)*abs(norminv(1-gama/2,0,1));pdh0_beta=normpdf(beta,0,sigma/sqrt(N));pdh1_beta=normpdf(beta,a,sqrt((sigma^2+asigma^2)/N));num1=floor((2-b

5、eta)/0.01);line([mx(num1),mx(num1)],[0,pdh1(num1)],'Color','k','linestyle',':');fill([mx(1:num1),mx(num1)],[pdh0(1:num1),0],'g')num2=ceil((2+beta)/0.01);line([mx(num2),mx(num2)],[0,pdh1(num2)],'Color','k','linestyle',':');fill([mx(num2),mx(num2:end)],[0,pdh0(num2:e

6、nd)],'g')判決性能仿真代碼及仿真圖fori=1:51gama(i)=0.02*(i-1);err(i)=errper(gama(i));endplot(gama,err)xlabel('虛警概率gama')ylabel('虛警概率為gama,仿真1000次的漏警概率');errper：計算仿真1000次的漏警概率functionerr=errper(gama)%虛警概率約束為gama的條件下，在仿真1000次的漏警概率mx=[-2:0.01:2];a=0;%隨機變量A取一確定值sigma=1;%發(fā)射隨機信號x的標準差

7、asigma=2;%隨機變量A的標準差N=20;%接收隨機信號采樣點數(shù)beta=sqrt(sigma^2/N)*abs(norminv(1-gama/2,0,1));%接收信號采樣均值的判決門限值errnum=0;fori=1:1000ifround(rand())==1xx=normrnd(0,sqrt(sigma^2+asigma^2),1,N);ifabs(mean(xx))

8、概率