How to code equation

조회 수: 13 (최근 30일)
Cesar Cardenas
Cesar Cardenas 2023년 1월 24일
댓글: Cesar Cardenas 2023년 2월 3일
Hello, just wondering if this equation could be coded like this:
Any feedback will be greatly appreciated. Thanks
m = 10e-3
a = 0.05
c = 0.1
x(t) = m/a e^-(c/a)t

채택된 답변

Sulaymon Eshkabilov
Sulaymon Eshkabilov 2023년 1월 24일
A couple of small errors have been corrected:
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t)(m/a)*exp(-(c/a)*t);
t = linspace(0,1.3);
plot(t, x(t))
xlabel('t')
ylabel('x(t)')
grid on
  댓글 수: 1
Cesar Cardenas
Cesar Cardenas 2023년 1월 27일
Thanks much for your help. I added another equation and would like to know how to plot both on the same graph? not sure about it. Thanks
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t)(m/a)*exp(-(c/a)*t);
x = @(t)(m/a)*1-exp(-(a/m)*t);
t = linspace(0,1.5);
plot(t, x(t))
xlabel('t')
ylabel('x(t)')
grid on

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

추가 답변 (2개)

Torsten
Torsten 2023년 1월 24일
이동: Torsten 2023년 1월 24일
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t) m/a*exp(-c/a*t);
t = 0:0.01:2;
plot(t,x(t))
  댓글 수: 1
Cesar Cardenas
Cesar Cardenas 2023년 1월 27일
Thanks much for your help. I added another equation and would like to know how to plot both on the same graph? not sure about it. Thanks
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t)(m/a)*exp(-(c/a)*t);
x = @(t)(m/a)*1-exp(-(a/m)*t);
t = linspace(0,1.5);
plot(t, x(t))
xlabel('t')
ylabel('x(t)')
grid on

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


Sulaymon Eshkabilov
Sulaymon Eshkabilov 2023년 1월 27일
Put both functions in plot() command or use hold on:
m = 10e-3;
a = 0.05;
c = 0.1;
x1 = @(t)(m/a)*exp(-(c/a)*t);
x2 = @(t)(m/a)*1-exp(-(a/m)*t);
t = linspace(0,1.5);
plot(t, x1(t), 'r-', t, x2(t), 'b-.', 'linewidth', 2)
xlabel('$t$', 'Interpreter', 'latex')
ylabel('$x(t)$', 'Interpreter', 'latex')
legend('$x_1(t)$', '$x_2(t)$','Interpreter', 'latex', 'location', 'SE')
grid on
%% Alternative solution:
m = 10e-3;
a = 0.05;
c = 0.1;
x1 = @(t)(m/a)*exp(-(c/a)*t);
x2 = @(t)(m/a)*1-exp(-(a/m)*t);
figure
fplot(x1, [0, 1.5], 'r-', 'linewidth', 2)
hold on
fplot(x2, [0, 1.5], 'b-.', 'linewidth', 2)
xlabel('$t$', 'Interpreter', 'latex')
ylabel('$x(t)$', 'Interpreter', 'latex')
legend('$x_1(t)$', '$x_2(t)$','Interpreter', 'latex', 'location', 'SE')
grid on
  댓글 수: 1
Cesar Cardenas
Cesar Cardenas 2023년 2월 3일
Hello, thanks much for your response, really useful. Now, I'm trying to code these equations. This is my attempt. I would like to plot Cf vs x/L, but x/L as I have it shows me an error..the same with delta/x, not sure why? any help will be greatly appreciated. Thanks much.
Re = 0:500:1000;
x/L = 0:0.2:1;
Cf = 0.664*(Re)^-0.5;
delta/x = 5.0*(Re)^-0.5;

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

카테고리

Help CenterFile Exchange에서 Specifying Target for Graphics Output에 대해 자세히 알아보기

태그

Community Treasure Hunt

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

Start Hunting!

Translated by