Scribd Logo
Upload

Welcome to Scribd!

  • 50 views

    Name Balcha Areda ID NO pgr/18233/11: y X For y X y X U y X For y X y X U

    Uploaded by

    hundasa chala
    The document describes solving a Poisson problem using a finite difference method in MATLAB. It provides the Poisson equation, boundary conditions, and MATLAB code to discretize and solve the equation on a uniform grid. The code uses centered finite differences to approximate the Laplace operator and solves the resulting linear system to obtain the solution u on the grid. It then plots the numerical solution.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd

    Name Balcha Areda ID NO pgr/18233/11: y X For y X y X U y X For y X y X U

    Uploaded by

    hundasa chala
    50 views3 pages
    The document describes solving a Poisson problem using a finite difference method in MATLAB. It provides the Poisson equation, boundary conditions, and MATLAB code to discretize and solve the equation on a uniform grid. The code uses centered finite differences to approximate the Laplace operator and solves the resulting linear system to obtain the solution u on the grid. It then plots the numerical solution.

    Original Description:

    ggg

    Original Title

    Math Piosson

    Copyright

    © © All Rights Reserved

    Available Formats

    DOCX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    The document describes solving a Poisson problem using a finite difference method in MATLAB. It provides the Poisson equation, boundary conditions, and MATLAB code to discretize and solve the equation on a uniform grid. The code uses centered finite differences to approximate the Laplace operator and solves the resulting linear system to obtain the solution u on the grid. It then plots the numerical solution.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    50 views3 pages

    Name Balcha Areda ID NO pgr/18233/11: y X For y X y X U y X For y X y X U

    Uploaded by

    hundasa chala
    The document describes solving a Poisson problem using a finite difference method in MATLAB. It provides the Poisson equation, boundary conditions, and MATLAB code to discretize and solve the equation on a uniform grid. The code uses centered finite differences to approximate the Laplace operator and solves the resulting linear system to obtain the solution u on the grid. It then plots the numerical solution.

    Copyright:

    © All Rights Reserved
    Download as DOCX, PDF, TXT or read online from Scribd
    Download as docx, pdf, or txt
    of 3

    Name Balcha Areda ID NO pgr/18233/11

    Mathlab practice

    Use the mathlab code to solve the Poisson problem

     ux, y   5 
    2
    sin( x) cos( 2y ) for( x, y )    (0,1)  (0,1)
    2

    u ( x, y )  sin( x) cos( 2y ) for( x, y )  

    function [u,x,y] = fd2poisson(pfunc,bfunc,a,b,n)


    h = (b-a)/(n+1); % Mesh spacing
    [x,y] = meshgrid(a:h:b); % Uniform mesh, including boundary
    points.
    % Compute u on the boundary from the Dirichlet boundary condition
    ub = zeros(n,n);
    idx = 2:n+1;
    idy = 2:n+1;
    % West and East boundaries need special attention
    ub(:,1) = feval(bfunc,x(idx,1),y(idy,1)); % West Boundary
    ub(:,n) = feval(bfunc,x(idx,n+2),y(idy,n+2)); % East Boundary
    % Now the North and South boundaries
    ub(1,1:n) = ub(1,1:n) + feval(bfunc,x(1,idx),y(1,idy));
    ub(n,1:n) = ub(n,1:n) + feval(bfunc,x(n+2,idx),y(n+2,idy));
    % Convert ub to a vector using column reordering
    ub = (1/h^2)*reshape(ub,n*n,1);
    % Evaluate the RHS of Poisson’s equation at the interior points.
    f = feval(pfunc,x(idx,idy),y(idx,idy));
    % Convert f to a vector using column reordering
    f = reshape(f,n*n,1);
    % Create the D2x and D2y matrices
    z = [-2;1;zeros(n-2,1)];
    D2x = 1/h^2*kron(toeplitz(z,z),eye(n));
    D2y = 1/h^2*kron(eye(n),toeplitz(z,z));
    % Solve the system
    u = (D2x + D2y)\(f-ub);
    % Convert u from a column vector to a matrix to make it easier
    to work with
    % for plotting.
    u = reshape(u,n,n);
    % Append on to u the boundary values from the Dirichlet
    condition.
    u = [[feval(bfunc,x(1,1:n+2),y(1,1:n+2))];...
    [[feval(bfunc,x(2:n+1,1),y(2:n+1,1))] u ...
    [feval(bfunc,x(2:n+1,n+2),y(2:n+1,n+2))]];...
    [feval(bfunc,x(n+2,1:n+2),y(n+2,1:n+2))]];

    function testFdPoisson
    f = inline('-5*pi^2*sin(pi*x).*cos(2*pi*y)');
    g = inline('sin(pi*x).*cos(2*pi*y)');
    [u,x,y] = fd2poisson(f,g,0,1,5)
    h = x(1,2) - x(1,1);
    % Plot solution
    figure, set(gcf,'DefaultAxesFontSize',8,'PaperPosition', [0
    0 3.5 3.5]),
    mesh(x,y,u), ([0 0 0]), xlabel('x'), ylabel('y'),
    zlabel('u(x,y)'),
    title(strcat('Numerical Solun to Poisson Equation' ));
    >> testFdPoisson

    u =

    0 0.5000 0.8660 1.0000 0.8660 0.5000 0.0000

    0 0.2441 0.4227 0.4881 0.4227 0.2441 0.0000

    0 -0.2892 -0.5009 -0.5784 -0.5009 -0.2892 -0.0000

    0 -0.5572 -0.9652 -1.1145 -0.9652 -0.5572 -0.0000

    0 -0.2892 -0.5009 -0.5784 -0.5009 -0.2892 -0.0000

    0 0.2441 0.4227 0.4881 0.4227 0.2441 0.0000

    0 0.5000 0.8660 1.0000 0.8660 0.5000 0.0000

    You might also like