Метод конечных разностей

@Archij · Регистрация: 06.09.2011

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

Здравствуйте! Нужна помощь с методом конечных разностей. На этом сайте http://www.home.uni-osnabrueck... rences.pdf есть хороший код для этого метода. Сам код:

Matlab M

function [x,f]=od_fdm(xlo,flo,xhi,fhi,NX)
% solve f"+f=0 on x=linspace(xlo,xhi,NX)
% with boundary values flo and fhi, resp.
% NX must be 3 or larger
N=NX-2; % number of unknowns
x=linspace(xlo,xhi,NX); %NX точек от xlo до xhi
h=x(2)-x(1); %шаг
main=(-2/h^2+1)*ones(1,N);
next=(1/h^2)*ones(1,N-1);
DE=diag(next,-1)+diag(main,0)+diag(next,1);
BV=zeros(N,1);
BV(1)=-flo/h^2;
BV(N)=-fhi/h^2;
sol=DE\BV;
f=[flo,sol',fhi];
end % od_fdm

Проблема в том, что эта функция для решения конкретного примера, f"+f=0. Где надо менять код, что-бы решить любое другое дифференциальное уравнение? Скажем, f''+1=0. Вопрос, эти строки всегда будут такие:

Matlab M

1 2	main=(-2/h^2+1)ones(1,N); next=(1/h^2)ones(1,N-1);

и эти

Matlab M

1 2	BV(1)=-flo/h^2; BV(N)=-fhi/h^2;

Эти по моему надо менять для другой функции. Но я не знаю как именно. Спасибо за помощь!

	09.01.2013, 17:52