1 function [u,x,s]=HMM_REKONSTRUKTION(U,xfront,X,dx,n_x)
3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 %%% stueckweise konstante Rekonstruktion %%%
5 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
7 %
function [u,x,s]=HMM_REKONSTRUKTION(U,xfront,X,dx,n_x)
9 % Input: U derzeitige makroskopische Loesung
10 % xfront Stelle der Unstetigkeit
11 % X Vektor, makroskopische Gitterpunkte
12 % dx, n_x mikroskopische Ortsschrittweite bzw. Anzahl
13 % mikroskopische Ortsschritte
14 % Output: u stueckweise konstante Rekonstruktion der L?sung
15 % x Vektor, mikroskopische Gitterpunkte
16 % s Gitterpunkt an dem sich die Unstetigkeit befindet
20 x=(xfront:dx:xfront+(n_x-1)*dx);
22 % Berechnung des Intervalls I=(A,B], A,B \in X, so dass xfront \in I
24 s=find(a==1, 1,
'last' );
28 % Konstruiere u=RU durch stueckweise konstante Rekonstruktion
29 u=zeros(1,n_x); % initialisieren