Производная (разност. ур-ие и МНК)

@incognito · Регистрация: 08.04.2010

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

Здравствуйте!

Не могу понять, почему результаты 2х методов аналогичны? Ошибку тоже не могу найти...
ребят, посморите, пожалуйста, ..

Matlab M
1
2
3
4
5
6
7
8
 
clear all
clc
X=[1 2 3 4 5];
Y=[4 34 2 6 8];
N=3;
fd=1;
Tdiscr=1;

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
 
function diff_poly=Interpolation_by_eqpolyfit(N,X,Y, fd)
    for  i=1:(N+1);
        for j=1:Lx
B(i,j)=X(j)^(N+1-i);
        end
    end
F=B';
A=inv(F'*F)*F'*Y'%Calculation of matrix A of coefficients polynomials
poly=A';
p_diff=polyder(poly);
poly_diff_curr_X=polyval(p_diff, X);
diff_poly=poly_diff_curr_X(3);

Matlab M
1
2
3
 
function diff_eq=difference_equation(Y, Tdiscr)
diff_eq=(Y(1)-8*Y(2)+8*Y(4)-Y(5))/(12*Tdiscr);

Добавлено через 2 минуты
что важно, при предыдущих испытаниях результаты отличались и все было норм...
по-моему, реузльтаты должны отличаться...

Добавлено через 37 минут
ух, ребят, закипаю

эта проблема достала

Зосима · 22.08.2013, 09:51

А почему ты вместо этого:

Matlab M
1
2
3
4
5
6
7
8
 for  i=1:(N+1);
        for j=1:Lx
          B(i,j)=X(j)^(N+1-i);
        end
    end
F=B';
A=inv(F'*F)*F'*Y'%Calculation of matrix A of coefficients polynomials
poly=A';

не используешь:

Matlab M
1
p = polyfit(X,Y,N);

алгоритм совершенно тот-же самый, только короче, быстрее и понятней.
И вместо этого:

Matlab M
1
2
poly_diff_curr_X=polyval(p_diff, X);
diff_poly=poly_diff_curr_X(3);

надо написать:

Matlab M
1
diff_poly = polyval(p_diff, X(3));

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

**обрати внимание на мою подпись!

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .