判決式仿真代碼及仿真圖

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

1、判決式仿真代碼及仿真圖mx=[-2:0.01:2];a=0;%隨機(jī)變量A取一確定值sigma=1;%發(fā)射隨機(jī)信號(hào)x的標(biāo)準(zhǔn)差asigma=2;%隨機(jī)變量A的標(biāo)準(zhǔn)差gama=0.1;%虛警概率約束為10%N=20;%接收隨機(jī)信號(hào)采樣點(diǎn)數(shù)%H0為真的接收信號(hào)采樣均值的概率密度函數(shù)pdh0=normpdf(mx,0,sigma/sqrt(N));plot(mx,pdh0);holdon;%H1為真的接收信號(hào)采樣均值的概率密度函數(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)');%接收信號(hào)采樣均值的判決門(mén)限值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：計(jì)算仿真1000次的漏警概率functionerr=errper(gama)%虛警概率約束為gama的條件下，在仿真1000次的漏警概率mx=[-2:0.01:2];a=0;%隨機(jī)變量A取一確定值sigma=1;%發(fā)射隨機(jī)信號(hào)x的標(biāo)準(zhǔn)差

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

8、概率