Аппроксимация в matlab

@Tetya_Motya · Регистрация: 24.10.2010

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

Здравствуйте! Не могли бы помочь, нужно выполнить аппроксимацию в matlab(simulink)!
Есть кривая разгона (график) и функция описывающая его без коэффициентов, нужно найти эти коэффициенты по минимальному отклонению от данного графика.

@Sledge5 · 23.04.2011, 17:45

Смотрите функцию polyfit

@Tetya_Motya · 23.04.2011, 21:37 **[ТС]**

Я про эту функцию читала, только я такой профан в матлабе, что ничего понять не могу((
Может у кого-нибудь есть готовый пример, поделитесь плз!

@Sledge5 · 23.04.2011, 22:12

Легко!

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
clear x y p0 p1 p2, close, clc 
x = 1:10
y = sort(randi([0 100], 1, 10))
p0 = polyfit(x,y,0)
p1 = polyfit(x,y,1)
p2 = polyfit(x,y,2)
figure('Units','normalized','OuterPosition',[0 0 1 1])
plot(x, y, 'rx')
hold on
fplot(@(x)polyval(p0,x), [1 10], 'b')
fplot(@(x)polyval(p1,x), [1 10], 'k')
fplot(@(x)polyval(p2,x), [1 10], 'g')
legend('табличные значения', 'полином 0-й степени', 'полином 1-й степени', 'полином 2-й степени',-1)

@Tetya_Motya · 23.04.2011, 22:30 **[ТС]**

Тут как я поняла он находит коэффициент при х^3,х^2 и т.д.
А если у меня функция: k1/(k2*x+1)^3-(1/k3*x), как быть тогда?

@Sledge5 · 24.04.2011, 00:12

Тогда смотреть документацию по Curve Fitting Toolbox. Это дополнение к MATLABу, кажется, может такое.

@Tetya_Motya · 29.04.2011, 15:17 **[ТС]**

Спасибо, она подходит вроде все сделала правильно, но почему-то график получается прямая и выводит ошибку:

Matlab M
1
2
3
4
5
6
7
8
9
time = [0 6.75 12.5 25 50 75 100 150];
level = [0 0.75 2 3 3.6 3.3 2.5 1];
plot(time,level,'o');
xlabel('Time'); ylabel('Level');
>> modelFun = @(p,x) (p(1) ./(p(2) .* x-1).^3-1 ./ p(3) .*x);
startingVals = [1 2 40];
coefEsts = nlinfit(time, level, modelFun, startingVals);
xgrid = linspace(0,150,1000);
line(xgrid, modelFun(coefEsts, xgrid), 'Color','r');

Ошибка:
Rank deficient, rank = 2, tol = 2.4437e-015.
> In nlinfit>LMfit at 312
In nlinfit at 166
Warning: Rank deficient, rank = 2, tol = 2.4426e-015.
> In nlinfit>LMfit at 312
In nlinfit at 166
Warning: Rank deficient, rank = 2, tol = 2.4425e-015.
> In nlinfit>LMfit at 312
In nlinfit at 166
Warning: Rank deficient, rank = 2, tol = 2.4425e-015.
> In nlinfit>LMfit at 312
In nlinfit at 166
Warning: Rank deficient, rank = 1, tol = 2.4425e-015.
> In nlinfit>LMfit at 312
In nlinfit at 166
Warning: Some columns of the Jacobian are effectively zero at the solution, indicating that the model is insensitive to some of its
parameters. That may be because those parameters are not present in the model, or otherwise do not affect the predicted values. It
may also be due to numerical underflow in the model function, which can sometimes be avoided by choosing better initial parameter
values, or by rescaling or recentering. Parameter estimates may be unreliable.
> In nlinfit at 231

Не подскажите , что не правильно?

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Работа со звуком через SDL3_mixer 8Observer8 08.02.2026 Содержание блога Пошагово создадим проект для загрузки звукового файла и воспроизведения звука с помощью библиотеки SDL3_mixer. Звук будет воспроизводиться по клику мышки по холсту на Desktop и по. . .	SDL3 для Web (WebAssembly): Основы отладки веб-приложений на SDL3 по USB и Wi-Fi, запущенных в браузере мобильных устройств 8Observer8 07.02.2026 Содержание блога Браузер Chrome имеет средства для отладки мобильных веб-приложений по USB. В этой пошаговой инструкции ограничимся работой с консолью. Вывод в консоль - это часть процесса. . .	SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .
SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .	SDL3 для Web (WebAssembly): Сборка C/C++ проекта из консоли 8Observer8 30.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .	SDL3 для Web (WebAssembly): Установка Emscripten SDK (emsdk) и CMake для сборки C и C++ приложений в Wasm 8Observer8 30.01.2026 Содержание блога Для того чтобы скачать Emscripten SDK (emsdk) необходимо сначало скачать и уставить Git: Install for Windows. Следуйте стандартной процедуре установки Git через установщик. . . .	SDL3 для Android: Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 29.01.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами. Версия v3 была полностью переписана на Си, в. . .

@Tetya_Motya 0 / 0 / 0 Регистрация: 24.10.2010 Сообщений: 9

	Аппроксимация в matlab 23.04.2011, 17:25. Показов 21920. Ответов 6 Метки нет (Все метки) Здравствуйте! Не могли бы помочь, нужно выполнить аппроксимацию в matlab(simulink)! Есть кривая разгона (график) и функция описывающая его без коэффициентов, нужно найти эти коэффициенты по минимальному отклонению от данного графика. 0

@Sledge5 536 / 523 / 38 Регистрация: 13.03.2011 Сообщений: 728
	23.04.2011, 17:45
	Смотрите функцию `polyfit` 0

@Tetya_Motya 0 / 0 / 0 Регистрация: 24.10.2010 Сообщений: 9
	23.04.2011, 21:37 [ТС]
	Я про эту функцию читала, только я такой профан в матлабе, что ничего понять не могу(( Может у кого-нибудь есть готовый пример, поделитесь плз! 0

@Tetya_Motya 0 / 0 / 0 Регистрация: 24.10.2010 Сообщений: 9
	23.04.2011, 22:30 [ТС]
	Тут как я поняла он находит коэффициент при х^3,х^2 и т.д. А если у меня функция: k1/(k2x+1)^3-(1/k3x), как быть тогда? 0

@Sledge5 536 / 523 / 38 Регистрация: 13.03.2011 Сообщений: 728
	24.04.2011, 00:12
	Тогда смотреть документацию по Curve Fitting Toolbox. Это дополнение к MATLABу, кажется, может такое. 1