zhenmeizi 发表于 2007-7-6 17:43

关于坐标转换的问题

想请教一下高人:
这个循环怎么表述啊,这样算出来时矩阵,还不匹配,其实设计计算的是一些点,就像最后一项用它来求几何特征的,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 发表于 2007-7-6 17:53

就按你写出的式子
先给定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 发表于 2007-7-6 18:10

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

lxq 发表于 2007-7-6 19:28

回复 #3 zhenmeizi 的帖子

b(i)不就是得到具体的数?
页: [1]
查看完整版本: 关于坐标转换的问题