finite difference method code

Question

Ronald Aono 2019-11-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/489050-finite-difference-method-code

回答： Walter Roberson 2019-11-4

% set domains limits and boundary conditions
xo = pi/2; xf = pi; yxo = 1; yxf = 1; N = 10;
% compute interval size and discrete x vector 
  dx = (xf-xo)/N;  dx2 = dx*dx; x = (xo+dx):dx:xf;
  
  % analytica solution (exact)
    xe = linspace(xo,xf,N);
    ye = (pi./(2*xe)).*(sin(xe) - 2*cos(xe));
    
    % arranging the matrix a
    %node 1
      a(1,1)=dx2-2; a(1,2)=1+(dx/(xo+dx));  b(1)= ((yxo*dx) /(xo*dx))-yxo;  
    
    
    for i = 2:N-1
        a(i,i-1) = (1-(dx/x(i)));
        a(i,i) = dx2-2;
        a(i,i+1) = (1+(dx/x(i)));
        b(i)=0;
    end
        
    a(N,N-1)=(2*xf+2*dx)/xf;  a(N,N-2)=-1;  b(N)=yxf*dx2+yxf+((2*yxf*dx)/xf);
    
    yi=a\b;

i keep getting the following error code

finite_1

Error using \

Matrix dimensions must agree.

Error in finite_1 (line 26)

yi=a\b

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2019-11-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/489050-finite-difference-method-code#answer_399623

https://www.mathworks.com/help/matlab/ref/mldivide.html

If A is a rectangular m-by-n matrix with m ~= n, and B is a matrix with m rows, then A\B returns a least-squares solution to the system of equations A*x= B.