|
|
|
Scilab - polynomial fit und ableitung davon |
|
| xquadxq |
Forum-Newbie
|
 |
Beiträge: 8
|
|
|
|
Anmeldedatum: 13.03.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 04.05.2020, 17:46
Titel: Scilab - polynomial fit und ableitung davon
|
 |
Hallo,
ich habe ein .csv-file mit engauge digitizer erstellt wie hier:
x,D
100.1,200.3
101.3,200.4
...
Ich lade es in Scilab mit:
clear
Array=csvRead('C:\Users\Karl\Desktop\SCILAB\D_Werte.csv');
col1 = Array(:, 1);
col2 = Array(:, 2);
plot(col1, col2)
function y=poly2dfit.sci
Wie kann lässt sich ein polynomial-fit ausüben und diesen dann ableiten?
|
|
|
|
|
|
| xquadxq |
Themenstarter
Forum-Newbie
|
|
Beiträge: 8
|
|
|
|
Anmeldedatum: 13.03.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 04.05.2020, 22:21
Titel:
|
|
Habe etwas gefunden was funktioniert:
x = [-3 -1 0 1 3 5 7];
y = [50 74 62 40 19 35 52];
plot(x,y,"."); // plot sample points only
pcoeff = polyfit(x,y,4); // calculate polynominal coefficients (4th-degree)
xp = linspace(-3,7,100); // generate a little more x-values for a smoother curve fitting
yp = polyval(pcoeff,xp); // calculate the y-values for the curve fitting
plot(xp, yp,"k"); // plot the curve fitting in black
Habe polyfit und polyval installiert.
Wie kann ich nun die Vektoren x resp. y aus einer Matrix aulesen?
Bzw. wie erhalte ich diese aus einer zweispaltigen .csv-Datei?
|
|
|
|
| xquadxq |
Themenstarter
Forum-Newbie
|
|
Beiträge: 8
|
|
|
|
Anmeldedatum: 13.03.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 04.05.2020, 22:24
Titel:
|
|
so hier???
clear
Array=csvRead('C:\Users\Karl\Desktop\SCILAB\D_Werte.csv');
col1 = Array(:, 1);
col2 = Array(:, 2);
plot(x,y,"."); // plot sample points only
pcoeff = polyfit(x,y,4); // calculate polynominal coefficients (4th-degree)
xp = linspace(-3,7,100); // generate a little more x-values for a smoother curve fitting
yp = polyval(pcoeff,xp); // calculate the y-values for the curve fitting
plot(xp, yp,"k"); // plot the curve fitting in black
|
|
|
|
| xquadxq |
Themenstarter
Forum-Newbie
|
|
Beiträge: 8
|
|
|
|
Anmeldedatum: 13.03.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 04.05.2020, 22:38
Titel:
|
|
Also nochmal die Frage:
Wie kann ich aus der .csv-Datei spalte 1 als x und spalte 2 als y deklarieren?
|
|
|
|
| xquadxq |
Themenstarter
Forum-Newbie
|
|
Beiträge: 8
|
|
|
|
Anmeldedatum: 13.03.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 05.05.2020, 07:09
Titel:
|
|
clear
x_lst = [0 1 2 3 4 5]
y_lst = [2.1 7.7 13.6 27.2 40.9 61.1]
p=polyfit(x_lst, y_lst, 3)
f=polyval(p,x_lst)
scf()
plot(x_lst,f,"b-")
plot(x_lst,y_lst,"b+")
Wie erhalte ich das Polynom und wie die Ableitung davon?
|
|
|
|
| xquadxq |
Themenstarter
Forum-Newbie
|
|
Beiträge: 8
|
|
|
|
Anmeldedatum: 13.03.19
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 05.05.2020, 07:20
Titel:
|
|
clear
M=csvRead('e:\Users\Karl\Desktop\Simulation\test1.csv')
col1 = M(:, 1)
col2 = M(:, 2)
p=polyfit(col1, col2, 4)
f=polyval(p,col1)
scf()
plot(col1,f,"b-")
plot(col1,col2,"b+")
disp(f)
--> exec('E:\Users\Karl\Desktop\scilab\test1.sce', -1)
at line 68 of function polyfit ( C:\Users\Karl\AppData\Roaming\Scilab\scilab-6.1.0\matlab_equvivalent_fcn.sci line 358 )
at line 5 of executed file E:\Users\Karl\Desktop\scilab\test1.sce
With NaN or Inf a left division by scalar expected.
... was ist hier falsch???
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
|  RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|
