How to integrate the Rice distribution and build a graph?

12 次查看(过去 30 天)
Hello. I need to integrate the Rice distribution to calculate the probability of optical signal error with mathematical expectation v=1, variance s=0.033. There should be a double integration of the Rice distribution, where the inner integral is from 0.9 to infinity, and the outer integral is from 0 to infinity.
Please help me write the working code for calculating and plotting. I've been trying to write for 3 weeks, nothing works, every time there's some kind of mistake, probably due to the fact that I have little experience on matlab.

回答(1 个)

Torsten
Torsten 2024-3-25
编辑:Torsten 2024-3-25
As written, H is an integration variable - you cannot prescribe a value for it since it runs from 0 to Inf.
v = 1;
sigma = 0.033;
H = 0.9;
p = @(y)sqrt(y/v).*1/(sqrt(2*pi)*sigma).*exp(-(y-v).^2/(2*sigma^2));
I = integral2(@(H,x)p(x).*p(H),0,Inf,@(H)H,Inf)
I = 0.4999
  3 个评论
Torsten
Torsten 2024-3-25
v = 1;
sigma = 0.033;
p = @(y)sqrt(y/v).*1/(sqrt(2*pi)*sigma).*exp(-(y-v).^2/(2*sigma^2));
y = 0.9:0.001:1.1;
plot(y,p(y))

请先登录,再进行评论。

Community Treasure Hunt

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

Start Hunting!

Translated by