How do I add a curve fitted plot to this code?
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
clc
clear
close all
xls=xlsread("Exam2_Data.xls");
t=xls(:,1);
y=xls(:,2);
SST=sum((y-mean(y)).^2);
for i=1:3
fprintf('For %d order fit\n\n',i);
disp('Polynomial coefficients are');
C=polyfit(t,y,i);
SSR=sum((y-polyval(C,t)).^2);
R2=1-SSR/SST;
fprintf('R^2 is %f\n\n',R2);
end
for i=2:length(t)-1
T=[t(i-1) t(i) t(i+1)];
Y=[y(i-1) y(i) y(i+1)];
if(t(i)<33.75&&t(i+1)>33.75)
C=polyfit(T,Y,2);
disp('Value at 33.75 is ');
polyval(C,33.75)
return
end
end
2 commentaires
Rik
le 9 Avr 2019
Do you mean like in a live script? Also, are you sure about the return? Don't you mean break? Return causes the current running function to return to the caller, while break causes the current (inner) loop to exit.
Réponses (1)
Sambit Senapati
le 17 Avr 2019
you can use the polyval and plot functions in tandem to do the job here:
instead of calling polyval(C,33.75) at a single data point, do the following:
>> Y_new=polyval(C,T);
>> plot(T,Y_new)
You could do this for your original data set as well.
0 commentaires
Voir également
Catégories
En savoir plus sur Interpolation dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!