k=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2),3^(1/2)];%射线的斜率，前17条射线的斜率为0，后17条射线的斜率为 。b=[-80,-70,-60,-50,-40,-30,-20,-10,0,10,20,30,40,50,60,70,80,-160,-140,-120,-100 ,-80,-60,-40,-20,0,20,40,60,80,100,120,140,160];%34条射线的截距，以上为所需的参数，这些参数可以根据射线条数和射线的间距以及两次测量的射线的角度计算得出，现在为了计算更加的直观，减少变量的数量。先以数据计算,程序编好以后再修改增加参数的计算syms xfor i=1:34a(1,i)=1a(2,i)=xa(3,i) =k(i)*x+b(i)a(4,i) =x^2a(5,i) =x*( k(i)*x+b(i))a(6,i) =( k(i)*x+b(i))^2a(7,i) =x^3a(8,i) =x^2*( k(i)*x+b(i))a(9,i) =x*( k(i)*x+b(i))^2a(10,i) =( k(i)*x+b(i))^3a(11,i) =x^4'a(12,i) =x^3*( k(i)*x+b(i))a(14,i) =x^2*( k(i)*x+b(i))^2a(15,i) =x*( k(i)*x+b(i))^3a(16,i) =( k(i)*x+b(i))^4a(17,i) =x^5a(18,i) =x^4*( k(i)*x+b(i))a(19,i) =x^3*( k(i)*x+b(i))^2a(20,i)=x^2*( k(i)*x+b(i))^3a(21,i) =x*( k(i)*x+b(i))^4a(22,i) =( k(i)*x+b(i))^5a(23,i) =x^6a(24,i) =x^5*( k(i)*x+b(i))a(25,i) =x^4*( k(i)*x+b(i))^2a(26,i) =x^3*( k(i)*x+b(i))^3a(27,i) =x^2*( k(i)*x+b(i))^4a(28,i) =x*( k(i)*x+b(i))^5a(29,i) =( k(i)*x+b(i))^6endfor i=1:34for j=1:29A(j,i)=int(a(j,i),X2(i),X1(i))/(X2(i)-X1(i))endend
for i=1:34;X1(i)=-(sqrt(R^2*(1+k(i)^2)-b(i)^2)+k(i)*b(i))/(1+k(i)^2);% X2(i)=(sqrt(R^2*(1+k(i)^2)-b(i)^2)-k(i)*b(i))/(1+k(i)^2);% end;R=130对于X1（i）和X2（i）以及Rs
|
