坐标转换

matlab计算过程

>> y1=[1.5035820007324219 2.1732330322265625 1.54083251953125 2.4599075317382812 1.8644142150878904 2.3998260498046875 1.8347167968750002 2.2776031494140625 2.7443504333496094 2.9309463500976562]';

>> y2=[1.2998800834337652 1.1871251376595038 1.6065399343116196 1.5844042017650257 1.877124587985416 1.8102113116972178 2.1531566635242307 2.3963804941359395 2.548845329495116 2.232749202125121]';

>> x1=[114.20813262462616 114.21124398708342 114.20827478170395 114.21253681182861 114.20974999666214 114.21223908662796 114.20958906412125 114.21162217855453 114.21377062797546 114.21466380357742]';

>> x2=[22.684188520923495 22.683531479720816 22.6860594006558 22.685960412949996 22.6877223834242 22.687336335651327 22.689405143066505 22.690899812244602 22.691855006194086 22.689932238528215]';

>> x=[ones(10,1),x1,x2];

>> [b,bint,r,rint,s]=regress(y1,x,0.0000001);

>> [b1,bint1,r1,rint1,s1]=regress(y2,x,0.0000001);

>> b=vpa(b,30);

>> b1=vpa(b1,30);

>> b

b =

-24707.2949757538917765486985445

  215.726452049523999221491976641

   3.1336279148976076136534629768

>> b1

b1 =

-3554.63217025167932661133818328

-1.45023396372361923489791024622

  164.059668462394597554521169513

>> y=[ones(10,1),y1,y2];

>> [c,cint,cr,cs]=regress(x1,y,0.0000001);

>> [c1,cint1,cr1,cs1]=regress(x2,y,0.0000001);

>> c=vpa(c,30);

>> c1=vpa(c1,30);

>> c

c =

       114.201276439319201472244458273

    0.00463490282764736881426781422988

-0.0000885273324555299037565830877483

>> c1

c1 =

      22.6762054032649906787355575943

0.0000409743154115068265296722982427

   0.00609455561743045635264559578559

>>