Heat transient equation not executing

2 ビュー (過去 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 件のコメント
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 日
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))

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

カテゴリ

Help Center および File ExchangeCalculus についてさらに検索

タグ

製品


リリース

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by