Eigenvalues и Eigenmodes квадрата

В этом примере показано, как вычислить собственные значения и eigenmodes квадратной области.

Проблема УЧП собственного значения -Δu=λu. Этот пример находит собственные значения меньшими, чем 10 и соответствующий eigenmodes.

Создайте модель. Импортируйте и постройте геометрию. Файл описания геометрии для этой проблемы называется squareg.m.

model = createpde();
geometryFromEdges(model,@squareg);

pdegplot(model,'EdgeLabels','on')
ylim([-1.5,1.5])
axis equal

Задайте граничное условие Дирихле u=0 для левого контура.

applyBoundaryCondition(model,'dirichlet','Edge',4,'u',0);

Задайте нулевое Нейманово граничное условие для верхней и нижней границы.

applyBoundaryCondition(model,'neumann','Edge',[1,3],'g',0,'q',0);

Задайте обобщенное Нейманово условие un-34u=0 для правильного контура.

applyBoundaryCondition(model,'neumann','Edge',2,'g',0,'q',-3/4);

Коэффициенты УЧП собственного значения для этой проблемы являются c = 1, = 0, и d = 1. Можно ввести область значений собственного значения r как векторный [-Inf 10].

specifyCoefficients(model,'m',0,'d',1,'c',1,'a',0,'f',0);
r = [-Inf,10];

Создайте mesh и решите задачу.

generateMesh(model,'Hmax',0.05);
results = solvepdeeig(model,r);
              Basis= 10,  Time=   0.87,  New conv eig=  0
              Basis= 11,  Time=   0.93,  New conv eig=  0
              Basis= 12,  Time=   0.96,  New conv eig=  1
              Basis= 13,  Time=   0.99,  New conv eig=  1
              Basis= 14,  Time=   1.02,  New conv eig=  1
              Basis= 15,  Time=   1.05,  New conv eig=  1
              Basis= 16,  Time=   1.08,  New conv eig=  1
              Basis= 17,  Time=   1.12,  New conv eig=  1
              Basis= 18,  Time=   1.17,  New conv eig=  2
              Basis= 19,  Time=   1.20,  New conv eig=  2
              Basis= 20,  Time=   1.26,  New conv eig=  2
              Basis= 21,  Time=   1.32,  New conv eig=  3
              Basis= 22,  Time=   1.37,  New conv eig=  3
              Basis= 23,  Time=   1.43,  New conv eig=  4
              Basis= 24,  Time=   1.63,  New conv eig=  6
End of sweep: Basis= 24,  Time=   1.65,  New conv eig=  3
              Basis= 13,  Time=   1.99,  New conv eig=  0
              Basis= 14,  Time=   2.02,  New conv eig=  0
              Basis= 15,  Time=   2.06,  New conv eig=  0
              Basis= 16,  Time=   2.09,  New conv eig=  0
              Basis= 17,  Time=   2.12,  New conv eig=  0
              Basis= 18,  Time=   2.15,  New conv eig=  0
              Basis= 19,  Time=   2.23,  New conv eig=  0
              Basis= 20,  Time=   2.29,  New conv eig=  0
              Basis= 21,  Time=   2.33,  New conv eig=  0
              Basis= 22,  Time=   2.40,  New conv eig=  1
              Basis= 23,  Time=   2.46,  New conv eig=  2
End of sweep: Basis= 23,  Time=   2.46,  New conv eig=  0
              Basis= 13,  Time=   2.76,  New conv eig=  1
              Basis= 14,  Time=   2.81,  New conv eig=  1
              Basis= 15,  Time=   2.85,  New conv eig=  1
              Basis= 16,  Time=   2.87,  New conv eig=  3
              Basis= 17,  Time=   2.92,  New conv eig=  3
              Basis= 18,  Time=   2.94,  New conv eig=  3
              Basis= 19,  Time=   2.97,  New conv eig=  3
              Basis= 20,  Time=   3.02,  New conv eig=  3
              Basis= 21,  Time=   3.09,  New conv eig=  4
End of sweep: Basis= 21,  Time=   3.09,  New conv eig=  4
              Basis= 17,  Time=   3.44,  New conv eig=  0
End of sweep: Basis= 17,  Time=   3.47,  New conv eig=  0

Существует шесть собственных значений, меньших, чем 10 для этой проблемы.

l = results.Eigenvalues
l = 5×1

   -0.4146
    2.0528
    4.8019
    7.2693
    9.4550

Постройте первые и последние собственные функции в заданной области.

u = results.Eigenvectors;
pdeplot(model,'XYData',u(:,1));

pdeplot(model,'XYData',u(:,length(l)));

Эта проблема отделима, означая

u(x,y)=f(x)g(y).

Функции f и g являются собственными функциями в направлениях X и Y, соответственно. В направлении X первый eigenmode является медленно увеличивающейся показательной функцией. Более высокие режимы включают синусоиды. В направлении Y первый eigenmode является (постоянной) прямой линией, вторым является половина косинуса, третьим является полный косинус, четвертыми являются полтора полных косинуса и т.д. Эти eigenmodes в направлении Y сопоставлены с собственными значениями

0,π24,4π24,9π24,...

Возможно проследить предыдущие собственные значения в собственных значениях решения. Смотря на график первого eigenmode, вы видите, что он составлен из первого eigenmodes в направлениях X и Y. Второй eigenmode составлен из первого eigenmode в направлении X и второго eigenmode в направлении Y.

Посмотрите на различие между первым и вторым собственным значением по сравнению с π2/4:

l(2) - l(1) - pi^2/4
ans = 1.6751e-07

Аналогично, пятый eigenmode составлен из первого eigenmode в направлении X и третьего eigenmode в направлении Y. Как ожидалось, l(5)-l(1) приблизительно равно π2:

l(5) - l(1) - pi^2
ans = 6.2135e-06

Можно исследовать более высокие режимы путем увеличения поисковой области значений, чтобы включать собственные значения, больше, чем 10.

Для просмотра документации необходимо авторизоваться на сайте