Scilab. Методом наименьших квадратов подобрать зависимость

@serega_ivanov · Регистрация: 15.03.2020

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

Методом наименьших квадратов подобрать зависимость. Для оценки корректности вида зависимости пользуйтесь индексом корреляции. Ошибка получилась. Помогите исправить, пожалуйста.

Данные:
G(s)=As^b

x=[0.5 1.5 2 2.5 3 3.5 4 4.5 5];

y=[3.99 5.65 6.41 6.71 7.215 7.611 7.83 8.19 8.3];

Вот как я начал:

Matlab M

function [ind]=IndexCorrel(y,Y)
  ind =sqrt(1-sum((y - Y)^2)/sum((y-mean(y))^2));
endfunction
 
function [zr]=G2(c,z)
zr=(2)-c(1)*z(1).^c(2)
endfunction
function [zr]=G3(c,z)
zr=(2)-c(1)*z(1).^c(2)
endfunction
 
x=1 :1:0;
y=[-1  5 -27 -163 -445 -891 -1495 -2227-3033 -3835  -4531 -4995 -5077 -4603 -3375 -1171 2255 7173];
 
z=[x;y];
 
c=zeros(3,1)
[a,err]=datafit(G2,z,c);
Ind=IndexCorrel(y, (a(1)*x.^a(2))
 
plot2d(x,y,-2);
t=1:1:0;
Ptc=a(1)*t.^a(2);
plot2d(t,Ptc);
 
c=zeros(4,1)
[a,err]=datafit(G3,z,c);
Ind(2)=IndexCorrel(y, (a(1)*x.^a(2))
 
plot2d(x,y,2);
t=1:1:22;
Ptc=a(1)*t.^a(2);
plot2d(t,Ptc);

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

Правила форума, пункт 4.9. Используйте тэги форматирования текста и редактор формул для удобства восприятия ваших сообщений другими пользователями.

@Krasme · 23.12.2020, 23:05

g2,g3 одинаковы, какой смысл во второй раз писать?
внутри функции пропущен символ, смотрите внимательно
вектор х неверно задан

@serega_ivanov · 23.12.2020, 23:43 **[ТС]**

Krasme, Исправил. Все получилось. Я сомневаюсь нужно ли искать вектор начальных приближений кубическая функция?

Matlab M

function [ind]=IndexCorrel(y,Y)
  ind =sqrt(1-sum((y - Y)^2)/sum((y-mean(y))^2));
   endfunction
 
  function [zr]=G2(c,z)
  zr= z(2)-c(1)*z(1).^c(2)
  endfunction
 
 x=1 :1:18;
 y=[-1. 5. -27. -163. -445. -891. -1495. -2227. -3033. -3835. -4531. -4995. -5077. -4603. -3375. -1171. 2255. 7173.];
z=[x;y];
c=zeros(3,1)
 c  = 
   0.
   0.
   0.
 [a,err]=datafit(G2,z,c);
Ind=IndexCorrel(y, (a(1)*x.^a(2)))
 Ind  = 
   0.1070840
 plot2d(x,y,-4);
t=1:1:18;
 Ptc= a(1)*t.^a(2);
plot2d(t,Ptc);

@Krasme · 24.12.2020, 00:30

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

подобрать зависимость.

вид зависимости самостоятельно надо подбирать.. вот и подумайте, что туда лучше подойдет..

Добавлено через 36 минут
вот приближенная функция

Matlab M

1
2
3

function [zr]=G2(c,z)
  zr=z(2)-(c(1)*z(1)^3+c(2)*z(1)^2+c(3)*z(1)+c(4))
endfunction

по какой формуле вы ищете корреляцию? источник покажите..

@serega_ivanov · 24.12.2020, 00:31 **[ТС]**

Плохо в этом соображаю. Если верно, то по этой

function [ind]=IndexCorrel(y,Y)
ind =sqrt(1-sum((y - Y)^2)/sum((y-mean(y))^2));
endfunction

@Krasme · 24.12.2020, 00:34

код я и так вижу, мне нужно, на основе чего вы эту формулу прописали..

@serega_ivanov · 24.12.2020, 00:42 **[ТС]**

G(s)=As^b

Добавлено через 3 минуты
Krasme, G(s)=As^b

@Krasme · 24.12.2020, 00:54

это функция аппроксимации из прошлой задачи.

@serega_ivanov · 24.12.2020, 00:55 **[ТС]**

Krasme, это все в одном задании. Вот именно по этой функции и в 1 задании и 2 тоже. А так больше никакая не прописывается(((

@Krasme · 24.12.2020, 01:36

ну, по заданию вам виднее..

а по аппроксимации, второй набор данных аппроксимируется по кубической функции, вид ее я уже дала (пост №4)

Добавлено через 25 минут
полный код по второму набору данных

Кликните здесь для просмотра всего текста

Matlab M

clc; clear; clf;
function [ind]=IndexCorrel(y,Y)
  ind =sqrt(1-sum((y - Y)^2)/sum((y-mean(y))^2));
endfunction
 
function y = f(x,c)    
    y = c(1)*x^3+c(2)*x^2+c(3)*x+c(4)
endfunction
 
function [zr]=G2(c,z)
  zr=z(2)-(c(1)*z(1)^3+c(2)*z(1)^2+c(3)*z(1)+c(4))
endfunction
 
x=1:18;
y=[-1 5 -27 -163 -445 -891 -1495 -2227 -3033 -3835 -4531 -4995 -5077 -4603 -3375 -1171 2255 7173];
 
z=[x;y];
c=[1; 1; 1; 1]
[a,err]=datafit(G2,z,c);
Ind=IndexCorrel(y, f(x,a))
disp(Ind)
 
plot2d(x,y,-2);
t=1:0.1:18;
Ptc=f(t,a)
plot2d(t,Ptc,5);

@serega_ivanov · 24.12.2020, 10:49 **[ТС]**

Krasme, Да, я все переделал. Получилось. Спасибо огромное за помощь.

Новые блоги и статьи
Это работает. Скорость асинхронной логики велика. Вопрос видимо останется в стабильности. Плата - огонь! Hrethgir 13.01.2025 По прошлому проекту в Logisim Evolution https:/ / www. cyberforum. ru/ blogs/ 223907/ blog8781. html прилагаю файл архива проекта Gowin Eda и снимок. Восьмибитный счётчик из сумматора+ генератор сигнала. . .	UserScript для подсветки кнопок языков программирования в зависимости от текущего раздела volvo 13.01.2025 В результате работы этого скрипта подсвечиваются нужные кнопки не только в форме быстрого ответа, но и при редактировании сообщения: / / ==UserScript== / / @name CF_DefaultLangSelect / / . . .	Введение в модели и алгоритмы машинного обучения InfoMaster 12.01.2025 Машинное обучение представляет собой одну из наиболее динамично развивающихся областей искусственного интеллекта, которая фокусируется на разработке алгоритмов и методов, позволяющих компьютерам. . .	Как на Python создать нейросеть для решения задач InfoMaster 12.01.2025 В контексте стремительного развития современных технологий особое внимание уделяется таким инструментам, как нейросети. Эти структуры, вдохновленные биологическими нейронными сетями, используются для. . .	Как создать нейросеть для генерации картинок на Python InfoMaster 12.01.2025 Генерация изображений с помощью искусственных нейронных сетей стала одним из наиболее захватывающих направлений в области компьютерного зрения и машинного обучения. В этой статье мы рассмотрим. . .	Создание нейросети для генерации текста на Python InfoMaster 12.01.2025 Нейросети, или искусственные нейронные сети, представляют собой модели машинного обучения, вдохновленные работой человеческого мозга. Они состоят из множества взаимосвязанных узлов, или "нейронов",. . .
Как создать нейросеть распознавания изображений на Python InfoMaster 12.01.2025 Введение в распознавание изображений с помощью нейросетей Распознавание изображений с помощью нейронных сетей стало одним из самых впечатляющих достижений в области искусственного интеллекта. Эта. . .	Основы искуственного интеллекта InfoMaster 12.01.2025 Искусственный интеллект (ИИ) представляет собой одну из наиболее динамично развивающихся областей современной науки и технологий. В широком смысле под искусственным интеллектом понимается способность. . .	Python и нейросети InfoMaster 12.01.2025 Искусственные нейронные сети стали неотъемлемой частью современных технологий, революционизировав множество областей - от медицинской диагностики до автономных транспортных средств. Python, благодаря. . .	Python в машинном обучении InfoMaster 12.01.2025 Python стал неотъемлемой частью современного машинного обучения, завоевав позицию ведущего языка программирования в этой области. Его популярность обусловлена несколькими ключевыми факторами, которые. . .	Создание UI на Python с TKinter InfoMaster 12.01.2025 TKinter — это одна из наиболее популярных библиотек для создания графических интерфейсов пользователей (GUI) в языке программирования Python. TKinter входит в стандартную библиотеку Python, что. . .	HTML5 в разработке мобильных приложений InfoMaster 12.01.2025 Введение: Обзор роли HTML5 в мобильной разработке В современном мире мобильных технологий HTML5 стал ключевым инструментом для разработки кроссплатформенных приложений. Эта технология произвела. . .

@Krasme 6961 / 4990 / 2108 Регистрация: 02.02.2014 Сообщений: 13,249
	23.12.2020, 23:05	2
	g2,g3 одинаковы, какой смысл во второй раз писать? внутри функции пропущен символ, смотрите внимательно вектор х неверно задан 0

@serega_ivanov 0 / 0 / 0 Регистрация: 15.03.2020 Сообщений: 26
	24.12.2020, 00:31 [ТС]	5
	Плохо в этом соображаю. Если верно, то по этой function [ind]=IndexCorrel(y,Y) ind =sqrt(1-sum((y - Y)^2)/sum((y-mean(y))^2)); endfunction 0

@Krasme 6961 / 4990 / 2108 Регистрация: 02.02.2014 Сообщений: 13,249
	24.12.2020, 00:34	6
	код я и так вижу, мне нужно, на основе чего вы эту формулу прописали.. 0

@serega_ivanov 0 / 0 / 0 Регистрация: 15.03.2020 Сообщений: 26
	24.12.2020, 00:42 [ТС]	7
	G(s)=As^b Добавлено через 3 минуты Krasme, G(s)=As^b 0

@Krasme 6961 / 4990 / 2108 Регистрация: 02.02.2014 Сообщений: 13,249
	24.12.2020, 00:54	8
	это функция аппроксимации из прошлой задачи. 0

@serega_ivanov 0 / 0 / 0 Регистрация: 15.03.2020 Сообщений: 26
	24.12.2020, 00:55 [ТС]	9
	Krasme, это все в одном задании. Вот именно по этой функции и в 1 задании и 2 тоже. А так больше никакая не прописывается((( 0

@serega_ivanov 0 / 0 / 0 Регистрация: 15.03.2020 Сообщений: 26
	24.12.2020, 10:49 [ТС]	11
	Krasme, Да, я все переделал. Получилось. Спасибо огромное за помощь. 0

Scilab. Методом наименьших квадратов подобрать зависимость

Решение