Mathematical modeling in MATLAB of solar activity cycles according to the growth-decline of the Wolf number

В этой статье проводится математическое моделирование динамики солнечной активности. Исследуются данные наблюдений по средне-ежемесячному числу солнечных пятен, называемых числом Вольфа, в период за 24.5 года с мая 1996 года по октябрь 2022 года. Исходя из результатов подобного исследования данных по этому процессу, с применением уравнения Риккати дробного постоянного порядка, о том, что подъем и падение числа Вольфа со временем происходит по кривой, очень близкой к обобщенной логистической кривой, в данной статье также предлагается математическая модель, основанная на уравнении Риккати. Так как уравнение Риккати хорошо описывает процессы, которые подчиняются логистическому закону. Однако уравнение обобщается до интегро-дифференциального уравнения Риккати, введением дробной производной типа Герасимова-Капуто переменного порядка, а дробная производная с переменным порядком, позволяет получить более уточную математическую модель циклов числа Вольфа с насыщением, и позволяет учесть эффект переменной памяти. Все расчёты моделей, обработка данных и визуализации проводятся в программе FDRE 3.0 разработанной в пакете MATLAB. Параметры моделирования уточняются аппроксимацией известных исследуемых данных, при помощи регрессионного анализа. В результате модельные кривые и графики известных за 24.5 года наблюдаемых данных, показывают между собой хорошее соответствие. С помощью уточнённой математической модели делается прогноз на следующие 9 лет, который визуально хорошо согласуется с известными модельными результатами солнечной активности. In this article, mathematical modeling of the dynamics of solar activity is carried out. Observational data on the average monthly number of sunspots, called the Wolf number, for the period of 24.5 years from May 1996 to October 2022 are studied. Based on the results of a similar study of data on this process, using the Riccati equation of a fractional constant order, that the rise and fall of the Wolf number over time occurs along a curve very close to the generalized logistic curve, this article also proposes a mathematical model based on the Riccati equation. Since the Riccati equation describes well the processes that obey the logistic law. However, the equation is generalized to the integro-differential Riccati equation by introducing a fractional derivative of the Gerasimov-Caputo type of variable order, and a fractional derivative with a variable order, allows you to get a more precise mathematical model of Wolf number cycles with saturation, and allows you to take into account the effect of variable memory. All model calculations, data processing and visualization are carried out in the FDRE 3.0 program developed in the MATLAB package. Modeling parameters are refined by approximation of known data under study, using regression analysis. As a result, the model curves and graphs of the observed data known for 24.5 years show good agreement with each other. With the help of a refined mathematical model, a forecast is made for the next 9 years, which visually agrees well with the known model results of solar activity.