<< Chapter < Page Chapter >> Page >

icalcf.m [x,y,t,u,px,py,p] = icalcf(X,Y,PX,PY) is a function version of icalc, which allows arbitrary naming of variables.

function [x,y,t,u,px,py,p] = icalcf(X,Y,PX,PY)% ICALCF [x,y,t,u,px,py,p] = icalcf(X,Y,PX,PY) Function version of icalc% Version of 5/3/95 % Allows arbitrary naming of variablesx = X; y = Y;px = PX; py = PY;if length(X) ~= length(PX) error(' X and PX of different lengths')elseif length(Y) ~= length(PY) error(' Y and PY of different lengths')end [a,b]= meshgrid(PX,fliplr(PY)); p = a.*b; % Matrix of joint independent probabilities[t,u] = meshgrid(X,fliplr(Y)); % t, u matrices for joint calculations
Got questions? Get instant answers now!

icalc3.m Calculation setup for an independent class of three simple random variables.

% ICALC3 file icalc3.m Setup for three independent rv % Version of 5/15/96% Sets up for calculations for three % independent simple random variables% Uses m-functions rep, elrep, kronf X = input('Enter row matrix of X-values ');Y = input('Enter row matrix of Y-values '); Z = input('Enter row matrix of Z-values ');PX = input('Enter X probabilities '); PY = input('Enter Y probabilities ');PZ = input('Enter Z probabilities '); n = length(X);m = length(Y); s = length(Z);[t,u] = meshgrid(X,Y);t = rep(t,1,s); u = rep(u,1,s);v = elrep(Z,m,n); % t,u,v matrices for joint calculations P = kronf(PZ,kronf(PX,PY'));disp('Use array operations on matrices X, Y, Z,') disp('PX, PY, PZ, t, u, v, and P')
Got questions? Get instant answers now!

icalc4.m Calculation setup for an independent class of four simple random variables.

% ICALC4 file icalc4.m Setup for four independent rv % Version of 5/15/96% Sets up for calculations for four % independent simple random variables% Uses m-functions rep, elrep, kronf X = input('Enter row matrix of X-values ');Y = input('Enter row matrix of Y-values '); Z = input('Enter row matrix of Z-values ');W = input('Enter row matrix of W-values '); PX = input('Enter X probabilities ');PY = input('Enter Y probabilities '); PZ = input('Enter Z probabilities ');PW = input('Enter W probabilities '); n = length(X);m = length(Y); s = length(Z);r = length(W); [t,u]= meshgrid(X,Y); t = rep(t,r,s);u = rep(u,r,s); [v,w]= meshgrid(Z,W); v = elrep(v,m,n); % t,u,v,w matrices for joint calculationsw = elrep(w,m,n); P = kronf(kronf(PZ,PW'),kronf(PX,PY'));disp('Use array operations on matrices X, Y, Z, W') disp('PX, PY, PZ, PW, t, u, v, w, and P')
Got questions? Get instant answers now!

Calculations for random variables

ddbn.m Uses the distribution of a simple random variable (or simple approximation) to plot a step graph for the distribution function F X .

% DDBN file ddbn.m Step graph of distribution function % Version of 10/25/95% Plots step graph of dbn function FX from % distribution of simple rv (or simple approximation)xc = input('Enter row matrix of VALUES '); pc = input('Enter row matrix of PROBABILITIES ');m = length(xc); FX = cumsum(pc);xt = [xc(1)-1-0.1*abs(xc(1)) xc xc(m)+1+0.1*abs(xc(m))];FX = [0 FX 1]; % Artificial extension of range and domainstairs(xt,FX) % Plot of stairstep graph hold onplot(xt,FX,'o') % Marks values at jump hold offgrid xlabel('t')ylabel('u = F(t)') title('Distribution Function')
Got questions? Get instant answers now!

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Applied probability. OpenStax CNX. Aug 31, 2009 Download for free at http://cnx.org/content/col10708/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Applied probability' conversation and receive update notifications?

Ask