Heat transient equation not executing

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 件のコメント

Torsten
Torsten 2022 年 9 月 17 日
Please show us your attempt.
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 日

0 投票

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 日
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))

サインインしてコメントする。

カテゴリ

ヘルプ センター および File ExchangeMathematics についてさらに検索

製品

リリース

R2019b

タグ

質問済み:

2022 年 9 月 17 日

コメント済み:

2022 年 9 月 17 日

Community Treasure Hunt

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

Start Hunting!

Translated by