Heat transient equation not executing

1 次查看（过去 30 天）
Reji G
Reji G 2022-9-17
评论： Torsten 2022-9-17
I want to implement this formula in matlab. But i couldn't. So it'd be great if anyone can help me. Where x&t ranges from 0.1 to 10.
  2 个评论
Reji G
Reji G 2022-9-17
移动：Star Strider 2022-9-17
clc;
close all;
clear all;
T1=300;
Ti = 350;
L = 10;
alpha = 100;
syms n
T = symsum((exp(((-n*pi)/(2*L))*(alpha*t)))*sin((n*pi*x)/(2*L)),n,1,9)
T(x,t) = T1+((4/pi)*(Ti-T1)*T);

请先登录，再进行评论。

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

回答（1 个）

Torsten
Torsten 2022-9-17
Ran in:
The terms in your series do not agree with the formula you posted.
You must decide which one is correct and which one is wrong.
clc;
close all;
clear all;
T1 = 300;
Ti = 350;
L = 10;
alpha = 100;
syms n t x
T = symsum((exp(((-n*pi)/(2*L))*(alpha*t)))*sin((n*pi*x)/(2*L)),n,1,20);
T = T1+4/pi*(Ti-T1)*T;
T = matlabFunction(T);
x = 0:0.1:L;
t = 0:0.1:1;
[X,TIME] = meshgrid(x,t);
surf(TIME,X,T(TIME,X))
  2 个评论
Reji G
Reji G 2022-9-17
I need Temperature value(T(x,t)). an x y plot is required. Not the graph you plotted
Torsten
Torsten 2022-9-17
Ran in:
clc;
close all;
clear all;
T1 = 300;
Ti = 350;
L = 10;
alpha = 100;
syms n t x
T = symsum((exp(((-n*pi)/(2*L))*(alpha*t)))*sin((n*pi*x)/(2*L)),n,1,20);
T = T1+4/pi*(Ti-T1)*T;
T = matlabFunction(T);
x = 0:0.01:L;
t = 0.1;
figure(1)
plot(x,T(t,x))
x = L/2;
t = 0:0.001:0.5;
figure(2)
plot(t,T(t,x))

请先登录，再进行评论。

类别

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

标签

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by