Disk (mathematics): Difference between revisions
en>Tobias Bergemann m Reverted 1 edit by 32.173.153.198 (talk) identified as vandalism to last revision by Mcapdevila. (TW) |
Disambiguated: centre → centre (geometry) |
||
Line 1: | Line 1: | ||
[[Image:Finite element solution.svg|right|thumb|A solution to a discretized partial differential equation, obtained with the [[finite element method]].]] | |||
In [[mathematics]], '''discretization''' concerns the process of transferring [[continuous function|continuous]] models and equations into [[wiktionary:discrete|discrete]] counterparts. This process is usually carried out as a first step toward making them suitable for numerical evaluation and implementation on digital computers. Processing on a digital computer requires another process called [[Quantization (signal processing)|quantization]]. | |||
* [[Euler–Maruyama method]] | |||
* [[Zero-order hold]] | |||
Discretization is also related to [[discrete mathematics]], and is an important component of [[granular computing]]. In this context, ''discretization'' may also refer to modification of variable of category ''granularity'', as when multiple discrete variables are aggregated or multiple discrete categories fused. | |||
== Discretization of linear state space models == | |||
Discretization is also concerned with the transformation of continuous [[differential equation]]s into discrete [[difference equations]], suitable for [[Numerical analysis|numerical computing]]. | |||
The following continuous-time [[State space (controls)|state space model]] | |||
:<math>\dot{\mathbf{x}}(t) = \mathbf A \mathbf{x}(t) + \mathbf B \mathbf{u}(t) + \mathbf{w}(t)</math> | |||
:<math>\mathbf{y}(t) = \mathbf C \mathbf{x}(t) + \mathbf D \mathbf{u}(t) + \mathbf{v}(t)</math> | |||
where ''v'' and ''w'' are continuous zero-mean [[white noise]] sources with [[covariance]]s | |||
:<math>\mathbf{w}(t) \sim N(0,\mathbf Q)</math> | |||
:<math>\mathbf{v}(t) \sim N(0,\mathbf R)</math> | |||
can be discretized, assuming [[zero-order hold]] for the input ''u'' and continuous integration for the noise ''v'', to | |||
:<math>\mathbf{x}[k+1] = \mathbf A_d \mathbf{x}[k] + \mathbf B_d \mathbf{u}[k] + \mathbf{w}[k]</math> | |||
:<math>\mathbf{y}[k] = \mathbf C_d \mathbf{x}[k] + \mathbf D_d \mathbf{u}[k] + \mathbf{v}[k]</math> | |||
with covariances | |||
:<math>\mathbf{w}[k] \sim N(0,\mathbf Q_d)</math> | |||
:<math>\mathbf{v}[k] \sim N(0,\mathbf R_d)</math> | |||
where | |||
:<math>\mathbf A_d = e^{\mathbf A T} = \mathcal{L}^{-1}\{(s\mathbf I - \mathbf A)^{-1}\}_{t=T} </math> | |||
:<math>\mathbf B_d = \left( \int_{\tau=0}^{T}e^{\mathbf A \tau}d\tau \right) \mathbf B = \mathbf A^{-1}(\mathbf A_d - I)\mathbf B </math>, if <math>\mathbf A</math> is [[Invertible matrix|nonsingular]] | |||
:<math>\mathbf C_d = \mathbf C </math> | |||
:<math>\mathbf D_d = \mathbf D </math> | |||
:<math>\mathbf Q_d = \int_{\tau=0}^{T} e^{\mathbf A \tau} \mathbf Q e^{\mathbf A^T \tau} d\tau </math> | |||
:<math>\mathbf R_d = \mathbf R </math> | |||
and <math>T</math> is the sample time, although <math>\mathbf A^T</math> is the transposed matrix of <math>\mathbf A</math>. | |||
A clever trick to compute ''Ad'' and ''Bd'' in one step is by utilizing the following property, p. 215:<ref>Raymond DeCarlo: ''Linear Systems: A State Variable Approach with Numerical Implementation'', Prentice Hall, NJ, 1989</ref> | |||
:<math>e^{\begin{bmatrix} \mathbf{A} & \mathbf{B} \\ | |||
\mathbf{0} & \mathbf{0} \end{bmatrix} T} = \begin{bmatrix} \mathbf{M_{11}} & \mathbf{M_{12}} \\ | |||
\mathbf{0} & \mathbf{I} \end{bmatrix}</math> | |||
and then having | |||
:<math>\mathbf A_d = \mathbf M_{11}</math> | |||
:<math>\mathbf B_d = \mathbf M_{12}</math> | |||
=== Discretization of process noise === | |||
Numerical evaluation of <math>\mathbf{Q}_d</math> is a bit trickier due to the matrix exponential integral. It can, however, be computed by first constructing a matrix, and computing the exponential of it (Van Loan, 1978): | |||
:<math> \mathbf{F} = | |||
\begin{bmatrix} -\mathbf{A} & \mathbf{Q} \\ | |||
\mathbf{0} & \mathbf{A}^T \end{bmatrix} T</math> | |||
:<math> \mathbf{G} = e^\mathbf{F} = | |||
\begin{bmatrix} \dots & \mathbf{A}_d^{-1}\mathbf{Q}_d \\ | |||
\mathbf{0} & \mathbf{A}_d^T \end{bmatrix}.</math> | |||
The discretized process noise is then evaluated by multiplying the transpose of the lower-right partition of '''G''' with the upper-right partition of '''G''': | |||
:<math>\mathbf{Q}_d = (\mathbf{A}_d^T)^T (\mathbf{A}_d^{-1}\mathbf{Q}_d). </math> | |||
=== Derivation === | |||
Starting with the continuous model | |||
:<math>\mathbf{\dot{x}}(t) = \mathbf A\mathbf x(t) + \mathbf B \mathbf u(t)</math> | |||
we know that the [[matrix exponential]] is | |||
:<math>\frac{d}{dt}e^{\mathbf At} = \mathbf A e^{\mathbf At} = e^{\mathbf At} \mathbf A</math> | |||
and by premultiplying the model we get | |||
:<math>e^{-\mathbf At} \mathbf{\dot{x}}(t) = e^{-\mathbf At} \mathbf A\mathbf x(t) + e^{-\mathbf At} \mathbf B\mathbf u(t)</math> | |||
which we recognize as | |||
:<math>\frac{d}{dt}(e^{-\mathbf At}\mathbf x(t)) = e^{-\mathbf At} \mathbf B\mathbf u(t)</math> | |||
and by integrating.. | |||
:<math>e^{-\mathbf At}\mathbf x(t) - e^0\mathbf x(0) = \int_0^t e^{-\mathbf A\tau}\mathbf B\mathbf u(\tau) d\tau</math> | |||
:<math>\mathbf x(t) = e^{\mathbf At}\mathbf x(0) + \int_0^t e^{\mathbf A(t-\tau)} \mathbf B\mathbf u(\tau) d \tau</math> | |||
which is an analytical solution to the continuous model. | |||
Now we want to discretise the above expression. We assume that u is [[mathematical constant|constant]] during each timestep. | |||
:<math>\mathbf x[k] \ \stackrel{\mathrm{def}}{=}\ \mathbf x(kT)</math> | |||
:<math>\mathbf x[k] = e^{\mathbf AkT}\mathbf x(0) + \int_0^{kT} e^{\mathbf A(kT-\tau)} \mathbf B\mathbf u(\tau) d \tau</math> | |||
:<math>\mathbf x[k+1] = e^{\mathbf A(k+1)T}\mathbf x(0) + \int_0^{(k+1)T} e^{\mathbf A((k+1)T-\tau)} \mathbf B\mathbf u(\tau) d \tau</math> | |||
:<math>\mathbf x[k+1] = e^{\mathbf AT} \left[ e^{\mathbf AkT}\mathbf x(0) + \int_0^{kT} e^{\mathbf A(kT-\tau)} \mathbf B\mathbf u(\tau) d \tau \right]+ \int_{kT}^{(k+1)T} e^{\mathbf A(kT+T-\tau)} \mathbf B\mathbf u(\tau) d \tau</math> | |||
We recognize the bracketed expression as <math>\mathbf x[k]</math>, and the second term can be simplified by substituting <math>v = kT + T - \tau</math>. We also assume that <math>\mathbf u</math> is constant during the [[integral]], which in turn yields | |||
:<math>\mathbf x[k+1] = e^{\mathbf AT}\mathbf x[k] + \left( \int_0^T e^{\mathbf Av} dv \right) \mathbf B\mathbf u[k]=e^{\mathbf AT}\mathbf x[k] + \mathbf A^{-1}\left(e^{\mathbf AT}-\mathbf I \right) \mathbf B\mathbf u[k]</math> | |||
which is an exact solution to the discretization problem. | |||
=== Approximations === | |||
Exact discretization may sometimes be intractable due to the heavy matrix exponential and integral operations involved. It is much easier to calculate an approximate discrete model, based on that for small timesteps <math>e^{\mathbf AT} \approx \mathbf I + \mathbf A T</math>. The approximate solution then becomes: | |||
:<math>\mathbf x[k+1] \approx (\mathbf I + \mathbf AT) \mathbf x[k] + T\mathbf B \mathbf u[k] </math> | |||
Other possible approximations are <math>e^{\mathbf AT} \approx \left( \mathbf I - \mathbf A T \right)^{-1}</math> and <math>e^{\mathbf AT} \approx \left( \mathbf I +\frac{1}{2} \mathbf A T \right) \left( \mathbf I - \frac{1}{2} \mathbf A T \right)^{-1}</math>. Each of them have different stability properties. The last one is known as the bilinear transform, or Tustin transform, and preserves the (in)stability of the continuous-time system. | |||
== Discretization of continuous features == | |||
{{Main|Discretization of continuous features}} | |||
In [[statistics]] and machine learning, '''discretization''' refers to the process of converting continuous features or variables to discretized or nominal features. This can be useful when creating probability mass functions. | |||
==See also== | |||
*[[Discrete space]] | |||
*[[Time-scale calculus]] | |||
*[[Discrete event simulation]] | |||
*[[Stochastic simulation]] | |||
*[[Finite volume method for unsteady flow]] | |||
*[[Properties of discretization schemes]] | |||
== References == | |||
<references/> | |||
* {{cite book|author=Robert Grover Brown & Patrick Y. C. Hwang|title=Introduction to random signals and applied Kalman filtering|edition=3rd|isbn=978-0471128397}} | |||
* {{cite book|publisher=Saunders College Publishing|location=Philadelphia, PA, USA|year=1984|author=Chi-Tsong Chen|title=Linear System Theory and Design|isbn=0030716918}} | |||
* {{cite journal|author=C. Van Loan|title=Computing integrals involving the matrix exponential|doi=10.1109/TAC.1978.1101743|journal=IEEE Transactions on Automatic Control|volume=23|issue=3|pages=395–404|date=Jun 1978}} | |||
==External links== | |||
{{sisterlinks}} | |||
[[Category:Numerical analysis]] | |||
[[Category:Applied mathematics]] | |||
[[Category:Functional analysis]] | |||
[[Category:Iterative methods]] | |||
[[Category:Control theory]] | |||
[[de:Diskretisierung]] | |||
[[hr:Diskretizacija]] | |||
[[it:Discretizzazione]] | |||
[[pl:Dyskretyzacja (matematyka)]] | |||
[[zh:离散化]] |
Revision as of 12:28, 12 September 2013
![](https://upload.wikimedia.org/wikipedia/commons/thumb/a/a9/Finite_element_solution.svg/300px-Finite_element_solution.svg.png)
In mathematics, discretization concerns the process of transferring continuous models and equations into discrete counterparts. This process is usually carried out as a first step toward making them suitable for numerical evaluation and implementation on digital computers. Processing on a digital computer requires another process called quantization.
Discretization is also related to discrete mathematics, and is an important component of granular computing. In this context, discretization may also refer to modification of variable of category granularity, as when multiple discrete variables are aggregated or multiple discrete categories fused.
Discretization of linear state space models
Discretization is also concerned with the transformation of continuous differential equations into discrete difference equations, suitable for numerical computing.
The following continuous-time state space model
where v and w are continuous zero-mean white noise sources with covariances
can be discretized, assuming zero-order hold for the input u and continuous integration for the noise v, to
with covariances
where
- , if is nonsingular
and is the sample time, although is the transposed matrix of .
A clever trick to compute Ad and Bd in one step is by utilizing the following property, p. 215:[1]
and then having
Discretization of process noise
Numerical evaluation of is a bit trickier due to the matrix exponential integral. It can, however, be computed by first constructing a matrix, and computing the exponential of it (Van Loan, 1978):
The discretized process noise is then evaluated by multiplying the transpose of the lower-right partition of G with the upper-right partition of G:
Derivation
Starting with the continuous model
we know that the matrix exponential is
and by premultiplying the model we get
which we recognize as
and by integrating..
which is an analytical solution to the continuous model.
Now we want to discretise the above expression. We assume that u is constant during each timestep.
We recognize the bracketed expression as , and the second term can be simplified by substituting . We also assume that is constant during the integral, which in turn yields
which is an exact solution to the discretization problem.
Approximations
Exact discretization may sometimes be intractable due to the heavy matrix exponential and integral operations involved. It is much easier to calculate an approximate discrete model, based on that for small timesteps . The approximate solution then becomes:
Other possible approximations are and . Each of them have different stability properties. The last one is known as the bilinear transform, or Tustin transform, and preserves the (in)stability of the continuous-time system.
Discretization of continuous features
Mining Engineer (Excluding Oil ) Truman from Alma, loves to spend time knotting, largest property developers in singapore developers in singapore and stamp collecting. Recently had a family visit to Urnes Stave Church. In statistics and machine learning, discretization refers to the process of converting continuous features or variables to discretized or nominal features. This can be useful when creating probability mass functions.
See also
- Discrete space
- Time-scale calculus
- Discrete event simulation
- Stochastic simulation
- Finite volume method for unsteady flow
- Properties of discretization schemes
References
- ↑ Raymond DeCarlo: Linear Systems: A State Variable Approach with Numerical Implementation, Prentice Hall, NJ, 1989
- 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534 - 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534 - One of the biggest reasons investing in a Singapore new launch is an effective things is as a result of it is doable to be lent massive quantities of money at very low interest rates that you should utilize to purchase it. Then, if property values continue to go up, then you'll get a really high return on funding (ROI). Simply make sure you purchase one of the higher properties, reminiscent of the ones at Fernvale the Riverbank or any Singapore landed property Get Earnings by means of Renting
In its statement, the singapore property listing - website link, government claimed that the majority citizens buying their first residence won't be hurt by the new measures. Some concessions can even be prolonged to chose teams of consumers, similar to married couples with a minimum of one Singaporean partner who are purchasing their second property so long as they intend to promote their first residential property. Lower the LTV limit on housing loans granted by monetary establishments regulated by MAS from 70% to 60% for property purchasers who are individuals with a number of outstanding housing loans on the time of the brand new housing purchase. Singapore Property Measures - 30 August 2010 The most popular seek for the number of bedrooms in Singapore is 4, followed by 2 and three. Lush Acres EC @ Sengkang
Discover out more about real estate funding in the area, together with info on international funding incentives and property possession. Many Singaporeans have been investing in property across the causeway in recent years, attracted by comparatively low prices. However, those who need to exit their investments quickly are likely to face significant challenges when trying to sell their property – and could finally be stuck with a property they can't sell. Career improvement programmes, in-house valuation, auctions and administrative help, venture advertising and marketing, skilled talks and traisning are continuously planned for the sales associates to help them obtain better outcomes for his or her shoppers while at Knight Frank Singapore. No change Present Rules
Extending the tax exemption would help. The exemption, which may be as a lot as $2 million per family, covers individuals who negotiate a principal reduction on their existing mortgage, sell their house short (i.e., for lower than the excellent loans), or take part in a foreclosure course of. An extension of theexemption would seem like a common-sense means to assist stabilize the housing market, but the political turmoil around the fiscal-cliff negotiations means widespread sense could not win out. Home Minority Chief Nancy Pelosi (D-Calif.) believes that the mortgage relief provision will be on the table during the grand-cut price talks, in response to communications director Nadeam Elshami. Buying or promoting of blue mild bulbs is unlawful.
A vendor's stamp duty has been launched on industrial property for the primary time, at rates ranging from 5 per cent to 15 per cent. The Authorities might be trying to reassure the market that they aren't in opposition to foreigners and PRs investing in Singapore's property market. They imposed these measures because of extenuating components available in the market." The sale of new dual-key EC models will even be restricted to multi-generational households only. The models have two separate entrances, permitting grandparents, for example, to dwell separately. The vendor's stamp obligation takes effect right this moment and applies to industrial property and plots which might be offered inside three years of the date of buy. JLL named Best Performing Property Brand for second year running
The data offered is for normal info purposes only and isn't supposed to be personalised investment or monetary advice. Motley Fool Singapore contributor Stanley Lim would not personal shares in any corporations talked about. Singapore private home costs increased by 1.eight% within the fourth quarter of 2012, up from 0.6% within the earlier quarter. Resale prices of government-built HDB residences which are usually bought by Singaporeans, elevated by 2.5%, quarter on quarter, the quickest acquire in five quarters. And industrial property, prices are actually double the levels of three years ago. No withholding tax in the event you sell your property. All your local information regarding vital HDB policies, condominium launches, land growth, commercial property and more
There are various methods to go about discovering the precise property. Some local newspapers (together with the Straits Instances ) have categorised property sections and many local property brokers have websites. Now there are some specifics to consider when buying a 'new launch' rental. Intended use of the unit Every sale begins with 10 p.c low cost for finish of season sale; changes to 20 % discount storewide; follows by additional reduction of fiftyand ends with last discount of 70 % or extra. Typically there is even a warehouse sale or transferring out sale with huge mark-down of costs for stock clearance. Deborah Regulation from Expat Realtor shares her property market update, plus prime rental residences and houses at the moment available to lease Esparina EC @ Sengkang
External links
de:Diskretisierung hr:Diskretizacija it:Discretizzazione pl:Dyskretyzacja (matematyka) zh:离散化