适用于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公司- ...
最新文章
- 解决Http响应内容中文乱码问题
- JavaFX 之窗口大小自由拉伸(四)
- VUE el-input正则验证
- torch_geometric笔记:nn. graclus (图点分类)
- Mac : PS1变量不生效了
- [转载] QT中的connect用法总结
- iris数据集——决策树
- navicat 注册过程
- 公众号数据分析工具分析微信公众号的六大定位
- flash player 9 安装错误:您尝试安装的adobeflashplayer版本不是最新版本
- 截至9月28日聚合支付牌照知多少,附名单
- java spring boot 微信公众号 分享功能
- Thinkphp5 引入第三方类库
- IBM面试题 50条狗 有病狗 有几条病狗?
- Day03 BSP工程师小白的第一步-------今天开始正式学C语言
- 安装win11电脑必须支持TPM2.0和必须支持安全启动的解决方法
- 基因功能预测工具-HMMER的安装
- Unity入门03——Unity脚本
- 统计数据集图片通道的像素均值和方差
- java开发工具排名_干货:排名前16的Java工具类
热门文章
- [NOIP初赛复习]插空法与插板法解排列组合问题
- 用户划分——RMF方法
- 06-播放器-搭建媒体服务器测试videojs
- 迪拜“烧掉800亿”造了座烂尾岛,奇葩建筑惊呆网友:有钱人的世界,我不懂!
- Could not publish server configuration for Tomcat v8.0 Server at localhost.
- 分部积分法的快速运算 表格法
- 如何查看电脑端口占用
- php制作万年历的步骤_PHP制作万年历
- 【其他】神舟精盾K470-i3 D3更换支持5GHz的双频无线网卡
- WINDOWS更改无线网卡MAC地址的方法