N = 25;
sigmaratio = 0.01;
PD = zeros(1,1000);
SNR = zeros(1,1000);
PF = 0.00001;
count = 0;
for sigmaratio=0.01:0.01:10
count=count+1;
TH = sqrt(sigmaratio/N) * chi2inv(PF);
SNR(count) = 10*log10(sigmaratio);
PD(count) = chi2cdf(PF)*(TH/(sigmaratio+1));
end
figure;
plot(SNR, PD, 'b-');
xlabel('Signal to noise ratio(DB)');
ylabel('Probability of detection, PD');
hold on;
PF = PF*10;
count = 0;
for sigmaratio=0.01:0.01:10
count=count+1;
TH = sqrt(sigmaratio/N) * chi2inv(PF);
PD(count) = chi2cdf(PF)*(TH/(sigmaratio+1));
end
plot(SNR, PD, 'b-');
hold on;
PF = PF*10;
count = 0;
for sigmaratio=0.01:0.01:10
count=count+1;
TH = sqrt(sigmaratio/N) * chi2inv(PF);
PD(count) = chi2cdf(PF)*(TH/(sigmaratio+1));
end
plot(SNR, PD, 'b-');
hold on;
PF = PF*10;
count = 0;
for sigmaratio=0.01:0.01:10
count=count+1;
TH = sqrt(sigmaratio/N) * chi2inv(PF);
PD(count) = chi2cdf(PF)*(TH/(sigmaratio+1));
end
plot(SNR, PD, 'b-');
hold on;
PF = PF*10;
count = 0;
for sigmaratio=0.01:0.01:10
count=count+1;
TH = sqrt(sigmaratio/N) * chi2inv(PF);
PD(count) = chi2cdf(PF)*(TH/(sigmaratio+1));
end
plot(SNR, PD, 'b-');
hold on;

1 个评论

"Whats wrong with this code?"
No comments, no explanation, no introduction, no help, no examples, no idea what it should do...

请先登录,再进行评论。

 采纳的回答

Trying to run it, I get--
...
Error using chi2inv (line 22)
Requires two input arguments.
>>
Don't know what the intent is and with, as you note Stephen, no explanations, I'm not going to try to divine but need a DOF argument as the second argument to chi2inv
Whether that should be N or perhaps count OP will have to determine...

2 个评论

Moved OP's "Answer" followup here--dpb
I need to obtain this graph. Can you guide me through the coding?
Well, what else do you need besides that which I've already pointed out? I have no way of knowing the response function; what "Energy Detector" in particular???

请先登录,再进行评论。

更多回答(0 个)

类别

帮助中心File Exchange 中查找有关 Just for fun 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by