3.2. Пример программы:

function dirihle;

%решение задачи Дирихле по методу Галеркина

global be0 be1

clc

be0=0;

be1=1;

N=3;

M=40;

for i=1:N

F1 = @(x)f(x).*sin(i*pi*x)-g(x).*(be1-be0).*cos(i*pi*x);

d(i) = quad(F1,0,1);

for k=1:N

F2 = @(x)g(x).*cos(i*pi*x).*cos(k*pi*x)*i*k*pi^2;

G(i,k)=quad(F2,0,1);

end;

end

d

G

a=d/G

for i=1:M+1

xt(i)=(i-1)/M;

yt(i)=ut(xt(i));

y(i)=be0+(be1-be0)*xt(i);

for k=1:N

y(i)=y(i)+a(k)*sin(k*pi*xt(i));

end;

end;

hold off

plot(xt,y,xt,yt);

return

function z=f(x);

z=1; %1/(1+x/N);

return

function z=g(x);

z=1;%x./K;

return

function z=ut(x);%Точное решение

global be0 be1

z=-x*(x-1)/2+be0+(be1-be0)*x;

return

3.3. Индивидуальные задания

Требуется найти точное решение своего варианта задачи, после чего составить и написать программу решения задачи методом Галеркина. Сравнить точное решение с полученным, для чего вывести графики решений для нескольких значений K=3, 5, 10.

Для номера варианта N

,

Для N нечетного ⁰=0 ; ¹=(N+1)/N.

Для N четного ⁰=(N+1)/N, ¹=0.

При защите отчета ответить на вопросы по 4 лекции.

Задание 4 Нахождение минимума функционала

4.1. Пример программы

Изучить лекцию 5. В примере 1 для функционала , методом Ритца получено выражение для экстремали в виде ее разложения в ряд по синусам . В примере 3 для того же функционала для экстремали получено точное решение .

Запустите программу для получения графиков экстремали и убедитесь, что решения, полученные разными методами совпадают, если К>5.

function f=mygraf

%Графики экстремалей, полученных разными методами

clc

x=[];

y=[];

z=[];

hold off

for i=1:101

xi=(i-1)/100;

x(i)=xi;

y(i)=fy(xi);

z(i)=fz(xi);

end

plin1=plot(x,y,'b-',x,z,'k-');

set(plin1(1),'LineWidth',2);

set(plin1(2),'LineWidth',2);

grid on

return

function f=fy(r) %точное выражение для экстремали

f=(r^2-r)/2;

return

function f=fz(r)%выражение для экстремали полученное методом Ритца

s=0;

for k=1:2:3

s=s+sin(k*pi*r)/(k*pi)^3;

end;

f=4*s;

return

4.2. Индивидуальное задание

Решить эту же задачу при граничном условии , используя метод Ритца и метод Эйлера. После чего проверить правильность решения, построив графики полученных экстремалей с использованием выше приведенной программы.

Для номера варианта N, выбрать =N/(N+10).

Решение по методу Ритца искать в виде .

10