# Why do i get a blank graph?

조회 수: 32(최근 30일)
jpmontoya 2019년 6월 1일
댓글: jpmontoya 2019년 6월 2일
Hello, as you can see in the code below I am trying to make two plots but the problem is that the subplot (2,1,2) goes blank. I got nothing. I think maybe there is something wrong with the variables k1 and v1 since they depend on X(1) and X(3) which come from the solution of an ODE. Can somebody help me? Thanks in advance. function [T,X] = call_osc()
tspan = [0 2000];
x1_0=0.05;
x2_0=0.05;
x3_0=298.15;
x4_0=0.1;
odeset('RelTol', 1e-10, 'AbsTol', 1e-12);
[T,X]=ode15s(@osc,tspan,[x1_0 x2_0 x3_0 x4_0]);
R=8.314;
Ea1=80000;
alpha=11;
A=1;
k1= 0.8*exp(-Ea1/(alpha*R*X(3)));
v1=k1*A*X(1);
subplot(2,1,1)
plot(T,X(:,3))
grid on
xlabel('tiempo(s)','Fontsize',13,'FontWeight','Bold')
ylabel('Temperatura (K)','Fontsize',13,'FontWeight','Bold')
subplot(2,1,2)
plot(X(:,1),v1)
grid on
end
##### 댓글 수: 1표시숨기기 없음
Walter Roberson 2019년 6월 2일
It would help to have your osc function

댓글을 달려면 로그인하십시오.

### 채택된 답변

gonzalo Mier 2019년 6월 2일
편집: gonzalo Mier 2019년 6월 2일
You are computing v1 as
k1= 0.8*exp(-Ea1/(alpha*R*X(3)));
v1=k1*A*X(1);
v1 is a matriz of size 1x1. So when you do plot( "vector of size 1xn" , "cte" ) the output is an empty plot.
Check if the way to compute v1 is correct or you have to use X(:,1).
PD: You also can check this doing:
plot(X(:,1),v1,'*')
##### 댓글 수: 3표시숨기기 이전 댓글 수: 2
jpmontoya 2019년 6월 2일
Thank you so much. I did a lot of things considering the answers of Gonzalo and Walter and finally everything works perfect.

댓글을 달려면 로그인하십시오.

### 범주

Find more on 2-D and 3-D Plots in Help Center and File Exchange

### Community Treasure Hunt

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

Start Hunting!