Prog--oktave

I_UfE_Ag_oktave

octave:1> load ixi.data //in matrix laden
octave:4> plot(ixi(:,2),ixi(:,1),"r@"); //matrix plotten 2 spalte als x 1 spalte als y
octave:5> p,s=polyfit(ixi(:,2),ixi(:,1),1); //polynom 1 ordnung (gerade anfitten)
octave:6> p //steigung und offset ausgeben
p = 2.0479   192.6550
octave:7> hold on //nächsten plot auf den letzten plotten
octave:9> t=polyval(p,0:15000); //gerade aus p in matrix berechnen
octave:10> plot(t); // gerade plotten
octave:11> diff=ixi(:,1)-s.yf; 
octave:12> min(diff)
ans = -160.93
octave:13> max(diff)
ans =  281.70
octave:14> 281/30000*100
ans =  0.93667
octave:14>print bildname.png // druckt aktuelles bild als png
octave:14>plot(ixi(:,2),diff);
octave:14>plot(ixi(:,2),diff/30000*100);

load ixi.data
plot(ixi(:,2),ixi(:,1),"r@");
p,s=polyfit(ixi(:,2),ixi(:,1),1); //polynom 1 ordnung (gerade anfitten)
hold on
t=polyval(p,0:15000);
plot(t);
diff=ixi(:,1)-s.yf; 
print gerade.png
print gerade.svg
hold off
plot(ixi(:,2),diff/30000*100);
print abweichung.png
print abweichung.svg
\n\n