Finite Difference Method to find Steady State

2 次查看（过去 30 天）

显示更早的评论

Jake 2020-4-8

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/516366-finite-difference-method-to-find-steady-state

评论： Ragul Kumar 2020-11-6

采纳的回答： Torsten

Hello experts,

I have to solve the following equation to find a converged solution.

At i=0, T=10 and at i=5, T=50. (n=5)

I want to get the converged solution of T at i=2, 3, 4.

I'm fairly comfortable with MATLAB but this is the first time that I'm using MATLAB for these kind of mathematical approach and I'm learning completely alone, relying on online materials.

How should I approach this? I'd appreciate your recommendations, reading materials, suggestions and/or pointers rather than a complete code at this point.

I know the question is vague but I'd welcome any suggestions for starting points :)

TIA!

3 个评论
显示 1更早的评论隐藏 1更早的评论

Jake 2020-4-8

Apologies, I've missed including that info properly.

When n=0, T_i^(n+1)=10, 73, 100, 85, 50 respectively at i=1, 2, 3, 4, 5.

Ragul Kumar 2020-11-6

Hello experts,

I am trying to solve the finite difference methof for crank nicolson scheme to 2d heat equation. please let me know if you have any MATLAB CODE for this

Boundary codition are

If you can kindly send me the matlab code, it will be very useful for my research work . thank you very much.

请先登录，再进行评论。

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

采纳的回答

Torsten 2020-4-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/516366-finite-difference-method-to-find-steady-state#answer_424813

编辑：Torsten 2020-4-9

在 MATLAB Online 中打开

T = zeros(6,5);
T(1:6,1) = 10;
T(1:6,5) = 50;
T(1,2) = 73;
T(1,3) = 100;
T(1,4) = 85;
for n=1:5
  for i =2:4
     T(n+1,i) = T(n,i) + 0.3*(T(n,i+1)-2*T(n,i)+T(n,i-1));
  end
end
T