Plotting temperature distribution versus time.

9 次查看（过去 30 天）

Hasan Al Tarify 2021-11-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1584564-plotting-temperature-distribution-versus-time

编辑： 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);