Аппроксимация функцией Y=b+(a/X)

@Анирам2 · Регистрация: 07.01.2018

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

Добрый день! Подскажите где ошибка?

Аппроксимация функцией Y=b+(a/X)

Аппроксимация заключается в отыскании коэффициентов a и b уравнения таких, чтобы все экспериментальные точки лежали наиболее близко к аппроксимирующей прямой.

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
x=[0.659 1.364 2.070 2.775 3.388 4.001 4.660 5.366 5.979 6.638];
y=[3.371 1.062 0.321 -0.45 -0.22 -0.38 -0.45 -0.53 -0.64 -0.69];
u=y; z=x;
cd=polyfit(z,u,1);
a=cd(1); b=exp(cd(2));
yt=b.+x.^-a;
e=y-yt;
disp('x y yt e');
disp('-------------------------');
disp([x' y' yt' e']);
disp('-------------------------');
B=[z;u]'; D=corrcoef(B); R=D(1,2);
disp([a b R]);
xz=3.65; yz=b.+xz.^-a;
disp('xz yz'); disp([xz yz]);
plot(x,y,'ko',x,yt,'b',xz,yz,'R*'); grid on;
axis([0 7 -1 3.5]);

И как записать Y=1/(b+aX)

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

@ida2 · 07.01.2018, 14:36

Функция polyfit аппроксимирует МНК, полиномом положительной степени.

Функция lsqcurvefit аппроксимирует МНК, функцией заданной пользователем:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
x=[0.659 1.364 2.070 2.775 3.388 4.001 4.660 5.366 5.979 6.638];
y=[3.371 1.062 0.321 -0.45 -0.22 -0.38 -0.45 -0.53 -0.64 -0.69];
 
fun = @(a,xdata)a(1)+a(2)./xdata;
a = lsqcurvefit(fun,[100,-1],x,y);
yt = fun(a,x);
 
e=y-yt;
disp('x y yt e');
disp('-------------------------');
disp([x' y' yt' e']);
disp('-------------------------');
B=[z;u]'; D=corrcoef(B); R=D(1,2);
disp([a b R]);
xz=3.65; yz=fun(a,xz);
disp('xz yz'); disp([xz yz]);
plot(x,y,'ko',x,yt,'b',xz,yz,'R*'); grid on;
axis([0 7 -1 3.5]);

@Nick07 · 07.01.2018, 23:56

Cm. CFTOOL

Новые блоги и статьи Все статьи Все блоги /
Автозаполнение реквизита при выборе элемента справочника 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. Перед реализацией необходимо выполнить настройку системной учетной. . .	модель ЗдравоСохранения 5. Меньше увольнений- больше дохода! anaschu 24.03.2026 Теперь система здравосохранения уменьшает количество увольнений. 9TO2GP2bpX4 a42b81fb172ffc12ca589c7898261ccb/ https:/ / rutube. ru/ video/ a42b81fb172ffc12ca589c7898261ccb/ Слева синяя линия -. . .	Midnight Chicago Blues kumehtar 24.03.2026 Такой Midnight Chicago Blues, знаешь?. . Когда вечерние улицы становятся ночными, а ты не можешь уснуть. Ты идёшь в любимый старый бар, и бармен наливает тебе виски. Ты смотришь на пролетающие. . .

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	07.01.2018, 23:56
	Cm. CFTOOL 0