# Normpdf - but WITHOUT ANY WARRANTY without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU General

This preview shows page 1. Sign up to view the full content.

function p = normpdf(x,m,s); % NORMPDF returns normal probability density % % pdf = normpdf(x,m,s); % % Computes the PDF of a the normal distribution % with mean m and standard deviation s % default: m=0; s=1; % x,m,s must be matrices of same size, or any one can be a scalar. % % see also: NORMCDF, NORMINV % % Reference(s): % % \$Revision: 1.7 \$ % \$Id: normpdf.m,v 1.7 2005/05/25 02:35:36 pkienzle Exp \$ % Version 1.28 Date: 13 Mar 2003 % Copyright (C) 2000-2003 by Alois Schloegl <[email protected]> C % This program is free software; you can redistribute it and/or modify % it under the terms of the GNU General Public License as published by % the Free Software Foundation; either version 2 of the License, or % (at your option) any later version. % % This program is distributed in the hope that it will be useful, %
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: but WITHOUT ANY WARRANTY; without even the implied warranty of % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the % GNU General Public License for more details. % % You should have received a copy of the GNU General Public License % along with this program; if not, write to the Free Software % Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA % if nargin==1, m=0;s=1; elseif nargin==2, s=1; end; e % allocate output memory and check size of argument z = (x-m)./s; % if this line causes an error, input arguments do not fit. %p = ((2*pi)^(-1/2))*exp(-z.^2/2)./s; SQ2PI = 2.5066282746310005024157652848110; p = exp(-z.^2/2)./(s*SQ2PI); p p((x==m) & (s==0)) = inf; p p(isinf(z)~=0) = 0; p p(isnan(x) | isnan(m) | isnan(s) | (s<0)) = nan;...
View Full Document

## This note was uploaded on 07/08/2011 for the course EE 440 taught by Professor Jenq-nenghwang during the Spring '11 term at University of Washington.

Ask a homework question - tutors are online