1、用机算(MATLAB)和建模实践改造工科线性代数（项目介绍）,西安电子科技大学 陈怀琛,1. 立项依据,根据高教司函2008192号文件，建立了“使用信息技术工具改造课程项目”，资助20个利用计算机软硬件改造课程的项目。其中工科基础课8项（数、理、化、电各2项），经管类6项，艺术类6项。经过申报评审，教育部下达了教高厅函20091号文，我校被批准作为 “用机算和建模实践改造工科线性代数课程” 项目的牵头单位。要求我校组织联合1520个参与单位，推广试点两年，总结经验。将来再考虑向全国推广。,2. 课程改造的理念,1)实行“需求牵引”，根据机械、电子专业后续课对矩阵运算的需求分析，把线性代数解题
2、的目标定位为：“6阶以上、复数和超定问题”。 2)实行“技术推动”：引入MATLAB软件，把传统的“只会笔算”改造为“笔算与机算相结合”。 本次会议的论文1“科学计算能力培养与线性代数改革”就是分析这个问题的。 3)把矩阵建模和机算引导到后续课中去，实现无缝连接，推动相关后续课程的现代化。 本次会议的论文2“复杂线性系统建模的矩阵方法”给出了四门后续课程（电路、信号与系统、数字信号处理和自动控制原理）利用线性代数改革推动其内容现代化的例子。,3. 项目的实施进程,项目于2009年1月启动。 一、“机算”层次：在校内，主要是对全校5000名本科生的线性代数都加入了机算的内容，作为课程的要求；在校
3、外，联合了16所参与大学，他们将在今年进行几百名学生的机算教学的试点，明年要超过1000人。这是第一个层次“机算”，关键是要培训教师。要使线性代数老师在思想上接受这个概念，需要做艰苦的工作，（其实初步掌握软件并不十分困难，）同时要准备相应的教材，我们准备了三本教材。 二、“建模”层次，主要针对后续课的改革，要争取后续课教师的参与，这个工作更加艰难。因为课程很多，教师分散，加上过去他们学的线性代数多时不用，当时也没有学软件，所以很难与线性代数课程改革挂上钩。,教材一、线性代数实践及MATLAB,本书的标志性特征： 线性代数抽象吗？看了本书后，你会知道它的概念都基于空间形象。 线性代数冗繁吗？学了
4、本书后，你会懂得它的计算全可有简明程序。 线性代数枯燥吗？读了本书后，你会发现它的应用极其广泛又精彩。 这本书是针对已学过线性代数理论的老师、工程师、研究生的“补丁书”,教材二、工程线性代数(MATLAB版),四个主要特点： 1.基本概念都从几何图形引入，做到抽象与形象的结合； 2.一切繁琐计算都有简明程序，推动笔算与机算的结合； 3.大量实例诠释了课程的价值，实现理论与实践的结合； 4.能与后续课的需要无缝衔接，体现课程的辐射效应。,关于抽象与形象的结合, 三元方程组解的几何意义（适定、超定与欠定）； 两个向量的行列式是它们组成的平行四边形的面积； 三个向量的行列式是它们组成的平行六面体的体
5、积； 二维、三维向量线性相关和线性无关的几何意义； 平面（二维）线性变换的几何特征及其意义； 二维特征值和特征向量的几何意义； 三元齐次方程基础解系的几何特征； 二元超定方程最小二乘解的几何表述； 二次型化为标准型的不同方法的几何解释,关于笔算与机算的结合, 矩阵的赋值和其加、减、乘、除（求逆）命令； 矩阵化为最简行阶梯型的计算命令；U0,ip=rref(A) 多元线性方程组MATLAB求解的几种方法；x=inv(A)*b, U=rref(A) 行列式的几种计算机求解方法；D=det(A),L,U=lu(A);D=prod(diag(L) n个m维向量组的相关性及其秩的计算方法和命令; r=r
6、ank(A),U=rref(A) 求欠定线性方程组的基础解系及超定方程解的MATLAB命令；xb=null(A) 矩阵的特征方程、特征根和特征向量的计算命令；f=poly(A);P,D=eig(A) 化二次型为标准型的MATLAB命令；yTDy=xTAx; 其中y=P-1x,关于实现理论与实践的结合, 多项式插值系数的计算 平板稳态温度的计算 交通流量的分析 成本核算问题 图及其矩阵表述 网络的矩阵分割和连接 弹性梁的柔度矩阵 用行列式计算面积,关于实现理论与实践的结合(续), 化学方程的配平 减肥配方的实现 刚体平面运动的计算和绘图 混凝土配料中的应用 圆锥截面二次型方程插值问题 人口迁徙模
7、型 物料混合问题（用到微分方程）,关于与后续课应用的衔接, 用线性代数解直流电路举例 用线性代数解交流电路举例 用线性代数解线性系统中常微分方程的举例 用线性代数解线性系统中信号流图的举例 用线性代数求数字信号处理中的系统函数举例 用线性代数解静力学问题的举例 用线性代数解运动学问题的举例 用线性代数解机械测量学问题的举例 用线性代数解文献管理问题的举例 用线性代数解经济管理问题的举例,4. 用系统工程的概念来推动,（1）教育现代化是党中央的一贯方针，体现在邓小平“三个面向”的题词和胡锦涛十七大报告中。过去教育部却没有给予足够重视，这次的立项可以说是一个好的开端。现代化就是要创新，就要和传统观
8、念发生冲突，创新开始必定是少数，维护传统的一定是多数，而且有全套的体制和书籍文件支持。因此领导部门必须特别支持创新。尽管道路不会平坦，但这个历史潮流是绝对不会改变的。 （2）谁来推动课程的现代化？教育属于上层建筑，首先要由经济基础的四个现代化对人才能力的需求反映到教育上来，首先反映到对专业课和学生的能力的现代化要求，对基础课的现代化要求则是要由专业课对基础课提出需求。这一环不抓，基础课很容易安于现状，几十年一贯制。,续用系统工程的概念来推动,（3）制定和改革教学计划的关键在于把握课程与课程之间的衔接，如果这种衔接是建立在现代化的水平之上，那么这个计划就会是先进的。这种衔接要能充分体现后续课对基
9、础课的需求，正是有这种明确不断增长的需求才能促使基础课的现代化；另一方面，基础课的改革给后续课提供了更高的平台，后续的专业课必须利用新的平台，提高本课程的起点，为此我们提出了“教学计划矩阵”的概念 3，希望教学行政部门采纳，使以后的教学计划修订建立在更加科学的基础上。 （4）参与教学计划和教学大纲制订的领导和教师，必须具备全局的观念，不仅要熟悉自己教的课程，而且要对本课程的先修课和后续课都有较好的了解。特别是，对数学基础课的大纲不能由清一色的数学专家来规定。因为从线性代数可以看出，数学专家通常并不了解现代化对工科的最新需求。所以应该由工科专家为主确定课程的外延（目标），由数学专家为主确定它的内
10、涵（实现方法）。,5、总目标实现钱学森的预见,“我想今天已是二十世纪后期，我们正面临世纪之交，所以要考虑二十一世纪会需要什么样的工科教育；保持五十年代的模式不行，保持八十年代的模式也不行。我想现在已经可以看到电子计算机对工程技术工作的影响；今后对一个问题求解可以全部让电子计算机去干，不需要人去一点一点算。而直到今天，工科理科大学一二年级的数学课是构筑在人自己去算这一要求上的。从解析几何、微积分、微分方程、复变函数论、偏微分方程等，无不如此。将来全部可以用电子计算机了，这套课就失去意义。所以理工科的数学课必须改革，数学课不是为了学生学会自己去求解，而是为了学生学会让电子计算机去求解，学会理解电子
11、计算机给出的答案，知其所以然，这就是工科教学改革的部分内容。”,自动控制中应用实例之一,右图的控制系统 求输出输入关系,1. 用代入法的计算(30分钟）,2. 用矩阵法笔算(90分钟),为了求这个4*4矩阵的逆，推导就用了两页纸，整个推导用了90分钟。因为它不能只求一个输出，必须把四个输出全部求出来,3. 矩阵+计算机(5分钟),写出程序（用5分钟） syms W1 W2 W3 W4 u1 u2 Q(1,4)=-W3; Q(2,1)=W1; Q(3,2)=W2; Q(4,3)=1; P(2,1)=W4;P(1,1)=1;P(4,2)=1; W=inv(eye(4)-Q)*P X=W*u1;u2
12、 运行此程序即得到结果,用三种方法解题的比较,1、用中学的代入法，30分钟； 2、用矩阵+笔算，90分钟 3、用矩阵+软件，5分钟 就解题的效率而言，线性代数比代数方法还要差！终于明白了，原来线性代数不是为手算设计的，它的优点必须在使用了计算机后才能显现。学线性代数而不用计算机，在学术上就是一个谬误，就好像学载重汽车的构造而不教开车、不给汽油一样，一点没有用，还不如手推车有用。,滤波器计算的实例,x1 u k3x4; x2x1; x3k3x2x4; x4qx7; x5 x2 k2x8; x6x5; x7k2x6x8; x8qx11; x9 x6 k1x12; x10x9; x11k1x10x1
13、2;x12qx10; x13y C0x12 C1x11 C2x7 C3x3 13元13个线性代数方程，怎么解？用笔算还是计算机？,本题的矩阵模型是,本题的MATLAB程序,clear, syms q % 对符号变量进行定义 k1=1/4; k2=1/2; k3=1/3; k0=1; % 给网络中的常数 C0=-0.2; C1=0.8; C2=1.5; C3=1; Q(4,7)=q; % Q的第一个赋值元素为符号变量 Q(5,2)=1; Q(1,4)=-k3; Q(2,1)=1;% 为连接矩阵的非零元素赋值 Q(3,2)=k3; Q(3,4)=1; Q(5,8)=-k2; Q(6,5)=1; Q
14、(7,6)=k2; Q(7,8)=1; Q(8,11)=q; Q(9,6)=1; Q(9,12)=-k1; Q(10,9)=1; Q(11,10)=k1; Q(11,12)=1; Q(12,10)=q; Q(13,3)=C3; Q(13,7)=C2; Q(13,11)=C1; Q(13,12)=C0; Q(13,13)=0; % 给右下角元素赋值，使未赋值元素全部为零 P(1,1)=k0; P(13,1)=0; % 给P赋值， W=inv(eye(size(Q)-Q)*P % 信号流图方程解 pretty(W(13) % 以美观的形式显示W中第13行,程序运行的结果,手工将q换成z-1，得到,
15、大大显示了线性代数的威力,1. 线性代数对几乎所有后续课都非常有用； 2. 线性代数是组织海量数据进行科学计算的最好工具；（13*13=169个元素可变） 3. 算得最多的课程变成算得最少的课程。得益于它的规范化。 我们终于看到，钱学森“指挥计算机算”的话在矩阵建模问题上得到了证实，这使我们有更大的信心沿着这个方向去努力！历史的潮流说明了这是正确而可行的。但需要基础课和后续课老师共同努力。推动自己的课程与计算机相结合的改革！这样中国教育才能迈向现代化，才能与国际接轨！,获取改革资料的两个网址,教育部“使用信息技术工具改造课程项目”的官方网站http:/ 我们“用机算和建模实践改造工科线性代数课程” 项目的网站http:/,谢谢大家！,
