Блог

Главная - Вопросы по программированию - Как я могу выполнить расчет вероятности с помощью функции собственных волн, собранной функцией eigs() в MATLAB?

Как я могу выполнить расчет вероятности с помощью функции собственных волн, собранной функцией eigs() в MATLAB?

#matlab #eigenvalue

Вопрос:

     h = 0.0454973;
    n = 10;
    x = chebfun('x',[0.5,3]);
    V = 545.0362942355468 - 1998.817146282548*x   2980.608037604876*x^2 
        - 2314.110473743709*x^3   989.3591697175247*x^4 
        - 221.48586861701781*x^5   20.33254549765507*x^6;
    %% Eigenvalue computation:
    [xmin, xmax] = domain(V);                          % domain of problem
    % Create a CHEBOP with Dirichlet BCs
    L = chebop([xmin, xmax]);
    L.lbc = 0; L.rbc = 0;
    L.op = @(x,u) -h^2*diff(u,2)   V.*u;               % Schroedinger operator
    [U, D] = eigs(L, n, 'sr');                         % compute evals/efuns
    d = diag(D);                                       % vector of evals
    [d, ii] = sort(d);                                 % sort them
    U = U(:,ii);
    plot(U(:,8))

Вот код, который я использую для поиска собственной функции и собственного значения. Я могу построить график собственной функции, используя последнюю строку кода. Однако я не знаю, как преобразовать U во что-то, что можно интегрировать, чтобы вычислить вероятность соответствующей собственной функции. Кто-нибудь может мне помочь?

Метки: ,