PSCAD求方程组

crosssky · 发表于 2016-12-11 19:58:32

crosssky 发表于 2016-12-11 19:502 |1 D; X T* }& U% B4 W6 l' h
楼主脑洞很大，哈哈。。
; Z( u8 X" t0 _. O E8 @& L但是我要说这不是不可能的，世界属于有梦想的人。

怎么发不出来

crosssky · 发表于 2016-12-11 20:06:18

感觉这论文网站有bug发了一个图，整个页面都再加载。

crosssky · 发表于 2016-12-11 20:07:47

ZGGLLJT 发表于 2016-12-11 19:01
' m; X9 R8 v, w, F7 {, _我只是路过打酱油的。第一次听说用PSCAD求解方程组的的

借你的宝地回复一下楼主，我自己发的帖子好像遇到了网站的bug,显示不出来了

楼主脑洞很大，哈哈。。

但是我要说这不是不可能的，世界属于有梦想的人。

接下来让我我来分析一下楼主的方程，楼主没有交代beta是否就是导通角。我假装楼主想告诉我们f()函数是无法求出解析式的，而楼主求的变量是beta, 那么方程变成 beta - f(U,P) =0 ，已知U和P，方程最终化为一个变量beta的方程，但是没有解析形式的解，无法用pscad计算模块直接计算，所以需要解方程，类似 beta^2+sin(beta)=0这种。

那么怎么用pscad求解？？我给楼主提供一下思路，这种方程无法解析求，只能通过非线性方程的数值解法去求解，所以要在pscad里面求解，就需要编写非线性方程求解函数，那么pscad中可以实现非线性方程的求解么，答案是肯定的。首先pscad中有fortran程序自定义功能，我们可以自定义一个模块，输入是U和P的值，仿真中U和P的值自动生成，输入到我们的自定义模块中，我们编写程序求解beta的值，然后作为自定义模块的输出，输出到仿真界面来。这是肯定能实现的，但是好像很麻烦。不用怕，我来提供第二种思路。

思路2就是，我们不用自己编写非线性数值求解函数，因为Matlab中有现成的，我们只需要实现将pscad中的U和P输出到maltab中，利用maltab中数值求解函数解出来beta,然后再返回pscad即可。那怎么实现这两个软件之间的数据传递过程呢？我再提供给楼主一个思路，首先利用pscad中fortran自定义功能,做一个写文件的模块，只需要每个步长将仿真产生的U和P写入一个txt文件，然后Matlab读入这个文件中U和P，利用fsovle等非线性方程求解功能求解出beta,然后再写入另外一个txt文件，pscad再次读入这个文件，用到仿真中去即可。。当然我说的很简单，但是实际上需要两个程序之间配合，pscad再进入下次求解之前要等待maltab返回求解结果，这些配合细节需要去设计。

crosssky · 发表于 2016-12-11 20:08:50

ZGGLLJT 发表于 2016-12-11 19:01* a3 i8 D+ [. Z% H# I3 C
我只是路过打酱油的。第一次听说用PSCAD求解方程组的的

借你的宝地回复一下楼主，我自己的回复好像遇到了网站的bug,显示不出来了

cuicui · 发表于 2016-12-12 10:47:35

crosssky 发表于 2016-12-11 19:58
5 b; g2 n$ d" G0 [( s% W6 }怎么发不出来

我想的也是用MATLAB算

ZGGLLJT · 发表于 2016-12-12 12:57:09

crosssky 发表于 2016-12-11 20:06% h! V8 n# ^1 @1 m9 [
感觉这论文网站有bug发了一个图，整个页面都再加载。

你在自我陶醉吗

cuicui · 发表于 2016-12-20 12:33:48

crosssky 发表于 2016-12-11 19:58' \4 D6 n& K0 s, H
怎么发不出来

我用的MATLAB和pscad交互，不知道这个.m文件对不对,第二张图是函数,I是pscad给MATLAB的输入,α是MATLAB计算出导入到PSCAD中.能帮忙看看么,太感谢了. QQ图片20161220121746.png

ZGGLLJT · 发表于 2016-12-20 13:16:36

淡定，淡定，淡定……

cuicui · 发表于 2016-12-20 15:16:12

cuicui 发表于 2016-12-20 12:333 `+ p: x/ V3 B0 F$ n1 ^
我用的MATLAB和pscad交互，不知道这个.m文件对不对,第二张图是函数,I是pscad给MATLAB的输入,α是MATLAB计 ...

多传了一张图

crosssky · 发表于 2016-12-21 20:57:47

cuicui 发表于 2016-12-20 12:33
& q J' h/ \& t9 E7 D* x0 ]我用的MATLAB和pscad交互，不知道这个.m文件对不对,第二张图是函数,I是pscad给MATLAB的输入,α是MATLAB计 ...

你写的这是什么鬼，，没有这么求方程的。。。。你好好学习下，matlab怎么解方程。。网上那么多例子

账号		自动登录	找回密码
密码			立即加入

PSCAD求方程组

点评

点评

点评