Квадратная и кубическая аппроксимация функции

@Roozevelt · Регистрация: 22.10.2013

Студворк — интернет-сервис помощи студентам

Добрый день есть задача и решение к ней, помогите пожалуйста написать код к другой похожей задаче.
Вложение 448074

C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
 
namespace laba1
{
    class Program
    {
        public static void output(double[] mas, double[] mas1, double[] maznamen, double[] mas3)
        {
            double f, p1, p2;
            int j = 0;
            Console.WriteLine("1.по формулi");
            Console.WriteLine("2.Лагранж");
            Console.WriteLine("3.кубическая апроксимация");
            Console.WriteLine("4.1-а производная по формулi");
            Console.WriteLine("5.1-а производная апроксимация");
            Console.WriteLine("6.2-а производная по формулi");
            Console.WriteLine("7.2-а производная апроксимация");
            for (double i = 0.12; i < 1.44; i += 0.06)
            {
                f = Math.Log(i * i + i + 1); //наша функция
                p1 = (2 * i + 1) / (i * i + i + 1);//1-a производная
                p2 = -((2 * i * i + 2 * i - 1) / (i * i * i * i + 2 * i * i * i + 3 * i * i + 2 * i + 1));//2-a производная
                Console.Write("{0:F3} \t", f);
                Console.Write("{0:F3}   \t", mas[j]);//за Лагранжом            
                Console.Write("{0:F3}   \t", mas1[j]);//кубическая апроксимация
                Console.Write("{0:F3}   \t", p1);//1 призводная
                Console.Write("{0:F3}  \t", maznamen[j]);
                Console.Write("{0:F3}  \t", p2);
                Console.WriteLine("{0:F3}", mas3[j]);
                j++;
            }
        }
        public static double[] koef_lagr(double x, double[] dot_s, double[] f_x)//возвращает значения многочлена Лагранжа
        {
            double[] koef = new double[10];
            for (int i = 0; i < 10; i++)
            {
                double chiseln = 1, znamen = 1;
                for (int j = 0; j < 10; j++)
                {
                    if (i != j)
                    {
                        chiseln = chiseln * (x - dot_s[j]);
                        znamen = znamen * (dot_s[i] - dot_s[j]);
                    }
                }
                koef[i] = f_x[i] * chiseln / znamen;
            }
            return koef;
        }
        public static double[] Lagranj(double[] dot_s, double[] f_x)//апроксимация за Лагранжом
        {
 
            double[] za = new double[25];
            double[] koef = new double[10];
            int q = 0;
            for (double j = 0.12; j < 1.5; j += 0.06)
            {
                double s = 0;
                koef = koef_lagr(j, dot_s, f_x);
                for (int k = 0; k < 10; k++)
                    s += koef[k];
                za[q] = s;
                q++;
            }
            return za;
        }
        public static double[] koef_aprocs(double x, double[] dot_s, double[] f_x)//коефициенты кубической апроксимации
        {
            double[] koef = new double[4];
            int k = 0;
            for (int i = 0; i < dot_s.Length - 1; i++)
            {
                if ((x > dot_s[i]) && (x < dot_s[i + 1]))
                {
                    switch (i)
                    {
                        case 0: k = 0; break;
                        case 1: k = 0; break;
                        case 2: k = 1; break;
                        case 3: k = 2; break;
                        case 4: k = 3; break;
                        case 5: k = 4; break;
                        case 6: k = 5; break;
                        case 7: k = 6; break;
                        case 8: k = 6; break;
                        case 9: k = 6; break;
                    }
                    break;
                }
            }
            for (int i = 0; i < 4; i++)
            {
                double chiseln = 1, znamen = 1;
                for (int j = 0; j < 4; j++)
                {
                    if (i != j)
                    {
                        chiseln = chiseln * (x - dot_s[j + k]);
                        znamen = znamen * (dot_s[i + k] - dot_s[j + k]);
                    }
                }
                koef[i] = f_x[i + k] * chiseln / znamen;
            }
            return koef;
        }
        public static double[] Aprocsim(double[] dot_s, double[] f_x)//кубическая апроксимация
        {
            double[] za = new double[25];
            double[] koef = new double[4];
            int q = 0;
            for (double j = 0.12; j < 1.5; j += 0.06)
            {
                double s = 0;
                koef = koef_aprocs(j, dot_s, f_x);
                for (int k = 0; k < 4; k++)
                    s += koef[k];
                za[q] = s;
                q++;
            }
            return za;
        }
        public static double[] pohidna_1(double[] znach2)
        {
            double[] za1 = new double[25];
            for (int i = 0; i < znach2.Length - 1; i++)
            {
                if (i == 0)
                    za1[i] = (znach2[i + 1] - znach2[i]) / 0.06;//крайняя 1-а разничная производная
                else
                    za1[i] = (znach2[i + 1] - znach2[i - 1]) / 2 / 0.06;//центральная 1-я разничная производная
            }
            return za1;
        }
        public static double[] pohidna_2(double[] znach)
        {
            double[] z = new double[25];
            for (int i = 0; i < znach.Length - 1; i++)
            {
                if (i == 0)
                    z[i] = Math.Abs((pohidna_1(znach)[i + 1] - pohidna_1(znach)[i])) / 0.06;//левая 2-я разничная производная
                else
 
                    if (i + 1 == znach.Length)
                        z[i + 1] = Math.Abs((pohidna_1(znach)[i + 1] - pohidna_1(znach)[i])) / 0.06;//права 2-я разничная производная
                    else
                        z[i] = Math.Abs((znach[i + 1] + znach[i - 1] - 2 * znach[i])) / 0.06 / 0.06;//центральная 2-я разничная производная
            }
            return z;
        }
        static void Main(string[] args)
        {
            double[] dot_s = new double[10] { 0.5, 0.69, 0.78, 0.99, 1.21, 1.34, 1.51, 1.63, 1.71, 1.83 };
            double[] f_x = new double[10] { 0.55962, 0.77293, 0.87062, 1.0886, 1.30131, 1.41963, 1.566561, 1.66523, 1.72884, 1.82114 };
            double[] znach = Lagranj(dot_s, f_x);
            double[] znach2 = Aprocsim(dot_s, f_x);
            double[] znach3 = pohidna_1(znach2);
            double[] znach4 = pohidna_2(znach2);
            output(znach, znach2, znach3, znach4);
            Console.ReadKey();
        }
    }
}

Вот собственно и она Вложение 448077
Заранее благодарен!

@Roozevelt · 28.10.2014, 21:45 **[ТС]**

Подскажите пожалуйста как записать в C# квадратную и кубическую аппроксимацию функции

@tezaurismosis · 28.10.2014, 21:56

Комментарий модератора

Roozevelt, перепечатайте текст задания в тему,
задания в виде картинок недопустимы.
Языки форума - русский и английский.

@Roozevelt · 28.10.2014, 21:59 **[ТС]**

Задание: Подскажите пожалуйста как записать в C# квадратную и кубическую аппроксимацию функцию.

Новые блоги и статьи Все статьи Все блоги /
Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .

@Roozevelt 2 / 2 / 2 Регистрация: 22.10.2013 Сообщений: 109
	28.10.2014, 21:45 [ТС]
	Подскажите пожалуйста как записать в C# квадратную и кубическую аппроксимацию функции 0

@Roozevelt 2 / 2 / 2 Регистрация: 22.10.2013 Сообщений: 109
	28.10.2014, 21:59 [ТС]
	Задание: Подскажите пожалуйста как записать в C# квадратную и кубическую аппроксимацию функцию. 0