Plotting temperature distribution versus time.

4 vues (au cours des 30 derniers jours)
Hasan Al Tarify
Hasan Al Tarify le 11 Nov 2021
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);

Réponses (0)

Catégories

En savoir plus sur Mathematics dans Help Center et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by