适用于DSP的四阶矩阵求逆算法
若:A= [[a(1),a(5),a(9),a(13)];
[a(2),a(6),a(10),a(14)];
[a(3),a(7),a(11),a(15)];
[a(4),a(8),a(12),a(16)]];
则:invA为A的逆求法如下:
invA=[ [invA(1), invA(5), invA(9) ,invA(13)]
[invA(2), invA(6) ,invA(10) ,invA(14)]
[invA(3), invA(7), invA(11) ,invA(15)]
[invA(4), invA(8) ,invA(12) ,invA(16)]];
detA=a(1)*a(6)*a(11)*a(16)-a(1)*a(6)*a(15)*a(12)-a(1)*a(10)*a(7)*a(16)+a(1)*a(10)*a(15)*a(8)+a(1)*a(14)*a(7)*a(12)-a(1)*a(14)*a(11)*a(8)-a(5)*a(2)*a(11)*a(16)+a(5)*a(2)*a(15)*a(12)+a(5)*a(10)*a(3)*a(16)-a(5)*a(10)*a(15)*a(4)-a(5)*a(14)*a(3)*a(12)+a(5)*a(14)*a(11)*a(4)+a(9)*a(2)*a(7)*a(16)-a(9)*a(2)*a(15)*a(8)-a(9)*a(6)*a(3)*a(16)+a(9)*a(6)*a(15)*a(4)+a(9)*a(14)*a(3)*a(8)-a(9)*a(14)*a(7)*a(4)-a(13)*a(2)*a(7)*a(12)+a(13)*a(2)*a(11)*a(8)+a(13)*a(6)*a(3)*a(12)-a(13)*a(6)*a(11)*a(4)-a(13)*a(10)*a(3)*a(8)+a(13)*a(10)*a(7)*a(4);
adjA(1)=a(6)*a(11)*a(16)+a(7)*a(12)*a(14)+a(8)*a(10)*a(15)-a(8)*a(11)*a(14)-a(7)*a(10)*a(16)-a(6)*a(12)*a(15);
adjA(2)=a(2)*a(11)*a(16)+a(10)*a(15)*a(4)+a(14)*a(3)*a(12)-a(14)*a(11)*a(4)-a(10)*a(3)*a(16)-a(2)*a(15)*a(12);
adjA(3)=a(2)*a(7)*a(16)+a(6)*a(15)*a(4)+a(14)*a(3)*a(8)-a(14)*a(7)*a(4)-a(6)*a(3)*a(16)-a(2)*a(15)*a(8);
adjA(4)=a(2)*a(7)*a(12)+a(6)*a(11)*a(4)+a(10)*a(3)*a(8)-a(10)*a(7)*a(4)-a(6)*a(3)*a(12)-a(2)*a(11)*a(8);
adjA(5)=a(5)*a(11)*a(16)+a(9)*a(15)*a(8)+a(13)*a(7)*a(12)-a(13)*a(11)*a(8)-a(9)*a(7)*a(16)-a(5)*a(15)*a(12);
adjA(6)=a(1)*a(11)*a(16)+a(9)*a(15)*a(4)+a(13)*a(3)*a(12)-a(13)*a(11)*a(4)-a(9)*a(3)*a(16)-a(1)*a(15)*a(12);
adjA(7)=a(1)*a(7)*a(16)+a(5)*a(15)*a(4)+a(13)*a(3)*a(8)-a(13)*a(7)*a(4)-a(5)*a(3)*a(16)-a(1)*a(15)*a(8);
adjA(8)=a(1)*a(7)*a(12)+a(5)*a(11)*a(4)+a(9)*a(3)*a(8)-a(9)*a(7)*a(4)-a(5)*a(3)*a(12)-a(1)*a(11)*a(8);
adjA(9)=a(5)*a(10)*a(16)+a(9)*a(14)*a(8)+a(13)*a(6)*a(12)-a(13)*a(10)*a(8)-a(9)*a(6)*a(16)-a(5)*a(14)*a(12);
adjA(10)=a(1)*a(10)*a(16)+a(9)*a(14)*a(4)+a(13)*a(2)*a(12)-a(13)*a(10)*a(4)-a(9)*a(2)*a(16)-a(1)*a(14)*a(12);
adjA(11)=a(1)*a(6)*a(16)+a(5)*a(14)*a(4)+a(13)*a(2)*a(8)-a(13)*a(6)*a(4)-a(5)*a(2)*a(16)-a(1)*a(14)*a(8);
adjA(12)=a(1)*a(6)*a(12)+a(5)*a(10)*a(4)+a(9)*a(2)*a(8)-a(9)*a(6)*a(4)-a(5)*a(2)*a(12)-a(1)*a(10)*a(8);
adjA(13)=a(5)*a(10)*a(15)+a(9)*a(14)*a(7)+a(13)*a(6)*a(11)-a(13)*a(10)*a(7)-a(9)*a(6)*a(15)-a(5)*a(14)*a(11);
adjA(14)=a(1)*a(10)*a(15)+a(9)*a(14)*a(3)+a(13)*a(2)*a(11)-a(13)*a(10)*a(3)-a(9)*a(2)*a(15)-a(1)*a(14)*a(11);
adjA(15)=a(1)*a(6)*a(15)+a(5)*a(14)*a(3)+a(13)*a(2)*a(7)-a(13)*a(6)*a(3)-a(5)*a(2)*a(15)-a(1)*a(14)*a(7);
adjA(16)=a(1)*a(6)*a(11)+a(5)*a(10)*a(3)+a(9)*a(2)*a(7)-a(9)*a(6)*a(3)-a(5)*a(2)*a(11)-a(1)*a(10)*a(7);
invA(1)=adjA(1)/detA;
invA(2)=-adjA(2)/detA;
invA(3)=adjA(3)/detA;
invA(4)=-adjA(4)/detA;
invA(5)=-adjA(5)/detA;
invA(6)=adjA(6)/detA;
invA(7)=-adjA(7)/detA;
invA(8)=adjA(8)/detA;
invA(9)=adjA(9)/detA;
invA(10)=-adjA(10)/detA;
invA(11)=adjA(11)/detA;
invA(12)=-adjA(12)/detA;
invA(13)=-adjA(13)/detA;
invA(14)=adjA(14)/detA;
invA(15)=-adjA(15)/detA;
invA(16)=adjA(16)/detA;
适用于DSP的四阶矩阵求逆算法相关推荐
- 基于DSP的图像旋转算法数据调度策略
基于DSP的图像旋转算法数据调度策略 [日期:2008-8-14] 来源:电子技术应用 作者:李筱琳,冯 燕,何亦征 [字体:大 中 小] <script src="http://w ...
- 【算法】机器人走迷宫(适用于走迷宫、最短路径算法)-20200412
标题:机器人走迷宫(适用于走迷宫.最短路径算法) 问题描述: 一块矩形方格,含有障碍和可通行格子,求从某一点到另外一点的最短距离?N*M的矩阵: 其中,1代表障碍,0代表可通行:示例:给定二维矩阵 0 ...
- 【算法】机器人走迷宫破壁解法(适用于走迷宫、最短路径算法)-20200412
标题:机器人走迷宫破壁解法(适用于走迷宫.最短路径算法)-20200412 问题描述: 一块矩形方格,含有障碍和可通行格子,求从某一点到另外一点的最短距离?N*M的矩阵: 其中,1代表障碍,0代表可通 ...
- DSP平台实现图像识别算法总体来说不如嵌入式?
基于DSP的嵌入式视觉应用平台存在的主要问题为:运算速度较低,通常需要设计专用的多核DSP.开发相应的并行处理机制.并利用其提供的指令集进行算法优化,使得专注于图像处理算法的工程师"望而生畏 ...
- DSP 中的基础算法和模型的详细解析
查看全文 http://www.taodudu.cc/news/show-2989078.html 相关文章: 唐端荣|DSP投放引擎的设计与实现 dsp 精准投放_招商加盟行业如何精准获客 DSP投 ...
- DSP音频软件与算法开发
1加我V hezkz17 进同行数字音频系统设计开发交流答疑群\n\n在目前中国的学科体系中,自然辩证法专业的哲学名称是科学技术哲学.1990年之前,"科学技术哲学"的名字是&qu ...
- c 语言矩阵求逆算法,矩阵的逆 C 语言 算法一
#include #include voidmain() {inti,j;intdimension;double deterValue=1;double **array,**deterArray, * ...
- 《基于ADAU1452 DSP车速音量补偿算法实现方法 》
加我微信 hezkz17 进同行微信交流群答疑 1 功能需求曲线图 Volume= f(v) 2 不同等级的速度曲线 3 Sigma实现方法 4 直接在MCU 中的实现方法 5 CAN 周立功验证方 ...
- DSP基础算法与模型研究
DSP基础算法与模型研究 (转载请保留原文链接 http://www.techinads.com/archives/41 authored by 江申_Johnson) 美国有一家很优秀的DSP公司- ...
最新文章
- RDKit | 删除方差低的描述符
- TSNE算法是什么?该算法有什么缺陷?
- 词性标注,实体识别,ICTCLAS分析系统的学习
- 13道Python3实例
- 《噬血代码》:一款轻度的魂Like游戏
- LeetCode Valid Sudoku
- CSS+DIV定位分析(relative,absolute,static,fixed)
- byte数组添加数据_C#基于S7协议实现对PLC中DB块字节数据的获取及自定义textbox实现数据解析...
- HTML autofocus
- python人脸识别库_基于facenet的实时人脸识别系统
- 【BZOJ 2957】 2957: 楼房重建 (线段树)
- 聚类(2)-- Fuzzy k-means clustering
- PJzhang:kali linux安装网易云音乐、Visual Studio Code、skype
- 使用vue-cropper
- 练习题2 -和可被 K 整除的子数组
- C语言字谜游戏(函数嵌套、指针)
- 最小生成树问题:算法分析 Java 实现
- 浅谈一万小时定律的表与里
- SoftWare Engineering -- WEEK.3
- 张雨潇 计算机科学,三亚考生本科第二批录取名单[一]