【算法】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实现

待实现

相关推荐
晨非辰13 小时前
【数据结构入坑指南】--《层序分明:堆的实现、排序与TOP-K问题一站式攻克(源码实战)》
c语言·开发语言·数据结构·算法·面试
hansang_IR13 小时前
【题解】P2217 [HAOI2007] 分割矩阵 [记忆化搜索]
c++·数学·算法·记忆化搜索·深搜
fl17683113 小时前
基于matlab实现的DnCNN网络
开发语言·matlab
奋斗的牛马13 小时前
FPGA—ZYNQ学习Helloward(二)
单片机·嵌入式硬件·学习·fpga开发
Lee_yayayayaya14 小时前
《通信之道—从微积分到5G》阅读笔记
开发语言·matlab
Voyager_415 小时前
算法学习记录03——二叉树学习笔记:从两道题看透后序位置的关键作用
笔记·学习·算法
我搞slam20 小时前
快乐数--leetcode
算法·leetcode·哈希算法
WWZZ202521 小时前
快速上手大模型:机器学习3(多元线性回归及梯度、向量化、正规方程)
人工智能·算法·机器学习·机器人·slam·具身感知
东方佑1 天前
从字符串中提取重复子串的Python算法解析
windows·python·算法
西阳未落1 天前
LeetCode——二分(进阶)
算法·leetcode·职场和发展