Расчет коэффициента полинома в матлаб

Большую часть времени у океанолога занимает общение с различными программами и написание небольших программ для обработки и визуализации данных. Обсуждению этого и посвящён данный раздел.

Модераторы: Magik, koldunov aleksey, kostek

Расчет коэффициента полинома в матлаб

Сообщение Саймон » Пт окт 13, 2017 12:19 pm

Нужно рассчитать коэффициент полинома в Матлаб.
Есть массив данных. Нужно по ним, сначала вычесть линейный тренд, (отбросив из данных 9999), а затем рассчитать коэффициент полинома второй степени. И если второй коэфициент <0, то ветви параболы направлены вверх, ищем минимум, в противном случае ищем максимум.
что то среднее между
d=importdata('TXT file (2).txt')';

for i=1:size(d,2)
D=d(:,i);
V=D(D~=9999);
v=V(5:end);
% plot(v,'*'); hold on;
fprintf('%5d %5d %5d %10.3f %10.3f\n',D(1), D(3), D(4), min(v), max(v))
end
grid on;
и.....

di=importdata('TXT file (2).txt')';

nt=100; % число точек для построения графика полинома
d=di(5:end,:);
for i=1:size(d,2)
D(:,1)=d(:,i);
D(:,2)=1:size(d,1);
V=D(D(:,1)~=9999,:);

p3=polyfit(V(:,2),V(:,1),2); % 2 - степень полинома
x=V(1,2):(V(end,2)-V(1,2))/nt:V(end,2);
yP3=polyval(p3,x);
figure
plot(V(:,2),V(:,1),'r*'); hold on;
plot(x,yP3,'b-');
grid on;
end



di=importdata('F:\argo\DATA_GRID\.txt')';

nt=100; % число точек для построения графика полинома
d=di(5:end,:);
for i=1:size(d,2)
D(:,1)=d(:,i);
D(:,2)=1:size(d,1);
V=D(D(:,1)~=9999,:);

p3=polyfit(V(:,2),V(:,1),2); % 2 - степень полинома
x=V(1,2):(V(end,2)-V(1,2))/nt:V(end,2);
yP3=polyval(p3,x);
figure
plot(V(:,2),V(:,1),'r*'); hold on;
plot(x,yP3,'b-');
grid on;
end
Саймон
Гость
 
Сообщения: 8
Зарегистрирован: Пн сен 25, 2017 1:47 pm

Re: Расчет коэффициента полинома в матлаб

Сообщение Саймон » Пт окт 13, 2017 12:21 pm

данные текстовые
Вложения
111(3).xls
(27.5 Кб) Скачиваний: 16
Саймон
Гость
 
Сообщения: 8
Зарегистрирован: Пн сен 25, 2017 1:47 pm

Re: Расчет коэффициента полинома в матлаб

Сообщение Саймон » Пн окт 23, 2017 10:57 am

Есть соображения как доделать код?
Саймон
Гость
 
Сообщения: 8
Зарегистрирован: Пн сен 25, 2017 1:47 pm


Вернуться в Программы и программирование

cron