Why is my plot plotting blank?
显示 更早的评论
clear all
close all
%DEFINING THE GRID (MESH) AND SIZE OF THE DOMAIN
L = 50; %Length of the vertical rock section (m)(use meter, kelvin, second always)
dz = 0.1; % size of the elements (m)
N = L/dz+1; % Number of nodes in the rock column
%DEFINING THE PHYSICAL PROPERTIES
K = 1.8; %Thermal conductivity of the rocks
Cp = 840 %specific heat of the rocks
rho = 3000; %density of the rocks
kappa = K/(Cp*rho); %Diffusivity of the rocks
%INITIALIZE THE VARIABLES
time = 0; %I define the initial time
Tbot = 140;
Tbotnew = 140;
Ttop = 122; %Initial temperature, in degrees centigrade (C)
Ttopnew = linspace(Ttop,Tbot,N); %Initialize the new temperature array
%Plot initial T conditions
x= 1:N:50
y= Ttopnew
plot(x,y,'r');
xlabel('Depth of rock column (m)')
ylabel('Temperature (C)')
1 个评论
Stephen23
2015-3-6
Please edit your question and format the code correctly using the {} Code button that you will find above the textbox.
采纳的回答
Jan
2015-3-7
0 个投票
You N is 501. Then x= 1:N:50 let x be the scalar 1. Inconsequence the plot use one x value only and when I run your code I see a vertical line, as expected.
When you omit the evil clear all you can use the debugger to step through the code line by line.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Animation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
