Hộp công cụ Matlab Pde: Vẽ giải pháp trên một dòng hoặc trên một thư mục con

Tôi đang sử dụng hộp công cụ pde Matlab để giải một phương trình elip nhất định trong 2D.

Giải pháp là tốt, mặc dù tôi cần phải vẽ nó dọc theo một đường nhất định, tức là cắt một lát phẳng từ lưới 3D đại diện cho giải pháp.

Tôi không thể tìm ra một cách thông minh liên quan đến các chức năng của hộp công cụ (nghĩa là không liên quan đến phép nội suy mức thấp trên lưới hình tam giác).

Bất kỳ trợ giúp đánh giá cao.

pde matlab

— Acorbe
nguồn

Tôi đã đăng câu hỏi tương tự trên trung tâm MATLAB và Bill Greene vui lòng cung cấp câu trả lời .

Tôi báo cáo ngắn gọn và mở rộng nó để nó có thể hữu ích cho người khác.

Giải pháp của Bill:

Đây là một cách để tạo ra một cốt truyện như vậy. Giả sử bạn có ma trận điểm được tạo bởi mesher PDE Toolbox, p và vectơ giải pháp, u. Hàm bên dưới sẽ tạo một biểu đồ của giải pháp đó dọc theo một đường được xác định bởi vị trí x và y của hai điểm cuối. Ví dụ của tôi là cho một giải pháp trên một hình vuông đơn vị và tôi muốn một âm mưu dọc theo dòng (0, .5) đến (1, .5). Tôi muốn bao gồm 25 điểm trong cốt truyện. Như bạn có thể thấy, công việc thực sự đang được thực hiện bởi TriScatteredInterpchức năng từ MATLAB cốt lõi.
plotAlongLine(p, u, [0,.5], [1,.5], 25);

function plotAlongLine(p, u, xy1, xy2, numpts)
    x = linspace(xy1(1),xy2(1),numpts);
    y = linspace(xy1(2),xy2(2),numpts);
    F = TriScatteredInterp(p(1,:)', p(2,:)', u);
    uxy = F(x,y);
    figure; plot(x, uxy);  //REM: x is chosen here as a curvilinear coordinate
end

Tôi muốn nhận xét thêm rằng hàm trước cho phép một biểu đồ giải pháp uhoặc các hàm của nó f(u), được cung cấp uđược xác định trên các nút lưới (như thường xảy ra đối với các giải pháp gần đúng FEM).

Nếu một nhu cầu cho các phần âm mưu của các chức năng xác định trên các trung tâm lưới (ví dụ như chức năng của grad u), cô có thể sử dụng trước các chức năng pdeprtnitrong đó sản xuất nút có giá trị chức năng ra khỏi trung tâm có giá trị chức năng .

— Acorbe
nguồn