Peter Berthold

32 Reputation

2 Badges

17 years, 205 days

MaplePrimes Activity


These are answers submitted by Peter Berthold

Sorry, i wasn't aware that you could upload here, and i opened this thread a bit premature and have therefore moved it to here: http://www.mapleprimes.com/forum/help-needed-for-3d-plot-fitting Liblip seems to be exactly what i need, even though i have some problems with the programming, but thank you for the advice. Will try to make it work, and will post in the thread specified above from now on.
Didn't know you could upload here, but here it is: Download: http://www.mapleprimes.com/files/4416_All_Math%20dep%20FORUM.mw Full view: http://maplenet.maplesoft.com/maplenet/worksheet/mapleprimes/4416_All_Math%20dep%20FORUM.mw Thanks again to everybody that wants to help! Main view:

>

My results as 3d coordinates:

> Maple EquationPx:=pointplot3d([[-165.196852270254,0,899.3]....]],shading=zhue, symbol=point):

> Maple Equation

My attempt on fitting a function to the data:

> f:=y->y+200*sin((1/100)*y)+(1/2)*y*cos((1/100)*y)-(1/3)*y-((y/5000)*(y+200*sin((1/100)*y)+(1/2)*y*cos((1/100)*y)-(1/3)*y));

Maple Equation

> with(plots):

Warning, the name changecoords has been redefined

> plot(f,0..3000);

Maple Plot

> g:=(x,y)->y+100*sin((1/100)*y)+(1/2)*y*cos((1/100)*y)-(1/3)*y-(1/70)*x^2-((0.75)*10^(-5))*y*x^2-((y/5000)*(y+200*sin((1/100)*y)+(1/2)*y*cos((1/100)*y)-(1/3)*y))+1000;

Maple Equation

> J:=plot3d(g(x,y), x=-300..300, y=0..3000, axes=boxed, grid=[80,80], style=wireframe ,color=gray):

> Eo:=plot3d(0,x=-300..350, y=0..3000,color=black, axes=boxed):

> Ey:=plot3d([x,3000,z],x=-300..350,z=0..3500, style=wireframe, color=gray):

> Ex:=plot3d([-300,y,z],y=0..3000,z=0..3500,color=black):

> display({Px,Eo,Ey,Ex},view=[-350..380,0..3000,0..3500],axes=normal);

Maple Plot

> display({J,Px},view=[-400..400,0..3000,0..3500]);

Maple Plot

> display({Px,Eo,Ey,Ex},view=[-300..300,1300..2200,0..3500]);

Maple Plot

> Restart;

Maple Equation

> with(plots):

> c2989:=pointplot3d([[-219.146070008111,2801.123,2000], [-170.205757834452,2801.123,2249.6], [-141.760361173355,2801.123,2005.0], [-69.0796641566822,2801.123,2511.5], [-42.8018691180654,2801.123,2733.5], [-24.0416305603426,2801.123,2376.2], [29.0688837074973,2801.123,2736.6], [55.9016994374947,2801.123,2478.9], [88.6002257333467,2801.123,2940.2], [133.150291024842,2801.123,2857.4], [269.066906177627,2801.123,1904.4], [169.59363195592,2801.123,2705.4], [207.060377667964,2801.123,2510.6], [-286.048946860498,2801.123,1580.9], [324.237567225021,2801.123,1381]], axes=boxed, shading=zhue, symbol=point):

> display({c2989});

Maple Plot

Finding Least-Squares Fit Curve and Plotting in 3D

Modelling to one time frame at a time

http://www.maplesoft.com/applications/app_center_view.aspx?AID=978&CID=10&SCID=148

> with(stats):

>

> Kx:=[-219.146070008111,-170.205757834452,-141.760361173355,-69.0796641566822,-42.8018691180654,-24.0416305603426,29.0688837074973,55.9016994374947,88.6002257333467,133.150291024842,269.066906177627,169.59363195592,207.060377667964,-286.048946860498,324.237567225021]:

> Ky:=[2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123,2801.123]:

> Kz:=[2000,2249.6,2005.0,2511.5,2733.5,2376.2,2736.6,2478.9,2940.2,2857.4,1904.4,2705.4,2510.6,1580.9,1381]:

> leastsquare[[x,y,z]]([Kx,Ky,Kz]):

> fit[leastsquare[[x,y,z], y=a*x^2+b*x+c]]([Kx,Ky,Kz]);

Maple Equation

> s:=x->0.3729830443e-10*x^2-0.4041446601e-8*x+2801.122999;

Maple Equation

> S:=plot3d(s(x), x=-300..300, y=2800..2900, axes=boxed, style=wireframe , color=gray):

> display({c2989,S}, scaling=constrained, labeldirections=[horizontal,horizontal,horizontal],tickmarks=[3,3,3]);

Maple Plot

> fit[leastsquare[[x,y], y=a*x^2+b*x+c]]( [[-219.146070008111,2801.123,2000], [-170.205757834452,2801.123,2249.6], [-141.760361173355,2801.123,2005.0], [-69.0796641566822,2801.123,2511.5], [-42.8018691180654,2801.123,2733.5], [-24.0416305603426,2801.123,2376.2], [29.0688837074973,2801.123,2736.6], [55.9016994374947,2801.123,2478.9], [88.6002257333467,2801.123,2940.2], [133.150291024842,2801.123,2857.4], [269.066906177627,2801.123,1904.4], [169.59363195592,2801.123,2705.4], [207.060377667964,2801.123,2510.6], [-286.048946860498,2801.123,1580.9], [324.237567225021,2801.123,1381]]);

Maple Equation

> h:=x->-0.1330956511e-3*x^2+1.327445943*x+127.0907208;

Maple Equation

> H:=plot3d(h(x), x=-300..320, y=2800..2810, axes=boxed, style=wireframe ,shading=zhue):

> display({c2989,H});

Maple Plot

>

This post was generated using the MaplePrimes File Manager

View 4416_All_Math dep FORUM.mw on MapleNet or Download 4416_All_Math dep FORUM.mw
View file details

Link to data: http://www.megaupload.com/?d=JK70IPS5 Thank you to everybody who are willing just to have a look at it.
Page 1 of 1