【算法】Cordic算法的原理及matlab/verilog应用

一、前言

单片机或者FPGA等计算能力弱的嵌入式设备进行加减运算还是容易实现,但是想要计算三角函数(sin、cos、tan),甚至双曲线、指数、对数这样复杂的函数,那就需要费些力了。通常这些函数的计算需要通者查找表或近似计算(如泰勒级数逼近)等技术来转换为硬件易于实现的方式。

CORDIC(Coordinate Rotation Digital Computer, 坐标旋转数字计算方法)算法就是一种化繁为简的算法,通过基本的加减和移位运算代替乘法运算,逐渐逼近目标值,得出函数的数值解。

二、Cordic算法理论推导

理论推导参考:CORDIC算法理论详解_cordic算法详解-CSDN博客,这篇博客的推导仔细而全面。

Cordic算法的基石在于一个规律:从tan45°开始,角度每减半,tan值也接近减半。这一规律直接将三角函数运算变成2的幂运算,而这在数字电路中可直接用移位运算来实现。

三、Cordic算法 matlab实现

3.1 已知相位(角度)求坐标(正弦余弦)

Matlab 复制代码
function [sin_theta,cos_theta] = cordic_sincos(theta,n)
% n:iterations
% theta: -180~180
tan_table = 2.^-(0 : n-1);
angle_rad_lut = atan(tan_table);
%angle_deg_lut = rad2deg( atan(tan_table) );
 
k = 1;
for i = 0 : n-1
    k = k*(1/sqrt(1 + 2^(-2*i)));
end
 

x = k;
y = 0;
theta_tar = theta*pi/180;  % to be rad
z=theta_tar;
 
% preprocess
if (theta_tar > pi/2)  
    theta_tar = theta_tar - pi/2;
elseif (theta_tar < -pi/2)
    theta_tar = theta_tar + pi/2;
end
     
 for i = 0 : n-1   
         if (z > 0)  
             d =1;
         else
             d=-1;
         end
        x_temp = x;
        y_temp =y;
        z_temp = z;
        x= x_temp - d*y_temp*2^(-i);
        y = y_temp + d*x_temp*2^(-i);
        z = z_temp - d*angle_rad_lut(i+1);
        
end
    
if (theta_tar > pi/2)  
    sin_theta = -x;
    cos_theta = y;
elseif (theta_tar < -pi/2)
    sin_theta = x;
    cos_theta = -y;
else
    sin_theta = y;
    cos_theta = x;
end
end

3.2 已知坐标求相位

待实现

四、cordic算法的verilog实现

待实现

相关推荐
Swift社区9 小时前
LeetCode 394. 字符串解码(Decode String)
算法·leetcode·职场和发展
tt55555555555510 小时前
LeetCode进阶算法题解详解
算法·leetcode·职场和发展
让我们一起加油好吗10 小时前
【基础算法】DFS中的剪枝与优化
算法·深度优先·剪枝
Q741_14710 小时前
C++ 模拟题 力扣495. 提莫攻击 题解 每日一题
c++·算法·leetcode·模拟
璞致电子11 小时前
fpga开发板ZYNQ 璞致 PZ7010/7020 邮票孔核心板简介-ZYNQ7000系列小系统学习板
linux·嵌入式硬件·学习·fpga开发·fpga·fpga开发板·xilinx开发板
Felven11 小时前
A. Be Positive
算法
小O的算法实验室11 小时前
2026年COR SCI2区,自适应K-means和强化学习RL算法+有效疫苗分配问题,深度解析+性能实测,深度解析+性能实测
算法·论文复现·智能算法·智能算法改进
青岛少儿编程-王老师12 小时前
CCF编程能力等级认证GESP—C++7级—20250927
数据结构·c++·算法
夏鹏今天学习了吗12 小时前
【LeetCode热题100(39/100)】对称二叉树
算法·leetcode·职场和发展
鲸大鱼的自我修养13 小时前
Matlab/Simulink技巧
matlab