ALGORITMO CORDIC
Publicado por loubna (1 intervención) el 06/06/2014 13:47:32
hola:
tengo que programar el algoritmo cordic para seno y coseno hiperbólicos. el tema es que el código funciona solo en el rango de convergencia y cuando lo expando no . a ver si alguien me puede echar una mano . el código esta implementado en scilab
tengo que programar el algoritmo cordic para seno y coseno hiperbólicos. el tema es que el código funciona solo en el rango de convergencia y cuando lo expando no . a ver si alguien me puede echar una mano . el código esta implementado en scilab
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
function [xn,yn,zn,k]= cos(a)
k=1;k1=1;k2=1;y(1)=0;
z(1)=a;
l=4;//x1=0;
//y1=0;
//z1=0;
for j=-5:-1:0
k1 =k1*(sqrt(1-(1-(2^(2* j)/4))));
end
for j=1:1:100
k2=k2*(sqrt(1-(1/(2^(2* j)))));
end
k=k1*k2;
x(1)=1/k;
for i=-5:-1:0 // iteraciones adicionales
if z(i)<0
d=-1;
else
d=1;end
x(i+1)=x(i)+y(i)*(d*(1-(2^i/4)));
y(i+1)=y(i)+x(i)*(d*(1-(2^i/4)));
z(i+1)=z(i)-d*atanh((1-(2^i/4)));
x1=x(i+1);
y1=y(i+1);
z1=z(i+1);
end
i=1;x(1)=x1;
y(1)=y1;
z(1)=z1;
for i=1:1:100 // iteraciones normales
if z(i)<0
d=-1;
else
d=1;end
if i+1==l
x(i+1)=x(i)+y(i)*(d/2^(i));
y(i+1)=y(i)+x(i)*(d/2^(i));
z(i+1)=z(i)-d*atanh(1/2^(i));
l=3*l+1;
end
x(i+1)=x(i)+y(i)*(d/2^(i));
y(i+1)=y(i)+x(i)*(d/2^(i));
z(i+1)=z(i)-d*atanh(1/2^(i));
end
gracias :)
if z==0.0001
break;
end
plot(x,y,'-',x,y,'b-')
n=i;xn=x(i+1);
yn=y(i+1);
zn=z(i+1);
endfunction
Valora esta pregunta
0