QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 15213|回复: 30
打印 上一主题 下一主题

CORDIC算法的matlab实现

[复制链接]
字体大小: 正常 放大
大笨象 实名认证       

42

主题

11

听众

2119

积分

di_dar

  • TA的每日心情
    无聊
    2015-1-15 22:05
  • 签到天数: 79 天

    [LV.6]常住居民II

    自我介绍
    隐秘盛开

    优秀斑竹奖 新人进步奖 发帖功臣

    群组Matlab讨论组

    群组数学趣味、游戏、IQ等

    群组数学建模

    群组SIMULINK

    群组LINGO

    跳转到指定楼层
    1#
    发表于 2011-6-7 11:12 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    本帖最后由 大笨象 于 2011-6-7 12:26 编辑
    % N1 L6 e9 f# ?& [
    0 Z8 z0 `* [4 w" r+ A不知道有木有人研究这个,比较偏向于硬件描述。+ s# t* g% k/ p
    今天先发一个sine cose函数的实现。8 q) Y9 c, c3 W, U
    有兴趣的可以看一下,m文件见附件。6 s; H$ I3 L# V* i5 ~# }( ]
    function v=f_cordic(beta,n); N: p+ a) |' u1 [
    if beta < -pi/2 || beta > pi/2
    7 S4 P$ v- ~; i7 E% ^: t$ W7 j    if beta < 0
    ) @$ T; h" \, y  r        v=f_cordic(beta + pi,n);
    - e9 A; P3 O2 x' l7 N9 V1 N    else4 r/ R. h  m! n% I; m
            v=f_cordic(beta - pi,n);' T1 `, a/ B6 u( l) [
        end
    6 G2 `: I" t# \: P    v = -v;% flip the sign for second or third quadrant
    " |- p3 k& `: j1 B5 I3 w( L    return
    5 \( {# b' P- {9 u" Y. I4 n0 G, rend& {, Z& E% G  {8 H9 S0 f
    % Initialization of tables of constants used by CORDIC
    ) e+ B) G9 U. E% l5 }% need a table of arctangents of negative powers of two, in radians:
    ' S4 u8 u5 B8 Z9 I' ^; o$ ?; R% angles = atan(2.^-(0:27));% c' h  E8 m, ?. J$ c) `
    angles =  [  ...* Y5 T  s3 s: c8 G5 h8 h' c* g
        0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    / }9 X. N/ p* p& u/ S2 T+ ~    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    ( A( R5 W0 M0 P  y    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ..., K' `& J: }) j. x, u/ ]: @' v& M
        0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    & d: ]5 d, Z2 y    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...6 b* J, x2 _/ Y6 r
        0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    6 T7 t, q! Z9 V9 I; X    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
    2 }1 h  z& ?5 @/ J+ V3 P! \  l( T2 M% and a table of products of reciprocal lengths of vectors [1, 2^-j]:, J' {3 I, J6 G* d
    Kvalues = [ ...
    0 l4 `: V6 b- b* l    0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...# W0 ]( H( ~& _( h9 b  C2 k* r; U
        0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...  g+ n; q' g/ Z% W8 a8 G& A
        0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    1 g' m" X# O' D    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    3 [' m) s" L0 u) ~3 D7 ^9 n3 Y    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    / Q8 S( T, r6 e; Z# t* G* t0 `/ O    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
    : w& {1 R9 X7 L' wKn = Kvalues(min(n, length(Kvalues)));
    . l* B- i- V0 h4 |: i6 s6 Y! M
    6 l9 _4 M1 J% I) v% Initialize loop variables:
    4 P: v/ \( X2 p! |v = [1;0]; % start with 2-vector cosine and sine of zero
    % K1 U, O& O3 i# |) ^% g6 }poweroftwo = 1;
    6 a8 t1 J3 V# X- i% Iterations
    2 S3 h7 u& k& c4 zfor j = 1:n;
    0 n* s9 S* q: f) ~! F1 T$ u! G    if beta < 0
    / w8 C  M) q, c+ X& U5 l        sigma = -1;' c1 L7 x7 ?2 L0 S& f
        else. Q: A1 n" V; |: z+ V  x
            sigma = 1;& I$ X- P0 d4 ^3 ^" y& ?
        end
      [% E( v, Z# {7 f! [    % update the angle from table, or eventually by just dividing by two5 b* j3 A& W: A( X1 D8 y
        if j <= length(angles)* B& @7 q, K+ z
            angle = angles(j);
    & n7 [3 S. W6 e: R6 m3 H- L6 k    else
    . ^9 H0 Y5 W+ \5 d$ Y! G5 W3 `        angle = angle/2;
    / Y% W. I& H- t# N+ H    end
    9 `7 ?$ Q- _6 ?0 M! }    factor = sigma * poweroftwo;* i- S3 S1 z+ o/ _
        R = [1, -factor; factor, 1];- _! }3 l  r. t
        v = R * v; % 2-by-2 mtrix multiply5 r2 D0 a2 G) j
        beta = beta - sigma * angle; % update the remaining angle+ \; m4 \. k0 f) f
        poweroftwo = poweroftwo / 2;$ K+ V+ F2 k: E
    end$ h8 D! B& L2 N! o
    % Adjust length of output vector to be [cos(beta), sin(beta)]:
    7 f0 W  a, m# F$ U& Z& f, t  nv = v * Kn;" m3 \8 ]& U3 m0 `0 u# H
    return
    ) |- p/ m, \, M+ Y) a6 N
    8 x3 ~' [+ f/ D; @% V使用20位寄存器,仿真波形如下:* I  Q9 d  j, f( v
    仿真波形
    1 K- ?# l& M3 _2 E

    f_cordic.m

    2.25 KB, 下载次数: 21, 下载积分: 体力 -2 点

    CORDIC算法

    点评

    Ling619400  给力  发表于 2015-2-1 22:31
    zan
    已有 2 人评分体力 收起 理由
    sdccumcm + 5 很给力!
    厚积薄发 + 5 赞一个!

    总评分: 体力 + 10   查看全部评分

    转播转播0 分享淘帖0 分享分享1 收藏收藏2 支持支持3 反对反对0 微信微信
    我正在用金山快盘。http://www.kuaipan.cn/?channel=6ohsg7
    杨帆 实名认证       

    4

    主题

    4

    听众

    626

    积分

    升级  6.5%

  • TA的每日心情
    奋斗
    2012-12-4 13:38
  • 签到天数: 314 天

    [LV.8]以坛为家I

    新人进步奖 最具活力勋章 发帖功臣

    群组数学趣味、游戏、IQ等

    群组数学建模

    群组我们一定会赢

    回复

    使用道具 举报

    1

    主题

    3

    听众

    453

    积分

    升级  51%

  • TA的每日心情
    开心
    2015-7-7 09:53
  • 签到天数: 50 天

    [LV.5]常住居民I

    群组华南理工大学

    群组Matlab讨论组

    群组狂热数模爱好者

    群组数学建摸协会

    群组学术交流A

    回复

    使用道具 举报

    赛才        

    0

    主题

    2

    听众

    46

    积分

    升级  43.16%

    该用户从未签到

    回复

    使用道具 举报

    大笨象 实名认证       

    42

    主题

    11

    听众

    2119

    积分

    di_dar

  • TA的每日心情
    无聊
    2015-1-15 22:05
  • 签到天数: 79 天

    [LV.6]常住居民II

    自我介绍
    隐秘盛开

    优秀斑竹奖 新人进步奖 发帖功臣

    群组Matlab讨论组

    群组数学趣味、游戏、IQ等

    群组数学建模

    群组SIMULINK

    群组LINGO

    zhangjichang 发表于 2011-6-8 22:31
    5 Y: O1 j) ~* L8 Q5 [) ]  J+ k这是什么算法?第一次听说

    , J$ O+ ~  e5 b7 U( h2 A. O" u这个,就是那个。感兴趣的话你还是自己搜搜吧。
    我正在用金山快盘。http://www.kuaipan.cn/?channel=6ohsg7
    回复

    使用道具 举报

    大笨象 实名认证       

    42

    主题

    11

    听众

    2119

    积分

    di_dar

  • TA的每日心情
    无聊
    2015-1-15 22:05
  • 签到天数: 79 天

    [LV.6]常住居民II

    自我介绍
    隐秘盛开

    优秀斑竹奖 新人进步奖 发帖功臣

    群组Matlab讨论组

    群组数学趣味、游戏、IQ等

    群组数学建模

    群组SIMULINK

    群组LINGO

    赛才 发表于 2011-6-9 15:01
    0 J* b- O, @% N9 ^5 H1 Y3 S我表示很给力

    2 e5 p- V, C. w: S8 P) G给力就好。呵呵
    我正在用金山快盘。http://www.kuaipan.cn/?channel=6ohsg7
    回复

    使用道具 举报

    1

    主题

    2

    听众

    48

    积分

    升级  45.26%

  • TA的每日心情

    2011-7-1 22:07
  • 签到天数: 10 天

    [LV.3]偶尔看看II

    自我介绍
    最强的人往往是含泪奔跑的人!!!

    群组C 语言讨论组

    群组东北三省联盟

    群组数学专业考研加油站

    群组数学建模

    回复

    使用道具 举报

    yingzhen        

    1

    主题

    3

    听众

    297

    积分

    升级  98.5%

  • TA的每日心情

    2011-12-18 20:51
  • 签到天数: 30 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    羅雲琦        

    5

    主题

    2

    听众

    603

    积分

    升级  0.75%

  • TA的每日心情
    开心
    2014-11-2 16:11
  • 签到天数: 237 天

    [LV.7]常住居民III

    自我介绍
    我是來自四川省合江縣的男孩,樂觀自信

    社区QQ达人

    群组科技写作基础培训

    群组2011年第一期数学建模

    群组第四届cumcm国赛实训

    群组学术交流A

    群组学术交流B

    回复

    使用道具 举报

    jt202010 实名认证    中国数模人才认证  会长俱乐部认证 

    109

    主题

    165

    听众

    1万

    积分

    升级  0%

  • TA的每日心情
    奋斗
    2026-7-7 15:03
  • 签到天数: 3619 天

    [LV.Master]伴坛终老

    社区QQ达人 邮箱绑定达人 最具活力勋章 发帖功臣 风雨历程奖 新人进步奖

    群组数学建模

    群组自然数狂想曲

    群组2013年数学建模国赛备

    群组第三届数模基础实训

    群组第四届数学中国美赛实

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-7-9 08:53 , Processed in 1.163521 second(s), 108 queries .

    回顶部