matpower运行是loadcse错误

matpower运行是loadcse错误,第1张

matpower运行是loadcse错误,表示安装matpower时有错误。

完全卸载干净后重新去官网再进行下载,将下载的压缩包解压,复制放在matlab安装目录中的toolbox文件夹下。打开matlab,将工作路径切换到上一步的matpower文件夹,输入命令:install_matpower,回车。选择Polyspace,回车。测试安装的工具,输入命令:test_matpower,回车。等待1分钟左右,可以看到测试结果为所有工具都安装成功。在matlab主页的“设置路径”中,添加matpower文件夹到路径中,这里一定要选择“添加文件夹”,而不要选择“添加并包含子文件夹”,否则后面使用过程中会报错。

任何直流电动机都有可逆性、即可输入电流是电动机,外力带着转动是发电机。减速比大的电机转速高,完全可做手摇发电机用。什么型号好可现场选用,如欲照明用,选12v电机。注意、理论上可行,实际输出3w电力,人力都不能持久。搞试验可以。

哈哈 还是专门找我的啊 我受宠若惊 好怀念啊 你也搞这个啊,是毕业论文么?

好闲话少说,下面是步骤,当初求助未果,自己瞎搞搞出来的:

以下是我当初论文所需要的命令,希望对你有所帮助。主要是靠设置断点,找出那个雅克比矩阵J,其实就一部而已,具体我记不清了,电脑的matlab被我卸载了,太大而且机器慢,我记得是打开其m文件,然后设置断点,貌似是在左边点一下吧会由“-” 变成“○”。然后在运行runpf(case'9'),会在中间停下来,这时那个“J”就是你要的矩阵。之后的那些命令可能对你没用,我当时做的是对其雅克比矩阵做预处理,从而减少它的迭代次数。如果你也是做这个那爽了,干脆直接借你抄.....哈哈哈哈

1.求标准系统IEEE9节点系统刚开始迭代的雅可比矩阵的条件数

(1)首先要在matlab内的matpower中的m文件设置断点

(2)输入命令runpf(case'9')对其进行牛顿法潮流计算

(3)求该系统的矩阵的条件数,输入命令cond(J)

(其他节点的求解方法与之相同,所以省略,以下各程序命令都以IEEE9节点为例)

2.求标准系统IEEE9节点系统雅可比矩阵的谱图

(1)首先将稀疏矩阵J还原full(J)

(2)求其特征值im=eig(ans)

(3)对其特征根求谱图h=plot(im,'*')

3.对IEEE9节点系统运用矩阵的平衡的预处理方法

(1)首先要在matlab内的matpower中的m文件设置断点

(2)输入命令runpf(case'9')对其进行牛顿法潮流计算

(3)对矩阵进行计算

diag(diag(J))

ans*max(det(J))/diag(det(J))

inv(ans)

ans*J*eye(14,14)

(4)求其条件数cond(ans)

4.对IEEE9节点系统运用不完全LU分解的预处理方法

(1)首先要在matlab内的matpower中的m文件设置断点

(2)输入命令runpf(case'9')对其进行潮流计算

(3)对矩阵进行计算

[L,U]=luinc(J,'0')

A=inv(L)*J*inv(U)

(4)求其条件数cond(A)

5.对IEEE9节点系统运用J的分块对角阵的预处理方法

(1)首先要在matlab内的matpower中的m文件设置断点

(2)输入命令runpf(case'9')对其进行潮流计算

(3)对矩阵进行计算

A=J,再对其分块矩阵J置零

(4)求其条件数cond(A)

6.对IEEE9节点系统运用

(1)首先要在matlab内的matpower中的m文件设置断点

(2)输入命令mpopt=mpoption('PF_ALG',2)

runpf('case9',mpopt)用快速解耦法对其进行潮流计算

(3)对矩阵进行计算

A=Bp

B=Bpp

C=[A,zeros(8,6)zeros(6,8),B]

D=inv(C)

runpf('case9')

J*D*eye(14,14)

A=ans

(4)求其条件数

cond(ans)

7.各种预处理法的作图程序

(1)求计算后矩阵的特征值

im=eig(J’)

(2)对其特征值作图

h=plot(im,'*')

有其他需要请留下qq等联系方式


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/bake/11270750.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-14
下一篇 2023-05-14

发表评论

登录后才能评论

评论列表(0条)

保存