Аппроксимация. Метод наименьших квадратов

@Klau · Регистрация: 04.02.2022

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

Добрый вечер! НЕобходимо выполнить аппроксимацию точек, методом наименьших квадратов. Не хочет запускать команду, помогите найти ошибку

Visual Basic
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
Sub lb3()
Dim mx(5) As Variant
Dim my(5) As Variant
Mas2 = Array(14, 6, 2, 20, 3, 27, 15, 8, 10, 99)
For i = LBound(Mas2) To UBound(Mas2)
If i Mod 2 = 0 Then
mx(j) = Mas2(i)
j = j + 1
Else
my(jj) = Mas2(i)
jj = jj + 1
End If
Next i
'Debug.Print "x: " & Join(mx, vbTab)
'Debug.Print "y: " & Join(my, vbTab)
Call lab4(mx, my, a, b)
MsgBox "y=" & Round(a, 3) & "*x+" & Round(b, 3)
End Sub
 
Sub lab4(x As Variant, y As Variant, a As Double, b As Double)
n% = UBound(x, 1)
Sx# = 0
Sy# = 0
Sxy# = 0
Sx2# = 0
SY2# = 0
For i% = 1 To n%
    Sx# = Sx# + x(i%)
    Sy# = Sy# + x(i%)
    Sxy# = Sxy# + x(i%) * y(i%)
    Sx2# = Sx2# + x(i%) ^ 2
    SY2# = SY2# + y(i%) ^ 2
Next i%
    D# = CDb1(n%) * Sx2# - Sx# ^ 2
    Da# = CDb1(n%) * Sxy# - Sx# * Sy#
    Db# = Sy# * Sx2# - Sx# * Sxy#
    a = Da# / D#
    b = Db# / D#
End Sub

I can · 24.02.2022, 17:59

Сообщение от Klau

помогите найти ошибку

Где код взял? Копипастил или переписывал?

Добавлено через 3 минуты
1. Переменные a и b надо объявить.
2. Исправить функцию CDb1 на CDbl

@Klau · 24.02.2022, 18:04 **[ТС]**

Сообщение от I can

Где код взял? Копипастил или переписывал?

переписывал

I can · 24.02.2022, 18:06

Visual Basic
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
Sub lb3()
Dim mx(5) As Variant, a As Double, b As Double
Dim my(5) As Variant
Mas2 = Array(14, 6, 2, 20, 3, 27, 15, 8, 10, 99)
For i = LBound(Mas2) To UBound(Mas2)
If i Mod 2 = 0 Then
mx(j) = Mas2(i)
j = j + 1
Else
my(jj) = Mas2(i)
jj = jj + 1
End If
Next i
'Debug.Print "x: " & Join(mx, vbTab)
'Debug.Print "y: " & Join(my, vbTab)
Call lab4(mx, my, a, b)
MsgBox "y=" & Round(a, 3) & "*x+" & Round(b, 3)
End Sub
 
Sub lab4(x As Variant, y As Variant, a As Double, b As Double)
n% = UBound(x, 1)
Sx# = 0
Sy# = 0
Sxy# = 0
Sx2# = 0
SY2# = 0
For i% = 1 To n%
    Sx# = Sx# + x(i%)
    Sy# = Sy# + x(i%)
    Sxy# = Sxy# + x(i%) * y(i%)
    Sx2# = Sx2# + x(i%) ^ 2
    SY2# = SY2# + y(i%) ^ 2
Next i%
    D# = CDbl(n%) * Sx2# - Sx# ^ 2
    Da# = CDbl(n%) * Sxy# - Sx# * Sy#
    Db# = Sy# * Sx2# - Sx# * Sxy#
    a = Da# / D#
    b = Db# / D#
End Sub

@Klau · 24.02.2022, 18:10 **[ТС]**

Спасибо! помогло)

@Argus19 · 26.02.2022, 04:34

Сообщение от Klau

НЕобходимо выполнить аппроксимацию точек, методом наименьших квадратов.

Что именно получается в результате работы программы?

I can · 26.02.2022, 05:39

Сообщение от Argus19

Что именно получается в результате работы программы?

Код посмотри, и поймёшь

@Argus19 · 26.02.2022, 06:14

Сообщение от I can

Код посмотри, и поймёшь

Мне нужен метод МНК и способ решения, полученных в результате его работы, системы уравнений, для получения, в свою очередь, коэффициентов полинома нужной степени. Для чего я и задал вопрос.
Уважаемый Catstail пытался помочь, но пока не получилось.

@Catstail · 26.02.2022, 07:52

Argus19, да, тема забылась... Но она подвешена.

@Argus19 · 26.02.2022, 10:05

Сообщение от Catstail

Но она подвешена.

Спасибо.

Новые блоги и статьи Все статьи Все блоги /
Автоматическое создание документа при проведении другого документа Maks 29.03.2026 Реализация из решения ниже выполнена на нетиповых документах, разработанных в конфигурации КА2. Есть нетиповой документ "ЗаявкаНаРемонтСпецтехники" и нетиповой документ "ПланированиеСпецтехники". В. . .	Настройка движения справочника по регистру сведений Maks 29.03.2026 Решение ниже реализовано на примере нетипового справочника "ТарифыМобильнойСвязи" разработанного в конфигурации КА2, с целью учета корпоративной мобильной связи в коммерческом предприятии. . . .	Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .
Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .	Команды формы и диалоговое окно Maks 26.03.2026 1. Команда формы "ЗаполнитьЗапчасти". Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. В качестве источника данных. . .	Кому нужен AOT? DevAlt 26.03.2026 Решил сделать простой ланчер Написал заготовку: dotnet new console --aot -o UrlHandler var items = args. Split(":"); var tag = items; var id = items; var executable = args;. . .	Отправка уведомления на почту при создании или изменении элементов справочника Maks 24.03.2026 Программная отправка письма электронной почты на примере типового справочника "Склады" в конфигурации БП3. Перед реализацией необходимо выполнить настройку системной учетной записи электронной. . .

@Klau 0 / 0 / 0 Регистрация: 04.02.2022 Сообщений: 6
	24.02.2022, 18:10 [ТС]
	Спасибо! помогло) 0

@Catstail Супер-модератор 38195 / 21128 / 4309 Регистрация: 12.02.2012 Сообщений: 34,733 Записей в блоге: 14
	26.02.2022, 07:52
	Argus19, да, тема забылась... Но она подвешена. 1

Опции темы