用Excel进行矩阵计算 一、Excel的数组、数组名和矩阵函数的设置 1矩阵不是一个数，而是一个数组。在Excel里，数组占用一片单元域，单元域用大括号表示，例如{A1：C3}，以便和普通单元域A1：C3相区别。设置时先选定单元域，同时按Shift＋Ctrl＋Enter键，大括弧即自动产生，数组域得以确认。 2Excel的一个单元格就是一个变量，一片单元域也可以视为一组变量。为了计算上的方便，一组变量最好给一个数组名。例如A={A1：C3}、B={E1:G3}等。数组名的设置步骤是：选定数组域，点“插入”菜单下的“名称”，然后选择“定义”，输入数组名如A或B等，单击“确定”即可。 3矩阵函数是Excel进行矩阵计算的专用模块。常用的矩阵函数有MDETERM(计算一个矩阵的行列式)、MINVERSE(计算一个矩阵的逆矩阵)、MMULT(计算两个矩阵的乘积)、SUMPRODUCT(计算所有矩阵对应元素乘积之和)……函数可以通过点击“=”号，然后用键盘输入，可以通过点击“插入”菜单下的“函数”，或点击fx图标，然后选择“粘贴函数”中相应的函数输入。 二、矩阵的基本计算 数组计算和矩阵计算有很大的区别，我们用具体例子说明。 已知A={3 -2 5，6 0 3，1 5 4}，B={2 3 -1，4 1 0，5 2 -1}，将这些数据输入Excel相应的单元格，可设置成图1的形状，并作好数组的命名，即第一个数组命名为A，第二个数组命名为B。计算时先选定矩阵计算结果的输出域，3×3的矩阵，输出仍是3×3个单元格，然后输入公式，公式前必须加上=号，例如=A ＋B、=A－B、=A＊B等。A＋B、A－B数组运算和矩阵运算没有区别，“=A＊B”是数组相乘计算公式，而“=MMULT(A，B)”则是矩阵相乘计算公式，“=A/B”是数组A除数组B的计算公式，而矩阵相除是矩阵A 乘B的逆矩阵，所以计算公式是“=MMULT(A,MINVERSE(B))”。公式输入后，同时按Shift＋Ctrl＋Enter 键得到计算结果。图1中的数组乘除写作A＊B、A/B，矩阵乘除写作A·B、A÷B，以示区别。 三、矩阵计算的应用 下面让我们来计算一个灰色预测模型。 灰色预测是华中理工大学邓聚龙教授创立的理论，其中关键的计算公式是计算微分方程＋B1x=B2的解，{B1，B2}=(XTX)－1(XTY)，式中：XT是矩阵X的转置。 作为例子，已知X={-45.5 1，-79 1，-113.5 1，-149.5 1}Y={33，34，35，37} 在Excel表格中，{B2：C5}输入X，{E2：H3}输入X的转置。处理转置的方法是：选定原数组{B2：C5}，点“编辑”菜单的“复制”，再选定数组转置区域{E2：H3}，点“编辑”菜单的“选择性粘贴”，再点“转置”即可。{J2：J5}输入Y，然后选取{L2：L3}为B1、B2的输出区域，然后输入公式： =MMULT(MINVERSE(MMULT(E2:H3,B2:C5)),MMULT(E2:H3,J2:J5)) 公式输入完毕，同时按Shift＋Ctrl＋Enter键，B1、B2的答案就出来了，如图2。 如果计算的矩阵更复杂一些，就必须分步计算。不过，使用Excel也是很方便的。(江苏陈岁松) ==== POWERPOINT 演示文档https://www.doczj.com/doc/7a16611565.html,.tw/~ccw/manage_math/array.ppt EXCEL矩陣運算（繁体中文）参考文献：https://www.doczj.com/doc/7a16611565.html,/4/wenzi/wz042.htm
C语言课程设计 题目 XXXXXXX 二级学院电子与信息工程学院 班级 1X级XXXXX 专业X班 姓名 XXX 学号 XXX 指导教师梁建勇 2016年07月01日
通讯录管理系统 电子与信息工程学院C语言课程设计任务书
2014022061 刘言松 通讯录管理系统 刘言松安康学院 2014级电子与信息工程系电子信息工程2班 摘要：通讯录管理系统是生活中不可缺少的部分编写了一个通讯录管理系统是十分必要的。本程序设计具有通讯信息的输入，进行保存、查找、删除等功能，操作界面简洁美观，易于操作。程序用了条件、循环、结构体等知识点。可运用于小的管理软件。软件运用菜单实现交互式管理，用户轻松的按键既可实现对软件的操作，操作简单易懂，功能丰富，可以很好的满足的需要。通过对通讯录管理系统的设计，进一步理解和掌握C语言这门课程的知识点，能够熟练的调用各种函数，把各种C语句有机的结合起来，提高自己C语言程序设计的能力，为今后解决实际问题打下良好基础。 关键词：通讯录管理系统；结构体；条件语句；循环结构；指针 1引言 学习了C语言程序设计这门课程之后，使我对程序设计有了一定的了解，为了熟练掌握C 语言的语法特点、及其数据结构，提高自身编写程序的能力。通过课程设计实现理论知识与实际情况的结合，使所学知识能够在现实社会中起到一定的作用，防止所编写的程序脱离实际，让程序可以更好的发挥作用，以便减轻实际工作中所遇到的繁琐步骤，提高现实工作中的效率。此外通过课程设计提高逻辑思考能力和解决实际问题的能力。程序设计中，语句的语法和常量变量的定义都有严格的要求，有时输了一个中文标点、打错了一个字母，编译就不通过，程序无法正常运行。因此，程序设计初学阶段，学生经常会犯这样的错误，可能要通过几次乃至十多次的反复修改、调试，才能成功，但这种现象会随着学习的深入而慢慢改观。这当中就有一个严谨治学、一丝不苟的科学精神的培养，又有一个不怕失败、百折不挠品格的锻炼。通讯录信息管理系统是生活中不可缺少的部分编写了一个通讯录信息管理系统是十分必要的。通过程序设计，达到理论与实际应用相结合，提高对信息管理系统的分析能力，能够根据实际应用，初步实现系统功能模块及算法的设计，通过编程基本实现信息管理系统。 2设计方案 2.1设计思路 建立一个函数，用来添加、显示、删除、查询、修改等信息，完成通讯录的功能。主函数可以调用五个子函数，分别完成添加记录、显示记录、删除记录、查询记录、修改记录、退出系统等功能。在主函数中可以以1、2、3、4、5、数字键分别可以执行某个功能模块 3程序设计与实施 3.1根据该流程图定义数据类型建立函数 typedef struct { char score; //编号 char name[10]; //姓名 char num[15]; //电话号码
第3章 矩阵及其运算 3.1 基本要求、重点难点 基本要求： 1．1．掌握矩阵的定义. 2．2．掌握矩阵的运算法则. 3．3．掌握伴随矩阵的概念及利用伴随矩阵求逆矩阵的方法. 4．4．掌握矩阵秩的概念及求矩阵秩的方法. 5．5． 掌握初等变换和初等矩阵的概念，能够利用初等变换计算矩阵的秩，求可逆矩阵的逆矩阵. 6．6．掌握线形方程组有解得判定定理及其初等变换解线形方程组的方法. 重点难点：重点是矩阵定义，矩阵乘法运算，逆矩阵的求法，矩阵的秩，初等 变换及线性方程组的解. 难点是矩阵乘法，求逆矩阵的伴随矩阵方法. 3.2 基本内容 3.2.1 3.2.1 重要定义 定义3.1 由n m ?个数)2,1;,2,1(n j m i a ij ==组成的m 行n 列的数表成为一个m 行n 列矩阵，记为 ????????????mn m m n n a a a a a a a a a 2122221 11211 简记为A n m ij a ?=)(，或A )(ij a =,n m A ?,mn A 注意行列式与矩阵的区别： （1） （1） 行列式是一个数，而矩阵是一个数表. （2） （2） 行列式的行数、列数一定相同，但矩阵的行数、列数不一定相 同. （3） （3） 一个数乘以行列式，等于这个数乘以行列式的某行（或列）的所有元素，而一个数乘以矩阵等于这个数乘以矩阵的所有元素. （4） （4） 两个行列式相等只要它们表示的数值相等即可，而两个矩阵相等则要求两个矩阵对应元素相等. （5） （5） 当0||≠A 时，||1A 有意义，而A 1 无意义.
n m =的矩阵叫做阶方阵或m 阶方阵.一阶方阵在书写时不写括号，它在 运算中可看做一个数. 对角线以下（上）元素都是0的矩阵叫上（下）三角矩阵，既是上三角阵， 又是下三角的矩阵，也就是除对角线以外的元素全是0的矩阵叫对角矩阵.在对角矩阵中，对角线上元素全一样的矩阵叫数量矩阵；数量矩阵中，对角线元素全是1的n 阶矩阵叫n 阶单位矩阵，常记为n E （或n I ），简记为E （或I ），元素都是0的矩阵叫零矩阵，记为n m 0?，或简记为0. 行和列分别相等的两个矩阵叫做同型矩阵，两个同型矩阵的且对应位置上的 元素分别相等的矩阵叫做相等矩阵. 设有矩阵A =n m ij a ?)(，则A -n m ij a ?-=)(称为A 的负矩阵. 若A 是方阵，则保持相对元素不变而得到的行列式称为方针A 的行列式，记 为||A 或A Det . 将矩阵A 的行列式互换所得到的矩阵为A 的转置矩阵，记为T A 或A '. 若方阵A 满足A A T =，则称A 为对称矩阵，若方阵A 满足A A T -=，则称A 为反对称矩阵. 若矩阵的元素都是实数，则矩阵称为实矩阵.若矩阵的元素含有复数，则称矩 阵为复矩阵，若A =n m ij a ?)(是复矩阵，则称矩阵n m ij a ?)(（其中ij a 为ij a 的共轭矩阵，记为A n m ij a ?=)(. 定义3.2 对于n 阶矩阵A ，如果存在n 阶矩阵B ，使得E BA AB ==，则 称方阵A 可逆，B 称为A 的逆矩阵，记做1-=A B . 对于方阵A n m ij a ?=)(，设ij a 的代数余子式为ij A ，则矩阵 *A ????????????=nm n n n n A A A A A A A A A 2122212 12111 称为A 的伴随矩阵，要注意伴随矩阵中元素的位置. 定义3.3 设有矩阵A ，如果： （1） （1） 在A 中有一个r 阶子式D 不为零.
struct jdlb *insert1(struct jdlb *tp,struct jdlb *z) //节点导纳矩阵插入指针数据{ struct jdlb *p0,*p111,*p112; double r,r1,x,x1; kk=0; p111=tp; p0=z; if(p0==null) return(tp); if(tp==null) { tp=p0; p0->next=null; return(tp); } if(p0->lnxt
} return(tp); } //线路部分形成节点导纳矩阵 p1=(struct jdlb *)malloc(len); p2=(struct jdlb *)malloc(len); p3=headlij; //线路部分 while(p3!=null) //形成节点导纳矩阵,可为双边的{ r=p3->fu; x=p3->jd; bb=p3->bb; i=p3->i; j=p3->j; gij=r/(r*r+x*x); bij=-x/(r*r+x*x); r=-gij; x=-bij; tmp=sqrt(r*r+x*x); if(tmp!=0) { p1->irow=i; p1->lnxt=j; p1->fu=tmp; p1->jd=atan2(x,r); a[i]=insert1(a[i],p1); if(kk==0) p1=(struct jdlb *)malloc(len); p2->irow=j; p2->lnxt=i; p2->fu=tmp; p2->jd=atan2(x,r); a[j]=insert1(a[j],p2); if(kk==0) p2=(struct jdlb *)malloc(len); } p3->fu=tmp; p3->jd=atan2(x,r); p3->bb=bb; g[i]=g[i]+gij; b[i]=b[i]+bij+bb; g[j]=g[j]+gij; b[j]=b[j]+bij+bb;
C语言课程设计报告 班级：16 物联网 一、题目说明 （1）读取输入信息，能正确进行数字运算
（2）拥有记忆功能，可查询最近几次计算结果并继续进行数字运算。 二、基本功能实现 （1）能实现两个数之间的加减，乘，除，乘方，开方等运算。 （2）实现记忆功能，可保留近十位结果，可用#键查看。（3）*键作为清0键，清除当前屏幕上的值为0表示重新开三、设计流程（可以画流程图，可以直接写说明）
四、代码 #include
fgets(str,40,stdin); //将从键盘的一整行标准输入保存在字符串str中 if(sscanf(str,"%lf", &a) == 1){ //将str以长浮点型的格式放入a中，如果返回值为1，说明输入的是数值， scanf("%c %lf", &opt, &b); //获取操作符的值和另外一个操作数的值 res = yunsuan(a,opt,b); //调用yunsuan（）函数计算，并将结果给res printf("%lf %c %lf = %lf\n\n",a, opt, b, res); PutRes(history, res, &i); //调用函数PutRes()将本次结果保存到数组history中（i用引用传递） // printf("%d\n", i); }else if(sscanf(str,"%c",&opt) == 1){ //将str以字符型的格式放入opt中，如果返回值为1，说明输入的是字符 if(opt == '#'){ //如果opt为‘#’，则显示history 中的数值 int j = 0; for(j = 0; j < i; ++j){ //i为history中存放的数值个数 printf(" %lf ", history[j]); } printf("\n");
计算实习报告 一 实习目的 （1）了解矩阵特征值与相应特征向量求解的意义，理解幂法和反幂法的原理， 能编制此算法的程序，并能求解实际问题。 （2）通过对比非线性方程的迭代法，理解线性方程组迭代解法的原理，学会编 写Jacobi 迭代法程序，并能求解中小型非线性方程组。初始点对收敛性质及收 敛速度的影响。 （3）理解 QR 法计算矩阵特征值与特征向量的原理，能编制此算法的程序，并 用于实际问题的求解。 二 问题定义及题目分析 1. 分别用幂法和幂法加速技术求矩阵 2.5 2.5 3.00.50.0 5.0 2.0 2.00.50.5 4.0 2.52.5 2.5 5.0 3.5-?? ?- ?= ?-- ?--?? A 的主特征值和特征向量． 2. 对于实对称矩阵n n ?∈A R ，用Jacobi 方法编写其程序，并用所编程序求下列矩阵的全部 特征值． 1515 4 1141144114114?-?? ?-- ? ?- ?= ? ?- ?-- ? ?-??A 3. 对于实矩阵n n ?∈A R ，用QR 方法编写其程序，并用所编程序求下列矩阵的全部特征值: 111 21 113,4,5,62311111n n n n n n ? ???? ?????==+? ????? ??+??A 三 概要设计 （1） 幂法用于求按模最大的特征值及其对应的特征向量的一种数值算法，
它要求矩阵 A 的特征值有如下关系： 12n ...λλλ>≥≥ ，对于相应 的特征向量。其算法如下： Step 0：初始化数据0,, 1.A z k = Step 1：计算1k k y A z +=。 Step 2：令 k k m y ∞=。 Step 3：令 k k k z y m = ；如果1k k m m +≈或1k k z z +≈，则 goto Step 4；否则 ， k = k + 1 ，goto Step 1。 Step 4：输出结果 算法说明与要求 输入参数为实数矩阵、初始向量、误差限与最大迭代次数。输出 参数为特征值及相对应的特征向量。注意初始向量不能为“0”向量。 （2） 迭代法的原理 如果能将方程 Ax =b 改写成等价形式：x=Bx+f 。如果B 满足：ρ(B )<1，则对于任意初始向量 x (0) ，由迭代 x ( k + 1) = Bx (k ) + f 产生的序列均收敛到方程组的精确解。迭代法中两种最有名的迭代法就是Jacobi 迭代法，它的迭代矩阵 B 为： 1()J D L U -=-+,1 f D b -= 其中，D 为系数矩阵 A 的对角元所组成对角矩阵，L 为系数矩阵 A 的对角元下方所有元素所组成的下三角矩阵，U 为系数矩阵 A 的对角元上方所有元素所组成的上三角矩阵。 算法如下： Step 0：初始化数据 00,,,,k A b x δ=和ε。 Step 1：计算D,L,U,J 或G, 得到迭代矩阵B. Step 2：:1k k =+ 0x B x f * =+ 0x x = 如果0x x δ-<或()f x ε≤，goto Step 3?否则 goto Step ２。 Step 3：输出结果。 程序说明与要求
电力网节点导纳矩阵计算例题与程序 佘名寰 编写 用计算机解算电力网潮流电压和短路电流问题首先需确定电力网的节点导纳矩阵或节点阻抗矩阵。本文通过例题介绍用网络拓扑法计算节点导纳矩阵的方法和程序，程序考虑了线路并联电容和变压器支路标么变比不为1时的影响。程序用MATLAB 语言编写，线路参数均采用标么值。本文稿用office word 2007 版编写，可供电气专业人员计算相关问题时参考。 1.用网络拓扑计算节点导纳矩阵 1.1网络拓扑矩阵： 【例1.1】 例图1-1是有5 个节点和5条支路的网络，节点5作为基准参考点，1 ,2, 3, 4为独立节点，支路编号和方向图中已标识。 例图1-1 对于具有n 个节点b 条支路的有向图，它的关联矩阵为一个N ×B 的矩阵A a ： A a =[a ij ] 若支路j 与节点i 相关，且箭头背离节点i ，则a ij =1，若箭头指向节点则a ij =-1，若支路j 与节点i 无关，则a ij =0， 图1-1所示的有向图的关联矩阵为 ① ② ③ ④ ⑤ 支路编号 A ij =[ ?10100110?100?1?10?100.01000001] 行编号从上到下为1 2 3 4 5节点编号（5为参考节点） 去掉第5行即为独立节点的关联矩阵。 以下介绍生成网络关联矩阵的M 函数文件 ffm.m ： % M FUNCTION ffm.m
% Np is number of node point,Nb is number of braches % nstart--the start point of branches ,nend -- the end point, % A -- network incidence matrix function [A]=ffm(nstart,nend) global Np Nb n=length(nstart); A=zeros(Np,Nb); for i=1:n A(nstart(i),i)=1; A(nend(i),i)=-1; end 以例图1-1网络为例调用ffm.m 文件求其关联矩阵 运算以上程序可得关联矩阵 mm ij 如下： mm = -1 0 1 0 0 1 1 0 -1 0 0 -1 -1 0 -1 0 0 0 1 0 0 0 0 0 1 Mm ij 明显与A ij 是相同的。 1.2生成节点导纳矩阵程序： ⑴ 由网络原始矩阵计算节点导纳矩阵公式 Y=AY s0A t (1-1) Y ----节点导纳矩阵 A------网络关联矩阵 A t -----A 的转置矩阵 Y S0----网络原始导纳矩阵 若网络各支路阻抗为 Z b =[z b1,z b2,……,z bn ] 则Z S0=[ z b100000z b20000 0:0000.0:00 00 0z bn ] Y s0=Z s0-1 (1-2) Y=A Z s0-1A t (1-3) ⑵ 节点导纳计算程序 以例1-1网络为例，在不计对地电容和变压器变比假定为1条件下，节点导纳矩阵计算程序如下： clear global Np Nb % Np is number of node point,Nb is number of braches, Np=5;Nb=5;
Excel矩阵计算 目前有很多软件可进行矩阵运算，特别是Matlab，其矩阵运算功能尤为强大。但这些专业软件所占空间很大，价格昂贵。其实Excel就有矩阵运算功能，虽然比不上专业软件，但不比一些数学小软件差多少。下面把从网上搜集到的一些有关利用Excel进行矩阵计算的资料整理如下： 资料一： （1）数组和矩阵的定义 矩阵不是一个数，而是一个数组。在Excel里，数组占用一片单元域，单元域用大括号表示，例如{A1：C3}，以便和普通单元域A1：C3相区别。设置时先选定单元域，同时按Shift＋Ctrl＋Enter键，大括弧即自动产生，数组域得以确认。 一个单元格就是一个变量，一片单元域也可以视为一组变量。为了计算上的方便，一组变量最好给一个数组名。例如A={A1：C3}、 B={E1:G3}等。数组名的设置步骤是：选定数组域，单击“插入”菜单，选择“名称”项中的“定义”命令，输入数组名，单击“确定”按钮即可。更简单的命名办法为：选择数组域，单击名称框，直接输入名称就行了。 矩阵函数是Excel进行矩阵计算的专用模块。用“插入”－“函数”命令打开“粘贴函数”对话框（如图11），选中函数分类栏中的“数学与三角函数”，在右边栏常用的矩阵函数有：MDETERM--计算一个矩阵的行列式；MINVERSE--计算一个矩阵的逆矩阵；MMULT--计算两个矩阵的乘积； SUMPRODUCT--计算所有矩阵对应元素乘积之和。 （2）矩阵的基本计算 数组计算和矩阵计算有很大的区别，比如下面这个例子中，A和B都是定义好的数组，因为这两个数组都是3×3的，输出结果也是3×3个单元格。计算时先选定矩阵计算结果的输出域，为3×3的单元格区域，然后输入公式。如果输入“=A＋B”或“=A-B”，计算结果是数组对应项相加或相减，输入“=A*B”表示数组A和B相乘，输入“=A/B”表示数组A 除数组B。如果要进行矩阵计算，就要用到相应的矩阵函数。矩阵相加、相减与数组的加减表达形式是一样的，也是“=A＋B”和“=A-B”，表示矩阵相乘可以输入“=MMULT(A，B)”，而矩阵相除是矩阵A乘B 的逆矩阵，所以计算公式是“=MMULT(A，MINVERSE（B）)”。公
山西工程技术学院《C语言程序设计》 课程设计说明书 姓名张雅男 学号 150721035 班级 15级计应1班 指导老师刘红梅 2016年1月6日
目录 摘要------------------------------------------------------ 2 1 引言------------------------------------------------------ 3 1.1设计目的------------------------------------------------ 3 1.2设计要求------------------------------------------------ 3 2 系统分析------------------------------------------------- 4 2.1功能需求------------------------------------------------ 4 2.2数据需求------------------------------------------------ 5 2.3程序所需函数-------------------------------------------- 5 3 详细设计与实现-------------------------------------------- 6 3.1功能模块的设计与实现------------------------------------ 6 3.2 程序的详细设计过程-------------------------------------- 7 4 系统测试-------------------------------------------------- 14 4.1设计测试数据-------------------------------------------- 14 4.2 测试结果及分析------------------------------------------ 14 总结-------------------------------------------------------- 22 参考文献---------------------------------------------------- 23
电力系统分析作业——电网节点导纳矩阵的计算机形成 编程软件：matlab R2010b 程序说明： 1.如果已经输入i-j支路的信息，则不可再输入j-i支路的信息。 2.变压器支路的第一个节点编号默认为变压器一次侧，即变压器的等值电路中的阻抗归算侧，亦即变压器非标准变比的1:k中的‘1’。 3.标幺值等值电路中，如果变比为1:1,则默认为线路，因此，变压器的非标准变比不可以是1:1。 5.如果变压器支路也有导纳B不为零，则说明此导纳就是励磁导纳，与线路的导纳B/2不同含义，只算作变压器原边的自导纳。 4.由于程序执行的是复数运算，所以即使实部为零时，也会输出实部‘0’。 程序代码: a=load('');%从’’中读入数据 [m,n]=size(a); w=1i; u=1; while (u<=m) hnode=a(u,1); enode=a(u,2); z=a(u,3)+a(u,4)*w; b=a(u,5)*w; k=a(u,6); y(hnode,enode)=-1/(k*z); y(enode,hnode)=-1/(k*z); y(hnode,hnode)=y(hnode,hnode)+1/(k*z)+(k-1)/(k*z); y(enode,enode)=y(enode,enode)+1/(k*z)+(1-k)/(k*k*z); if (abs(k-1)<%如果为线路 y(hnode,hnode)=y(hnode,hnode)+b; y(enode,enode)=y(enode,enode)+b; end
if (abs(k-1)>%如果为变压器 y(hnode,hnode)= y(hnode,hnode)-b; end u=u+1; end [m,n]=size(y); disp(‘Y=’); disp(y(1:m,1:n)); clear; 算例 输入数据： 首端编号末端编号电阻电抗电纳/2 变比 2 3 1 4 2 0 0 5 3 0 0 1 2 1 1 3 0 1 输出数据： Y= - + + 0 0 + + 0 + 0 + + 0 0 + 0 0 + 0 0 0 0 0 0 + 0 0 经手算校验，程序结果准确。
/************************ Matrix Chain Multiplication ***************************/ /************************ 作者：Hugo ***************************/ /************************ 最后修改日期：2015.09.10 ***************************/ /************************ 最后修改人：Hugo ***************************/ using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.Threading.Tasks; using System.Text.RegularExpressions; using System.Collections; namespace Matrix { class Program { public static int nummulti = 0; static ArrayList list1 = new ArrayList();//定义计算式存储列表 static ArrayList listrow = new ArrayList();//定义矩阵行数存储列表 static ArrayList listcolumn = new ArrayList();//定义矩阵列数存储列表 static void Main(string[] args) { /****************************************************************************** *****************/ //从键盘上获取矩阵 int nummatrix = Int32.Parse(Console.ReadLine()); int countmat = 0; for (countmat = 0; countmat < nummatrix; countmat++) { string s = Console.ReadLine(); string[] str = s.Split(' ');//把输入的一行字符按空格拆分 listrow.Add(Int32.Parse(str[1]));//行数存储到矩阵行数存储列表 listcolumn.Add(Int32.Parse(str[2]));//列数存储到矩阵列数存储列表
Matlab形成节点导纳矩阵 学号：0214393 姓名：侯成滨 引言：电力网的运行状态可用节点方程或回路方程来描述。节点导纳矩阵是以系统元件的等值导纳为基础所建立的、描述电力网络各节点电压和注入电流之间关系的线性方程。导纳矩阵计算是电力系统分析最基本的计算。除它自身的重要作用之外，还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础。本次任务是用计MATLAB语言编写程序求出潮流计算中要用到的导纳矩阵。为了确定结果是否正确，与一个手工计算比较运算结果，验证程序是否正确。 一、分析网络等效电路 此电力系统是一个6节点，7支路的电力网络。可以把系统等值网络画出来，如图1-1。 图1-1 某电力系统的等值网络 在计算电力系统网络的潮流分布时，我们需要把变压器转化成变压器的∏型等值电路来进行计算器等效导纳，根据等效电路的等效原则，可以把上图等效成如图1-2导纳等值网络图。对导纳等值网络图简化电路图，可以得到图1-3简化导纳等值电路图，方便潮流计算中导纳矩阵的计算。
图1-2电力系统网络的导纳等值电路 图1-3电力系统简化等值电路图 二、MATLAB程序形成导纳矩阵 导纳矩阵的计算总结如下： 1）导纳矩阵的阶数等于电力系统网络的节点数； 2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连的不接地支路数；
3）导纳矩阵的对角元素，即各节点的自导纳等于相应节点所连之路的导纳之和： Y ij=y ij j∈i 其中，y ij为节点i与节点支路阻抗Z ij的倒数，符号j∈i表示j属于i或与i相连的j，即∑内只包括与节点i直接相连的节点j。当节点i有接地支路时，还应包括j=0的情况。 4）导纳矩阵非对角元素等于节点i与节点j之间的导纳的负数。 2.1 MATLAB程序及其运行 节点导纳程序如下： N=input('请输入节点数: N='); L=input('请输入支路数: L='); B=input('请输入支路信息: B='); X=input('请输入由节点号及其对地阻抗形成的矩阵：X='); Y=zeros(N); for n=1:N; if X(n,2)~=0; p=X(n,1); Y(p,p)=1./X(n,2); end end for n=1:L if B(n,6)==0 p=B(n,1);q=B(n,2); else p=B(n,2);q=B(n,1); end Y(p,q)=Y(p,q)-1./(B(n,3)*B(n,5)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B(n,3)*B(n,5)^2)+B(n,4)./2; Y(p,p)=Y(p,p)+1./B(n,3)+B(n,4)./2; end disp('导纳矩阵Y='); disp(Y) 运行结果如下：
Excel中矩阵的计算一、求逆矩阵 （1）打开一个新的空工作簿，如图所示 （2）输入数据。 （3）在另外的活动单元格中拉黑同行列，如图
（4）点击“公式”-fx，点击“数学与三角函数” （5）点击求逆矩阵的函数键MINVERSE (6)点击确定,输入原数据所在块的第一个数据的行列 ,加“:”,输入最后一个数据的行列
(7) 点击确定后，计算后会返回一个值 (8)按F2，然后CRTL+SHIFT+ENTER，就会显示出一个三行三列的矩阵，即原矩阵的逆矩阵 二、其他：矩阵法解方程组 Excel的数组、数组名和矩阵函数的设置 矩阵不是一个数，而是一个数组。在Excel里，数组占用一片单元域，单元域用大括号表示，例如{A1：C3}，以便和普通单元域A1：C3相区别。设置时先选定单元域，同时按Shift＋Ctrl＋Enter键，大括弧即自动产生，数组域得以确认。 Excel的一个单元格就是一个变量，一片单元域也可以视为一组变量。为了计算上的方便，一组变量最好给一个数组名。例如A={A1：C3}、B={E1:G3}等。 具体操作 1、数组名的设置。
选定数组域，点“插入”菜单下的“名称”，然后选择“定义”，输入数组名如A或B等，单击“确定”即可。 或是： 选定要命名的单元格，点右键——“命名单元格区域”出现下图即可“命名A”； 选定要命名的单元格，点右键——“命名单元格区域”出现下图即可“命名B”； 如：已知A={3 -2 5，6 0 3，1 5 4}，B={2 3 -1，4 1 0，5 2 -1}，将这些数据输入Excel 相应的单元格，可设置成图1的形状，并作好数组的命名，即第一个数组命名为A，第二个数组命名为B。 2、矩阵函数是Excel进行矩阵计算的专用模块。 常用的矩阵函数有： MDETERM(计算一个矩阵的行列式) MINVERSE(计算一个矩阵的逆矩阵) MMULT(计算两个矩阵的乘积) SUMPRODUCT(计算所有矩阵对应元素乘积之和) TRANSPOSE(计算矩阵的转置矩阵)…… 函数可以通过点击“=”号，然后用键盘输入，可以通过点击“插入”菜单下的“函数”； 或点击fx图标，然后选择“粘贴函数”中相应的函数输入。 3、计算时先选定矩阵计算结果的输出域，3×3的矩阵，输出仍是3×3个单元格，然后输入公式，公式前必须加上=号，例如=A＋B、=A－B、=A＊B等。 A＋B、A－B数组运算和矩阵运算没有区别 =A+B 按enter，F2，Shift＋Ctrl＋Enter 然后选定预选区域，按F2，Shift＋Ctrl＋Enter，即可得结果 “=A＊B”是数组相乘计算公式，而“=MMULT(A，B)”则是矩阵相乘计算公式， “=A/B”是数组A除数组B的计算公式，而矩阵相除是矩阵A乘B的逆矩阵，所以计算公式是“=MMULT(A,MINVERSE(B))”。 公式输入后，同时按F2, 然后Shift＋Ctrl＋Enter键得到计算结果。 图中的数组乘除写作A＊B、A/B，矩阵乘除写作A·B、A÷B，以示区别。
C语言课程设计报告 设计名称C语言课程设计 专业班级 姓名 学号 指引教师 实习时间 成绩
目录 一、设计目................................................................. 错误!未定义书签。 二、设计内容............................................................. 错误!未定义书签。 三、详细设计阐明..................................................... 错误!未定义书签。 四、软件使用阐明..................................................... 错误!未定义书签。 五、设计总结（课程设计心得与体会） ................ 错误!未定义书签。 六、附录1：参照文献 ............................................. 错误!未定义书签。 七、附录2：某些程序清单（带有较详细注释） . 错误!未定义书签。 八、教师评语............................................................. 错误!未定义书签。
一、设计目 通过课程设计实践环节教学，可以加深学生对课堂所学基本知识掌握与理解，提高学生对所学内容综合运用能力；同步也可以通过查询有关资料，培养学生自学能力、接受新知识能力，提高学习兴趣；增强学生程序设计能力，掌握编程技巧，并可培养学生事实上机调试程序能力。“理论与实践”相结合，使学生得到较好锻炼，为后来学习、工作打下坚实基本。每个学生需要独立完毕教师安排设计题目或由学生自拟实际设计题目，但自拟自拟题目须经指引教师批准方可生效。
E x c e l中矩阵的运算
nxn方阵对应行列式的值 第二步，选中A4单元格，在“插入”菜单中选中“函数”菜单项： 第三步，在打开的“函数”对话框中，选中“MDETERM”函数如图2，并按“确定”按钮： 第四步，在弹出的对话框中输入矩阵所在的地址，按确定即得到行列式的值。 矩阵求和 已知 第二步，在A5单元格中输入公式：=A1+El，按回车，这时A5中显示数字7； 第三步，选中A5单元格，移动鼠标至其右下角，鼠标形状变为黑色十字时，按下鼠标左键往右拖至C5，B5和C5中分别显示一3．3。同样的方法选中A5：C5，往下拖至A7：C7，便得到A+B的值。 矩阵求逆 第一步，在A1：C3中输入矩阵A； 第二步。选中A5：C7，“插入”→“函数”→“MINVERSE”→“确定”： 第三步，在“array”项中输入A1：C3，按F2，同时按CTRL+SHIFF+ENTER 即可如图6。 5矩阵转置
第一步，在Al：C3中输入矩阵A，并选中； 第二步，“编辑”→“复制”； 第三步，选中A5，“编辑”→“选择性粘贴”→“转置”→确定”。 矩阵求秩 6．1矩阵秩的概念 定义设A是mxn矩阵，从A中任取k行k列(k≤min(m，n))，由这些行、列相交处的元素按原来的次序所构成的阶行列式，称为矩阵A的一个k阶子行列式，简称k阶子式。 定义矩阵A的所有不为零的子式的最高阶数r称为矩阵A的秩，记作 r(A)，即r(A)=r。 6．2矩阵秩的数学求法 6．2．1行列式法：即定义从矩阵的最高阶子式算起，计算出不等于零的子式的最高阶数r，此r即为该矩阵的秩。 6．2．2行初等变换法：用初等行变换化矩阵为阶梯形矩阵，此阶梯形矩阵非零行的行数r就是该矩阵的秩。 6．3利用EXCEL求矩阵秩 方法一，根据矩阵秩的定义，可以求所有不为零子式的最高阶数。 求矩阵A的秩． 显然A是4x4矩阵，4为其所有子式的最高阶数。先求IAI的值，若｜A｜不为零，则矩阵A的秩为4。若｜A｜为零，求所有阶数为3的子式的值。若存
HUBEI UNIVERSITY OF AUTOMOTIVE TECHNOLOGY C语言程序设计 课程设计报告 课设题目：学生信息管理系统 专业：电气工程及其自动化 班级: 学生信息管理系统 姓名： 成绩： 完成日期: 2012年6月18日－2012年6月22日 指导教师:
目录 1.课程设计目的 (2) 2.需求分析 (2) 2.1要求 (2) 2.2任务描述 (2) 2.3运行环境 (3) 2.4开发工具 (3) 3.总体设计 (3) 3.1系统流设计 (3) 3.2系统模块设计 (4) 4.详细设计 (5) 4.1数据结构设计 (5) 4.2模块接口设计 (6) 4.3 以下是几个主要模块的流程图： (7) 5.运行结果 (10) 5.1主菜单界面： (10) 5.2信息录入界面： (10) 5.3信息浏览界面： (11) 5.4 信息查询界面： (11) 5.5信息修改界面： (12) 6.程序调试与测试 (12) 6.1程序调试 (12) 6.2程序测试 (12) 7.问题与不足 (12) 8.建议 (12) 9.编程体会 (12) 10.程序源代码： (12) 11.参考文献 (13) 12.用户使用手册 (13)
1.课程设计目的 1)通过课程设计，加深对结构化程序设计思想的理解，能对系统功能进行分析， 并设计合理的模块化结构。 2)通过课程设计，进一步掌握C语言中重要的数据结构。掌握并熟练运用指针、 链表、结构体、文件等数据结构。 3)通过课程设计，提高程序开发能力，能运用合理的控制流编写清晰高效的程序。 4)通过课程设计，培养C语言知识运用和自学能力，真正的把课堂上获得的知识 运用起来，培养自己对程序编写的兴趣，并能独立设计和实现一个小型系统。 5)通过课程设计，加强了自己的程序分析能力。 2.需求分析 2.1要求 1)遵循软件开发过程的基本规范；能灵活应用，初步具备运用C程序解决实际问 题的能力。 2)运用结构化程序设计方法，按照课程设计题目要求，完成包括需求分析、数据 结构设计、概要设计、详细设计、编写代码、调试和测试C程序等过程。 3)学会利用数组、结构体、链表进行相关信息表示。 4)功能模块划分清晰，系统结构图符合功能操作要求，且有详尽的功能描述 5)各主要功能模块有详尽的流程图及描述，用函数的形式实现，且调用入口、功 能、出口等关系正确。 6)所设计的C程序较好的可读性和可靠性；程序书写风格规范，注释清晰。 7)界面友好（良好的人机互交），凡需要输入信息时，都有较详尽的提示信息，并 对输入信息有较完善的值域分析。 8)程序经上机调试通过，不应出现语法型、功能型错误，并获得满意的结果。 9)学会通过各功能性测试用例测试程序的能力，尽可能避免原则性错误的出现。 10)掌握编写、调试和测试C语言程序的基本技巧，充分理解结构化程序设计的基 本方法。 11)C语言课程设计报告的内容反应实际设计过程，且符合格式要求，无抄袭现象。 2.2任务描述 （根据具体的题目来写这部分，不要照抄）