3 矩阵运算_【科研心得】基于MATLAB的矩阵运算与OD矩阵处理的简易教程

大家好！我是凯森赛。终于又跟大家见面了！我真是想死你们了！这一期我将介绍如何用MATLAB玩转矩阵运算以及OD矩阵处理。其实矩阵运算还可以用其他很多工具能够实现，例如Excel，R等。笔者正是由于MATLAB对矩阵运算具有易操作性才彻底喜欢上用它。有兴趣的小伙伴可以一起学习，多提意见和评论。亦可以大胆提出想法，说出诉求，碰到疑难杂症可以相互沟通，实现学术共鸣。只要能帮到大家，笔者和编辑组都会十分欣慰的(?)。阅读本文需要30-40分钟时间，理解本文需要有一定的专业基础。有条件的读者朋友可以亲自上机实践一下好过任何教程。另外上一节期 《【科研心得】基于MATLAB的曲线拟合》之中，笔者已经介绍了MATLAB的基础知识，相关的PPT课件在百度文库也有海量，大家可以尝试下载看看，在此不赘述了。U•ェ•*U  

01 【基础篇】

好的！下面开始了，快上车吧(￣)￣)↗　。矩阵，相信大家都不会陌生。矩阵的运算穿插了整个大学课程——《线性代数》。在这里笔者想再龟毛几句。记得大学那会儿《线代》也没怎么认真听，写作业就喜欢偷懒，会耍小聪明，自己就琢磨出来的歪门邪道一有空就拿Excel去算题。作业呢也基本上一口气写完不费劲 ( ‵▽′)ψ，考试前突击了一下90分飘过了。(编辑注：老凡尔赛了)至于很复杂的矩阵运算，其实笔者(擅自地)认为只要不是真心想在考研上下功夫提分或者本身就是数学专业人士(在此希望不要被我无意冒犯)，站在十几年后本人的角度去想，当时在大学中花费太多的精力去深入学习高数其实并没有那么地需要，可能会触及很多高深的我们无法触及的高数世界而受到打击导致对高数世界的探寻感到望而却步(那种求知欲被碾压的感觉你懂的)。相反如何活用高数的知识给将来的科研之路打基础，如何让高数的知识跟编程语言做结合，如何让高数思想跟算法技巧结合，如何遇到问题能够融会贯通以及如何学以致用才更重要。我们只需要记得曾经学过矩阵运算，将来记得线性代数的老师那句话就中：几乎全部的矩阵运算都可以用Excel运算，所以能来坚持听课，能好好写作业就一定不会挂科的？！

图1中演示了MATLAB文法如何定义矩阵。在变量A的后面用一个等号后接表示赋值；被中括号括住的都是矩阵内的元素；3×3矩阵里，分号表示换行，空格或者逗号表示行方向的分割。

图1

图2中演示了元素的更替和指定。只需按照图中所示范的command去执行就可以了。没有什么技术含量，在此就不赘述了。

图 2

图3演示了子矩阵的抽出和指定。顺便补充一句，冒号在矩阵之中元素(子矩阵)声明和引用的时候十分常用。假如矩阵的尺寸很大例如1000×1000的时候，可以在格式上进行如此的活用：例如A(2:end,:)就表示第2行到最后1行的所有元素所构成的子矩阵；A(2:end-1,:)就表示第2行到倒数第2行的所有元素构成的矩阵(当然前提是A这个矩阵得足够大，不然这里的例子就会返回一个空矩阵)。有兴趣的小伙伴不妨用很大的矩阵来定义A，再来体验一下这个细节，在此不赘述。