PDEModel 객체를 사용하여 문제 풀기

  1. 문제를 Partial Differential Equation Toolbox™ 솔버에 맞는 올바른 형식으로 구성합니다. 자세한 내용은 Partial Differential Equation Toolbox를 사용하여 풀 수 있는 방정식 항목을 참조하십시오. 문제를 발산 형식으로 변환해야 하는 경우 Put Equations in Divergence Form 항목을 참조하십시오.

  2. PDEModel 모델 컨테이너를 만듭니다. 스칼라 PDE의 경우 인수 없이 createpde를 사용합니다.

    model = createpde();

    시스템(연립방정식)에 N개의 방정식이 있는 경우 입력 인수 N과 함께 createpde를 사용합니다.

    model = createpde(N);

  3. 지오메트리를 가져오거나 만듭니다. 자세한 내용은 지오메트리 및 메시 항목을 참조하십시오.

    importGeometry(model,"geometry.stl"); % importGeometry for 3-D
geometryFromEdges(model,g); % geometryFromEdges for 2-D

  4. 경계의 레이블을 알 수 있도록 지오메트리를 표시합니다.

    pdegplot(model,FaceLabels="on") % FaceLabels for 3-D
pdegplot(model,EdgeLabels="on") % EdgeLabels for 2-D

    3차원 모델의 레이블을 보려면 모델을 회전하거나, 투명하게 만들거나, 확대해야 할 수 있습니다. STL 파일 가져오기 항목을 참조하십시오.

  5. 경계 조건을 만듭니다. 자세한 내용은 Specify Boundary Conditions 항목을 참조하십시오.

    % Face for 3-D
applyBoundaryCondition(model,"dirichlet",Face=[2,3,5],u=[0,0]);
% Edge for 2-D
applyBoundaryCondition(model,"neumann",Edge=[1,4],g=1,q=eye(2));

  6. PDE 계수를 만듭니다.

    f = [1;2];
a = 0;
c = [1;3;5];
specifyCoefficients(model,m=0,d=0,c=c,a=a,f=f);

  7. 시간 종속 방정식에 대해, 또는 선택적으로 비선형 정상(stationary) 방정식에 대해 초기 조건을 만듭니다. Set Initial Conditions 항목을 참조하십시오.

  8. 메시를 만듭니다.

    generateMesh(model);

  9. 적절한 솔버를 호출합니다. 고유값 문제를 제외한 모든 문제에 대해서는 solvepde를 호출합니다.

    result = solvepde(model); % for stationary problems
result = solvepde(model,tlist); % for time-dependent problems

    고유값 문제에 대해서는 solvepdeeig를 호출합니다.

    result = solvepdeeig(model);

  10. 해를 검토합니다. 시각화 항목을 참조하십시오.

참고 항목

| | | | | | |

Teaching Resources