Типовые задачи VBA

Помогите пожалуйста решить задачи VBA
Вот текст задач... функция второй задачи в приложенной картинке + текст задач в приложенном документе doc.
Задача 1. Построить линейную регрессионную модель вида
у = aо + a1*х
формализации процесса функционирования некоторой системы, представленного следующей таблицей экспериментальных данных
№ опыта/ Значение x/ Значение y
1/ 0,2/ 2,92
2/ 0,4/ 3,25
3/ 0,6/ 3,74
4/ 0.8/ 4,15
5/ 1,0/ 4,45
6/ 1,2/ 4,88
7/ 1,4/ 5,33
8/ 1,6/ 5,72
используя метод наименьших квадратов

Задача 2. Разработать алгоритм исследования системы, описывающейся следующей системой обыкновенных дифференциальных уравнений

см. рис 1


где А1, А2, B1, B2 заданные постоянные коэффициенты;
x(t), y(t), – искомые функции времени;
cx, cy , kx, ky – нелинейные коэффициенты, зависящие от x и y;
qx, – воздействие внешней среды на систему и ее скорость;
при начальных условиях x(tо)=0,05 м, y(tо)= 0,025 м
на интервале времени t = 0…20 с.

Буду очень благодарен. В силу моих знаний vba решение не представляется возможным.))) Предполагаю, что для знающего человека эти задачи пустяк.

Миниатюры

 

Вложения

Задачи ОИТ.docx (19.4 Кб, 16 просмотров)

0

Сообщение от winecs Задача 1.

Реализовать программу для апроксимации методом найменьших квадратов данной таблицы (только тренд там не не y, а z, он такой цвета хаки — увидите на графике)

Разумеется, что, открыв код*, вы подставите туда ряды x и y из своей задачи.
________________
* Alt-F8 → Макрос7 → Войти (либо Alt-F11)

1

Огромное спасибо!!!
Однако, по причине отсутствия решения задачи 2... да и совершенного не представления его пути... Взял другую задачу... Аналогичную первой, но функция нелинейная: у = а*х^b, вместо у = aо + a1*х,
Поскажите, что изменить в коде

Option Explicit
Option Base 1

Sub Макрос7() ' Строит прямой тренд по методу наименьших квадратов.
' Сочетание клавиш: Ctrl+й
Dim x, y 'массивы значений для иксов и игреков
Dim a, b, D 'искомые коэффициенты уравнения y = ax + b и детерминант матрицы
Dim Da, Db 'детерминанты при замене столбцов соответственно при неизвестных a и b
Dim СуммаИксов, СуммаИгрек, СуммаКвадратовИкс, СуммаПроизведений, ЧислоТочек
Dim z(), i 'переменные для вывода аппроксимирующего массива точек

x = Array(0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4, 1.6)
y = Array(2.92, 3.25, 3.74, 4.15, 4.45, 4.88, 5.33, 5.72)

If UBound(y) < UBound(x) Then MsgBox "В массиве y мало элементов.": Exit Sub
ЧислоТочек = UBound(x)
ReDim z(1 To ЧислоТочек)

With Application.WorksheetFunction
СуммаИксов = .Sum(x): СуммаИгрек = .Sum(y)
СуммаКвадратовИкс = .SumSq(x): СуммаПроизведений = .SumProduct(x, y)

'СуммаКвадратовИкс * a + СуммаИксов * b = СуммаПроизведений
' СуммаИксов * a + ЧислоТочек * b = СуммаИгрек

D = .MDeterm(Array(Array(СуммаКвадратовИкс, СуммаИксов), _
Array(СуммаИксов, ЧислоТочек)))
Da = .MDeterm(Array(Array(СуммаПроизведений, СуммаИксов), _
Array(СуммаИгрек, ЧислоТочек)))
Db = .MDeterm(Array(Array(СуммаКвадратовИкс, СуммаПроизведений), _
Array(СуммаИксов, СуммаИгрек)))
a = Da / D
b = Db / D
End With

For i = 1 To UBound(x)
z(i) = CDbl(FormatNumber(a * x(i) + b, 3)) 'заодно форматируем
Next
'Теперь проведём прямую y = ax + b (аналитическую) и исходную (эмпирическую) кривую.
Cells.Clear 'ячейки чистим
[r1c1] = "x": [r2c1] = "y": [r3c1] = "z" 'здесь первые 3 строки получают имена

[r1c2].Resize(1, ЧислоТочек) = x 'здесь 1-я строка заполняется из массива x
[r2c2].Resize(1, ЧислоТочек) = y 'здесь 2-я строка заполняется из массива y
[r3c2].Resize(1, ЧислоТочек) = z 'здесь 3-я строка заполняется из массива z

ActiveSheet.Shapes.AddChart.Select 'добавили окно диаграмм
ActiveChart.SetSourceData Source:=Range("$1:$3") 'выделили первые 3 строки
ActiveChart.ChartType = xlXYScatterLinesNoMarkers 'провели линию, без маркеров
ActiveChart.SeriesCollection(1).Select: Selection.Delete 'удалили ряд 1 (x)
End Sub

0

А классно вы подставили свои данные (думал, не разберётесь)!

Сообщение от winecs Поскажите, что изменить в коде

Линейную зависимость на экспоненциальную. (Но как — ума не приложу.)

Сообщение от winecs Предполагаю, что для знающего человека эти задачи пустяк.

Вот если б это было так, я б не сидел всю ночь над предыдущей (линейной)! Это вам не дважды два.

0