Plotting temperature distribution versus time.

4 ビュー (過去 30 日間)
Hasan Al Tarify
Hasan Al Tarify 2021 年 11 月 11 日
編集済み: Hasan Al Tarify 2021 年 11 月 11 日
Hello everyone,
I have this temperature distribution function shown in the picture and wanted to plot it versus time. I would like to mention the (x) domain to be from 0 to 1. I have tried the following code but it did not work at all. Can someone please solve my problem? I would really appreciate your help and time spent on solving the problem. Thanks in advance.
clc
close all;
clear;
workspace;
format long g;
format compact;
fontSize = 18;
fprintf('Beginning to run %s.m ...\n', mfilename);
numPoints = 64;
t = 100;
L = 1;
x = linspace(0, L, numPoints);
t1 = linspace(0, t, numPoints);
n = 1;
T = zeros(numPoints);
Theta = zeros(numPoints);
for it = 1 : numPoints
this_t = t(it);
for jx = 1:numPoints
this_x = x(jx);
theSum = 0;
for n = 1 : 100
term1 = -100*(-1)^n/(pi*n);
term2 = exp(-(n*pi)^2)*this_t;
term3 = sin(n * pi * this_x);
theSum = theSum + term1 * term2 * term3;
end
T(it, jx) = -100*this_x + 100 + theSum;
end
end
imshow(T, [], 'XData', x, 'YData', t, 'InitialMagnification', 1600);
fplot(t,T);
axis('on', 'image');
title('T(t, x)', 'FontSize', 15);
xlabel('x', 'FontSize', 15);
ylabel('t', 'FontSize', 15);
fprintf('Done running %s.m ...\n', mfilename);

回答 (0 件)

カテゴリ

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

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by