QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 14945|回复: 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 编辑
    ; g7 c# [, t2 D5 n4 m  u3 p# B  g
    ; ]* ~5 ]7 _4 o1 S不知道有木有人研究这个,比较偏向于硬件描述。2 _1 m/ {& d5 Y1 s; I9 O
    今天先发一个sine cose函数的实现。! ^8 A5 Z' w2 F' d& I9 b+ \9 p; j
    有兴趣的可以看一下,m文件见附件。' `" \$ J) O" z6 [( F5 F/ R
    function v=f_cordic(beta,n)) K0 A0 T# Y& y; `* V. [
    if beta < -pi/2 || beta > pi/2  P4 J; u( _( h% ~$ Z# Q4 ^) j
        if beta < 0
    1 c% ~* ~# k( J, i4 z& s        v=f_cordic(beta + pi,n);
    6 D* v0 X6 O( T2 z) g$ S$ B, f2 E( q9 T    else
    : h, d( s8 C9 ~, M, q* q5 ]        v=f_cordic(beta - pi,n);
    8 E+ B: s1 ?. X+ J( v( `+ a7 R    end: g, U) I7 M+ d& u5 d1 M+ P
        v = -v;% flip the sign for second or third quadrant
    , {2 N( a+ |$ [    return$ H" F3 d% U' [8 P2 g, P( }8 G2 v7 D
    end
    " \7 H8 w  @' T0 L8 D% ]% Initialization of tables of constants used by CORDIC4 u  N' S. q4 b$ v* `5 J& y
    % need a table of arctangents of negative powers of two, in radians:1 V" G& W% C1 u  J9 G. v
    % angles = atan(2.^-(0:27));" t; F# g: b  y
    angles =  [  ...) n  A" l9 N9 f; R3 d4 \- F
        0.78539816339745   0.46364760900081   0.24497866312686   0.12435499454676 ...
    5 ^$ u4 _+ l$ Q  u" i8 x4 |. ]% t    0.06241880999596   0.03123983343027   0.01562372862048   0.00781234106010 ...
    7 {4 l' I  f% `) j, t8 ^    0.00390623013197   0.00195312251648   0.00097656218956   0.00048828121119 ...4 X" Y/ K* S6 b! w& m
        0.00024414062015   0.00012207031189   0.00006103515617   0.00003051757812 ...2 t' {* [0 `# x) f* `! k( f  u
        0.00001525878906   0.00000762939453   0.00000381469727   0.00000190734863 ...5 c; o2 K7 J$ [
        0.00000095367432   0.00000047683716   0.00000023841858   0.00000011920929 ...
    4 i  ]9 @/ _2 r    0.00000005960464   0.00000002980232   0.00000001490116   0.00000000745058 ];  p% C% f8 ]' v
    % and a table of products of reciprocal lengths of vectors [1, 2^-j]:0 ]1 V! I! S' T5 s
    Kvalues = [ ...( S2 |7 T2 P0 x- v# q" M
        0.70710678118655   0.63245553203368   0.61357199107790   0.60883391251775 ...' J! B. n4 V7 X) \
        0.60764825625617   0.60735177014130   0.60727764409353   0.60725911229889 ...8 ]. s3 f7 T  I3 _% C1 k! M8 V
        0.60725447933256   0.60725332108988   0.60725303152913   0.60725295913894 ...9 U& ^2 r9 b; I7 x+ k* R
        0.60725294104140   0.60725293651701   0.60725293538591   0.60725293510314 ...
    ) ^+ D/ q8 e' J4 I( m! i    0.60725293503245   0.60725293501477   0.60725293501035   0.60725293500925 ...+ t' S  U  F0 r3 @  f4 T* o
        0.60725293500897   0.60725293500890   0.60725293500889   0.60725293500888 ];
    9 j& y3 n; o) g9 Y; T$ {Kn = Kvalues(min(n, length(Kvalues)));  I6 A9 J4 y. j8 w0 H0 [" h; o
    / Y: X' I+ e7 w7 e
    % Initialize loop variables:4 J% h, v; z% x7 ^0 _
    v = [1;0]; % start with 2-vector cosine and sine of zero1 `( V/ _1 r$ F0 _9 l% O
    poweroftwo = 1; 7 @2 m7 X0 J: [
    % Iterations/ ?8 ]+ B) I5 O! {
    for j = 1:n;
    8 r3 a0 i# K$ g$ n    if beta < 0" R( f: B  a! c
            sigma = -1;9 Q, c" V: t6 w. e# g
        else- {$ K2 _: D) V% Q
            sigma = 1;" E( S, Q- W4 a4 L: K
        end
    " l% e- \( [+ r9 x& {    % update the angle from table, or eventually by just dividing by two4 i/ N" C. ~' l* \# d* K
        if j <= length(angles)
    # |7 T6 w& m+ t1 d! B( m7 Q        angle = angles(j);8 c5 N  A2 [9 Z% _/ _1 K3 t
        else
    5 e/ g4 r$ R/ S( G9 n) e# @$ C; I        angle = angle/2;
    8 L' K; y" M+ x    end0 v" P! I! q; b3 q8 }
        factor = sigma * poweroftwo;) _+ T! N- f& a2 T4 a9 n
        R = [1, -factor; factor, 1];  ~9 B/ Y( a8 T" O. l
        v = R * v; % 2-by-2 mtrix multiply
      T, b* O' k. a    beta = beta - sigma * angle; % update the remaining angle4 I. P3 v. t8 g& b; F, s" _4 W
        poweroftwo = poweroftwo / 2;
    ! V/ H" ]- X, W+ _5 K, f1 Iend
    5 H2 L) Q6 O+ T$ L7 @: e+ E& t% Adjust length of output vector to be [cos(beta), sin(beta)]:
    ) ~, N6 x8 ~, m; R% T( o9 M- u& @v = v * Kn;
    ( o' W0 G8 ?/ T: creturn3 e# N9 K. u0 y4 e

    : e  a) u+ {/ e' {% Z使用20位寄存器,仿真波形如下:; W' r/ S4 m( f( A5 k$ `* v
    仿真波形 2 A$ p) |) ]2 D! e6 n

    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 ( W0 x* g! ]. z2 \2 I
    这是什么算法?第一次听说

    ' \" N7 t  e5 I$ H& g* k这个,就是那个。感兴趣的话你还是自己搜搜吧。
    我正在用金山快盘。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 9 v3 F4 c& u4 c$ o) a
    我表示很给力

    1 g) a1 I& R) j: ^+ N9 b给力就好。呵呵
    我正在用金山快盘。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的每日心情
    奋斗
    2025-12-13 17:58
  • 签到天数: 3589 天

    [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-4-9 20:16 , Processed in 0.795287 second(s), 106 queries .

    回顶部