12345678910111213141516 |
- % Setting up data;
- Y = hist(degMat,4);
- X = 1:length(Y);
- % Plotting the data
- plot(X,Y,'+r'), hold on
- % Polynomial curve fitting of log values so that we have a linear equation
- % Simplifying, log(Y) = log(a*X.^k) = log(X)*k + log(a)
- p = polyfit(log(X),log(Y),1);
- % Evaluating Coefficients
- k = p(1);
- % Accounting for the log transformation
- a = exp(p(2));
- % Final plot
- fplot(@(X) a*X.^k,[X(1) X(end)])
|