Heat transient equation not executing
1 回表示 (過去 30 日間)
古いコメントを表示

2 件のコメント
回答 (1 件)
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 件のコメント
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))
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!