MATLAB简单潮流计算程序？

MATLAB简单潮流计算程序如下：

function lianxuchaoliu

clear

clc

n=9%节点数；

nl=9%支路数；

isb=1%平衡节点号；

pr=0.00001%误猛桐局差精度；

MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中。

MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首轮歼屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的枝让程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

首先你要明白以下几点：

1、.m文件是一个程序计算过程，还是一个function。

2、里面有哪些变量与你的数据有关。

知樱升道上边的脊凳老两个，就按部就班的修改一下，就能用了。

如果是个计算过程，直接运行就能用，在这之前，要确粗枣保变量名都一致，如果不一致，可以修改一下。

如果是个函数，要调用。调用参数也要正确。

希望对你有帮助~~

1．1 基本方法和模型的建立

在小变形条件下，根据力的独立作用原理，无

论载荷多么复杂，都可以将其分解为若干简单载

1 复杂载荷作用F的简支梁

然后应用叠加法得到复杂载荷下梁内的弯矩：如

图1所示为受任意载荷的简支梁。在集中力偶、

集中力和分布载荷单独作用下的弯矩方程分别为

MM( )= RMA-丁+M( 一＆)

R ：一 M (1)

Mp( )= R^P —P( —b)

= 一 P

MfJ( )=R 一 1 q( —c)生二【】： ± 2 L (3)

则在集中力偶M 、集中力P和分布载荷q共同作

用下的弯矩方程为

M ( )= MM( )+MfJ( )+M ( ) (4)

以上各式中应用了跳跃函数，其意义如下

f 0 (I『≤ )

一L『 1( 。， ) ( >)

1．2 计算机分析的实现过程

根据上面介绍的计算模型可应用Matlab编

制如下的计算程序

clear；

L=input( L(In)= )．

M=input( M(KNm)= )．

a input( a(In)= )．

P=input( P(KNm)= )．

b=input( b(In)= )．

q=input( q(KN／m)= )；

c=input( c(In)= )．

d=input( d(In)= )．

nd= 3000；

nf=nd+l：

x=linspace(0，L，nf)；

dx=L／nd；

％ * * * * * * * * * * * * * * * * *

RMA= M／L；nl=a／dx+l；

MM1=RMA X(1：n1)：

MⅣI2=RMA X(nl+l：nf)+M ：

MM=[MMl，MM2 J；

％ * * * * * * * * * * * * * * * * * *

nl=b／dx+1；bb=L—b；

RPA=bb／L*P：

M[Pl=RPA*X(1：n1)；

M口f)2=RPA X(nl+l：nf)

P*(x(nl+l：nf)一b)；

MP=[MP1，MP2 J；

％ * * * * * * * * * * * * * * * * * *

nl=c／dx+l：

n2=d／dx+1；

RqA (L一0．5*(c+d))／L q*(d—c)；

Mql=RqA*X(1：n1)；

Mq2=RqA X(nl+l：n2)一

0．5 q (X(nl+l：n2)一禅祥c)． 2；

Mq3=RqA X(n2+l：nf)一

0．5 q*(X(n2+l：nf) c)．陆世 2

+0．5 q*(X(n2+l：nf)一d)． 2；

Mq=[Mql，Mq2，Mq3]；

‘J／n* * * * * * * * * * * * * * * * *

M = MM +MP+Mq：

subplot(2，l，1)；Mmax=max(M)，Mmin=

rain(M)

plot(X，M)，

title(’复杂早袭肢载荷作用下的弯矩图’)

grid

当L =3 In、M =15 kNm、q=30 kN／m、“

= 0．5 nl、6= l In、f=1．5 In、d =2．5 In时，运

行程序时，得到如图2所示的弯矩图，最大和最小

弯矩分别为Mmax=33．333 kN、Mmin=0=

---