clear;close all;clc; %Working condition intializing
m=-1; % hyperbolic coordinate system
u=1; % rotating direction
K=0.82816; %修正因子
seq=[1 2 3 4 4 5 6 7 8 9 10 11 12 13 13 14 15];%旋转序列
x0=1/K;
y0=0;
z0=0.2;
for l=1:length(seq)
x=x0-m*u*y0*2^(-(seq(l)));
y=y0+u*x0*2^(-(seq(l)));
z=z0-u*atanh(2^(-(seq(l))));
x0=x;
y0=y;
z0=z;
u=sign(z0);
end
x0
y0
z0
文章评论(0条评论)
登录后参与讨论