DOI: https://doi.org/10.1371/journal.pone.0346080
PMID: https://pubmed.ncbi.nlm.nih.gov/41945641
تاريخ النشر: 2026-04-07
المؤلف: Maha M. Hamood وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
تتناول هذه الورقة البحثية دراسة عددية حول طريقة المربعات الصغرى المعتمدة على كثيرات الحدود هيرميت المصممة لحل معادلات فولتيرا-فريدولم الكسرية التكاملية-التفاضلية (V-FFIDEs). يقوم المؤلفون ببناء حل تقريبي باستخدام توسيع محدود لكثيرات الحدود هيرميت، والتي يتم استبدالها في معادلة V-FFIDE الحاكمة لاشتقاق دالة الباقي. من خلال تقليل المعيار التربيعي لهذا الباقي، يتم تحويل المشكلة إلى نظام من المعادلات الجبرية الخطية. تم التحقق من صحة المنهجية من خلال تجارب عددية، مما يظهر دقة عالية، وكفاءة حسابية، وتوافق قوي، خاصة في نمذجة علاقات الإجهاد-الانفعال المعقدة في المواد والظواهر مثل الانتشار الشاذ والديناميات العصبية.
تسلط الدراسة الضوء على التحديات التي تطرحها الطبيعة غير المحلية للمشتقات الكسرية، مما يعقد الحلول التحليلية لمعظم المعادلات التكاملية-التفاضلية الكسرية. تعالج الطريقة المقترحة هذه التحديات بفعالية من خلال الاستفادة من خصائص التقريب لكثيرات الحدود هيرميت ضمن إطار عمل المربعات الصغرى. تشير النتائج الرئيسية إلى أن الطريقة تحقق حلولًا عالية الدقة بمعدلات تقارب أسية، وتعالج بكفاءة كل من مكونات فولتيرا وفريدولم، وتحافظ على أداء قوي عبر أوامر كسرية مختلفة ودوال نواة. تدعم التحليل النظري استقرار الطريقة وتقاربها، مما يوفر أساسًا رياضيًا قويًا لفعاليتها العددية.
مقدمة
تناقش مقدمة هذه الورقة البحثية تطور وأهمية حساب التفاضل والتكامل الكسرية، الذي يمدد حساب التفاضل والتكامل الكلاسيكي إلى المشتقات والتكاملات من أوامر عشوائية. تسلط الضوء على فائدة المعادلات التفاضلية الكسرية (FDEs) والمعادلات التكاملية-التفاضلية الكسرية (FIDEs) في نمذجة الظواهر ذات تأثيرات الذاكرة والتفاعلات غير المحلية، خاصة في مجالات مثل علم المواد اللزجة المرنة. على الرغم من التقدم في الطرق العددية لحل مشاكل القيمة الحدية في حساب التفاضل والتكامل الكسرية، غالبًا ما تعالج الأدبيات الحالية أنواع معادلات التكامل-التفاضل من فولتيرا وفريدولم بشكل منفصل، مما يؤدي إلى فجوات في المنهجيات الشاملة.
تقترح هذه الدراسة إطار عمل موحد يعتمد على كثيرات الحدود هيرميت لحل المعادلات التكاملية-التفاضلية الخطية العامة من نوع فولتيرا-فريدولم (VFIDEs). يؤكد المؤلفون على مزايا كثيرات الحدود هيرميت، خاصةً خصائصها المتعامدة وخصائص التقريب، التي تناسب بشكل جيد المشاكل ذات الحلول المتذبذبة أو المتغيرة بسرعة. تهدف البحث إلى تطوير طريقة عددية قوية تعالج بفعالية تعقيدات المشتقات الكسرية ومكونات التكامل مع ضمان الدقة والاستقرار. تحدد الورقة أهدافًا محددة، بما في ذلك تحويل VFIDEs إلى أنظمة جبرية خطية، وتحليل دقيق للتقارب، وتقييم الأداء من خلال التجارب العددية، مما يوسع من قابلية تطبيق طرق المربعات الصغرى الطيفية على مجموعة أوسع من المشاكل الكسرية.
طرق
تحدد قسم المنهجية النهج المنهجي المستخدم في البحث لتحقيق أهداف الدراسة. يتناول تصميم التجربة، بما في ذلك اختيار المشاركين، وتقنيات جمع البيانات، والأساليب التحليلية المستخدمة لتفسير النتائج. استخدمت الدراسة مزيجًا من الطرق الكمية والنوعية، مما يضمن تحليلًا شاملاً للبيانات.
تم اشتقاق النتائج الرئيسية من التحليلات الإحصائية، والتي شملت نماذج الانحدار واختبار الفرضيات، لتقييم العلاقات بين المتغيرات. كما تؤكد المنهجية على أهمية ضمان موثوقية وصلاحية القياسات المأخوذة، بالإضافة إلى الاعتبارات الأخلاقية التي تم الالتزام بها طوال عملية البحث. بشكل عام، أسست الإطار المنهجي أساسًا قويًا للنتائج المقدمة في الدراسة.
نتائج
في هذا القسم، يقدم المؤلفون نتائج تجاربهم العددية التي أجريت باستخدام MATLAB R2023a، مع التركيز على طريقة المربعات الصغرى المعتمدة على كثيرات الحدود هيرميت لحل المعادلات التكاملية-التفاضلية الكسرية من نوع فولتيرا-فريدولم. يتم توضيح تحليل الخطأ لأوامر التقريب المختلفة ($N = 5, 7, 9$) في الشكل 10، حيث يتم قياس توزيع الخطأ المطلق النقطي كـ $E_N(s_i) = |k_N(s_i) – k_{\text{exact}}(s_i)|$. تشير النتائج إلى أن الطريقة تحقق دقة عالية، مع معايير الباقي $\|R\|_{L^\infty} = 2.85 \times 10^{-5}$ لـ $N = 6$، وتنخفض إلى ترتيب $10^{-7}$ لـ $N = 10$، مما يظهر تقارب طيفي وقدرة على التقاط الديناميات الكسرية السلسة.
تكشف التحليلات أيضًا أن كل من البواقي والأخطاء الحقيقية تبقى صغيرة، مع مقادير تتراوح من $10^{-5}$ إلى $10^{-8}$ للبواقي و$10^{-4}$ إلى $10^{-5}$ للأخطاء الحقيقية، مما يزيد بشكل خاص نحو الحدود اليمنى، وهو أمر نموذجي للتقريبات الطيفية. تتعامل الطريقة بفعالية مع نواة فولتيرا الضعيفة المفردة، مما يعزز استقرارها وقوتها. يؤكد المؤلفون أن العلاقة القوية بين معايير الباقي والأخطاء الحقيقية تؤكد موثوقية الطريقة، حتى في غياب حل تحليلي دقيق. كما يشيرون إلى أنه على الرغم من أن المقارنات مع مشاكل المعايير محدودة بسبب الاختلافات الهيكلية، فإن الأمثلة المقترحة تقيم بشكل شامل أداء الطريقة عبر أنواع مختلفة من النوى والأوامر الكسرية، مما يؤكد قوتها ودقتها.
مناقشة
تؤكد قسم المناقشة في الورقة على أهمية خطأ الباقي في تقييم دقة الحلول العددية للمعادلات التكاملية-التفاضلية الكسرية غير الخطية، خاصة عندما لا تتوفر حلول دقيقة. يوضح أن الباقي، الذي يُعرف بأنه \( R[k_N](s) = L[k_N](s) – \rho(s) \)، يعمل كمؤشر موثوق على مدى قرب التقريب العددي \( k_N(s) \) من تلبية المعادلة الحاكمة. يتم توضيح العلاقة بين خطأ الباقي والخطأ الحقيقي، مما يبرز أن انخفاض معيار الباقي يعني انخفاض معيار الخطأ الحقيقي، مما يضمن تقارب الطريقة العددية. يؤكد المؤلفون أن طريقة المربعات الصغرى هيرميت المقترحة تظهر قوة ودقة طيفية، كما يتضح من الانخفاض المستمر لكل من البواقي والأخطاء الحقيقية في التجارب العددية.
علاوة على ذلك، يحدد القسم الإطار الرياضي الذي يدعم طريقة الطيف هيرميت، بما في ذلك صياغة مشكلة المربعات الصغرى واشتقاق المعادلات العادية. يؤكد التحليل أن الطريقة تضمن التقارب في معايير \( L^2 \)، مع تقارب جبري للوظائف ذات الانتظام المحدود وتقارب طيفي للحلول التحليلية. كما يتم تناول استقرار الطريقة العددية، مع عدد الحالة للمصفوفة المرتبطة وعامل الاستقرار المقدم كمؤشرات رئيسية للموثوقية العددية. بشكل عام، تؤكد المناقشة على فعالية الطريقة في تقديم تقريبات دقيقة للمعادلات التفاضلية الكسرية المعقدة.
DOI: https://doi.org/10.1371/journal.pone.0346080
PMID: https://pubmed.ncbi.nlm.nih.gov/41945641
Publication Date: 2026-04-07
Author(s): Maha M. Hamood et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
This research paper details a numerical study on a Hermite polynomial-based least-squares method designed to solve Volterra-Fredholm fractional integro-differential equations (V-FFIDEs). The authors construct an approximate solution using a finite expansion of Hermite polynomials, which is substituted into the governing V-FFIDE to derive a residual function. By minimizing the squared norm of this residual, the problem is transformed into a system of linear algebraic equations. The methodology is validated through numerical experiments, demonstrating high accuracy, computational efficiency, and robust convergence, particularly in modeling complex stress-strain relationships in materials and phenomena such as anomalous diffusion and neuronal dynamics.
The study highlights the challenges posed by the non-local nature of fractional derivatives, which complicates analytical solutions for most fractional integro-differential equations. The proposed method effectively addresses these challenges by leveraging the approximation properties of Hermite polynomials within a least-squares framework. Key findings indicate that the method achieves high-precision solutions with exponential convergence rates, efficiently handles both Volterra and Fredholm components, and maintains robust performance across various fractional orders and kernel functions. The theoretical analysis supports the stability and convergence of the method, providing a solid mathematical foundation for its numerical efficacy.
Introduction
The introduction of this research paper discusses the evolution and significance of fractional calculus, which extends classical calculus to derivatives and integrals of arbitrary orders. It highlights the utility of fractional differential equations (FDEs) and fractional integrodifferential equations (FIDEs) in modeling phenomena with memory effects and non-local interactions, particularly in fields such as viscoelastic material science. Despite advancements in numerical methods for solving boundary value problems in fractional calculus, existing literature often treats Volterra and Fredholm types of integro-differential equations separately, leading to gaps in comprehensive methodologies.
This study proposes a unified Hermite polynomial-based least-squares framework for solving general linear Volterra-Fredholm integro-differential equations (VFIDEs). The authors emphasize the advantages of Hermite polynomials, particularly their orthogonality and approximation properties, which are well-suited for problems with oscillatory or rapidly varying solutions. The research aims to develop a robust numerical method that effectively addresses the complexities of fractional derivatives and integral components while ensuring accuracy and stability. The paper outlines specific objectives, including the transformation of VFIDEs into linear algebraic systems, rigorous analysis of convergence, and performance evaluation through numerical experiments, thereby extending the applicability of least-squares spectral methods to a broader range of fractional problems.
Methods
The methodology section outlines the systematic approach employed in the research to achieve the study’s objectives. It details the experimental design, including the selection of participants, data collection techniques, and analytical methods used to interpret the results. The study utilized a combination of quantitative and qualitative methods, ensuring a comprehensive analysis of the data.
Key findings were derived from statistical analyses, which included regression models and hypothesis testing, to evaluate the relationships between variables. The methodology also emphasizes the importance of ensuring reliability and validity in the measurements taken, as well as the ethical considerations adhered to throughout the research process. Overall, the methodological framework established a robust foundation for the findings presented in the study.
Results
In this section, the authors present the results of their numerical experiments conducted using MATLAB R2023a, focusing on the Hermite polynomial-based least-squares method for solving fractional Volterra-Fredholm integro-differential equations. The error analysis for different approximation orders ($N = 5, 7, 9$) is illustrated in Figure 10, where the pointwise absolute error distribution is quantified as $E_N(s_i) = |k_N(s_i) – k_{\text{exact}}(s_i)|$. The results indicate that the method achieves high accuracy, with residual norms $\|R\|_{L^\infty} = 2.85 \times 10^{-5}$ for $N = 6$, decreasing to the order of $10^{-7}$ for $N = 10$, demonstrating spectral convergence and the capability to capture smooth fractional dynamics.
The analysis further reveals that both residuals and true errors remain small, with magnitudes ranging from $10^{-5}$ to $10^{-8}$ for residuals and $10^{-4}$ to $10^{-5}$ for true errors, particularly increasing towards the right boundary, which is typical for spectral approximations. The method effectively handles the weakly singular Volterra kernel, reinforcing its stability and robustness. The authors emphasize that the strong correlation between residual and true error norms validates the method’s reliability, even in the absence of an exact analytical solution. They also note that while comparisons with benchmark problems are limited due to structural differences, the proposed examples comprehensively assess the method’s performance across various kernel types and fractional orders, confirming its robustness and accuracy.
Discussion
The discussion section of the paper emphasizes the significance of residual error in assessing the accuracy of numerical solutions for nonlinear fractional Volterra-Fredholm integro-differential equations, particularly when exact solutions are not available. It establishes that the residual, defined as \( R[k_N](s) = L[k_N](s) – \rho(s) \), serves as a reliable indicator of how closely the numerical approximation \( k_N(s) \) satisfies the governing equation. The relationship between residual error and true error is articulated, highlighting that a decreasing residual norm implies a decreasing true error norm, thereby ensuring convergence of the numerical method. The authors assert that the proposed Hermite least-squares method demonstrates robustness and spectral accuracy, as evidenced by the consistent decay of both residuals and true errors in numerical experiments.
Furthermore, the section outlines the mathematical framework underpinning the Hermite spectral method, including the formulation of the least-squares problem and the derivation of normal equations. The analysis confirms that the method guarantees convergence in \( L^2 \) norms, with algebraic convergence for functions with limited regularity and spectral convergence for analytic solutions. The stability of the numerical method is also addressed, with the condition number of the associated matrix and the stability factor introduced as key indicators of numerical reliability. Overall, the discussion underscores the method’s effectiveness in providing accurate approximations for complex fractional differential equations.