实验12-中成药的配制问题.doc

实验12 中成药的配制问题 12.1实验指导 本实验给出了一个中成药的配制问题，其中涉及到了向量组的线性相关性、向量组的最大无关组、向量的线性表示及向量空间等线性代数知识。并利用MATLAB软件完成具体计算。 12.2实验内容 例12.1 某中药厂用9种中草药（A、B、……、I），根据不同的比例制成了7种特效药。表12.1给出了每种特效药每包所需各种成分的质量（单位：克）。 表12.1 特效药的成分含量（单位：克） 1号成药 2号成药 3号成药 4号成药 5号成药 6号成药 7号成药 A 10 2 14 12 20 38 100 B 12 0 12 25 35 60 55 C 5 3 11 0 5 14 0 D 7 9 25 5 15 47 35 E 0 1 2 25 5 33 6 F 25 5 35 5 35 55 50 G 9 4 17 25 2 39 25 H 6 5 16 10 10 35 10 I 8 2 12 0 0 6 20 （1）某医院要购买这7种特效药，但药厂的第3号和第6号特效药已经卖完，请问能否用其它特效药配制出这两种脱销的药品。 （2）现在该医院想用这9中种草药配制三种新的特效药，表12.2给出新药所需的成分质量（单位：克）。请问如何配制。 表12.2 三种新药的成分含量（单位：克） 1号新药 2号新药 3号新药 A 40 162 88 B 62 141 67 C 14 27 8 D 44 102 51 E 53 60 7 F 50 155 80 G 71 118 38 H 41 68 21 I 14 52 30 解：（1）把每一种特效药都看成一个9维列向量，分析向量组的线性相关性。若向量组线性无关，则无法配制脱销的特效药；若向量组相关，并且能找到不含和的一个最大无关组，则可以利用现有的特效药来配制第3号和第6号药品。 在MATLAB的命令窗口输入： u1=[10;12;5;7;0;25;9;6;8]; u2=[2;0;3;9;1;5;4;5;2]; u3=[14;12;11;25;2;35;17;16;12]; u4=[12;25;0;5;25;5;25;10;0]; u5=[20;35;5;15;5;35;2;10;0]; u6=[38;60;14;47;33;55;39;35;6]; u7=[100;55;0;35;6;50;25;10;20]; U=[u1,u2,u3,u4,u5,u6,u7] [U0,r]=rref(U) 计算结果为： U = 10 2 14 12 20 38 100 12 0 12 25 35 60 55 5 3 11 0 5 14 0 7 9 25 5 15 47 35 0 1 2 25 5 33 6 25 5 35 5 35 55 50 9 4 17 25 2 39 25 6 5 16 10 10 35 10 8 2 12 0 0 6 20 U0 = 1 0 1 0 0 0 0 0 1 2 0 0 3 0 0 0 0 1 0 1 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 r = 1 2 4 5 7 从最简行阶梯矩阵U0中可以看出：特效药向量组是线性相关的，它的秩是5，它的一个最大无关组是：，则可以用现有特效药来配制出3号和6号药品：，。 （2）三种新药分别用向量和表示，把特效药向量组和新药向量组放入同一个矩