基于MATLAB的断开式转向梯形的优化设计
华北水利水电大学
《汽车设计》课程设计任务书
设计题目：乘用车整车设计
转向系统——转型梯形的优化设计
专业：机械设计制造及其自动化
班级学号：201108207
姓名：刘鹏飞
指导教师：郭朋彦
设计期限：2014年12 月29日开始
2015年1 月9日结束
机械学院
2014年12月26日
一．设计的目的和意义
课程设计题目——乘用车整车设计是针对2011级汽车方向《汽车设计》课程设计而设置的。设置本选题具有以下的目的和意义：
1.通过对轻型乘用车的设计，可以使我们的理论知识更扎实，加深我们对于《汽车构造》、《汽车理论》、《汽车设计》等专业知识的理解，同时使我们学到的理论知识得以应用。
2.在设计的过程中，需要对参考车型的零部件进行了解、分析、设计、建模与装配、验证等这个过程，可以使我们了解产品的研发过程，位我们步入工作岗位，快速适应工作打下良好的基础。
3.本次设计运用三维设计软件CATIA、UG、Pro-E、Solidworks、Solidedge进行建模和仿真，使我们有机会学习和应用目前三维软件领域最为领先的软件的具体操作，了解行业最前沿，同时使用三维软件进行设计可以缩短产品开发周期，提高设计质量。
二、设计参数
1. 加速时间(0—100 km/h)：11.8s；
2. 最小转弯半径：5.3m；
3. 整备质量：1457kg；
4. 满载质量：1940kg；
5. 最高车速：190km/h；
6.外形尺寸（长X宽X高）：4850mmX 1795mmX 1475mm；
7 轴距：2775mm；
8．前轮距：1560mm；
9．后轮距：1560mm；
10. 最小离地间隙：135mm；
11. 行李箱容积：506L；
12. 燃油箱容积：70L；
13．驱动方式：前置前驱，发动机横置；
14：供油方式：多点电喷；
15．发动机排量、燃油、气缸排列型式、进气型式：2000mL、汽油93号(北京92号)、L型、自然吸气式；
16．压缩比、环保标准、缸体材料：10:1、国4、铝合金；
17. 最大功率/转速：100kW/ 5600r/min(rpm)；
18.最大扭矩/转速：190Nm/ 4400r/min(rpm)；
19.百公里油耗：市区工况：9.4L，综合工况：7.8L，市郊工况：6.1L；
20. 转向助力：电子液压；
21. 前制动类型、后制动类型、手刹类型：通风盘、盘式、机械驻车制动；
22．制动距离（(100—0km/h)）：42.27m；
23．前悬挂类型、前麦弗逊式独立悬架；
24.后悬挂类型：后多连杆式独立悬架；
25．轮胎规格：205/65 R16。
三．课程设计内容
1.分组进行、每组课设内容不相同：依据总体设计、离合器设计、变速器设计、转向驱动桥设计、转向从动桥（前桥）、前麦弗逊式悬架设计、后多连杆式悬架设计、制动系统设计、转向系统设计共分成9组进行，要求各组分别进行计算、设计和三维建模，并独立设计出各自的内容，最后能够形成完整的装配图。
2.各组的设计包括以下几部分：
（1）各总成组成部分的形式分析，参数选择；
（2）各总成组成部件的结构设计与计算。
（3）各总成组成零部件三维实体模型的建立与整体装配。
（4）有余力同学建议进行各零件运动学仿真、有限元分析和工程图的生成。
四．设计方法
1.零部件的选型：即根据使用情况，初步确定零部件的形式。
2.参数计算与设计：根据总体设计要求，选择各参数，并进行相关的强度、刚度、疲劳等校核，最终确定出零部件的各个合理参数。该步骤是设计的关键步骤。
3.计算机三维建模：根据理论计算的主要参数，对所设计部件各零件和总成进行三维建模和装配，要遵循三维建模的原则，注意造型细部规划，并按照软件设计小组的要求进行相关格式的转变。该步骤是设计的重要组成部分。
五．设计成果要求
1.设计组成果的提交
（1）设计计算说明书：各组提交本组课程设计说明书一份。说明书要求内容简洁完整、排版整齐、条理清楚、文字通顺、书写规范。含目录、摘要、正文、总结等。要求正文字数应在7000字以上，要求将三维模型抓图后放到说明书中。提交纸质版和电子版。
（2）三维模型：各组需要提交各自总成的三维模型和装配模型（并确保在整车上装配没有问题），总体设计组应提交整车三维模型。提交电子版。
2.每个同学成果的提交
（1）自己的设计说明书部分（要求将三维模型抓图后放到说明书中），提交纸质版和电子版。
（2）自己做的三维模型，提交电子版。
（3）将各自的说明书纸质版打印、装入档案袋；档案袋上写明学号、姓名及所装内容。
五．考核
提交成果时，对个人设计内容及小组设计内容进行现场提问和答辩，包括参数选择，计算方法及建模过程等。
六．课程设计进程表（时间：第18—19周）
目录
一、转向系统结构模型 (6)
二、理想的左右转向轮转角关系 (7)
三、用解析法求实际的内外轮转角关系 (8)
四、转向传动机构的优化设计 (9)
五、三维实体建模 (11)
六、结论与感想 (14)
基于MATLAB的断开式转向梯形机构的优化设计摘要：以与某汽车齿轮齿条式转向器配用的转向传动机构为实例，在传统的转向系设计中引入Solidworks软件建立断开式式转向梯形机构的三维运动实体分析模型，并利用MATLAB中的最小二乘法进行转向梯形的优化设计。通过对比了优化前后的特性曲线，表明优化之后的转向梯形使车轮在转向时左右车轮转角更加符合理论转角关系，从而降低了轮胎磨损，提高的行车平顺性和安全性并更好的保证汽车转弯时车轮作无滑动的纯滚动运动。
引言：转向梯形有整体式和断开式两种，选择整体式还是断开式转向梯形方案与悬架采用何种方案有关。无论采用哪一种方案，都必须正确选择转向梯形参数，做到汽车转弯时，保证全部车轮绕一个瞬时转向中心行驶，使在不同圆周上运动的车轮，作无滑动的纯滚动运动。同时，为达到总体布置要求的最小转弯直径值，转向轮应有足够大转角一、转向系统结构模型
所选用的某汽车转向系统采用齿轮齿条式转向器与麦弗逊独立悬架相匹配的断开式转向梯形机构，其结构原理，如图一转向梯形机构所示。其中梯形臂O1A=OB 的长为L
，横
1
，E 为齿条滑块，长M=637.5mm，齿条E 的许用行程[S]=100mm，轮距拉杆AD=BC 的长为L
2
K=1560mm，轴距L=2775mm，车轮的滚动半径r=334.15mm，主销后倾角准Ψ=2.5°（图中未示出），考虑此角后的计T=L+rtanΨ=2789.58mm，根据最小转弯半径的要求，最大外轮转角α=31.57°
图一 CAD绘制转向梯形机构
通过初步计算得到了转向机构的各项参数后，采用目前在机械产品设计中应用较为广泛的三维设计软件Solidworks建立该机构的运动实体模型。本部分由该组其他成员完成，如图二所示，为梯形机构的运动实体模型。在该模型中，以分段式转向梯形机构取代，如图1 所示的整体式转向梯形机构。如图2 所示的运动实体模型考虑了前轮定位参数，通过建立各部件之间的相互约束和主动件的运动形式，可以模拟分段式转向梯形机构的运动，其性能与真实系统几乎相同，更直观反映转向梯形结构特点。
图二三维实体运动模型
二、理想的左右转向轮转角关系
由转向基本要求可知，在不计轮胎侧偏时，实现转向轮纯滚动、无侧滑转向的条件是内、外轮转角具有如图所示的理想关系，以左转弯为例，即cotα-cotβ=K/L，如图三内轮和外轮转角的关系视图所示。将理想的内轮转角β表示为α的函数即：
β
0=arccot(cotα-
L
K
)=arctan
α
α
tan
1
tan
L
K
-
图三内轮和外轮转角的关系视图
三、用解析法求实际的内外轮转角关系
由转向梯形机构所提供的内、外轮实际转角关系可以根据几何关系来求解。当转动转向盘时，齿条便向左或右移动，使左右两边的杆系产生不同的运动，从而使左右车轮车轮分别获得一个转角。以汽车左转弯为例，此时右轮为外轮，外轮一侧的杆系的杆系运动，如图四左转外轮转角几何关系视图所示。设齿条向右移过某一行程S，通过右横拉杆推动右梯形臂，使之转过α。
图四 左转外轮转角几何关系视图
如图四左转外轮转角几何关系视图所示，取梯形右底角顶点O 为坐标原点，X ，Y 轴如图四 所示，则可导出齿条行程s 与外轮转角a 的关系:
S =
2
M K --L 1cos (α?+)-2
122])sin([e L L -+-α? ?ψγα++=
上式中：=γarctan
S
M K e
22--
ψ=arccos 2
2122
222
1)2
(2)2(e S M K L L e S M K L +---+--+ 而内轮转向的杆系运动与外轮类似。齿条右移了相同的行程s ，通过左横拉杆拉动左梯形臂转过β。
取梯形左底角顶点O ，为坐标原点，X ，Y 轴如图示，则同样可求出齿条行程s 与内轮转角β 的关系，即:
S=L 1cos(β?-)-2
1
22])sin([e L L ---β?-2
M
K - ''γψ?α--=
='
ψarccos
2
2122
2221)2(2)2(
e S M K L L e S M
K L ++--++-+ ='γarctan
S
M K e
22+- 将上述公式结合起来便可将β 表示为α 的函数。 四、转向传动机构的优化设计
1.目标函数的建立
转向机构优化设计的目标就是要在规定的转角范围内使实际的内轮转角尽量的接近对应的理想的内轮转角，采用离散化方法，给出了优化设计目标函数为：
2i 01
i )](αβ)β(α[-=∑=n
i F
当F 取得最小值时，即车轮实际转角与理想值最为接近，优化结果最为理想。由解析法中的公式可设1L 为优化变量x(1)，2L 为优化变量x(2),e 为优化变量x(3)，φ为优化变量x(4)。则，上述优化目标的M 文件可写成：
arf=[0:0.00551:0.551]；
K=1560；%input （‘输入主销中心线间距（mm ）’)；%给与几个变量赋值 M=640；
L=2775；%input(‘输入轴距（mm ）’)；
thetamax=31.57；%input(‘输入外转向轮最大转角（度）’)； thetamax=thetamax*pi/180；%以下将各公式单位转换，并代入公式； B=acot(cot(arf)-K/L)；
x(2)=sqrt((K-M)/2-x(1)*cos(x(4))^2+(x(1)*sinx(4)-x(3))^2)；
s=(K-M)/2-x(1)*cos(x(4)+arf)-sqrt(x(2)^2-(x(1)*sin(x(4)+arf)-x(3))^2； b=x(4)-arccos(x(1)^2+((K-M)/2+s)^2+x(3)^2-x(2)^2)/(2*x(1)*sqrt(((K-M)/2+s )^2+x(3)^2)))-atan(2*x(3)/(K-M+2*s))；
arf=x(4)+arccos(x(1)^2+((K-M)/2-s)^2+x(3)^2-x(2)^2)/2*x(1)*sqrt(((K-M)/2-s)^2)))+atan(2*x(3)/(K-M+2*s))；
s=x(1)*cos(x(4)-b)+sqrt(x(2)^2-(x(1)*sin(x(4)-b)-x(3))^2)-(K-M)/2； fun=abs(B-b)；
优化问题就是求解f=abs(β-0β)的极小值。
五、三维实体建模
1，调用solidworks中的toolbox.
2、插入齿轮
3、与齿轮配合的轴
4、将齿轮与轴配合，配合出装配体
5、另一零件图
首先先进行边线法兰操作，然后生成边线法兰。
最后，生成转向前桥
六、结论与感想
推导出某汽车断开式转向梯形机构的运动学关系式，利用最小二乘法对该转向梯形进行优化设计。研究表明：优化设计后的转向梯形机构明显地提高左右轮转角跟踪理想阿克曼转角的能力。在转向过程中，传动比变化更为合理，保证了车辆做纯滚动，提高了转向性能和轮胎寿命。
通过此次课程设计，我已经基本掌握转向梯形结构的优化过程，和一些基本的三维建模方法，不过还是没有深入了解，还有在优化过程中由于不会使用matlab导致程序只编写了一半，没有完成，关于Solidworks的使用不是很精通，只是画了一些简单的零件，并将其简单的装配。总之，经过这次的课设后，了解到自己的不足之处，我相信在以后的学习中会多加的注意，课下以后，应该多学习一些软件的，为以后走上工作岗位做好准备，以备将来不时之需。
matlab电力系统潮流计算
华中科技大学 信息工程学院课程设计报告书题目: 电力系统潮流计算 专业：电气工程及其自动化 班级： 学号： 学生姓名： 指导教师： 2015年 11 月 10 日
2015年11月12日
信息工程学院课程设计成绩评定表
摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。 关键词：电力系统潮流计算 MATLAB仿真
Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation
复化梯形公式及复化辛普森公式的精度比较
实验四、复化梯形公式和复化Simpson公式的精度比较 （2学时） 一、实验目的与要求 1、熟悉复化Simpson公式和复化梯形公式的构造原理； 2、熟悉并掌握二者的余项表达式； 3、分别求出准确值，复化梯形的近似值，复化Simpson的近似值，并比较后两 者的精度； 4、从余项表达式，即误差曲线，来观察二者的精度，看哪个更接近于准确值。 二、实验内容： 对于函数 sin () x f x x =，试利用下表计算积分1 sin x I dx x =?。 表格如下： 注：分别利用复化梯形公式和复化Simpson公式计算，比较哪个精度更好。其中：积分的准确值0.9460831 I=。 三、实验步骤
1、熟悉理论知识，并编写相应的程序； 2、上机操作，从误差图形上观察误差，并与准确值相比较，看哪个精度更好； 3、得出结论，并整理实验报告。 四、实验注意事项 1、复化梯形公式，程序主体部分： for n=2:10 T(n)=0.5*T(n-1) for i=1:2^(n-2) T(n)=T(n)+(sin((2*i-1)/2^(n-1))/((2*i-1)/2^(n-1)))/2^(n-1)； end end 2、复化Simpson公式，程序主体部分： for i=1:10 n=2.^i x=0:1/n:1 f=sin(x)./x f(1)=1 s=0 for j=1:n/2
s=s+f(2*j) end t=0 for j=1:(n/2-1) t=t+f(2*j-1) end S(i)=1/3/n*(f(1)+4*s+2*t+f(n+1)) end 五．实验内容 复化梯形公式和复化辛普森公式的引入 复化梯形公式： 1 10[(()]2 n n k k k h T f x f x -+==+∑； 复化辛普森公式： 1 1102 [(4()()]6n n k k k k h S f x f x f x -++ ==++∑； 根据题意和复化梯形公式、复化辛普森公式的原理编辑程序求解代码如下： Matlab 代码 clc s=quad('sin(x)./x',0,1) p1=zeros(10,1);
汽车转向梯形机构图解解析
轮式车辆转向梯形结构的图解解析 常州工业技术学院钨华芝常州市政工程管理处魏晓静 摘要介绍几种简单实用的车辆转向梯形结构的图解解析设计法。通过事先设定内、外转向轮实际特性曲线与理论特性的交点位置来控制转角偏差的方法，选择转向梯形机构参数，可以大大减少图次数，提高工作效率，减小转角误差。 关键词：转向梯形机构解析图解 1 引言 轮式车辆一般都是依靠转向车轮偏转一个角度来实现转弯或曲线行驶。转向是的基本要求是保证所有车轮滚动而不发生滑动，这一要求通常由平面四杆机构来达到。传统的设计都采用图解转向梯形的方法。这种方法需要按经验数据选择机构的几何参数，然后作图校核该梯形机构在运动过程中转向轮的转角偏差是否大于允许偏差，若大于允许偏差，则重新选择或调整几何参数，再校核图，直至转角偏转小于允许偏差为止。这实际上是一种试凑的方法，带有较大的盲目性，工作量大。随着计算机的发展，解析法得到了较好的应用，但是传统的图解法仍有它直观、方便的优点，因此仍然被工程设计人员广泛采用。本文介绍一种简单高效且实用的图解解析设计法，可以大大减少作图校核的次数，提高工作效率。 2 转向理论特性 机动车辆或装卸搬运车辆的转向大多采用双轴线式转向方式，见图1。为了满足纯滚动条件，转向时所有车轮必须以不同的半径围绕同一转向中心滚动，各个车轮的轴线交于瞬时转向中心O点。虽然两个转向轮偏转的角度不同，但是两个转角之间应满足下列几何关系： ctg?-ctga=M/L （1）式中?-外轮转角a-内轮转角M-转向轴两主销中心距L-车辆前后轴轴距 为了满足运动学上的这一几何关系，一般都是通过设计转向梯形机构来实现的。式（1）称为转向理论特性。
利用复化梯形公式、复化simpson 公式计算积分
实验 目 的 或 要 求1、利用复化梯形公式、复化simpson 公式计算积分 2、比较计算误差与实际误差 实 验 原 理 （ 算 法 流 程 图 或 者 含 注 释 的 源 代 码 ） 取n=2,3,…,10分别利用复化梯形公式、复化simpson 公式计算积分1 20I x dx =?，并与真值进行比较，并画出计算误差与实际误差之间的曲线。 利用复化梯形公式的程序代码如下： function f=fx(x) f=x.^2; ％首先建立被积函数，以便于计算真实值。 a=0; ％积分下线 b=1; ％积分上线 T=[]; ％用来装不同n 值所计算出的结果 for n=2:10; h=(b-a)/n; ％步长 x=zeros(1,n+1); ％给节点定初值 for i=1:n+1 x(i)=a+(i-1)*h; ％给节点赋值 end y=x.^2; ％给相应节点处的函数值赋值 t=0; for i=1:n t=t+h/2*(y(i)+y(i+1)); ％利用复化梯形公式求值 end T=[T,t]; ％把不同n 值所计算出的结果装入 T 中 end R=ones(1,9)*(-(b-a)/12*h.^ 2*2); ％积分余项（计算误差） true=quad(@fx,0,1); ％积分的真实值 A=T-true; ％计算的值与真实值之差（实际误差） x=linspace(0,1,9); plot(x,A,'r',x,R,'*') ％将计算误差与实际误差用图像画出来 注：由于被积函数是x.^2，它的二阶倒数为2，所以积分余项为：(-(b-a)/12*h.^ 2*2)
实 验 原 理 （ 算 法 流 程 图 或 者 含 注 释 的 源 代 码）利用复化simpson 公式的程序代码如下： 同样首先建立被积函数的函数文件： function f=fx1(x) f=x.^4; a=0; ％积分下线 b=1; ％积分上线 T=[]; ％用来装不同n值所计算出的结果 for n=2:10 h=(b-a)/(2*n); ％步长 x=zeros(1,2*n+1); ％给节点定初值 for i=1:2*n+1 x(i)=a+(i-1)*h; ％给节点赋值 end y=x.^4; ％给相应节点处的函数值赋值 t=0; for i=1:n t=t+h/3*(y(2*i-1)+4*y(2*i)+y(2*i+1)); ％利用复化simpson公式求值end T=[T,t] ; ％把不同n值所计算出的结果装入T中 end R=ones(1,9)*(-(b-a)/180*((b-a)/2).^4*24) ; ％积分余项（计算误差） true=quad(@fx1,0,1); ％积分的真实值 A=T-true; ％计算的值与真实值之差（实际误差） x=linspace(0,1,9); plot(x,A,'r',x,R,'*')
matlab潮流计算
附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵:1、支路首端号;2、末端号;3、支路阻抗;4、支路对地电纳 % 5、支路的变比;6、支路首端处于K侧为1,1侧为0 % B2矩阵:1、该节点发电机功率;2、该节点负荷功率;3、节点电压初始值 % 4、PV节点电压V的给定值;5、节点所接的无功补偿设备的容量 % 6、节点分类标号:1为平衡节点(应为1号节点);2为PQ节点;3为PV 节点; %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1、/(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1、/(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1、/B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵 disp('导纳矩阵Y='); disp(Y) %------------------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部与虚部 for i=1:n %给定各节点初始电压的实部与虚部 e(i)=real(B2(i,3)); f(i)=imag(B2(i,3));
matlab 无约束优化问题
实验八 无约束优化问题 一.实验目的 掌握应用matlab 求解无约束最优化问题的方法 二.实验原理及方法 1：标准形式： 元函数 为其中n R R f X f n R x n →∈:) (min 2．无约束优化问题的基本算法一．最速下降法（共轭梯度法）算法步骤：⑴ 给定初始点 n E X ∈0，允许误差0>ε,令k=0； ⑵ 计算() k X f ?； ⑶ 检验是否满足收敛性的判别准则： () ε≤?k X f ， 若满足，则停止迭代，得点k X X ≈*，否则进行⑷； ⑷ 令() k k X f S -?=，从k X 出发，沿k S 进行一维搜索， 即求k λ使得： ()() k k k k k S X f S X f λλλ+=+≥0 min ； ⑸ 令k k k k S X X λ+=+1，k=k+1返回⑵. 最速下降法是一种最基本的算法，它在最优化方法中占有重要地位.最速下降法的优点是工作量小，存储变量较少,初始点要求不高；缺点是收敛慢，最速下降法适用于寻优过程的前期迭代或作为间插步骤，当接近极值点时，宜选用别种收敛快的算法.．牛顿法算法步骤： (1) 选定初始点n E X ∈0，给定允许误差0>ε，令k=0； (2) 求()k X f ?,()() 1 2-?k X f ,检验：若() ε
基于MATLAB的电力系统潮流计算
基于MATLAB的电力系统潮流计算 %简单潮流计算的小程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号，将变压器的串联阻抗置于低压侧处理。 %第三列为支路的串列阻抗参数。 %第四列为支路的对地导纳参数。 %第五烈为含变压器支路的变压器的变比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，%“0”为不含有变压器。 %B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中 %“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。 %X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地%参数。 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n); Times=1; %置迭代次数为初始值 %创建节点导纳矩阵 for i=1:n1 if B1(i,6)==0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3);
基于MATLAB的潮流计算源程序代码(优.选)
%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;
matlab实现复化梯形公式,复化simpson公式以及romberg积分
（一） 实验目的 熟悉并掌握数值积分的方法，重要训练复化梯形公式，复化simpson 公式以及romberg 积分。 （二） 问题描述 问题三数值积分椭圆周长的计算。考虑椭圆22221x y a b +=,为计算其周长,只要计算其第一象限的长度即可. 用参数方程可以表示为cos (0/2)sin x a t t y b t π=?≤≤?=? , 计算公式为/0π? 为计算方便,我们可以令1a =,即计算下面的积分 / 0π?/0π=? (/0π?/0a π=?可以归结为上面的形式) 采用复化梯形公式，复化Simpson 公式以及Romberg 积分的方法计算积分 / 0()I b π=? 给出通用程序，该通用程序可以计算任何一个函数在任意一个区间在给定的精度下的数值积分。程序输出为计算出的数值积分值以及计算函数值的次数。 （三） 算法介绍 首先利用给出的各迭代公式，设计程序。在matlab 对话框中输入要计算的函数，给出区间和精度。
复化梯形的迭代公式为： ; 复化simpson迭代公式为： ; Romberg迭代公式为： 。 （四）程序 对于复化梯形公式和复化simpson公式，我们放在中。 （%标记后的程序可用来把b看为变量时的算法实现） %复化梯形公式 function y=jifenn(f,n,a,b) （说明：f表示任一函数，n精度，a，b为区间）fi=f(a)+f(b); h=(b-a)/n; d=1; %function f=jifen(n,a,b,c) %syms t %y=sqrt(1+(c^2-1)*cos(t)^2); %ya=subs(y,t,a); %yb=subs(y,t,b); %fi=ya+yb; for i=1:n-1 x=a+i*h; fi=fi+2*f(x); d=d+1; %yx=subs(y,t,x); %fi=fi+2*yx; end f4=h/2*fi,d %复化simposon公式 f1=0; f2=0; dd=1;
汽车转向梯形的优化设计
齿轮齿条式转向梯形的优化设计 学院：车辆与能源学院 专业：2012级车辆工程 学号：S12085234009 姓名：刘建霞 日期：2014年4月15日
齿轮齿条式转向器（如图1）具有结构简单紧凑，制造工艺简便等优点，不仅适用于整体式前轴也适用于前轮采用独立悬架的断开式前轴，目前被广泛地用于轿车、轻型客货车、微型汽车等车辆上。与该转向器相匹配的转向梯形机构与传统的整体式转向梯形机构相比有其特殊之处，下面举一实例加以说明。 图1 齿轮齿条式转向梯形机构运动实体模型 题目：已知某微型汽车（如图2所示）各参数如下：1274.24K mm =， 0()=2.5β主销后倾角，L(轴距)=2340mm ，=mm r （车轮滚动半径）266， =oy B y 梯形臂球头销中心的（）42坐标.12mm ，由最小转弯半径得最大外轮转角为 28o ，许用齿条行程[]62.3S mm =,选用参数624M mm =，试设计转向传动机构。 要求： （1）用优化方法设计此转向梯形传动机构。 （2）优化后校验，压力角40o α≤。 （3）计算出l 1长度，齿条左右移动最大距离。
图2 齿轮齿条转向梯形机构 一 建模 由转向基本要求可知，在不计轮胎侧偏时，实现转向轮纯滚动、无侧滑转向的条件是内、外轮转角符合Arckerman 理想转角关系：cot cot /O i k L θθ-=，如图3所示。 图3 理想的内外轮转角关系 （1）设计变量： 选取变量 1(,,) X l h γ=
图4 外轮一侧杆系运动情况 由图4内外轮转角的关系得： 221o 21o l cos(r )l [sin()h]2 K M S l r θθ-=-+-+- S M K h 22arctan +-=? (2) 2 212 2 2221)2 (2)2(arccos h S M K l l h S M K l ++--++-+=γ (3) i r θφγ=-- (4) 联立上式可得o ()i g θθ=的函数关系式。 对于给定的汽车和选定的转向器，转向梯形机构有横拉杆长l 1和梯形臂长m 两个设计变量。在计算过程中，以梯形底角r 代替横拉杆长l 1作为设计变量，再代入式(1)得到l 1。底角r 可按经验公式先选一个初始值 43r arctan()67.88L K ==，进行优化搜索。 12
基于matlab--psat软件的电力系统潮流计算课程设计
东北电力大学课程设计改革试用任务书: 电力系统潮流计算课程设计任务书 设计名称：电力系统潮流计算课程设计 设计性质：理论计算，计算机仿真与验证 计划学时：两周 一、设计目的 1.培养学生独立分析问题、解决问题的能力； 2.培养学生的工程意识，灵活运用所学知识分析工程问题的能力 3.编制程序或利用电力系统分析计算软件进行电力系统潮流分析。 二、原始资料 1、系统图：IEEE14节点。 2、原始资料：见IEEE14节点标准数据库 三、课程设计基本内容： 1.采用PSAT仿真工具中的潮流计算软件计算系统潮流； 1)熟悉PSAT仿真工具的功能； 2)掌握IEEE标准数据格式内容； 3)将IEEE标准数据转化为PSAT计算数据； 2.分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况， 分析两种方法计算潮流的优缺点； 3.分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。
4.选择以下内容之一进行分析： 1)找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施 的特点，并仿真验证； 2)找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证； 3)找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较 各种措施的特点，并仿真验证； 5.任选以下内容之一作为深入研究：（不做要求） 1)找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有 功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗； 2)找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功 功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗； 3)找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压， 并有效调整发电机端电压，提高该节点电压水平； 四、课程设计成品基本要求： 1.绘制系统潮流图，潮流图应包括： 1)系统网络参数 2)节点电压幅值及相角 3)线路和变压器的首末端有功功率和无功功率 2.撰写设计报告，报告内容应包括以下几点： 1)本次设计的目的和设计的任务； 2)电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点； 3)对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性； 4)找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施； 5)分析各种调整措施的特点并比较它们之间的差异； 6)结论部分以及设计心得； 五、考核形式 1.纪律考核：学生组织出勤情况和工作态度等； 2.书面考核：设计成品的完成质量、撰写水平等； 3.答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩； 4.采用五级评分制：优、良、中、及格、不及格五个等级。
转向梯形机构计算及优化案.doc
转向梯形机构确定、计算及优化 转向梯形有整体式和断开式两种，选择整体式或断开式转向梯形方案与悬架采用何种方案有联系。无论采用哪一种方案，必须正确选择转向梯形参数，做到汽车转弯时，保证全部车轮绕一个瞬时转向中心行驶，使在不同圆周上运动的车轮，作无滑动的纯滚动运动。同时，为达到总体布置要求的最小转弯直径值，转向轮应有足够大的转角。 5.5.1转向梯形结构方案分析 1.整体式转向梯形 图5－14 整体式转向梯形 1—转向横拉杆2—转向梯形臂3—前轴 整体式转向梯形是由转向横拉杆1，转向梯形臂2和汽车前轴3组成，如图5－14所示。其中梯形臂呈收缩状向后延伸。这种方案的优点是结构简单，调整前束容易，制造成本低；主要缺点是一侧转向轮上、下跳动时，会影响另一侧转向轮。 当汽车前悬架采用非独立悬架时，应当采用整体式转向梯形。整体式转向梯形的横拉杆可位于前轴后或前轴前(称为前置梯形)。对于发动机位置低或前轮驱动汽车，常采用前置 梯形。前置梯形的梯形臂必须向前外侧方向延伸，因而会与车轮或制动底板发生干涉，所以在布置上有困难。为了保护横拉杆免遭路面不平物的损伤，横拉杆的位置应尽可能布置得高些，至少不低于前轴高度。 2.断开式转向梯形 转向梯形的横拉杆做成断开的，称之为断开式转向梯形。断开式转向梯形方案之一如图5－15所示。断开式转向梯形的主要优点是它与前轮采用独立悬架相配合，能够保证一侧车轮上、下跳动时，不会影响另一侧车轮；与整体式转向梯形比较，由于杆系、球头增多，所以结构复杂，制造成本高，并且调整前束比较困难。
图5－15断开式转向梯形 横拉杆上断开点的位置与独立悬架形式有关。采用双横臂独立悬架，常用图解法(基于三心定理)确定断开点的位置。其求法如下(图5－16b)： 1)延长B K B 与A K A ，交于立柱AB 的瞬心P 点，由P 点作直线PS 。S 点为转向节臂球销中心在悬架杆件(双横臂)所在平面上的投影。当悬架摇臂的轴线斜置时，应以垂直于摇臂轴的平面作为当量平面进行投影和运动分析。 2)延长直线AB 与B A K K ，交于AB Q 点，连AB PQ 直线。 3)连接S 和B 点，延长直线SB 。 4)作直线BS PQ ，使直线AB PQ 与BS PQ 间夹角等于直线A PK 与PS 间的夹角。当S 点低于A 点时，BS PQ 线应低于AB PQ 线。 5)延长PS 与B BS K Q ，相交于D 点，此D 点便是横拉杆铰接点(断开点)的理想的位 置。
MATLAB下的潮流计算实现-稀疏技术毕业设计
毕业设计（论文）MATLAB下的潮流计算实现-稀疏技术
毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作者签名：日期： 指导教师签名：日期： 使用授权说明 本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：日期：
摘要 电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。因此潮流计算是研究电力系统的一种很重要和很基础的计算。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流问题的方程式阶数越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。 本文旨在于研究潮流计算的牛顿—拉夫逊法的基本原理，在Matlab环境中实现牛顿—拉夫逊法潮流计算的数学模型，程序流程以及编制相应程序，并在程序中融合了节点优化编号和稀疏技术，以提高计算效率。最后用IEEE-3O节点标准测试系统验证所编程序。 关键词：潮流计算Newtom-Raphson法节点优化稀疏技术Matlab ABSTRACT Power flow calculation is fundanmental of analysis. Network reconfiguration,fault management,state estimator etc also need the data of electrial system power flow.There is important significance to develop power flow calculation in allusion to traits of distribution network. This paper introduces the principle of Newtom-Raphson algorithm, which is developed for calculation of power flow calculation ,where zero sequence network is open.With this algorithm,the three-phase load is resolved into positive/negative sequence power and coupling power,thus,decoupling three phase power flow into sequencet component power flow.The power flow can be obtained by just finding the positive sequence power flow and then finding the negative sequent component from the coupling https://www.360docs.net/doc/3c5077347.html,pared with the existing methods,the jacobian matrix with the proposed algorithm is of much lower order,thus substantially reducing the computation burden.The proposed algorithm,together with a reference algorithm,has been simulated on an actual IEEE-30 system using statistic load date.And then it will
最优化方法的Matlab实现(公式(完整版))
第九章最优化方法的MatIab实现 在生活和工作中，人们对于同一个问题往往会提出多个解决方案，并通过各方面的论证从中提取最佳方案。最优化方法就是专门研究如何从多个方案中科学合理地提取出最佳方案的科学。由于优化问题无所不在，目前最优化方法的应用和研究已经深入到了生产和科研的各个领域，如土木工程、机械工程、化学工程、运输调度、生产控制、经济规划、经济管理等，并取得了显著的经济效益和社会效益。 用最优化方法解决最优化问题的技术称为最优化技术，它包含两个方面的内容： 1）建立数学模型即用数学语言来描述最优化问题。模型中的数学关系式反映了最优化问题所要达到的目标和各种约束条件。 2）数学求解数学模型建好以后，选择合理的最优化方法进行求解。 最优化方法的发展很快，现在已经包含有多个分支，如线性规划、整数规划、非线性规划、动态规划、多目标规划等。 9.1 概述 利用Matlab的优化工具箱，可以求解线性规划、非线性规划和多目标规划问题。 具体而言，包括线性、非线性最小化，最大最小化，二次规划，半无限问题，线性、非线性方程（组）的求解，线性、非线性的最小二乘问题。另外，该工具箱还提供了线性、非线性最小化，方程求解，曲线拟合，二次规划等问题中大型课题的求解方法，为优化方法在工程中的实际应用提供了更方便快捷的途径。 9.1.1优化工具箱中的函数 优化工具箱中的函数包括下面几类： 1 ?最小化函数
2.方程求解函数 3.最小—乘（曲线拟合）函数
4?实用函数 5 ?大型方法的演示函数 6.中型方法的演示函数 9.1.3参数设置 利用OPtimSet函数，可以创建和编辑参数结构；利用OPtimget函数，可以获得o PtiOns优化参数。 ? OPtimget 函数 功能：获得OPtiOns优化参数。 语法：
汽车转向梯形机构设计
设计题目：汽车转向梯形机构的设计 班级：机自 xx 姓名： xxx 指导老师： xx 2010年10月10日 西安交通大学
汽车转向梯形机构设计 机自84班李亚敏 08011098 设计要求： (1)设计实现前轮转向梯形机构； (2)转向梯形机构在运动过程中有良好的传力性能。 原始数据： 车型：无菱兴旺，转向节跨距M：1022mm，前轮距D：1222mm，轴距L:1780mm，最小转弯半径R：4500mm。 前言： 汽车转向系统是用来改变或恢复其行驶方向的专设机构，由转向操纵机构、转向器和转向传动机构三部分组成。转向操纵机构主要由方向盘、转向轴、转向管柱等组成：转向器将方向盘的转动变为转向摇臂的摆动或齿条轴的往复运动，并对转向操纵力进行放大的机构：转向传动机构将转向器输出的力和运动传给车轮，并使左右车轮按一定关系进行偏转运动的机构。 设计过程： 一、设计原理简介 1采用转向梯形机构转向的机动车辆，左右转弯时应具有相同的特征，因此左右摇臂是等长的。 2内外侧转向轮偏转角满足无侧滑条件时的关系式为：
cotα?cotβ=M L (1) 3.转向过程中转向梯形机构应满足的方程为 cos(α+α0)=cos(β+β0)?a M cos(β+β0?α? α0)+2a2?b2+M2 2Ma (2) 且 b=M?2acosα0 (3) 代人整理得： cos(α+α0)=?cos(β?α0)+a M cos(β?α?2α0)+2cosα0? 2cos2α0 M +a M (4) 式中αβ为无侧滑状态下梯形臂转角的对应位置，可视为已知。由(1)式算出来，因此，方程中有两个独立的未知量需求解，要梯形臂转角的两个对应位置即两个方程来求解。 4梯形臂转角的两个对应位置的确定 由函数逼近理论确定梯形臂转角的两个对应位置的方程为：αi= qq 2[1?cos2i?1 4 π](i=1,2) (5)式中， qq为外偏转角的最佳范围值，由计算机逐步搜索获得。由汽车的最
Matlab牛拉法潮流计算程序
%本程序的功能是用牛顿——拉夫逊法进行潮流计算 % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点； % 3为PV节点； clear; n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); isb=input('请输入平衡母线节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); % % %--------------------------------------------------- for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧 end %求导纳矩阵 disp('导纳矩阵Y='); disp(Y) %---------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部 e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); %PV节点电压给定模值 end for i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i节点无功补偿量 end %=================================================================== P=real(S);Q=imag(S); %分解出各节点注入的有功和无功功率 ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数IT2
潮流计算(matlab)实例计算
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不
平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图