How to simulate the outage probability using Matlab?

2016 12 17
1 个回答
更新时间：2016 12 18
2 次查看（30 天）

0 个投票

Can any one help me ?? i want to simulate my CDF (or outage) equation and i have never done the simulation before, can anyone send me a sample code on Matlab to learn how to do it?

请先登录，再回答此问题。

回答（1 个）

the cyclist
the cyclist 2016-12-17

0 个投票

Does the answer to this question help you?

5 个评论
显示 3更早的评论 隐藏 3更早的评论

Deema42
Deema42 2016-12-17
编辑：Deema42 2016-12-17
Thank you for your reply. Unfortunately, no!. i use Matlab very well, but i'm still a beginner on the simulation, so i need a good written code to learn from.
the cyclist
the cyclist 2016-12-17
In that case I suggest you
  • Give significantly more detail about what you are trying to do
  • Show the code that you have tried to write, so we can help you
Deema42
Deema42 2016-12-18
编辑：Deema42 2016-12-18
N = 1000;
mu = -1;
sigma = 5;
X = mu + sigma*randn(N,1);
% Points for which CDF and PDF are to be evaluated
x = linspace(-10,10,1000); % Theoretical PDF and CDF
fx = (1/sqrt((2*pi*sigma*sigma)))*exp(-(((x - mu).^2)/(2*sigma*sigma))); Fx = 0.5*(1 + erf((x - mu)/(sqrt(2*sigma*sigma)))); % Estimate PDF and CDF
[f,F] = EstimateDistribution(X,x); % Plot Results
figure(2); plot(x,f,x,fx,x,F,x,Fx); xlabel('x'); ylabel('PDF & CDF'); str1 = strcat('Simulated PDF;','Area = ',num2str(trapz(x,f))); str2 = strcat('Theoretical PDF;','Area = ',num2str(trapz(x,fx))); legend(str1,str2,'Simulated CDF','Theoretical CDF','Location','northwest');% code
Deema42
Deema42 2016-12-18
function [f,F] = EstimateDistribution(X,x) % Checking Input Arguments if nargin<2||isempty(x), x = linspace(-10,10,1000);end if nargin<2||isempty(X) error('Missing Input Arguments: Please specify vector random variables'); end
% Impelementation Starts Here f = zeros(1,length(x)); % Preallocation of memory space F = f; % Preallocation of memory space h = 0.000000001; % Small value closer to zero for evaluating % numerical differentiation.
% Estimating CDF by its definition for m = 1:length(x) p = 0; % True Probability q = 0; % False Probability for n = 1:length(X) if X(n)<=x(m) % Definition of CDF p = p + 1; else q = q + 1; end end F(m) = p/(p + q); % Calulating Probability end
% Estimating PDF by differentiation of CDF for k = 1:length(x) fxph = interp1(x,F,x(k) + h,'spline'); % Interpolating value of F(x+h) fxmh = interp1(x,F,x(k) - h,'spline'); % Interpolating value of F(x-h) f(k) = (fxph - fxmh)/(2*h); % Two-point formula to compute end % Numerical differentiation f = smooth(f);
Deema42
Deema42 2016-12-18
编辑：Deema42 2016-12-18
sorry for the delay, these are two codes i found on Mathworks, they are fine, i want to do exactly like this but for a different PDF and CDF, mine is for cascaded Rician random variable, so the pdf is different than a normal rician. how can i modify this code?? if you like, i can send you my PDF and CDF equations.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

提问:

Deema42
Deema42
2016-12-17

编辑：

Deema42
Deema42
2016-12-18

Community Treasure Hunt

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

Start Hunting!

Translated by