How to solve this question? Help :'(

2020 6 30
1 个回答
回答已采纳
更新时间：2020 6 30
3 次查看（30 天）

1 个投票

Question 4 (20 marks)
Laplace’s equation: Determine the distribution of temperature in a rectangular plane section, subject to a temperature distribution around its edges as follows:
The section shape, the boundary temperature distribution section and two chosen nodes are shown in the figure.
The temperature distribution is described by Laplace’s equation. Solving this equation by the finite difference method to nodes 1 and 2 of the mesh shown in the figure. This gives
(233.33 +T2 + 0 100 - 4T1)/h^2 = 2
(333.33 + 250 + 0 + T1 - 4T2)/h^2 = 2
where and are the unknown temperatures at nodes 1 and 2, respectively, and
Rearranging these equations gives
[■(-4&1@1&-4)][■(T_1@T_2 )]=[■(-333.33@-583.33)][■(-4&1@1&-4)][■(T_1@T_2 )]=[■(-333.33@-583.33)]
Solving this equation, we have and .
If we require a more accurate solution of Laplace’s equation, then we must use more nodes and the computation burden increases rapidly. Write a MATLAB script to solve the Laplace’s equation for any number of nodes in a square domain only. Run your codes for number of nodes

4 个评论
显示 2更早的评论 隐藏 2更早的评论

Rik
Rik 2020-6-30
I would almost say good job on posting your homework question, but you forgot to fix the formatting.
You also forgot to ask a question. This is your homework. What did you try?
raizo dono
raizo dono 2020-6-30
编辑：raizo dono 2020-6-30
help me Rik, this is for my final assessment. I have attached the fixed formatting for the question and what I've tried for the question. Thank you for your time.
Rik
Rik 2020-6-30
What is your question? When I run your code I get an output.
You can find guidelines for posting homework on this forum here. Please use the tools explained on this page to make your question more readable and attach the ellipgen function as an m file.
raizo dono
raizo dono 2020-6-30
编辑：raizo dono 2020-6-30
The ellipgen and Untitled5 file that I attached earlier is the halfway answer for the question. What I am trying to find is how to run codes for number of nodes n=50, 100 and 150. I cannot find any solution for this. Attached is the m file for the ellipgen

请先登录，再进行评论。

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

 采纳的回答

Rik
Rik 2020-6-30

0 个投票

You should learn about the linspace function. I would also suggest you check if you have x and y correct. I kept your convention.
You were also limiting the axis so changing the number of nodes doesn't behave as expected.
Lx = 3; Ly = 2;
nx = 120; ny = 12; hx = Lx/nx; hy = Ly/ny;
by0 = linspace(0,Lx,nx+1); %by0 = 0*[0:hx:Lx];
byn = 200+(100/3)*linspace(0,Lx,nx+1).^2; %byn = 200+(100/3)*[0:hx:Lx].^2;
bx0 = 100*linspace(0,Ly,ny+1); %bx0 = 100*[0:hy:Ly];
bxn = 250*linspace(0,Ly,ny+1); %bxn = 250*[0:hy:Ly];
F = zeros(nx+1,ny+1); G = F;
a = ellipgen(nx,hx,ny,hy,G,F,bx0,bxn,by0,byn);
aa = flipud(a); colormap(gray)
surfl(aa)
xlabel('x direction')
ylabel('y direction')
zlabel('Temperature')
axis([0 nx 0 ny 0 500]) %axis([0 12 0 12 0 500])

2 个评论
显示 无 隐藏 无

raizo dono
raizo dono 2020-6-30
So, how do I unlimit the axis to change the node numbers as asked in the question?
Rik
Rik 2020-6-30
I already edited that line for you. The only thing that is left for you is to figure out how you can run the same code with different inputs. Hint: use a for-loop, wrap your code in a function, or both.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Functions 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by