QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 15220|回复: 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 编辑
    % C8 _% a, _" Y3 x3 S8 h2 v: H) F1 N% p. U7 x* a
    不知道有木有人研究这个,比较偏向于硬件描述。" w. C3 ?8 H% b3 G
    今天先发一个sine cose函数的实现。; m* [& q! r6 T3 ^' |# |" W  j  H
    有兴趣的可以看一下,m文件见附件。
      X* V4 K% E  _7 E8 efunction v=f_cordic(beta,n). @# b; L) D: h8 W
    if beta < -pi/2 || beta > pi/29 Y9 F. L: S. K
        if beta < 05 Z9 ~. x) ]6 t2 z( b
            v=f_cordic(beta + pi,n);
      h; n# G$ V4 R4 s7 L4 b0 {! q) L    else: _! S( B6 K3 a5 {. C
            v=f_cordic(beta - pi,n);" O# l, v3 E' |! z. |, o# Y' d
        end
    ' J- J* {- e+ M% z    v = -v;% flip the sign for second or third quadrant" w$ C1 a) p6 B( M5 n  h1 @
        return
    : e' N2 G0 W) j3 g2 Gend
    7 a; S- [4 `$ A) W% Initialization of tables of constants used by CORDIC
    1 y3 `+ A, M' t; Z: ^, Z% need a table of arctangents of negative powers of two, in radians:" {# t; k6 I1 j" s$ L
    % angles = atan(2.^-(0:27));
    & B! J: q) r0 E: uangles =  [  ...$ h; q/ i. K/ y5 L
        0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...7 e1 y8 F' ^/ w+ I* J; W: s0 J& Z# Q
        0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    ( q  ~& S7 |9 l( Q) W    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...
    : p" K: ~; t- [# n    0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...
    % k! R% _6 W2 ~! ~8 l    0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...
    2 _8 m- E# d1 ~    0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...$ u  u% z( [3 y
        0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];
    5 l2 N3 \! O! \! \9 e3 z% and a table of products of reciprocal lengths of vectors [1, 2^-j]:" |$ I) Z0 K0 ^+ n
    Kvalues = [ ...7 ~  w- O$ v* Z( E) C
        0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...
    / p, v: r  J) f" d: `    0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...% T' |# j! j: \- [
        0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...9 U9 m) t( U2 t' I7 a# f7 h
        0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...) V  c  J- A9 f
        0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...
    9 G  r0 I/ d% r    0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];' B: ]- ?5 k5 [3 o- d/ `8 J9 r
    Kn = Kvalues(min(n, length(Kvalues)));6 v4 y# j# e1 R( r( X  {5 v

    * l' u0 Y9 u" d% Initialize loop variables:+ z: u" D, f' o5 y0 s7 H& B
    v = [1;0]; % start with 2-vector cosine and sine of zero0 j- a1 T0 ?+ q. D+ {3 u3 v
    poweroftwo = 1;
    $ c, B" {  F. l% l% Iterations
    : c3 k2 _: V  x7 w* Z+ ?+ pfor j = 1:n;3 m- q/ t$ I% |8 h
        if beta < 0
    & c0 F" s& A. U, T% e  z6 B3 r* z        sigma = -1;8 s; J; o* A6 m4 \9 K
        else
    # S& t  F0 g4 s" \        sigma = 1;$ v  U* D" q! x
        end
    5 W; c8 X. \, ~; M( a6 T4 g3 Q    % update the angle from table, or eventually by just dividing by two. t' ^! M6 M4 v& ?3 N1 _' o. {
        if j <= length(angles)
    * c2 p, t1 o2 L8 q4 g$ Z! ^        angle = angles(j);: B1 k2 U; w8 P0 `( `
        else
    5 {1 z9 U: z7 R9 ]+ _2 ]0 V        angle = angle/2;$ ~. o4 s" t  R! z1 _, ]! D
        end' i* N# @: r" C7 {
        factor = sigma * poweroftwo;
    * C% g: e: }+ O) k' d" B9 q1 _$ E    R = [1, -factor; factor, 1];$ s" i8 c% c  Y' Y$ a% q
        v = R * v; % 2-by-2 mtrix multiply- d: ?) n1 F! L+ }) x
        beta = beta - sigma * angle; % update the remaining angle* b+ L: X) [: ]& K( a7 T- {
        poweroftwo = poweroftwo / 2;/ p* a  _) D/ @6 U5 E9 X
    end
      Q; G! i3 [% ^/ \% Adjust length of output vector to be [cos(beta), sin(beta)]:$ Z) x5 |8 U+ v9 `
    v = v * Kn;
    " d" s3 I6 E+ P) P' o2 O: X2 hreturn; I9 ^# o, _2 k1 L) P: N; K

    + l' N5 y& S* R. W7 K6 W% ]" y) ]使用20位寄存器,仿真波形如下:: `* R9 z( p# W1 H
    仿真波形 4 r& _5 ?% F6 t+ |7 c8 r

    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 9 I0 y) H6 D: z8 f
    这是什么算法?第一次听说
    3 T- H: s: p$ Y7 B8 J# ]" J
    这个,就是那个。感兴趣的话你还是自己搜搜吧。
    我正在用金山快盘。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 5 D) ?9 v7 M# P% d$ o% n
    我表示很给力

    ; y5 d8 G0 t, A% `5 E  s8 ^给力就好。呵呵
    我正在用金山快盘。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 11:03 , Processed in 0.503290 second(s), 105 queries .

    回顶部