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
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 CenterFile Exchange에서 Mathematics에 대해 자세히 알아보기

태그

제품


릴리스

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by