QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 15215|回复: 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 编辑
    0 o+ B8 I& f: T/ f1 e2 c3 ?" R# F! y1 g/ X- G
    不知道有木有人研究这个,比较偏向于硬件描述。
    1 J( S0 v2 Z: S2 e; L今天先发一个sine cose函数的实现。
    1 @: R& K6 g1 _9 m有兴趣的可以看一下,m文件见附件。
    $ j# u* L+ X1 {. W4 ?$ rfunction v=f_cordic(beta,n)
    7 ?6 m9 w# i1 U# ?if beta < -pi/2 || beta > pi/22 d" l7 u+ z. a7 o
        if beta < 0
    " o- e2 v$ j% y6 X% m" D- ?        v=f_cordic(beta + pi,n);4 g8 m; L& s, P
        else
    4 m: n. |6 I" {* F" D6 g        v=f_cordic(beta - pi,n);
    + w, j- `, d; _& P    end5 \4 X. y: ~* @  Q7 B
        v = -v;% flip the sign for second or third quadrant6 B; _) B: i% @  u; P4 e- U
        return  H: s; z4 ]0 v0 K$ P6 M1 P
    end
    ' @) E1 X' }" o: T" T/ ~% Initialization of tables of constants used by CORDIC
    3 k2 V) g( R2 |" {7 F2 r. n* L$ O% need a table of arctangents of negative powers of two, in radians:
    : |8 j1 e# @/ g# `4 t% angles = atan(2.^-(0:27));
    - N) l+ z" b: @angles =  [  ...8 D/ ?- x; l- p2 B
        0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...) S) l: q. ^& I6 {6 o9 t" |* J
        0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    1 t6 z, X5 R! I# C  R    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...  ]! [+ ^* x; T$ K8 _
        0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...) V# O7 ~6 _  l0 }
        0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...* F5 y; L3 N5 P/ R
        0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    ! s3 j# ^, Z+ ^) `, r' z1 |( q( u' |    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];# M9 Z0 X2 ~8 i$ D3 L% S& t' x
    % and a table of products of reciprocal lengths of vectors [1, 2^-j]:" h  X6 c: l2 y- c2 z3 V2 Q
    Kvalues = [ ...3 ?2 [. ~7 x/ e- u1 G! b
        0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...3 K1 P  S& M, A9 R! d3 e! k
        0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...
    , Y7 u, n$ O( m% g3 |* g    0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...
    2 J  e9 x- c( H6 ?    0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 .../ m: O2 ?7 m0 ]+ H! E6 j' A
        0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    : v' E6 C" a4 O( [& M    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
    0 }# P' @6 ^' F& pKn = Kvalues(min(n, length(Kvalues)));
    2 c1 x' i# R* f+ ^- B/ j 8 V" y3 T4 ?0 l8 t
    % Initialize loop variables:
    4 T# w7 p. O" d, rv = [1;0]; % start with 2-vector cosine and sine of zero
    ( T) z6 L5 D. t  ^- xpoweroftwo = 1; . a' R4 z( A+ C
    % Iterations5 _7 Q! r5 s# P
    for j = 1:n;2 l' x. j' E3 _& P2 q8 H
        if beta < 0
    % T- w* l$ V, A# ?7 @        sigma = -1;
    ; o4 d, o% N' p- _/ u( O- l3 Q    else
    0 q* L, Y5 Z  ]* j        sigma = 1;3 m4 n+ [! }9 m# }5 g) J
        end' h$ C6 u$ m4 B! M7 Y7 g
        % update the angle from table, or eventually by just dividing by two
    & K; |; E: E, }    if j <= length(angles)
    4 e% z, H) u  D        angle = angles(j);
    " l, R$ [7 c  }    else8 W/ p2 b) B1 r. r. J, p; n) S) @
            angle = angle/2;
    % g" e4 Z" O9 x    end
    & F0 B" K0 j/ F    factor = sigma * poweroftwo;
    ) f5 {. B6 e& o0 A    R = [1, -factor; factor, 1];. v& D9 i! V! \! @/ @9 G
        v = R * v; % 2-by-2 mtrix multiply9 T+ @+ s! Q" e, y5 r4 ?3 e
        beta = beta - sigma * angle; % update the remaining angle8 U" ?# j1 ^- w$ Z, E
        poweroftwo = poweroftwo / 2;
    ( Z- I" v+ k6 K4 F' dend
    ' B2 u  N! p2 j# q1 {1 }5 _% Adjust length of output vector to be [cos(beta), sin(beta)]:
    ' ]* n* o7 o, ]0 e) k* \2 Uv = v * Kn;
    + m( e# ?) [! B& ]4 v5 breturn
    3 B, }1 T9 W$ M8 J# W
      c* M! G& U8 j9 [$ u使用20位寄存器,仿真波形如下:
    8 |+ h2 t; Y0 q% K$ s 仿真波形
    . q) J7 ~* Q6 M: e& F! f, `# [

    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 ; l  f0 ^) k0 l* a2 @9 e5 f
    这是什么算法?第一次听说
    & U- `" T% \% q  o  B
    这个,就是那个。感兴趣的话你还是自己搜搜吧。
    我正在用金山快盘。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 7 L3 b! U* r- v  M" B
    我表示很给力
    1 v0 J" W: B( t& m. |% _5 O
    给力就好。呵呵
    我正在用金山快盘。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 10:01 , Processed in 0.519452 second(s), 106 queries .

    回顶部