Matlab - ALGORITMO CORDIC

   
Vista:

ALGORITMO CORDIC

Publicado por loubna bobolouna@hotmail.com (1 intervención) el 05/06/2014 16:01:50
hola:
tengo que programar el algoritmo cordic para seno y coseno hiperbólicos. el tema es que el código función solo el 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
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder