Path: news.mathworks.com!newsfeed-00.mathworks.com!nlpi057.nbdc.sbc.com!prodigy.net!Xl.tags.giganews.com!border1.nntp.dca.giganews.com!nntp.giganews.com!local02.nntp.dca.giganews.com!nntp.web-ster.com!news.web-ster.com.POSTED!not-for-mail
NNTP-Posting-Date: Mon, 21 Jul 2008 20:11:18 -0500
Date: Mon, 21 Jul 2008 18:11:28 -0700
From: Tim Wescott <tim@seemywebsite.com>
User-Agent: Thunderbird 2.0.0.14 (Windows/20080421)
MIME-Version: 1.0
Newsgroups: comp.soft-sys.matlab,sci.stat.math,sci.engr.control,comp.dsp
Subject: Re: Kalman filtering with multiplicative noise
References: <1309743c-3d5c-4165-9cdb-1ab51892727d@z16g2000prn.googlegroups.com> 	<g61e76$f06$1@fred.mathworks.com> <1d3da6d4-ec1f-4e99-ae7e-c3478412f54c@l8g2000prm.googlegroups.com>
In-Reply-To: <1d3da6d4-ec1f-4e99-ae7e-c3478412f54c@l8g2000prm.googlegroups.com>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
Message-ID: <Ht-dna1mP6ArrhjVnZ2dnUVZ_oDinZ2d@web-ster.com>
Lines: 44
X-Usenet-Provider: http://www.giganews.com
X-Trace: sv3-NCuAZ3uLfrbisjsYkukkXZv3LDmIh5GxZcOk6f6jnd5yJ7u5pQ6as4oatQUEbppW/zz1ApEBAvcFP7O!lun81M+AnGML2jX4ipKmWGZp4sLMgoFbMOS82VXUd4BAdQRi2RoCCD28tuLIZ8/BIFwPXjeDH1vB!NbNoWkVEypOZy1c=
X-Abuse-and-DMCA-Info: Please be sure to forward a copy of ALL headers
X-Abuse-and-DMCA-Info: Otherwise we will be unable to process your complaint properly
X-Postfilter: 1.3.39
Bytes: 3125
Xref: news.mathworks.com comp.soft-sys.matlab:480817 sci.stat.math:83514 sci.engr.control:50820 comp.dsp:245400



dsp@myallit.com wrote:
> On Jul 21, 5:35 pm, "Bruno Luong" <b.lu...@fogale.fr> wrote:
>> d...@myallit.com wrote in message
>>
>> You might consider Extended Kalman filtering (EKF). Be aware
>> about the eventual non-stability of the scheme.
>>
> 
> What do you mean by the eventual non-stability? I did look at the EKF,
> there is some simple sample MATLAB code here:
> http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=18189
> 
> But in the first few lines of this script it says:
> 
> % for nonlinear dynamic system:
> %           x_k+1 = f(x_k) + w_k
> %           z_k   = h(x_k) + v_k
> % where w ~ N(0,Q) meaning w is gaussian noise with covariance Q
> %       v ~ N(0,R) meaning v is gaussian noise with covariance R
> 
> so the EKF looks appropriate for non-linear process models and
> measurement models that can be represented by any arbitrary functions
> f(x) and h(x), but the noise is still assumed to be additive.

He means that it's a nonlinear system, and therefore you can't assume 
that just because it's locally stable around some state that it's 
globally stable for any state.  In particular, there may be input 
vectors that will drive it into a limit cycle or off toward infinity, 
either permanently or temporarily.

Contrary to the beliefs of some, the MathWorks isn't on high, and MatLab 
isn't an extension of the Bible.  I would trust them for examples, but 
don't take their documentation for anything but a means of selling you 
copies of MatLab.

-- 

Tim Wescott
Wescott Design Services
http://www.wescottdesign.com

Do you need to implement control loops in software?
"Applied Control Theory for Embedded Systems" gives you just what it says.
See details at http://www.wescottdesign.com/actfes/actfes.html