En las matemáticas, la discretización se refiere al proceso de transferir modelos y ecuaciones continuos en contrapartes discretas . Este proceso se realiza generalmente en primer lugar hacia la fabricación de ellas convenientes para la evaluación y la puesta en práctica numéricas en las calculadoras numéricas. Para ser procesado en una calculadora numérica otra cuantificación nombrada de proceso es esencial.
discretización de Euler
asimiento de la Cero-orden
La discretización también se relaciona con las matemáticas discretas, y es un componente importante computacional granular. En este contexto, la discretización del puede también referir a la modificación de la variable de la granulosidad del de la categoría, como cuando se agregan las variables discretas múltiples o las categorías discretas múltiples están fundidas.
Discretización de los modelos de espacio de estado linear
La discretización también se refiere a la transformación de las ecuaciones diferenciales continuo en las ecuaciones de diferencia discretas, conveniente para el computacional numérico.
El del mathbf D del mathbf C del {y} (t) \ del mathbf {x} (t) \ del mathbf {u} (t) donde están continuos el v y el w cero-significar las fuentes de ruido blanco con el se pueden individualizar, el asimiento asumido de la Cero-orden para el u de la entrada y la integración continua para el v del ruido, a = \ + \ mathbf B_d \ mathbf {u} del mathbf A_d del \ del mathbf del {x} \ del mathbf {x} + \ = \ + \ + \ mathbf {w} del mathbf D_d del mathbf C_d del \ del mathbf mathbf {v} {y} \ del mathbf {x} \ del mathbf {u} con covariación
\ mathbf {v} \ sim N (0, \ mathbf Q_d)
\ mathbf {} \ sim N (0, \ mathbf R_d) de w donde
\ mathbf A_d = e^ {\ mathbf A T} = \ mathcal ^ {L} {- 1} \ {(s \ mathbf I - \ mathbf A)^ {- 1} \} = \ dejado (\ e^ del ^ del int_ {\ tau=0} {T} {\ mathbf A \ tau} d \ tau \ derecho) de B_d del \ del mathbf _ {t=T} \ = \ mathbf A^ del mathbf B {- 1} (\ mathbf A_d -) \ mathbf B de I, si el \ el mathbf A es = \ mathbf C de C_d no singular \ del mathbf del
de \ mathbf D_d = \ mathbf D \ mathbf Q_d = \ e^ del ^ del int_ {\ tau=0} {T} {\ mathbf A \ tau} \ e^ del mathbf Q {\ mathbf A^T \ tau} = \ mathbf R de R_d del \ del mathbf del
de d \ del tau y T es el tiempo de la muestra.
Derivación
El comenzar con el
del del modelo continuo \ = \ mathbf A \ mathbf x del mathbf {\ punto {x}} (t) (t) + \ mathbf B \ mathbf u (t) nosotros saben que matriz exponencial es
\ frac {d} {despegue} e^ {\ mathbf en} = \ mathbf A e^ {\ mathbf en} = e^ {\ mathbf en} \ mathbf A y por premultiplicando modelo nosotros consiguen
e^ {- \ mathbf en} \ mathbf {\ punto {x}} (t) = e^ {- \ mathbf en} \ mathbf A \ mathbf x (t) + e^ {- \ mathbf en} \ mathbf B \ mathbf u (t) cuál nosotros reconocen como
\ frac {d} {despegue} (e^ {- \ mathbf en} \ mathbf x (t)) = e^ {- \ mathbf en} \ mathbf B \ mathbf u (t) e integrando.
e^ {- \ mathbf en} \ mathbf x (t) - e^0 \ mathbf x (0) = \ int_0^t e^ {- \ mathbf A \} \ mathbf B \ mathbf u (\ tau) d del tau \ del
de tau \ mathbf x (t) = e^ {\ mathbf en} \ 0) e^ + \ int_0^t del mathbf x ({\ mathbf A (t \ tau)} \ mathbf B \ mathbf u (\ tau) d \ tau cuál es una solución analítica al modelo continuo.
Ahora queremos individualizar la expresión antedicha. Asumimos que u es el constante durante cada timestep.
\ mathbf x \ \ stackrel {\ mathrm {def}} {=} \ \ mathbf x (kT) \ mathbf x = e^ {\ mathbf AkT} \ 0) e^ + \ int_0^ {kT} del mathbf x ({\ mathbf A (kT- \ tau)} \ mathbf B \ mathbf u (\ tau) d \ del
de tau \ mathbf x = e^ {\ mathbf A (k+1) T} \ 0) e^ + \ int_0^ del mathbf x ({(k+1) T} {\ mathbf A ((k+1) t \ tau)} \ mathbf B \ mathbf u (\ tau) d \ del
de tau \ mathbf x = e^ {\ mathbf EN} \ e^ dejado {\ mathbf AkT} \ 0) e^ + \ int_0^ {kT} del mathbf x ({\ mathbf A (kT- \ tau)} \ mathbf B \ mathbf u (\ tau) d \ tau \ right+ \ e^ del ^ del int_ {kT} {(k+1) T} {\ mathbf A (kT+T- \ tau)} \ mathbf B \ mathbf u (\ tau) d \ tau Reconocemos la expresión acorchetada como el \ mathbf x, y el segundo término puede ser simplificado substituyendo el v = el kT + - \ tau de T. Nosotros también asumen que \ mathbf u es constante durante integral, que alternadamente rinde
\ mathbf x = e^ {\ mathbf EN} \ mathbf x + \ a la izquierda dv del e^ (\ int_0^T {\ mathbf sistema de pesos americano} \) derecho \ mathbf B \ mathbf u cuál es una solución exacta al problema de la discretización.
Aproximaciones
La discretización exacta puede a veces ser insuperable debido a
la matriz pesada exponencial y a las operaciones integrales implicadas. Él es
mucho fácil calcular aproximado discreto modelo, basado en eso para pequeño timesteps
e^ {\ mathbf EN} \ aproximadamente \ + \ mathbf A T del mathbf I. La solución aproximada entonces se convierte:
\ mathbf x \ aproximadamente (\ + \ mathbf del mathbf I EN) \ mathbf x + (\ mathbf I T + \ frac {1} {2} \ mathbf A T^2) \ mathbf B \ mathbf u cuál puede más lejos ser aproximado si
\ el frac {1} {2} \ mathbf A T^2 es pequeño; rindiendo
\ mathbf x \ aproximadamente (\ + \ mathbf del mathbf I EN) \ mathbf x + T \ mathbf B \ mathbf u Otro posible aproximación son e^ {\ mathbf EN} \ aproximadamente \ se fue (\ mathbf I - \ mathbf A T \ derecho) ^ {- 1} y e^ {\ mathbf EN} \ aproximadamente \ se fue (\ mathbf I + \ frac {1} {2} \ mathbf A T \) derecho \ a la izquierda (\ mathbf I - \ frac {1} {2} \ mathbf A T \ derecho) ^ {- 1} . Cada uno de ellos tiene diversas características de la estabilidad. El pasado se conoce como el bilineario transforma, o Tustin transforma, y los cotos la estabilidad (la O.U) del sistema del continuo-tiempo.
ZenithicSaint-Nazaire-le-Désert