How to plot a log log equation

2 views (last 30 days)
Okechukwu Anene
Okechukwu Anene on 7 Mar 2018
Commented: Akira Agata on 8 Mar 2018
I am trying to plot the following equation in matlab.. log ts = log c + nlog(v/a). I am required to find the value of c(which is the intercept) and n(which is the slope). The plot is ts against v/a.
Pls how do i proceed on this?
  2 Comments
Benjamin Kraus
Benjamin Kraus on 7 Mar 2018
Can you please clarify...
You said you want to plot an equation, which you can only do if you know the values of c and n. However, you then say you are required to find the values of c and n, which suggests you have data you want to plot, then fit a curve to the data.
Are you trying to plot an equation, or are you trying to fit a line to your data to extract values for c and n?
Okechukwu Anene
Okechukwu Anene on 7 Mar 2018
I apologize for the confusion. Yes I have data that i want to plot. i have 4 values for Ts : 3.48, 15.78, 10.17 and 8.13. Also, values for v/a are: 0.226, 0.536, 0.416, 0.367.
I am required to plot Ts against v/a on a log log graph and to find the slope n and intercept c using matlab.
The equation of the line is log ts = log c + nlog (v/a).

Sign in to comment.

Answers (1)

Akira Agata
Akira Agata on 8 Mar 2018
It seems that your data can be better fitted by 'Ts = b1 + b2*(v/a)' rather than 'Ts = b1 + b2*log(v/a)'. Here is an example to fit the data and plot in log-log graph.
Ts = [3.48, 15.78, 10.17, 8.13]';
va = [0.226, 0.536, 0.416, 0.367]';
X1 = [ones(size(va)) va];
b1 = X1\Ts;
X2 = [ones(size(va)) log(va)];
b2 = X2\Ts;
vaa = linspace(0.2,0.6)';
Ts1 = [ones(size(vaa)), vaa]*b1;
Ts2 = [ones(size(vaa)), log(vaa)]*b2;
loglog(va,Ts,'o-')
hold on
loglog(vaa,Ts1)
loglog(vaa,Ts2)
legend({'Data','Ts = b1 + b2*va','Ts = b1 + b2*log(va)'})
xlabel('v/a','FontSize',12)
ylabel('Ts','FontSize',12)
grid on
  2 Comments
Okechukwu Anene
Okechukwu Anene on 8 Mar 2018
Thank you so much for your response. however am running into some issues applying the code. I attached a screenshot of the error. Pls what do you think?
Thanks so much
Akira Agata
Akira Agata on 8 Mar 2018
Hi Okechukwu-san,
Please insert space between 'ones(size(va))' and 'va', or just copy&paste my sample code.

Sign in to comment.

Tags

Community Treasure Hunt

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

Start Hunting!