MATLAB编程线性代数.pptVIP

t [0 .3 .8 1.1 1.6 2.3]; b [.82 .72 .63 .60 .55 .50]; A [ones size t exp -t ]; c inv A*A *A*b c 0.4760 0.3413 T 0:0.1:2.5 ; Y [ones size T exp -T ]*c; plot T,Y,-,t,y,o 解 系数矩阵为 clear format rat A [1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; b [1,4,0]’; A1 [A,b] %增广矩阵 [m,n] size A ; r1 rank A1 ;r2 rank A ; 综合题：解方程组 if r1 r2r2 n X inv A *b elseif r1 r2r2 n X null A,‘r’ % 求齐次方程组的通解 x0 A\b% 求非齐次的一个特解 else X inv A*A *A*b end 案例1：减肥配方的实现 设三种食物每100克中蛋白质、碳水化合物和脂肪的含量如下表，表中还给出了美国特别流行的剑桥大学医学院的简捷营养处方。现在的问题是：如果用这三种食物作为每天的主要食物，那么它们的用量应各取多少,才能全面准确地实现这个营养要求? 营 养 每100g食物所含营养 减肥期间所要求的每日营养量 脱脂牛奶 大豆面粉 鸡蛋 蛋白质 36 57 13 33 碳水 化合物 52 34 74 45 脂肪 0 7 1.1 3 设脱脂牛奶的用量为x1个单位（100g），大豆面粉的用量为x2个单位（100g），乳清的用量为x3个单位（100g），表中的三个营养成分列向量为： 则它们的组合所具有的营养为 使这个合成的营养与剑桥配方的要求相等，就可以得到以下的非齐次线性方程： 用MATLAB解这个问题非常方便： A [36,51,13;52,34,74;0,7,1.1] b [33;45;3] x A\b 程序执行的结果为： 即脱脂牛奶的用量为27.7g，大豆面粉的用量为39.2g，鸡蛋的用量为23.3g，就能保证所需的综合营养量。 案例2 希尔密码问题 在信息传递方法中，经常将每个字母与一个整数相对应，然后传输这串整数来代替原来的信息。 例如，信息SEND MONEY可以编码为 〔 19，5，14，4，13，15，14，5，25 〕 其中S表示为19， E表示为5，依此类推，这就是替代密码法 。 那么如何使密文不被攻破呢？用矩阵乘法可以对信息进行进一步的伪装。这就是希尔密码的思想。 首先，我们将替代密码中的数字放入方阵A，然后构造密钥矩阵B，B与A可乘，且B可逆，这样传输信息就为C＝AB，显然，对C的破译是困难的。而接收的一方只要利用逆矩阵就可以还原信息。 比如，要传出SEND MONEY的信息，可以将其作成一个3*3矩阵 然后我们构造一个保密的密钥矩阵 这样将A*B作为明文矩阵信息发出，对方收到后只需要用密钥矩阵的逆矩阵对其进行解密即可得知将要表达的真实信息。 思考：如何将I love you的信息准确无误而又保密的传达给你的心上人呢？ 练习 1.试构造矩阵,练习四则运算与幂运算，练习函数det,inv,trace,rref,rank,norm,fliplr,flipud, flipdim,transpose,reshape，以及矩阵的逻辑运算: ,|,~和xor。 2、设有向量a [1,2,3],b [4,5,6],c [7,8,9],…，试 利用rref找出它们的最大无关组。 MATLAB编程之“线性代数” 延边大学 信息管理专业 13级 崔基哲 矩阵的建立 矩阵分析 矩阵分解 线性方程组的求解 数学建模举例 MATLAB与线性代数 向量 创建数组向量的三种常用方法: 1.增量法（利用冒号生成向量） x i : j 如果 i j, 生成以i为初值, j为终值,1为步长的有限等差数列; x i: j: k 如果 i k, 生成以i为初值, k为终值,j为步长的有限等差数列; x i: -j: k 如果 i k, 生成以i为初值, k为终 值,-j为步长的有限等差数列; 只需要用空格或逗号间隔数组元素， 再用方括号括起来。如：X [0 2 3 6 7 8] x linspace 2,8 将生成100个数,以2开始,以8结束。 y linspace 2,8,10 将生成10个数,以2开始,以8结束。 3. 利用logspace a,b 生成等比向量 y logspace 1,3,3 2. 利用linspace a,b 生成等差向量 如：A 10:15 A 10 11 12 13 14 15 ——矩阵的创建/ 直接输入法—