Помогите переписать процедуру Spline с Pascal на C#

@Bar-bar1an · Регистрация: 12.12.2014

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

Помогите пожалуйста переписать процедуру Spline на C#
вот код паскаль:

Pascal
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
 Vector = Array [0..640] of Real;
      Vec   = Array [0..30]  of Real;
      Matr = Array [ 1..30,1..30] of Real;  Ne,Ngr :Integer;
   Xe,Ye,Xg,Yg:Vector;
   M :Integer;
   al,bl:Real;
       K:Vec; 
     Procedure Spline(Xe, Ye : Vector; Ne, Ngr : Integer; Var Xg, Yg : Vector);
  Var  a, b, c, d : Vec;
             h, R : Real;
          i, j, L : Integer;
          //cubspln
Procedure CubSpln(Xe,Ye:Vector;Ne:Integer;Var a,b,c,d:Vec);
 
Var      h, Y : Vec;
        Ma : Matr;
      i, j : Integer;
      R, Q : Real;
      //progonka
Procedure Progonka (A:Matr; b:Vec; N:Integer; Var X:Vec);
Var Alp, Bet : Vec;
           R : Real;
           i : Integer;
 Begin
  Alp[1] := - A[1, 2] / A[1, 1];
  Bet[1] := b[1] / A[1, 1];
  For i := 2 to N-1 do
   Begin
    R := A[i,i] + Alp[i - 1] * A[i, i - 1];
    Alp[i] := - A[i, i + 1] / R;
    Bet [i] := (b[i] - A[i, i - 1] * Bet[i - 1]) / R
   End;
   X[N] := (B[N] - A[N, N - 1] * Bet[N - 1]) / (A[N, N]
           + Alp[N - 1] * A[N, N - 1]);
   For i := N-1 downto 1 do
    X[i] := Alp [i] * X [i + 1] + Bet [i]
 End;
//progonka
Begin
  Ne := Ne - 1;
  For i := 1 to Ne do
   Begin
    H[i] := Xe[i] - Xe[i - 1];
    a[i] := Ye[i]
   End;
  For i := 1 to Ne - 1 do
   For j := 1 to Ne - 1 do
   If i=j then
    Begin
     R := H[i]; Q :=H[i + 1];
     Ma[i, i] := 2 * (R + Q);
     Y[i] := 6*((Ye[i+1]-Ye[i])/Q-(Ye[i]-Ye[i-1])/R);
    End   Else
    Begin
     Q := H[i + 1];
     Ma[j, i] := 0;
     Ma[i + 1, i] := Q;
     Ma[i, i + 1] := Q
    End;
Progonka (Ma,Y,Ne-1,c);
  c[0] := 0;  c[Ne] := 0;
  For i:=1 to Ne do
   Begin
    R := h[i];
    d[i] := (c[i] - c[i-1]) / R;
    b[i] := R*c[i]/2 - Sqr(R)*d[i]/6+(Ye[i] - Ye[i-1])/R
   End;
End;
  //cubspln
 Begin
  CubSpln (Xe, Ye, Ne, a, b, c, d);
  L := 1;
  h := (bl - al) / Ngr;
  Xg[0]:= al;
  For i := 0 to Ngr - 1 do
   Begin
    YG[i] := a[L] + b[L] * (Xg[i] - Xe[L])+ c[L] * Sqr(Xg[i] -
   Xe[L])/2 + d[L]*Sqr(Xg[i] - Xe[L])*(Xg[i] - Xe[L])/6;
    Xg[i + 1] := Xg[i] + h;
    If Xg[i] > Xe [L] then L := L + 1
   End
 END;
 
 
  //end spline

@Monty161 · 13.12.2015, 12:53

Bar-bar1an, лентяй

@HectorPrima · 13.12.2015, 12:58

Мне кажется что тут каша из процедур?
Что к чему относится не пойму.

@Bar-bar1an · 13.12.2015, 14:13 **[ТС]**

не лєнтяй) просто два дня не могу нормально переписать.
каша есть. там коментарии есть чтоб можно было разобраться

Добавлено через 6 минут
за кашу извеняйте)код не мой.

Konctantin · 13.12.2015, 14:52

о ужас, такое чувство что писал кто-то так, будто был пример полегче но докрутить функционал, в итоге получили эту "кашу". Переписать код не сложно, но вот только желание надо.

@Bar-bar1an · 13.12.2015, 23:56 **[ТС]**

ето мой код. помогите пожалуйста найти ошибку ( функция progonka заботает нормально)

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
 private void Progonka(Double[,] A, Double[] B, Int32 N, Double[] X)
        {
            Double[] Alp = new Double[N];
            Double[] Bet = new Double[N];
            Double R;
            Int32 i;
            Alp[0] = -A[0, 1] / A[0, 0];
            Bet[0] = B[0] / A[0, 0];
            for (i = 1; i < N; i++)
            {
                R = A[i, i] + Alp[i - 1] * A[i, i - 1];
                Alp[i] = -A[i, i + 1] / R;
                Bet[i] = (B[i] - A[i, i - 1] * Bet[i - 1]) / R;
            }
            X[N] = (B[N] - A[N, N - 1] * Bet[N - 1]) / (A[N, N] + Alp[N - 1] * A[N, N - 1]);
            for (i = N - 1; i >= 0; i--)
            {
                X[i] = Alp[i] * X[i + 1] + Bet[i];
            }
        }
        private void CubSpline(Double[] Xe, Double[] Ye, Int32 Ne, Double[] a, Double[] b, Double[] c, Double[] d) 
        {
            Double[,] Ma = new Double[50, 50];
            Double[] h = new Double[40];
            Double[] Y = new Double[40];
            Int32 i = 0, j = 0;
            Double R = 0, Q = 0;
            Ne--;
            for (i = 1; i <= Ne - 1; i++)
            {
                h[i] = Xe[i] = Xe[i - 1];
                a[i] = Ye[i];
            }
            for ( i = 1; i <= Ne - 1; i++)
            for ( j = 1; j <= Ne - 1; j++)
           {
                    if (i == j)
                    {
                        R = h[i];
                        Q = h[i + 1];
                        Ma[i, i] = 2 * (R + Q);
                        Y[i] = 6 * ((Ye[i + 1] - Ye[i]) / Q - (Ye[i] - Ye[i - 1]) / R);//Ye[i-1]???
                    }
                    else 
                    {
                        Q = h[i + 1];// 
                        Ma[j, i] = 0;
                        Ma[i + 1, i] = Q;
                        Ma[i, i + 1] = Q;
                    }
            }
                    Progonka(Ma, Y, Ne-1, c);
                    c[0] = 0;
                    c[Ne] = 0;
                 for (i = 1; i <= Ne; i++)
                   {
                        R = h[i];
                        d[i] = (c[i] - c[i-1]) / R;
                        b[i] = R * c[i] / 2 - R * R * d[i] / 6 + (Ye[i] - Ye[i - 1]) / R;
                   }
                
           }
        private void TabSpline(Double[] Xe, Double[] Ye, Int32 Ne, Int32 Ngr ,Double[]Xg, Double[]Yg) 
        {
            CubSpline(Xe, Ye, Ne, a, b, c, d);
            Int32 L = 1;
            Double h = (bl - al) / Ngr;
            Xg[0] = al;
            //Yg[0] = a[L] + b[L] * (Xg[0] - Xe[L]) + c[L] * Math.Pow(Xg[0] - Xe[L], 2) / 2 + d[L] * Math.Pow(Xg[0] - Xe[L], 2) * (Xg[0] - Xe[L]) / 6;
            for (int i = 0; i <= Ngr - 1 ; i++)
            {
                Yg[i] = a[L] + b[L] * (Xg[i] - Xe[L]) + c[L] * Math.Pow(Xg[i] - Xe[L], 2) / 2 + d[L] * Math.Pow(Xg[i] - Xe[L], 2) * (Xg[i] - Xe[L]) / 6;
                Xg[i + 1] = Xg[i] + h;
                if (Xg[i] > Xe[L]) L++;
            }
        }

Добавлено через 1 минуту
*работает

Добавлено через 8 часов 7 минут
уже не надто. сам сделал)

Новые блоги и статьи Все статьи Все блоги /
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 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .

@Monty161 6 / 6 / 6 Регистрация: 01.09.2015 Сообщений: 52
	13.12.2015, 12:53
	Bar-bar1an, лентяй 0

@HectorPrima 716 / 473 / 258 Регистрация: 19.12.2014 Сообщений: 1,801
	13.12.2015, 12:58
	Мне кажется что тут каша из процедур? Что к чему относится не пойму. 0

@Bar-bar1an 1 / 1 / 6 Регистрация: 12.12.2014 Сообщений: 50
	13.12.2015, 14:13 [ТС]
	не лєнтяй) просто два дня не могу нормально переписать. каша есть. там коментарии есть чтоб можно было разобраться Добавлено через 6 минут за кашу извеняйте)код не мой. 0

Konctantin 970 / 773 / 171 Регистрация: 12.04.2009 Сообщений: 1,700
	13.12.2015, 14:52
	о ужас, такое чувство что писал кто-то так, будто был пример полегче но докрутить функционал, в итоге получили эту "кашу". Переписать код не сложно, но вот только желание надо. 0