کد متلب روش های عددی، مثال های کتاب سادیکو، فصل 3، مثال 3-1

کد متلب روش های عددی، مثال های کتاب سادیکو، فصل 3، مثال 3-1

%Numerical Techiniques in EM

%Example 3.1 page 128

%rja 2/20/2007

clear all; format compact; tic

%(a) Explicit Method

    delx = 0.1; %change this to .01 for a pretty picture

    r = .5; 

    k = 1;

    delt = k*r*delx^2;

    Tsteps = round(.4/delt); %give about .4 sec

    X1 = zeros(Tsteps,1/delx+1);

    %Initial conditions

    X1(1,:) = 100;

    X1(1,[1,end])=50;

    for row = 2:size(X1,1)

        for col = 2:size(X1,2)-1

            X1(row,col) = .5*(X1(row-1,col-1)+X1(row-1,col+1));

        end

    end

    figure(1),imagesc(0:delx:1,(0:delt:Tsteps*delt),X1),colorbar

        ylabel('\leftarrow time (sec)')

        xlabel('x')

        title('Example 3.1: Explicit Method')

    if (delx==.1)&&(delt==.005)

        dispmat = [(0:delt:0.03)', X1(1:7,1:7)];        

        disp(sprintf('\nCompare to Table 3.3, page 130'))   

        disp(num2str(dispmat))

    end

%(b) Implicit Method

%clear all;

    delx = 0.2;

    r = 1;

    k = 1;

    delt = k*r*delx^2;

    Tsteps = round(.4/delt); %give about .4 sec

    X2 = zeros(Tsteps,1/delx+1);

    %Initial conditions

    X2(1,:) = 100;

    X2(1,[1,end])=50;

    %Create Crank-Nicholson structure

    %eqn (3.16)

    numX = 1/delx-1;

    PHI = diag(4*ones(1,numX),0)+diag(-ones(1,numX-1),-1)+diag(-ones(1,numX-1),1);

    

    for ii = 2:size(X2,1)

        b = (X2(ii-1,1:numX)+X2(ii-1,3:numX+2)).';

        X2(ii,2:end-1) = PHI\b;

    end

    

        figure(2),imagesc(0:delx:1,(0:delt:Tsteps*delt),X2),colorbar

        ylabel('\leftarrow time (sec)')

        xlabel('x')

        title('Example 3.1: Implicit Method')

toc

جزییات وام ازدواج

کد متلب روش های عددی، مثال های کتاب سادیکو، فصل 3، مثال 3-10

کد متلب روش های عددی، مثال های کتاب سادیکو، فصل 3، مثال 3-9