What is wrong in code of diffusion equation problem?

11 次查看（过去 30 天）
Attila Monteiro de Abreu
回答： Byron Castillo 2021-1-22
a=4; b=4; na=20; nb=16; k=5; to=80;
clear T
x=linspace(0,a,na);
y=linspace(0,b,nb);
[X,Y]=meshgrid(x,y);
for i=1:nb
for j=1:na
T(i,j)=0;
For n=1:k
Ns=2*n-1;
T(i,j)=T(i,j)+sin(ns*pi*X(i,j)/a).*sinh(ns*pi*Y(i,j)/a)/(sinh(ns*pi*b/a)*ns);
End
T(i,j)=T(i,j)*4*To/pi;
End
End
Mesh(X,Y,T)
Xlabel(‘x (m)’);ylabel(‘y (m)’);zlabel(‘T(^oC)’)

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

采纳的回答

Walter Roberson
Walter Roberson 2016-9-18
First of all you need to use for and end not For and End -- MATLAB is case sensitive.
Secondly, your line
Xlabel(‘x (m)’);ylabel(‘y (m)’);zlabel(‘T(^oC)’)
uses Unicode Character 'LEFT SINGLE QUOTATION MARK' (U+2018) and Unicode Character 'RIGHT SINGLE QUOTATION MARK' (U+2019) instead of apostrophe. U+2018 and U+2019 are also known as "smart quotes" because some editors automatically change apostrophes around words to these quotes because those quotes are used in publishing and writing text (but seldom in computer code.)
  1 个评论
Attila Monteiro de Abreu
Attila Monteiro de Abreu 2016-9-18
编辑：Attila Monteiro de Abreu 2016-9-18
Thank you so much brother!!! You help me in my work of conclusion of course in Engenering Mechanical! I'm really greatfull!!

请先登录，再进行评论。

更多回答（1 个）

Byron Castillo
Byron Castillo 2021-1-22
a=5; b=4; na=20; nb=16; k=5; T0=80;
clear T
x=linspace(0,a,na);
y=linspace(0,b,nb);
[X,Y]=meshgrid(x,y);
for i=1:nb
for j=1:na
T(i,j)=0;
for n=1:k
ns=2*n-1;
T(i,j)=T(i,j)+sin(ns*pi*X(i,j)/a).*sinh(ns*pi*Y(i,j)/
a)/(sinh(ns*pi*b/a)*ns);
end
T(i,j) = T(i,j)*4*T0/pi;
end
end
mesh(X,Y,T)
xlabel('x (m)'); ylabel('y (m)'); zlabel('T ( ^oC)')

类别

在 Help Center 和 File Exchange 中查找有关 Physics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by