how to remove matrix dimension error ?

2 次查看(过去 30 天)
Qudsiya  Irum
Qudsiya Irum 2017-7-27
评论: KSSV 2017-7-27
I want solve Non Linear 1-D Heat Equation using Galerkin Finite Element Method. As I entered its ODE to solve by backward euler method, it is showing matrix dimension error. AA and MM are element stiffness and mass matrix of order(2501*2501) double respectively.
MM= spdiags([dx/6*ones(nx+1,1) 2*dx/3*ones(nx+1,1) dx/6*ones(nx+1,1)],[-1 0 1],nx+1,nx+1)
AA = spdiags([-1/dx*ones(nx+1,1) 2/dx*ones(nx+1,1) -1/dx*ones(nx+1,1)],[-1 0 1],nx+1,nx+1);
for mn = 2:1:M
T_solid_new_t(mn+1) = {MM+ k.*AA.*dt + dt*MM.*ht}\(MM.*T_solid(mn)+(dt.*dx.*Hs)+(ht.*Tp_g.*dx.*dt))
end

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by