完全卸载干净后重新去官网再进行下载,将下载的压缩包解压,复制放在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等联系方式
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)