<p>The Kolmogorov <i>N</i>-width describes the best possible error one can achieve by elements of an <i>N</i>-dimensional linear space. Its decay has extensively been studied in approximation theory and for the solution of partial differential equations (PDEs). Particular interest has occurred within model order reduction (MOR) of parameterized PDEs, e.g., by the reduced basis method (RBM). While it is known that the <i>N</i>-width decays exponentially fast (and thus admits efficient MOR) for certain problems, there are examples of the linear transport and the wave equation, where the decay rate deteriorates to <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10444_2025_10224_Article_IEq1.gif" Format="GIF" Height="18" Rendition="HTML" Resolution="72" Type="Linedraw" Width="44" /> </InlineMediaObject> <EquationSource Format="TEX">\(N^{-1/2}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>N</mi> <mrow> <mo>-</mo> <mn>1</mn> <mo stretchy="false">/</mo> <mn>2</mn> </mrow> </msup> </math></EquationSource> </InlineEquation>. On the other hand, it is widely accepted that a smooth parameter dependence admits a fast decay of the <i>N</i>-width. However, a detailed analysis of the influence of properties of the data (such as regularity or slope) on the rate of the <i>N</i>-width seems to be lacking. In this paper, we state that the optimal linear space is a direct sum of shift-isometric eigenspaces corresponding to the largest eigenvalues, yielding an exact representation of the <i>N</i>-width as their sum. For the linear transport problem, which is modeled by half-wave symmetric initial and boundary conditions <i>g</i>, we obtain such an optimal decomposition by sorted trigonometric functions with eigenvalues that match the Fourier coefficients of <i>g</i>. Further, for normalized <i>g</i> in the Sobolev space <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10444_2025_10224_Article_IEq2.gif" Format="GIF" Height="14" Rendition="HTML" Resolution="72" Type="Linedraw" Width="24" /> </InlineMediaObject> <EquationSource Format="TEX">\(H^r\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>H</mi> <mi>r</mi> </msup> </math></EquationSource> </InlineEquation> of broken order <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10444_2025_10224_Article_IEq3.gif" Format="GIF" Height="13" Rendition="HTML" Resolution="72" Type="Linedraw" Width="40" /> </InlineMediaObject> <EquationSource Format="TEX">\(r&gt;0\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>r</mi> <mo>&gt;</mo> <mn>0</mn> </mrow> </math></EquationSource> </InlineEquation>, the sorted eigenfunctions give the sharp upper bound of the <i>N</i>-width, which is a reciprocal of a certain power sum. Yet, for ease, we also provide the decay <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10444_2025_10224_Article_IEq4.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="55" /> </InlineMediaObject> <EquationSource Format="TEX">\((\pi N)^{-r}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mrow> <mo stretchy="false">(</mo> <mi>π</mi> <mi>N</mi> <mo stretchy="false">)</mo> </mrow> <mrow> <mo>-</mo> <mi>r</mi> </mrow> </msup> </math></EquationSource> </InlineEquation>, obtained by the non-optimal space of ordering the trigonometric functions by frequency rather than by eigenvalue. Our theoretical investigations are complemented by numerical experiments which confirm the sharpness of our bounds and give additional quantitative insight.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

The Kolmogorov N-width for linear transport: exact representation and the influence of the data

  • Florian Arbes,
  • Constantin Greif,
  • Karsten Urban

摘要

The Kolmogorov N-width describes the best possible error one can achieve by elements of an N-dimensional linear space. Its decay has extensively been studied in approximation theory and for the solution of partial differential equations (PDEs). Particular interest has occurred within model order reduction (MOR) of parameterized PDEs, e.g., by the reduced basis method (RBM). While it is known that the N-width decays exponentially fast (and thus admits efficient MOR) for certain problems, there are examples of the linear transport and the wave equation, where the decay rate deteriorates to \(N^{-1/2}\) N - 1 / 2 . On the other hand, it is widely accepted that a smooth parameter dependence admits a fast decay of the N-width. However, a detailed analysis of the influence of properties of the data (such as regularity or slope) on the rate of the N-width seems to be lacking. In this paper, we state that the optimal linear space is a direct sum of shift-isometric eigenspaces corresponding to the largest eigenvalues, yielding an exact representation of the N-width as their sum. For the linear transport problem, which is modeled by half-wave symmetric initial and boundary conditions g, we obtain such an optimal decomposition by sorted trigonometric functions with eigenvalues that match the Fourier coefficients of g. Further, for normalized g in the Sobolev space \(H^r\) H r of broken order \(r>0\) r > 0 , the sorted eigenfunctions give the sharp upper bound of the N-width, which is a reciprocal of a certain power sum. Yet, for ease, we also provide the decay \((\pi N)^{-r}\) ( π N ) - r , obtained by the non-optimal space of ordering the trigonometric functions by frequency rather than by eigenvalue. Our theoretical investigations are complemented by numerical experiments which confirm the sharpness of our bounds and give additional quantitative insight.