MATLAB: Multiple graph plotting error

graph

ro=10;

ri=30;

t=0:0.1:2000;

ci=0.5

c_L=(ci*ri)/ro.*(1-exp(1).^((-ro.*t./((ri-ro).*t+1500))));

plot(t,c_L)

hold on

ro=15;

ri=25;

c_L=(ci*ri)/ro.*(1-exp(1).^((-ro.*t./((ri-ro).*t+1500))));

plot(t,c_L)

ro=25;

ri=15;

c_L=(ci*ri)/ro.*(1-exp(1).^((-ro.*t./((ri-ro).*t+1500))));

plot(t,c_L)

For the interval t it should be 0 to 2000, but when I run MATLAB with that interval graph looks very weired

However, with interval 0 to 100 looks fine.

What should I do?

Best Answer

See: multiple graphs plotting error

Related Solutions

MATLAB: Multiple graphs plotting error

The ‘t’ vector for the first series of curves is:

t=0:0.1:100;

and for the second:

t=0:0.1:2000;

That likely explains the difference.

To plot ‘t’ from 0 to 2000, define ‘t’ over the desired range.

MATLAB: How do determine the pressure? Please check me the code.

To do numerical integration, you need to use the integral function. Here, it is also necessary to use arrayfun to integrate over your ‘r’ vector. I changed ‘vth’ and also the plots, so ‘p_r’ is plotted as a function of ‘r’ (and ‘z’).

This runs, without errors or warnings. I will let your determine if it is correct:

t = 0:0.1:10*pi;
r = linspace (0, 1, numel (t));
z = linspace (0, 1, numel (t));
a = 0.009;
Gamm_inf = 2000;
nu = 3;
p0 = 1;
ro = 2;
vth = @(r) (Gamm_inf)./(2*pi*nu.*r) .* (1 - exp(-(a*r/(2*nu))));
p_r = p0 + ro.*arrayfun(@(r)integral(@(r)vth(r).*vth(r)./r,1E-8,r, 'ArrayValued',1),r)-(ro*a^2)/2.*(r.^2 + 4.*z.^2);
figure
plot(r, p_r);
grid
figure
plot3(r, z ,p_r);
grid on

Experiment to get the result you want.

Best Answer

Related Solutions

MATLAB: Multiple graphs plotting error

MATLAB: How do determine the pressure? Please check me the code.

Related Question