Погрешность аппроксимации

@Bryan1 · Регистрация: 23.11.2012

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

Помогите разобраться!

У меня задачи аппроксимации. Аппроксимировав функцию, нужно найти погрешность аппроксимации по методу наим квадратов.

у меня вот на этом шаге ошибка. Помогите разобраться!

Matlab M
1
2
3
4
5
6
7
8
9
10
11
clc;
clear all;
x=[22.0 24.0 27.0 30.0 31.0 35.0 40.0];
y=[1.24 1.37 1.46 1.26 1.66 1.84 1.99];
% Найдем коэффициэнты полинома
a=polyfit(x,y,5),
%  Найдем значения полиномов в интервале от 20 до 40
x1=20:0.5:40;
y1=polyval(a,x1);
% Погрешность аппроксимации по МНК
e=((y1(x1)-y(x))^2)/7,

Зосима · 03.12.2012, 11:00

Надо было сделать так:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
clc;
clear all;
x=[22.0 24.0 27.0 30.0 31.0 35.0 40.0];
y=[1.24 1.37 1.46 1.26 1.66 1.84 1.99];
% Найдем коэффициэнты полинома
a=polyfit(x,y,5),
%  Найдем значения полиномов в интервале от 20 до 40
x1=20:0.5:40;
y1=polyval(a,x1);
y0 = polyval(a,x); % расчет ф-ции в точках данных
% Погрешность аппроксимации по МНК
e=sum((y0-y).^2)/length(y)

@Vramo · 30.04.2016, 15:30

Прошу прощения, если например мне нужно показать на графике, локальные ошибки МНК(полинома 3-го порядка) функцией errorbar. То для вычисленения значений, локальных ошибок, мне необходимо делать примерно так ?

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
N = 3;
coeff1 = polyfit(x, y, N);%MNK
 
Y = polyval(coeff1,x);%значение полинома
 
%hold on; plot(x,Y,'b');
 
% Погрешность аппроксимации по МНК
for i=0:1:length(y)
    e = ((Y-y).^2)/length(y)
end
 
hold on;errorbar(x,Y,e,'r');

Зарание спасибо.

@Centurio · 01.05.2016, 08:48

Vramo, можно, только зачем вычисление е в цикл вводить? У вас каждую итерацию будет вычисляться один и тот же вектор. Выведите вычисление из цикла, цикл уберите

Новые блоги и статьи Все статьи Все блоги /
Почему дизайн решает? Neotwalker 09.01.2026 В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .	Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .
Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .	Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .

@Centurio Модератор 1765 / 1610 / 541 Регистрация: 13.09.2015 Сообщений: 5,629
	01.05.2016, 08:48
	Vramo, можно, только зачем вычисление е в цикл вводить? У вас каждую итерацию будет вычисляться один и тот же вектор. Выведите вычисление из цикла, цикл уберите 1