计算方法大作业_计算方法大作业班级X学号X任课老师贺力平(可编辑)doc下载

计算方式大作业开除XXXXXX学号XXXXXXX任课先生贺力平姓名XXX年代试验一幂法与矩阵固有值幂法求主固有值思绪幂法的次要思惟执意对认为的任性初始列矢径效能n次A矩阵(左乘A矩阵)后初始矢径就试图贿赂A矩阵的主固有值对应的表示特性的矢径。鉴于左乘n次A矩阵有能够会使成形计算量逃开因而每回都对列矢径作正态化处置。程序密码效能=MIFA(A,v)UNTITLEDSummaryofthisfunctiongoeshereDetailedexplanationgoeshereA=v=#u(:,)=v(:,)fori=:v(:,i)=A*u(:,i)ifabs(变憔悴(V),i))max(v(:,i)))^breakendu(:,i)=v(:,i)max(v(:,i))enddisp(u(:,i)显示(变憔悴(V),i)))k=u(:,i)#*A*u(:,i)(u(:,i)#*u(:,i))disp(k)x,c=eig(A)disp(c)disp(x)disp(i)end终结匹敌和意见初始矩阵A=初始矢径v=#幂法求得表示特性的矢径正态化后表示特性的矢径列矢径变憔悴相近主固有值Rayleigh商求出主固有值用eig()重大聚会算出的固有值和表示特性的矢径iiiiiiii影响的范围准确问所需次数:意见:可以看出初始列矢径表示方式屡次迭代后用幂法求出的固有值和用eig()重大聚会求出的A的固有值使满足计算准确在并且表示特性的矢径也具无数乘相干。对立的事物发现物取表示特性的矢径v的列变憔悴相近的固有值和用Rayleigh商求出的固有值相当的。试验二QR说明与固有值幂法求主固有值思绪源程序Household重大聚会密码functionQR=hshd(a)UNTITLEDSummaryofthisfunctiongoeshereDetailedexplanationgoesheren=length(a)h=eye(size(a))fori=:nif(a(i,i)sgn=elsesgn=endx=sgn*(总和(a(i),i)^))^u=a(i:n,i)x*眼(ni),)ro=norm(u,)^H=eye(ni)u*u#roH=eye(size(a))H(i:n,i:n)=Ha=H*ah=H*hendQ=inv(h)R=areturnendQR说明密码function=qrfenjie(a)UNTITLEDSummaryofthisfunctiongoeshereDetailedexplanationgoesherefori=:Q,R=hshd(a)a=R*Qenddisp(a)end终结与意见认为a矩阵为列举如下*矩阵a=想出a的固有值为ii三个用QR说明逆序相乘法后说服矩阵三不老实线上若非是实固有值若非是两两配合的虚固有值。在四周由实固有值和*矩阵结合的想象的固有值,t。即使经受住求得的矩阵为上三角矩阵这么固有值执意不老实线上元素并且固有值均为现实。用表示特性的重大聚会求矩阵的固有值:可以看出,这两种算法的固有值在ORD中是相当的。。试验三数值剖析在潮流计算旁边的请求潮流计算的计算图表算法简介潮流计算的计算图表算法是以电用网覆盖学说为根底的请求数值计算方式求解一组叙述供电系统定态表示特性的的方程。从算学上讲是一组多元的非线性相等的求解成绩这类方程的求解进程都离不开迭代。鉴于供电系统构成及决定并发症的一点点表示特性的同时跟随供电系统不竭增强潮流成绩的相等的阶数也越来越高非常的的非线性相等并挑剔若干算学方式都能公约塌下右键的。这种情况就变为使提前或突然发生供电系统计算权杖不竭追求新的且更塌实方式的一点钟要紧并发症。电网潮流计算的机能好丑普通本着的是打算塌实收敛计算职业的兴隆内存占领数量应用无论手巧的敏捷的整齐的和得体的无论轻易无论使满足工程需求等来判别在家以无论塌实收敛作为评价的次要基准。经用的剖析方式有高斯-赛德尔法。、牛顿-拉斐逊潮流算法、迅速解耦算法(PQ说明法)。潮流计算约束供电系统运转不可避免的使满足CERTA。这些问足以使满足P中已确定的变量的约束条件。:应使满足植物的节张力:(1)植物的节有效功率和有功功率:应使满足植物的节间张力的相差:植物的节导纳矩阵的使成形及植物的节张力方程的得体的、互导纳的身体检查意思自导纳在数值上当与该植物的节I坦率地衔接的拿分支导纳的总和。如。交互导纳在数值上当衔接植物的节、分支导纳i消极因素。如。植物的节导纳矩阵Yb为匀称的平方矩阵。植物的节导纳矩阵Yb是稀少矩阵。植物的节导纳矩阵具有不老实优势。坦率地使成形法扩展植物的节导纳矩阵Yb。植物的节导纳矩阵那就够了本着自导纳和互导纳的限制坦率地使成形也干分支mdashmdash植物的节关系矩阵计算。植物的节导纳矩阵的得体的()从原型用网覆盖精炼一分支同时增添一植物的节植物的节导纳矩阵将增添一阶。新的不老实线元素,新的非不老实线元素和不老实线元素。原始用网覆盖打中植物的节、在它们当中添加方式。原始用网覆盖打中植物的节当中切除一分支()原型用网覆盖的植物的节、导纳杂耍自:原用网覆盖植物的节一、j当中改革者的变比由改变为牛顿拉夫逊法(直角坐标)概述牛顿拉夫逊法的意思和衍生物进程把按泰勒系列节目在点大船上的小艇得体的方程 ().牛顿mdash拉夫逊法的表示特性的()牛顿拉夫逊法是迭代法逐步粗略估计的方式()得体的方程是校正方程它的校正进程表现在把非线性方程在按泰勒系列节目大船上的小艇并略去高阶少量()用牛顿mdash拉夫逊法解题时其初值问枯燥的(较试图贿赂真解)若非迭代不收敛。.多变的量非线性方程的解牛顿mdash拉夫逊法得体的方程:()缩写为()潮流计算时的得体的方程(直角坐标)PQ植物的节()PV植物的节()抵消植物的节只设一点钟,张力已知,无对迭代的援用,其张力为()得体的方程()雅可比较的矩阵元,雅可比较的矩阵打中非不老实元素是()。,雅可比较的矩阵中不老实元素为()图牛顿拉夫逊法计算走用MATLAB举行预调牛顿拉夫逊法(直角坐标)某电网接线图及确定的的决定并发症在家,,,PQ植物的节各分支阻抗为抵消植物的节:Z=Z=jZ=Z=jZ=Z=jZ=Z=jZ=Z=jZ=Z=jZ=Z=j各植物的节出力:j::j:j:潮流计算计算图表算法流程图小结数值剖析次要的思绪执意用数值解代表准确的解如下寻觅在必然准确内可以承受的数值算法。这是助长成绩的一点钟要紧方式。在工程上屡次地不问准确的解因而不需求那些的复杂不管到什么程度准确的的算法更多时辰能够也基本无准确的算法。跟随计算图表越来越难以对付的,数值解越来越多。在数值剖析进程中最要紧的思惟是多少设计。它如同与算学的思惟南辕北辙,但它的确表现了。世上无相对相当的的东西,也执意说,这是不能够的。。准确的解和现实境遇脱节请求代价珍奇地而数值剖析的思惟不但可以助长成绩并且也能公约准确问。

发表评论

电子邮件地址不会被公开。 必填项已用*标注