关于坐标转换的问题
想请教一下高人:这个循环怎么表述啊,这样算出来时矩阵,还不匹配,其实设计计算的是一些点,就像最后一项用它来求几何特征的,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)); 就按你写出的式子
先给定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
是不是这个意思? 对啊,但这样算出的是矩阵,我要得到的是具体的数,该怎样表述啊,谢谢了
回复 #3 zhenmeizi 的帖子
b(i)不就是得到具体的数?
页:
[1]