Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (52)

Search Parameters:
Keywords = non-linear conjugate gradient method

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 1797 KiB  
Article
Central Difference Variational Filtering Based on Conjugate Gradient Method for Distributed Imaging Application
by Wen Ye, Fubo Zhang and Hongmei Chen
Remote Sens. 2024, 16(18), 3541; https://doi.org/10.3390/rs16183541 - 23 Sep 2024
Abstract
The airborne distributed position and orientation system (ADPOS), which integrates multi-inertia measurement units (IMUs), a data-processing computer, and a Global Navigation Satellite System (GNSS), serves as a key sensor in new higher-resolution airborne remote sensing applications, such as array SAR and multi-node imaging [...] Read more.
The airborne distributed position and orientation system (ADPOS), which integrates multi-inertia measurement units (IMUs), a data-processing computer, and a Global Navigation Satellite System (GNSS), serves as a key sensor in new higher-resolution airborne remote sensing applications, such as array SAR and multi-node imaging loads. ADPOS can provide reliable, high-precision and high-frequency spatio-temporal reference information to realize multinode motion compensation with the various nonlinear filter estimation methods such as Central Difference Kalman Filtering (CDKF), and modified CDKF. Although these known nonlinear models demonstrate good performance, their noise estimation performance with its linear minimum variance estimation criterion is limited for ADPOS. For this reason, in this paper, Central Difference Variational Filtering (CDVF) based on the variational optimization process is presented. This method adopts the conjugate gradient algorithm to enhance the estimation performance for mean correction in the filtering update stage. On one hand, the proposed method achieves adaptability by estimating noise covariance through the variational optimization method. On the other hand, robustness is implemented under the minimum variance estimation criterion based on the conjugate gradient algorithm to suppress measurement noise. We conducted a real ADPOS flight test, and the experimental results show that the accuracy of the slave motion parameters has significantly improved compared to the current CDKF. Moreover, the compensation performance shows a clear enhancement. Full article
(This article belongs to the Special Issue Array and Signal Processing for Radar)
Show Figures

Figure 1

Figure 1
<p>The CDVF flowchart.</p>
Full article ">Figure 2
<p>Flight airplane.</p>
Full article ">Figure 3
<p>Flight experiment trajectory(the red lines reprent the imaging area).</p>
Full article ">Figure 4
<p>Attitude estimation comparison.</p>
Full article ">Figure 5
<p>Velocity estimation comparison.</p>
Full article ">Figure 6
<p>Position estimation comparison.</p>
Full article ">Figure 7
<p>2D imaging.</p>
Full article ">Figure 8
<p>3D imaging after compensation.</p>
Full article ">
21 pages, 6196 KiB  
Article
Unimodular Multi-Input Multi-Output Waveform and Mismatch Filter Design for Saturated Forward Jamming Suppression
by Xuan Fang, Dehua Zhao and Liang Zhang
Sensors 2024, 24(18), 5884; https://doi.org/10.3390/s24185884 - 10 Sep 2024
Abstract
Forward jammers replicate and retransmit radar signals back to generate coherent jamming signals and false targets, making anti-jamming an urgent issue in electronic warfare. Jamming transmitters work at saturation to maximize the retransmission power such that only the phase information of the angular [...] Read more.
Forward jammers replicate and retransmit radar signals back to generate coherent jamming signals and false targets, making anti-jamming an urgent issue in electronic warfare. Jamming transmitters work at saturation to maximize the retransmission power such that only the phase information of the angular waveform at the designated direction of arrival (DOA) is retained. Therefore, amplitude modulation of MIMO radar angular waveforms offers an advantage in combating forward jamming. We address both the design of unimodular MIMO waveforms and their associated mismatch filters to confront mainlobe jamming in this paper. Firstly, we design the MIMO waveforms to maximize the discrepancy between the retransmitted jamming and the spatially synthesized radar signal. We formulate the problem as unconstrained non-linear optimization and solve it using the conjugate gradient method. Particularly, we introduce fast Fourier transform (FFT) to accelerate the numeric calculation of both the objection function and its gradient. Secondly, we design a mismatch filter to further suppress the filtered jamming through convex optimization in polynomial time. The simulation results show that for an eight-element MIMO radar, we are able to reduce the correlation between the angular waveform and saturated forward jamming to −6.8 dB. Exploiting this difference, the mismatch filter can suppress the jamming peak by 19 dB at the cost of an SNR loss of less than 2 dB. Full article
(This article belongs to the Section Radar Sensors)
Show Figures

Figure 1

Figure 1
<p>Flowchart of waveform and mismatch filter design.</p>
Full article ">Figure 2
<p>Signal reception and processing model for co-located MIMO radar with spatial synthesis of the waveform and a saturated forward jamming signal in the mainlobe.</p>
Full article ">Figure 3
<p>Comparison of waveform performance of each method in 10° for 8 × 128 codes when <span class="html-italic">η</span> = 1. (<b>a</b>) Autocorrelation level; (<b>b</b>) jamming cross-correlation level.</p>
Full article ">Figure 4
<p>Comparison of waveform performance of each method in 10° for 8 × 128 codes when <span class="html-italic">η</span> = 10. (<b>a</b>) Autocorrelation level; (<b>b</b>) jamming cross-correlation level.</p>
Full article ">Figure 5
<p>Performance comparison for different wave numbers with code length Ns = 128. (<b>a</b>) JCSL of WF0, WF1 and WF2 at different jamming intensities; (<b>b</b>) JCL of WF0, WF1 and WF2 at different jamming intensities at zero lag.</p>
Full article ">Figure 6
<p>Performance comparison of the two filters for WF0. (<b>a</b>) The WF0 matched filter; (<b>b</b>) the WF0 mismatch filter.</p>
Full article ">Figure 7
<p>Performance comparison of the two filters for WF1. (<b>a</b>)The WF1 matched filter; (<b>b</b>) The WF1 mismatch filter.</p>
Full article ">Figure 8
<p>Performance comparison of the two filters for WF2. (<b>a</b>) The WF2 matched filter; (<b>b</b>) the WF2 mismatch filter.</p>
Full article ">Figure 9
<p>Comparison of pulse compression outputs with the MF and MMF, of which the input signal is composed of a WF2 angular echo at 100 with χ = 1 V, 5 saturated forward jamming bins at 50, 150, 200, 250 and 300 with the jamming intensity <span class="html-italic">η</span> = 10 V and Gaussian noise with σ2 = 10 dBw.</p>
Full article ">
20 pages, 2969 KiB  
Article
Numerical Simulation of Non-Darcy Flow in Naturally Fractured Tight Gas Reservoirs for Enhanced Gas Recovery
by João Gabriel Souza Debossam, Mayksoel Medeiros de Freitas, Grazione de Souza, Helio Pedro Amaral Souto and Adolfo Puime Pires
Gases 2024, 4(3), 253-272; https://doi.org/10.3390/gases4030015 - 20 Aug 2024
Viewed by 368
Abstract
In this work, we analyze non-Darcy two-component single-phase isothermal flow in naturally fractured tight gas reservoirs. The model is applied in a scenario of enhanced gas recovery (EGR) with the possibility of carbon dioxide storage. The properties of the gases are obtained via [...] Read more.
In this work, we analyze non-Darcy two-component single-phase isothermal flow in naturally fractured tight gas reservoirs. The model is applied in a scenario of enhanced gas recovery (EGR) with the possibility of carbon dioxide storage. The properties of the gases are obtained via the Peng–Robinson equation of state. The finite volume method is used to solve the governing partial differential equations. This process leads to two subsystems of algebraic equations, which, after linearization and use of an operator splitting method, are solved by the conjugate gradient (CG) and biconjugate gradient stabilized (BiCGSTAB) methods for determining the pressure and fraction molar, respectively. We include inertial effects using the Barree and Conway model and gas slippage via a more recent model than Klinkenberg’s, and we use a simplified model for the effects of effective stress. We also utilize a mesh refinement technique to represent the discrete fractures. Finally, several simulations show the influence of inertial, slippage and stress effects on production in fractured tight gas reservoirs. Full article
Show Figures

Figure 1

Figure 1
<p>Three−dimensional mesh.</p>
Full article ">Figure 2
<p>Carbon dioxide mole fraction for sugar cube configuration: classical Darcy’s law. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 2 Cont.
<p>Carbon dioxide mole fraction for sugar cube configuration: classical Darcy’s law. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 3
<p>Carbon dioxide mole fraction for sugar cube configuration: inertial effects. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 3 Cont.
<p>Carbon dioxide mole fraction for sugar cube configuration: inertial effects. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 4
<p>Carbon dioxide mole fraction for sugar cube configuration: effective stress. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 4 Cont.
<p>Carbon dioxide mole fraction for sugar cube configuration: effective stress. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 5
<p>Carbon dioxide mole fraction for sugar cube configuration: gas slippage. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 5 Cont.
<p>Carbon dioxide mole fraction for sugar cube configuration: gas slippage. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 6
<p>Carbon dioxide mole fraction for sugar cube configuration: all effects combined. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">Figure 6 Cont.
<p>Carbon dioxide mole fraction for sugar cube configuration: all effects combined. (<b>a</b>) t = 1500 days. (<b>b</b>) t = 3000 days. (<b>c</b>) t = 4500 days. (<b>d</b>) t = 6000 days.</p>
Full article ">
22 pages, 4890 KiB  
Article
An Improved Three-Term Conjugate Gradient Algorithm for Constrained Nonlinear Equations under Non-Lipschitz Conditions and Its Applications
by Dandan Li, Yong Li and Songhua Wang
Mathematics 2024, 12(16), 2556; https://doi.org/10.3390/math12162556 - 19 Aug 2024
Viewed by 357
Abstract
This paper proposes an improved three-term conjugate gradient algorithm designed to solve nonlinear equations with convex constraints. The key features of the proposed algorithm are as follows: (i) It only requires that nonlinear equations have continuous and monotone properties; (ii) The designed search [...] Read more.
This paper proposes an improved three-term conjugate gradient algorithm designed to solve nonlinear equations with convex constraints. The key features of the proposed algorithm are as follows: (i) It only requires that nonlinear equations have continuous and monotone properties; (ii) The designed search direction inherently ensures sufficient descent and trust-region properties, eliminating the need for line search formulas; (iii) Global convergence is established without the necessity of the Lipschitz continuity condition. Benchmark problem numerical results illustrate the proposed algorithm’s effectiveness and competitiveness relative to other three-term algorithms. Additionally, the algorithm is extended to effectively address the image denoising problem. Full article
(This article belongs to the Special Issue Advances in Computational Mathematics and Applied Mathematics)
Show Figures

Figure 1

Figure 1
<p>Performance profiles for time.</p>
Full article ">Figure 2
<p>Performance profiles for Nfunc.</p>
Full article ">Figure 3
<p>Performance profiles for Niter.</p>
Full article ">Figure 4
<p>The noise images for lighthouse, peppers, boat, Kiel, fruits, and brain with 30% salt and pepper noise (first column) and the images recovered by Algorithms ITTCG (second column), HTTCGP (third column), and ZYL (forth column).</p>
Full article ">Figure 5
<p>The noise images for clown, couple, trucks, baboon, Barbara, and cameraman with 30% salt and pepper noise (first column) and the images recovered by Algorithms ITTCG (second column), HTTCGP (third column), and ZYL (forth column).</p>
Full article ">
16 pages, 312 KiB  
Article
An Efficient Subspace Minimization Conjugate Gradient Method for Solving Nonlinear Monotone Equations with Convex Constraints
by Taiyong Song and Zexian Liu
Axioms 2024, 13(3), 170; https://doi.org/10.3390/axioms13030170 - 6 Mar 2024
Viewed by 1007
Abstract
The subspace minimization conjugate gradient (SMCG) methods proposed by Yuan and Store are efficient iterative methods for unconstrained optimization, where the search directions are generated by minimizing the quadratic approximate models of the objective function at the current iterative point. Although the SMCG [...] Read more.
The subspace minimization conjugate gradient (SMCG) methods proposed by Yuan and Store are efficient iterative methods for unconstrained optimization, where the search directions are generated by minimizing the quadratic approximate models of the objective function at the current iterative point. Although the SMCG methods have illustrated excellent numerical performance, they are only used to solve unconstrained optimization problems at present. In this paper, we extend the SMCG methods and present an efficient SMCG method for solving nonlinear monotone equations with convex constraints by combining it with the projection technique, where the search direction is sufficiently descent.Under mild conditions, we establish the global convergence and R-linear convergence rate of the proposed method. The numerical experiment indicates that the proposed method is very promising. Full article
(This article belongs to the Special Issue Numerical Analysis and Optimization)
Show Figures

Figure 1

Figure 1
<p>Performance profilesof the five algorithms with respect to number of iterations (Ni).</p>
Full article ">Figure 2
<p>Performance profiles of the five algorithms with respect to number of function evaluations (NF).</p>
Full article ">Figure 3
<p>Performance profiles of the five algorithms with respect to CPU time (T).</p>
Full article ">
29 pages, 1311 KiB  
Article
Preconditioning Technique for an Image Deblurring Problem with the Total Fractional-Order Variation Model
by Adel M. Al-Mahdi
Math. Comput. Appl. 2023, 28(5), 97; https://doi.org/10.3390/mca28050097 - 22 Sep 2023
Cited by 1 | Viewed by 1466
Abstract
Total fractional-order variation (TFOV) in image deblurring problems can reduce/remove the staircase problems observed with the image deblurring technique by using the standard total variation (TV) model. However, the discretization of the Euler–Lagrange equations associated with the TFOV model generates a saddle point [...] Read more.
Total fractional-order variation (TFOV) in image deblurring problems can reduce/remove the staircase problems observed with the image deblurring technique by using the standard total variation (TV) model. However, the discretization of the Euler–Lagrange equations associated with the TFOV model generates a saddle point system of equations where the coefficient matrix of this system is dense and ill conditioned (it has a huge condition number). The ill-conditioned property leads to slowing of the convergence of any iterative method, such as Krylov subspace methods. One treatment for the slowness property is to apply the preconditioning technique. In this paper, we propose a block triangular preconditioner because we know that using the exact triangular preconditioner leads to a preconditioned matrix with exactly two distinct eigenvalues. This means that we need at most two iterations to converge to the exact solution. However, we cannot use the exact preconditioner because the Shur complement of our system is of the form S=K*K+λLα which is a huge and dense matrix. The first matrix, K*K, comes from the blurred operator, while the second one is from the TFOV regularization model. To overcome this difficulty, we propose two preconditioners based on the circulant and standard TV matrices. In our algorithm, we use the flexible preconditioned GMRES method for the outer iterations, the preconditioned conjugate gradient (PCG) method for the inner iterations, and the fixed point iteration (FPI) method to handle the nonlinearity. Fast convergence was found in the numerical results by using the proposed preconditioners. Full article
Show Figures

Figure 1

Figure 1
<p>Cross sections.</p>
Full article ">Figure 2
<p>Right box.</p>
Full article ">Figure 3
<p>Middle box.</p>
Full article ">Figure 4
<p>Left box.</p>
Full article ">Figure 5
<p>TV-error.</p>
Full article ">Figure 6
<p>TFOV-error.</p>
Full article ">Figure 7
<p>Eigenvalues of <span class="html-italic">A</span>.</p>
Full article ">Figure 8
<p>Eigenvalues of <math display="inline"><semantics> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> <mi>A</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>Golden house image.</p>
Full article ">Figure 10
<p>Retinal image.</p>
Full article ">Figure 11
<p>Shape of the kernel.</p>
Full article ">Figure 12
<p>Residual versus iterations number when <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>=</mo> <mn>1</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </mrow> </semantics></math>.</p>
Full article ">Figure 13
<p>Residual versus iterations number when <math display="inline"><semantics> <mrow> <mi>λ</mi> <mo>=</mo> <mn>1</mn> <mo>×</mo> <msup> <mn>10</mn> <mrow> <mo>−</mo> <mn>5</mn> </mrow> </msup> </mrow> </semantics></math>.</p>
Full article ">Figure 14
<p>Golden house image (blurred).</p>
Full article ">Figure 15
<p>Retinal image (blurred).</p>
Full article ">Figure 16
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>1</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 17
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>2</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 18
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>1</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.8</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 19
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>2</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.8</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 20
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>1</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 21
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>2</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 22
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>1</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.8</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 23
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>2</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.8</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 24
<p>FGMRES vs. GMRES.</p>
Full article ">Figure 25
<p>Peppers image (exact).</p>
Full article ">Figure 26
<p>Peppers image (blurred).</p>
Full article ">Figure 27
<p>Using TV (<math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>).</p>
Full article ">Figure 28
<p>Using NP.</p>
Full article ">Figure 29
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>1</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.9</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 30
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>2</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.9</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 31
<p><math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>64</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 32
<p><math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>128</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 33
<p><math display="inline"><semantics> <mrow> <mi>N</mi> <mo>=</mo> <mn>256</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 34
<p>Satel image (blurred).</p>
Full article ">Figure 35
<p>Using NFOV.</p>
Full article ">Figure 36
<p>Using NP.</p>
Full article ">Figure 37
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>1</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.9</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 38
<p>Using <math display="inline"><semantics> <msub> <mi>P</mi> <mn>2</mn> </msub> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>=</mo> <mn>1.9</mn> </mrow> </semantics></math>.</p>
Full article ">
27 pages, 10099 KiB  
Article
Analytical and Numerical Results for the Diffusion-Reaction Equation When the Reaction Coefficient Depends on Simultaneously the Space and Time Coordinates
by Ali Habeeb Askar, Ádám Nagy, Imre Ferenc Barna and Endre Kovács
Computation 2023, 11(7), 127; https://doi.org/10.3390/computation11070127 - 29 Jun 2023
Cited by 5 | Viewed by 1776
Abstract
We utilize the travelling-wave Ansatz to obtain novel analytical solutions to the linear diffusion–reaction equation. The reaction term is a function of time and space simultaneously, firstly in a Lorentzian form and secondly in a cosine travelling-wave form. The new solutions contain the [...] Read more.
We utilize the travelling-wave Ansatz to obtain novel analytical solutions to the linear diffusion–reaction equation. The reaction term is a function of time and space simultaneously, firstly in a Lorentzian form and secondly in a cosine travelling-wave form. The new solutions contain the Heun functions in the first case and the Mathieu functions for the second case, and therefore are highly nontrivial. We use these solutions to test some non-conventional explicit and stable numerical methods against the standard explicit and implicit methods, where in the latter case the algebraic equation system is solved by the preconditioned conjugate gradient and the GMRES solvers. After this verification, we also calculate the transient temperature of a 2D surface subjected to the cooling effect of the wind, which is a function of space and time again. We obtain that the explicit stable methods can reach the accuracy of the implicit solvers in orders of magnitude shorter time. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The shape function of Equation (5) for the different parameter sets <math display="inline"><semantics><mrow><mfenced><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>,</mo><mo> </mo><msub><mi>c</mi><mn>2</mn></msub><mo>,</mo><mo> </mo><mi>c</mi><mo>,</mo><mo> </mo><mi>a</mi><mo>,</mo><mo> </mo><mi>D</mi></mrow></mfenced></mrow></semantics></math>; the black, red, blue, and green lines are for (0, 1, 0.5, 7.4, 0.3), (0, 1, 2, 7.4, 0.3), (0, 1, 0.5, 7.4, 1.8), and (0, 1, 0.5, 2.4, 4.1), respectively.</p>
Full article ">Figure 2
<p>The solution function of Equation (6) for the parameters of <span class="html-italic">D</span> = 0.3; <span class="html-italic">a</span> = 7.4; <span class="html-italic">c</span> = 1;<math display="inline"><semantics><mrow><mfenced><mrow><mi mathvariant="sans-serif">Ψ</mi><mo>=</mo><mn>5.572222</mn></mrow></mfenced></mrow></semantics></math>; <span class="html-italic">c</span><sub>1</sub> = 0; <span class="html-italic">c</span><sub>2</sub> = 1, respectively.</p>
Full article ">Figure 3
<p>The stability diagram of the shape function <math display="inline"><semantics><mrow><mi>f</mi><mo>(</mo><mi>η</mi><mo>)</mo></mrow></semantics></math> in the form of Equation (8) for <math display="inline"><semantics><mrow><mi>c</mi><mo>=</mo><mn>1</mn></mrow></semantics></math> for the Mathieu S part <math display="inline"><semantics><mrow><mfenced><mrow><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo> </mo><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow></mfenced></mrow></semantics></math>. The horizontal axis shows the range of the parameter <span class="html-italic">a</span>, which is responsible for the strength of the source term, while the vertical axis shows the range of the <span class="html-italic">D</span>, which is the diffusion coefficient. In the white regions, the function the integral of the function is infinite and therefore not stable. In the colored region, the integral is finite and has the value as the color shows.</p>
Full article ">Figure 4
<p>The solution <math display="inline"><semantics><mrow><mi>u</mi><mfenced><mrow><mi>x</mi><mo>,</mo><mi>t</mi></mrow></mfenced></mrow></semantics></math> of Equation (2) with form of Equation (9) for <math display="inline"><semantics><mrow><mi>c</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo> </mo><mi>a</mi><mo>=</mo><mn>2</mn><mo>,</mo><mo> </mo><mi>D</mi><mo>=</mo><mn>6</mn></mrow></semantics></math>. Only the Mathieu S function is presented with <math display="inline"><semantics><mrow><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mn>1</mn></mrow></semantics></math>; the Mathieu C function looks similar.</p>
Full article ">Figure 5
<p>Hopscotch-type space–time structures. The time elapses from the top <math display="inline"><semantics><mrow><mfenced><mrow><mi>t</mi><mo>=</mo><msup><mi>t</mi><mn>0</mn></msup></mrow></mfenced></mrow></semantics></math> to the bottom <math display="inline"><semantics><mrow><mfenced><mrow><mi>t</mi><mo>=</mo><msup><mi>t</mi><mrow><mi>fin</mi></mrow></msup></mrow></mfenced></mrow></semantics></math>.</p>
Full article ">Figure 6
<p>Errors as a function of the temporal step size for Case study 1.</p>
Full article ">Figure 7
<p>The concentration <span class="html-italic">u</span> as a function of <span class="html-italic">x</span> in the case of the initial function <math display="inline"><semantics><mrow><msup><mi>u</mi><mn>0</mn></msup></mrow></semantics></math>, the exact analytical solution at <math display="inline"><semantics><mrow><msup><mi>t</mi><mrow><mi>fin</mi></mrow></msup></mrow></semantics></math>, the implicit-PCG algorithm (with Tolerance 10<sup>−5</sup>), and the LH scheme for <math display="inline"><semantics><mrow><mi>h</mi><mo>=</mo><mn>0.002</mn></mrow></semantics></math> in the case of small value of <span class="html-italic">a</span> (Case study 1). The maximum errors (and running times) are 0.0015 (98.7 s) and 0.0022 (0.107 s) for the PCG and LH methods, respectively.</p>
Full article ">Figure 8
<p>Errors as a function of the temporal step size for Case study 2.</p>
Full article ">Figure 9
<p>The concentration <span class="html-italic">u</span> as a function of <span class="html-italic">x</span> in the case of the initial function <math display="inline"><semantics><mrow><msup><mi>u</mi><mn>0</mn></msup></mrow></semantics></math>, the exact analytical solution at <math display="inline"><semantics><mrow><msup><mi>t</mi><mrow><mi>fin</mi></mrow></msup></mrow></semantics></math>, the implicit–GMRES algorithm, and the LH scheme. For more data, please see the main text.</p>
Full article ">Figure 10
<p>Errors as a function of the temporal step size for Case study 3.</p>
Full article ">Figure 11
<p>The concentration <span class="html-italic">u</span> as a function of <span class="html-italic">x</span> in the case of the initial function <math display="inline"><semantics><mrow><msup><mi>u</mi><mn>0</mn></msup></mrow></semantics></math>, the exact analytical solution at <math display="inline"><semantics><mrow><msup><mi>t</mi><mrow><mi>fin</mi></mrow></msup></mrow></semantics></math>, the implicit-PCM algorithm, and the LH scheme.</p>
Full article ">Figure 12
<p>The surface of the modelled wall with two layers.</p>
Full article ">Figure 13
<p>The temperature distribution contour in Kelvin units for the surface area (<b>upper half</b>) constant convection and (<b>lower half</b>) the convection changes with time depending on weather data.</p>
Full article ">Figure 14
<p>The maximum errors as a function of the time step size <span class="html-italic">h</span> for the two implicit methods with two different tolerances and the two explicit methods in the case of the 50 by 50 system.</p>
Full article ">Figure 15
<p>The maximum errors as a function of the time step size <span class="html-italic">h</span> for the examined methods for the 100 by 100 system.</p>
Full article ">Figure 16
<p>The maximum errors as a function of the running time for the tested methods for the 50 by 50 system.</p>
Full article ">Figure 17
<p>The maximum errors as a function of the running time for the tested methods for the 100 by 100 system.</p>
Full article ">Figure 18
<p>The running time with total number of cells for <math display="inline"><semantics><mrow><mi>h</mi><mo>=</mo><mn>1</mn><mo> </mo><mi mathvariant="normal">s</mi></mrow></semantics></math>. The left-hand side vertical axis refers to the LH and DF methods, while it is on the right-hand side axis for the implicit methods.</p>
Full article ">Figure 19
<p>The running time in the same logarithmic scale for all the used solvers as a function of the total number of cells.</p>
Full article ">Figure 20
<p>The error of temperature distribution contour in Kelvin units (<b>A</b>) LH, (<b>B</b>) DF, (<b>C</b>) GMRES, and (<b>D</b>) PCG methods.</p>
Full article ">Figure 21
<p>The error of temperature distribution at <span class="html-italic">z</span> = 80 [cm] in Kelvin units for LH, DF, GMRES, and PCG methods along <span class="html-italic">x</span>-axis.</p>
Full article ">
22 pages, 8746 KiB  
Article
Optimization Analysis of Overlap Ratio in Wide and Heavy Plate Slitting Shear
by Yachu Liao, Wei Fang, Jiahui Li, Zhang Dang, Meng Li and Wenbin Shi
Machines 2023, 11(7), 683; https://doi.org/10.3390/machines11070683 - 28 Jun 2023
Viewed by 1024
Abstract
In studying planar multi-bar mechanisms with multiple degrees of freedom, mathematical modeling is undoubtedly a way to get closer to the expected trajectory. Compared with the analytical method, the optimization method has higher accuracy in solving nonlinear equations, and it can be searched [...] Read more.
In studying planar multi-bar mechanisms with multiple degrees of freedom, mathematical modeling is undoubtedly a way to get closer to the expected trajectory. Compared with the analytical method, the optimization method has higher accuracy in solving nonlinear equations, and it can be searched and iterated in an extensive range until it meets the real engineering solution. The research on the overlap ratio of the Slitting Shear essentially aims to study the kinematics of the two-DOF mechanism. In this study, a one-DOF mathematical model of mechanical motion was established for the Slitting Shear without considering the overlap ratio. The mathematical model was then verified through simulation. The overlap ratio was introduced based on the above mathematic model, and the three-DOF mathematic model of the shear mechanism was established. Finally, the setting of the overlap ratio was optimized using the conjugate gradient method and the global optimal NMinmize function. Under the condition of satisfying the overlap ratio, the recommended value for the overlap ratio setting of the two-DOF mechanism, which approximates the pure rolling cut, was determined. Full article
(This article belongs to the Section Advanced Manufacturing)
Show Figures

Figure 1

Figure 1
<p>Shearing process of Slitting Shear.</p>
Full article ">Figure 1 Cont.
<p>Shearing process of Slitting Shear.</p>
Full article ">Figure 2
<p>Vector diagram of linkage mechanism.</p>
Full article ">Figure 3
<p>Common origin and any plane rectangular coordinate system.</p>
Full article ">Figure 4
<p>Schematic diagram of Slitting Shear cutting mechanism.</p>
Full article ">Figure 5
<p>Curve change of <span class="html-italic">θ</span><sub>3</sub>, <span class="html-italic">θ</span><sub>4</sub>, <span class="html-italic">θ</span><sub>5</sub>, <span class="html-italic">θ</span><sub>6</sub>.</p>
Full article ">Figure 6
<p>Angular velocity and angular acceleration solve the curve.</p>
Full article ">Figure 7
<p>Trajectory curve of the hinge point of the tool table.</p>
Full article ">Figure 8
<p>Coordinates of any point on the top knife sledge.</p>
Full article ">Figure 9
<p>Curve of the angle of the top knife sledge <span class="html-italic">θ</span><sub>8</sub>.</p>
Full article ">Figure 10
<p>Trajectory curve of the midpoint K and the endpoint M of the top blade arc.</p>
Full article ">Figure 11
<p>Dynamic lowest point trajectories and overlap ratio curves.</p>
Full article ">Figure 12
<p>Solidworks trajectory check of arc dynamic lowest point on the top blade.</p>
Full article ">Figure 13
<p>Adams simulation model.</p>
Full article ">Figure 14
<p>Adams key point simulation trajectory diagram.</p>
Full article ">Figure 14 Cont.
<p>Adams key point simulation trajectory diagram.</p>
Full article ">Figure 15
<p>Longitudinal displacement curve of each marked point.</p>
Full article ">Figure 16
<p>Dynamic longitudinal displacement dynamic lowest point of top blade arc.</p>
Full article ">Figure 17
<p>Shearing process of steel plate.</p>
Full article ">Figure 17 Cont.
<p>Shearing process of steel plate.</p>
Full article ">Figure 18
<p>Relationship between shear force and percentage of cut depth.</p>
Full article ">Figure 19
<p>Curve of the relationship between steel plate thickness and penetration depth in actual production.</p>
Full article ">Figure 20
<p>One side overlap ratio adjustment model.</p>
Full article ">Figure 21
<p>Diagram of shear mechanism with overlap ratio adjustment.</p>
Full article ">Figure 22
<p>Flow chart of optimization design.</p>
Full article ">
35 pages, 686 KiB  
Article
Adaptive Hybrid Mixed Two-Point Step Size Gradient Algorithm for Solving Non-Linear Systems
by Eltiyeb Ali and Salem Mahdi
Mathematics 2023, 11(9), 2102; https://doi.org/10.3390/math11092102 - 28 Apr 2023
Cited by 2 | Viewed by 1292
Abstract
In this paper, a two-point step-size gradient technique is proposed by which the approximate solutions of a non-linear system are found. The two-point step-size includes two types of parameters deterministic and random. A new adaptive backtracking line search is presented and combined with [...] Read more.
In this paper, a two-point step-size gradient technique is proposed by which the approximate solutions of a non-linear system are found. The two-point step-size includes two types of parameters deterministic and random. A new adaptive backtracking line search is presented and combined with the two-point step-size gradient to make it globally convergent. The idea of the suggested method depends on imitating the forward difference method by using one point to estimate the values of the gradient vector per iteration where the number of the function evaluation is at most one for each iteration. The global convergence analysis of the proposed method is established under actual and limited conditions. The performance of the proposed method is examined by solving a set of non-linear systems containing high dimensions. The results of the proposed method is compared to the results of a derivative-free three-term conjugate gradient CG method that solves the same test problems. Fair, popular, and sensible evaluation criteria are used for comparisons. The numerical results show that the proposed method has merit and is competitive in all cases and superior in terms of efficiency, reliability, and effectiveness in finding the approximate solution of the non-linear systems. Full article
(This article belongs to the Special Issue Numerical Analysis and Optimization: Methods and Applications)
Show Figures

Figure 1

Figure 1
<p>The BItr, BEFs, and BTcpu obtained by the HRSG1, HRSG2 and Algo methods.</p>
Full article ">Figure 2
<p>The WItr, WFEs and WTcpu obtained by the HRSG1 and HRSG2 methods.</p>
Full article ">Figure 3
<p>The MEItr, MEFEs and METcpu obtained by the HRSG1 and HRSG2 methods.</p>
Full article ">Figure 3 Cont.
<p>The MEItr, MEFEs and METcpu obtained by the HRSG1 and HRSG2 methods.</p>
Full article ">
13 pages, 4886 KiB  
Article
A Radial Basis Scale Conjugate Gradient Deep Neural Network for the Monkeypox Transmission System
by Zulqurnain Sabir, Salem Ben Said and Juan L. G. Guirao
Mathematics 2023, 11(4), 975; https://doi.org/10.3390/math11040975 - 15 Feb 2023
Cited by 6 | Viewed by 1301
Abstract
The motive of this study is to provide the numerical performances of the monkeypox transmission system (MTS) by applying the novel stochastic procedure based on the radial basis scale conjugate gradient deep neural network (RB-SCGDNN). Twelve and twenty numbers of neurons were taken [...] Read more.
The motive of this study is to provide the numerical performances of the monkeypox transmission system (MTS) by applying the novel stochastic procedure based on the radial basis scale conjugate gradient deep neural network (RB-SCGDNN). Twelve and twenty numbers of neurons were taken in the deep neural network process in first and second hidden layers. The MTS dynamics were divided into rodent and human, the human was further categorized into susceptible, infectious, exposed, clinically ill, and recovered, whereas the rodent was classified into susceptible, infected, and exposed. The construction of dataset was provided through the Adams method that was refined further by using the training, validation, and testing process with the statics of 0.15, 0.13 and 0.72. The exactness of the RB-SCGDNN is presented by using the comparison of proposed and reference results, which was further updated through the negligible absolute error and different statistical performances to solve the nonlinear MTS. Full article
(This article belongs to the Section Mathematical Physics)
Show Figures

Figure 1

Figure 1
<p>A deep neural network process, mathematical formulations, and assessment of outcomes for the MTS.</p>
Full article ">Figure 2
<p>A multilayers procedure for the nonlinear MTS.</p>
Full article ">Figure 3
<p>A neural network training performance, an input, couple hidden and output layers for solving the nonlinear dynamics of the MTS. (<b>a</b>) A neural network training performance, (<b>b</b>) An input, couple hidden and output layers to solve the MTS.</p>
Full article ">Figure 3 Cont.
<p>A neural network training performance, an input, couple hidden and output layers for solving the nonlinear dynamics of the MTS. (<b>a</b>) A neural network training performance, (<b>b</b>) An input, couple hidden and output layers to solve the MTS.</p>
Full article ">Figure 4
<p>Optimal training and ToS for solving the nonlinear dynamics of MTS.</p>
Full article ">Figure 5
<p>Function fitness for solving the nonlinear dynamics of MTS.</p>
Full article ">Figure 5 Cont.
<p>Function fitness for solving the nonlinear dynamics of MTS.</p>
Full article ">Figure 6
<p>The performances of EHs for each case the nonlinear dynamics of MTS.</p>
Full article ">Figure 7
<p>Reg performances for each case of the MTS.</p>
Full article ">Figure 7 Cont.
<p>Reg performances for each case of the MTS.</p>
Full article ">Figure 8
<p>Result comparisons for each class of the nonlinear dynamics of the MTS.</p>
Full article ">
21 pages, 9635 KiB  
Article
Hydrogeological Study in Tongchuan City Using the Audio-Frequency Magnetotelluric Method
by Zhimin Xu, Huicui Xin, Yuren Weng and Guang Li
Magnetochemistry 2023, 9(1), 32; https://doi.org/10.3390/magnetochemistry9010032 - 14 Jan 2023
Cited by 1 | Viewed by 1618
Abstract
Tongchuan City, located in Shaanxi Province, northwest China, has limited groundwater resources. Rational planning and exploitation of groundwater are crucial to the sustainable development of the city, for which investigating the distribution of groundwater is the premise. Traditional resistivity sounding methods are often [...] Read more.
Tongchuan City, located in Shaanxi Province, northwest China, has limited groundwater resources. Rational planning and exploitation of groundwater are crucial to the sustainable development of the city, for which investigating the distribution of groundwater is the premise. Traditional resistivity sounding methods are often used to detect groundwater; however, these methods are not applicable in the study area where thick Quaternary loess is extensively distributed. In this study, we arranged five audio-frequency magnetotelluric (AMT) profiles to detect the deep clastic rock groundwater and carbonate karst fissure groundwater in Tongchuan. Firstly, we analyzed the electromagnetic interference (EMI) noises in Tongchuan City, revealing that the main EMI is power frequency interference (PFI). We used the dictionary learning processing technology to suppress the PFI. Secondly, the two-dimensional (2D) nonlinear conjugate gradient method was employed to invert a 2D electrical structure model for the area shallower than 1 km. We analyzed the characteristics of the electrical structure and its geological significance. Lastly, the three-dimensional (3D) electrical structure model of the study area was inverted using the 3D nonlinear conjugate gradient method, and the spatial distribution characteristics of the water-bearing strata were further analyzed. The results show that the PFI in urban environment can be suppressed by the dictionary learning processing technology. In Tongchuan city, the distribution of clastic rock fissure water is controlled by folds and faults, as well as the thickness of sandstone layers, and that of the carbonate karst fissure water is mainly controlled by faults. On this basis, we infer that the water-bearing areas are in the middle east and south of the study area. Full article
(This article belongs to the Special Issue Advances in Magnetotelluric Analysis)
Show Figures

Figure 1

Figure 1
<p>The regional geological map and approximate location of Tongchuan study area in China (upper). Geological map of the Tongchuan area (lower). The black triangles are the locations of audio-frequency magnetotelluric (AMT) sites. Geological profiles A–A’ and B–B’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>] are marked by black lines.</p>
Full article ">Figure 2
<p>Apparent resistivity and phase curve of the data at AMT sites (<b>a</b>) TC211 and (<b>b</b>) TC 402 in the study area.</p>
Full article ">Figure 3
<p>Time series (<b>a</b>) TS2, (<b>b</b>) TS3, and (<b>c</b>) TS4 of AMT site TC211.</p>
Full article ">Figure 4
<p>Time series (<b>a</b>) TS2, (<b>b</b>) TS3, and (<b>c</b>) TS4 of AMT site TC402.</p>
Full article ">Figure 5
<p>Apparent resistivity phase curve of AMT site TC211 (<b>a</b>) before and (<b>b</b>) after noise removal using the FFT–IOMP method.</p>
Full article ">Figure 6
<p>Apparent resistivity phase curve of AMT site TC402 (<b>a</b>) before and (<b>b</b>) after noise removal using the FFT–IOMP method.</p>
Full article ">Figure 7
<p>Inversion results of profile TC-4. (<b>a</b>) Occam 1D inversion model; (<b>b</b>) 2D inversion results with uniform half-space as the initial model; (<b>c</b>) 2D inversion results with Occam 1D inversion model as the initial model.</p>
Full article ">Figure 8
<p>The 2D inversion models of (<b>a</b>) profile TC1JM and (<b>b</b>) profile TC-1 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile B–B’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 9
<p>The 2D inversion models of (<b>a</b>) profile TC3JM and (<b>b</b>) profile TC-3 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile A–A’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 9 Cont.
<p>The 2D inversion models of (<b>a</b>) profile TC3JM and (<b>b</b>) profile TC-3 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile A–A’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 10
<p>The 2D inversion result of profile TC-4. The red line marks the locations of faults, and the blue dotted line marks clastic rock fissure water bearing area.</p>
Full article ">Figure 11
<p>The 2D inversion model of profile TC-5. The red line marks the locations of faults. The pink dotted line and blue dotted line delineate the karst fissure water bearing area and clastic rock fissure water-bearing area, respectively.</p>
Full article ">Figure 12
<p>The same as <a href="#magnetochemistry-09-00032-f012" class="html-fig">Figure 12</a>, but for profile TC-6.</p>
Full article ">Figure 13
<p>AMT sites selected for 3D inversion. A–A’ and B–B’ are geological profiles.</p>
Full article ">Figure 14
<p>(<b>a</b>) The 2D slice of the 3D inversion results along profile TC-1. (<b>b</b>) The 2D inversion model of profile TC-1 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile B–B’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 14 Cont.
<p>(<b>a</b>) The 2D slice of the 3D inversion results along profile TC-1. (<b>b</b>) The 2D inversion model of profile TC-1 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile B–B’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 15
<p>(<b>a</b>) The 2D slice of 3D inversion results along profile TC-2. (<b>b</b>) The 2D inversion results of profile TC-2 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile B–B’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 15 Cont.
<p>(<b>a</b>) The 2D slice of 3D inversion results along profile TC-2. (<b>b</b>) The 2D inversion results of profile TC-2 [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile B–B’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 16
<p>(<b>a</b>) The 2D slice of 3D inversion results along TC-3 profile. (<b>b</b>) The TC-3 2D inversion results [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile A–A’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 16 Cont.
<p>(<b>a</b>) The 2D slice of 3D inversion results along TC-3 profile. (<b>b</b>) The TC-3 2D inversion results [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>]. (<b>c</b>) Geological profile A–A’ [<a href="#B7-magnetochemistry-09-00032" class="html-bibr">7</a>].</p>
Full article ">Figure 17
<p>Spatial distribution of sandstone obtained through 3D inversion: (<b>a</b>) sandstone with high water content and a resistivity of 70 <math display="inline"><semantics> <mrow> <mrow> <mi mathvariant="sans-serif">Ω</mi> </mrow> <mo>.</mo> <mi mathvariant="normal">m</mi> </mrow> </semantics></math>; (<b>b</b>) sandstone with low water content and a resistivity of 200 <math display="inline"><semantics> <mrow> <mrow> <mi mathvariant="sans-serif">Ω</mi> </mrow> <mo>.</mo> <mi mathvariant="normal">m</mi> </mrow> </semantics></math>.</p>
Full article ">
15 pages, 505 KiB  
Article
An Algorithm for Solving Zero-Sum Differential Game Related to the Nonlinear H Control Problem
by Vladimir Milić, Josip Kasać and Marin Lukas
Algorithms 2023, 16(1), 48; https://doi.org/10.3390/a16010048 - 10 Jan 2023
Cited by 1 | Viewed by 2276
Abstract
This paper presents an approach for the solution of a zero-sum differential game associated with a nonlinear state-feedback H control problem. Instead of using the approximation methods for solving the corresponding Hamilton–Jacobi–Isaacs (HJI) partial differential equation, we propose an algorithm that calculates [...] Read more.
This paper presents an approach for the solution of a zero-sum differential game associated with a nonlinear state-feedback H control problem. Instead of using the approximation methods for solving the corresponding Hamilton–Jacobi–Isaacs (HJI) partial differential equation, we propose an algorithm that calculates the explicit inputs to the dynamic system by directly performing minimization with simultaneous maximization of the same objective function. In order to achieve numerical robustness and stability, the proposed algorithm uses: quasi-Newton method, conjugate gradient method, line search method with Wolfe conditions, Adams approximation method for time discretization and complex-step calculation of derivatives. The algorithm is evaluated in computer simulations on examples of first- and second-order nonlinear systems with analytical solutions of H control problem. Full article
(This article belongs to the Collection Feature Papers in Algorithms)
Show Figures

Figure 1

Figure 1
<p>Time dependence of the state variable (Example 4.1).</p>
Full article ">Figure 2
<p>Control and uncertainty variables in dependence on the state variable (Example 4.1).</p>
Full article ">Figure 3
<p>Time dependence of the control and uncertainty variables (Example 4.1).</p>
Full article ">Figure 4
<p>Time dependence of the state variables (Example 4.2).</p>
Full article ">Figure 5
<p>Time dependence of the control and uncertainty variables (Example 4.2).</p>
Full article ">
14 pages, 2230 KiB  
Article
Spectral Methods in Nonlinear Optics Equations for Non-Uniform Grids Using an Accelerated NFFT Scheme
by Pedro Rodríguez, Manuel Romero, Antonio Ortiz-Mora and Antonio M. Díaz-Soriano
Symmetry 2023, 15(1), 47; https://doi.org/10.3390/sym15010047 - 24 Dec 2022
Viewed by 1300
Abstract
In this work, we propose the use of non-homogeneous grids in 1D and 2D for the study of various nonlinear physical equations using spectral methods. As is well known, the use of spectral methods allow a faster resolution of the problem via the [...] Read more.
In this work, we propose the use of non-homogeneous grids in 1D and 2D for the study of various nonlinear physical equations using spectral methods. As is well known, the use of spectral methods allow a faster resolution of the problem via the application of the ubiquitous Fast Fourier Transform (FFT) algorithm. We will center our investigation on the search of fast and accurate schemes to solve the spectral operators in the Fourier space. In particular, we will use the Conjugate Gradient (CG) iterative method, with a preconditioning matrix to accelerate the inversion process of the non-uniform Fast Fourier Transform (NFFT). As it will be shown, the results obtained are in good agreement with the expected values. Full article
(This article belongs to the Special Issue Symmetry in Nonlinear Optics: Topics and Advances)
Show Figures

Figure 1

Figure 1
<p>Relative error for the different conjugate gradients algorithms used, with and without the preconditioner weight matrix <math display="inline"><semantics> <mi mathvariant="bold-italic">W</mi> </semantics></math>, measured as the relative L2 error obtained in the process of direct and inverse Fourier transformation of the initial state input for the 1D Raman case.</p>
Full article ">Figure 2
<p>Relative error for the two different conjugate gradients algorithms used, with and without the preconditioner weight matrix <math display="inline"><semantics> <mi mathvariant="bold-italic">W</mi> </semantics></math>, measured as the relative L2 error obtained in the process of direct and inverse Fourier transformation of the initial state input for the BEC asymmetric case. The Delaunay triangulation is the weight generator for <math display="inline"><semantics> <mi mathvariant="bold-italic">W</mi> </semantics></math>.</p>
Full article ">Figure 3
<p>Results obtained after filtering a order two solitonic pulse, <math display="inline"><semantics> <mrow> <mi>ψ</mi> <mo>(</mo> <mi>z</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> <mo>=</mo> <mn>2</mn> <mi>s</mi> <mi>e</mi> <mi>c</mi> <mi>h</mi> <mo>(</mo> <mi>z</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>, using the wavelet filter (DAUB4, threshold = 0.0002), details in <a href="#app2-symmetry-15-00047" class="html-app">Appendix B</a>. The values have a similar appearance for the different homogeneous grids (right up inset), their size is inversely proportional to the density of each one of them.The selected coordinates are those for which the result of applying the filter is above the dotted line (left up inset). The selection for every layer is shown on the bottom. The less dense grid is the basis to build the interpolation process when the pulse evolves.</p>
Full article ">Figure 4
<p>Intensity plots for different evolution lengths of a order two solitonic pulse input, i.e., <math display="inline"><semantics> <mrow> <mi>ψ</mi> <mo>(</mo> <mi>z</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> <mo>=</mo> <mn>2</mn> <mi>s</mi> <mi>e</mi> <mi>c</mi> <mi>h</mi> <mo>(</mo> <mi>z</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics></math>; a Raman term with value <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <mi>R</mi> </msub> <mo>=</mo> <mn>0.003</mn> </mrow> </semantics></math> has been added to the NLSE in this case. The fission process is clear and in the same way, it is appreciated how the grid is split following each one of the components. Considering a number of 2048 nodes available, base grid of 256 and four layers, the percentage of points used was about 20%.</p>
Full article ">Figure 5
<p>Contour plot of the position density for the BEC condensate in the asymmetric trap potential with the condition <math display="inline"><semantics> <mrow> <msub> <mi>γ</mi> <mi>y</mi> </msub> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math>, for different evolution times. In the lower panel, the standard deviation for the same times is shown. In this asymmetric case, the condensate contracts and expands alternately on both axes.</p>
Full article ">Figure 6
<p>Evolution of the wavelet-generated grid for the BEC asymmetric case, displaying shape fluctuations according to the position density changes as it is shown in <a href="#symmetry-15-00047-f005" class="html-fig">Figure 5</a>. Darker-colored areas indicate denser zones in a four-layer grid. Considering a number of 512 × 512 nodes available, a base grid of 64 × 64 and four layers, the percentage of points used was about 4%.</p>
Full article ">Figure A1
<p>In this figure we can see the flow diagram corresponding to the generation of the non-homogeneous grid. The application of the filter to the different areas of the homogeneous grid produces areas of greater or lesser density.</p>
Full article ">
18 pages, 569 KiB  
Article
Descent Derivative-Free Method Involving Symmetric Rank-One Update for Solving Convex Constrained Nonlinear Monotone Equations and Application to Image Recovery
by Aliyu Muhammed Awwal, Adamu Ishaku, Abubakar Sani Halilu, Predrag S. Stanimirović, Nuttapol Pakkaranang and Bancha Panyanak
Symmetry 2022, 14(11), 2375; https://doi.org/10.3390/sym14112375 - 10 Nov 2022
Cited by 2 | Viewed by 1366
Abstract
Many practical applications in applied sciences such as imaging, signal processing, and motion control can be reformulated into a system of nonlinear equations with or without constraints. In this paper, a new descent projection iterative algorithm for solving a nonlinear system of equations [...] Read more.
Many practical applications in applied sciences such as imaging, signal processing, and motion control can be reformulated into a system of nonlinear equations with or without constraints. In this paper, a new descent projection iterative algorithm for solving a nonlinear system of equations with convex constraints is proposed. The new approach is based on a modified symmetric rank-one updating formula. The search direction of the proposed algorithm mimics the behavior of a spectral conjugate gradient algorithm where the spectral parameter is determined so that the direction is sufficiently descent. Based on the assumption that the underlying function satisfies monotonicity and Lipschitz continuity, the convergence result of the proposed algorithm is discussed. Subsequently, the efficiency of the new method is revealed. As an application, the proposed algorithm is successfully implemented on image deblurring problem. Full article
(This article belongs to the Special Issue Nonlinear Analysis and Its Applications in Symmetry II)
Show Figures

Figure 1

Figure 1
<p>Performance profile for ITER.</p>
Full article ">Figure 2
<p>Performance profile for FVAL.</p>
Full article ">Figure 3
<p>Performance profile for TIME.</p>
Full article ">Figure 4
<p>First row (original images), second row (blurred images), third and last rows (recovered images by methods DFSR1 and CGD, respectively).</p>
Full article ">
18 pages, 668 KiB  
Article
Conjugate Gradient Algorithm for Least-Squares Solutions of a Generalized Sylvester-Transpose Matrix Equation
by Kanjanaporn Tansri and Pattrawut Chansangiam
Symmetry 2022, 14(9), 1868; https://doi.org/10.3390/sym14091868 - 7 Sep 2022
Cited by 5 | Viewed by 1708
Abstract
We derive a conjugate-gradient type algorithm to produce approximate least-squares (LS) solutions for an inconsistent generalized Sylvester-transpose matrix equation. The algorithm is always applicable for any given initial matrix and will arrive at an LS solution within finite steps. When the matrix equation [...] Read more.
We derive a conjugate-gradient type algorithm to produce approximate least-squares (LS) solutions for an inconsistent generalized Sylvester-transpose matrix equation. The algorithm is always applicable for any given initial matrix and will arrive at an LS solution within finite steps. When the matrix equation has many LS solutions, the algorithm can search for the one with minimal Frobenius-norm. Moreover, given a matrix Y, the algorithm can find a unique LS solution closest to Y. Numerical experiments show the relevance of the algorithm for square/non-square dense/sparse matrices of medium/large sizes. The algorithm works well in both the number of iterations and the computation time, compared to the direct Kronecker linearization and well-known iterative methods. Full article
(This article belongs to the Section Mathematics)
Show Figures

Figure 1

Figure 1
<p>The logarithm of the relative error <math display="inline"><semantics> <msub> <mrow> <mo>∥</mo> <msub> <mi>R</mi> <mi>r</mi> </msub> <mo>∥</mo> </mrow> <mi>F</mi> </msub> </semantics></math> for Example 1.</p>
Full article ">Figure 2
<p>The logarithm of the relative error for Example 2.</p>
Full article ">Figure 3
<p>The logarithm of the relative error for Example 3 with Y = 0.1 <math display="inline"><semantics> <mrow> <mo>×</mo> <mo form="prefix">ones</mo> <mo>(</mo> <mi>n</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The logarithm of the relative error for Example 3 with Y = <span class="html-italic">I</span>.</p>
Full article ">
Back to TopTop