关于坐标转换的问题

zhenmeizi

想请教一下高人：
这个循环怎么表述啊，这样算出来时矩阵，还不匹配，其实设计计算的是一些点，就像最后一项用它来求几何特征的，b只是一个边长，但和前一个角度有关系，不知道怎么表述，拜托帮一下忙，谢谢了！
   phi(n)=2*pi*(rand(1));
   theta(n)=acos(((1+g^2)/2*g)-((1-g^2)^2/2*g*(1-g+2*g*(rand(1)))^2));
   L(n)=log(rand(1));
   u(n)=sin(theta(n))*sin(phi(n));
   x(n)=x(n-1)+L(n)*u(n)
   b(n)=sin(theta(n-1))*d(n-1)/sin(pi-theta(n));

lxq

就按你写出的式子
先给定theta(1)=...
b(1)=...
for n=2:N
    theta(n)=acos(((1+g^2)/2*g)-((1-g^2)^2/2*g*(1-g+2*g*(rand(1)))^2));
    b(n)=sin(theta(n-1))*d(n-1)/sin(pi-theta(n));
end
是不是这个意思？

zhenmeizi

对啊，但这样算出的是矩阵，我要得到的是具体的数，该怎样表述啊，谢谢了

lxq

b(i)不就是得到具体的数？