螺旋线拟合 |
您所在的位置:网站首页 › 螺线图像和方程 › 螺旋线拟合 |
实验坐标是这样的 基本是一个阿基米德螺旋线,原数据是笛卡尔坐标系的,通过拟合半径与方位角之间的关系得到多项式 以下是多项式拟合 clc clear format long shuju=xlsread('jingpan12') x=flipud(shuju(:,1)); y=flipud(shuju(:,2)); r=sqrt(x.^2+y.^2); jiaodu=acosd(x./r); for i=1:1:19 jiaodu(i)=360-jiaodu(i); end for i=20:1:73 jiaodu(i)=360+jiaodu(i); end for i=74:1:147 jiaodu(i)=720-jiaodu(i); end for i=148:1:254 jiaodu(i)=720+jiaodu(i); end for i=255:1:257 jiaodu(i)=1080-jiaodu(i); end jiaodu=jiaodu*3.14./180; plot(jiaodu,r); p9=polyfit(jiaodu,r,9); n=3.2:0.1:16; y9=polyval(p9,n,'.'); hold on plot(n,y9)
将数据回代, clear all clc format long shuju=xlsread('jingpan12') p=[-5.00630375898308e-06,0.000460742305865134,-0.0182267190094672,0.404815073833663,-5.53480697324230,48.0789047375280,-264.371215185261,886.259566745230,-1645.21576472083,1297.95181501538] for n=3.2:0.1:16; p9=polyval(p,n); c=p9*cos(n); d=p9*sin(n); plot(c,d,'.') hold on end plot(shuju(:,1),shuju(:,2),'.')回代结果 实线是元数据,点是拟合结果
|
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |