The load at each step is a function of lambda, the continuation parameter, which is used to interpolate between the base case load (lambda = 0) and the target case load (lambda = 1). So you can use the following to find the load at the current value of lambda …
Pd_current = Pd_base * (1 - lambda) + Pd_target * lambda … where Pd_base is taken from the PD column of mpc_base and Pd_target from the PD column of mpc_target. The target and base cases and current value of lambda are available to CPF callback functions, if you really need it during the process. But the step-by-step history of lambda is also available after the CPF run is complete in results.cpf.lam. Hope this helps, Ray On Oct 4, 2020, at 12:45 AM, 走下去 <1411608...@qq.com<mailto:1411608...@qq.com>> wrote: Hi,everyone, I want to ask a question, which is that when I run the cpf function in MATPOWER, how can I get the Pload during the process? ------------------ 原始邮件 ------------------ 发件人: "MATPOWER discussion forum" <mirishtha...@gmail.com<mailto:mirishtha...@gmail.com>>; 发送时间: 2020年9月17日(星期四) 晚上9:29 收件人: "MATPOWER discussion forum"<matpowe...@list.cornell.edu<mailto:matpowe...@list.cornell.edu>>; 主题: Re: Optimal power flow with DG connected Hi, It's the file or casefile or M-file/MAT-file on which you are performing opf operations for analysis. e.g. runopf(case9) here case9 is the file he asked for. Thanks Mirish On Wed, Sep 16, 2020 at 6:43 AM Surbhi Aggarwal <surbhiaggarwa...@gmail.com<mailto:surbhiaggarwa...@gmail.com>> wrote: Which file do I need to send you? please specify it. Regards Surbhi Aggarwal On Tue, Sep 15, 2020 at 7:03 PM Ray Daniel Zimmerman <r...@cornell.edu<mailto:r...@cornell.edu>> wrote: Looks like the Jacobian is singular. Feel free to send me the case (off-list) and I’ll have a look. Ray On Sep 14, 2020, at 11:35 PM, Surbhi Aggarwal <surbhiaggarwa...@gmail.com<mailto:surbhiaggarwa...@gmail.com>> wrote: Even after doing everything, this error prompts out. Please help me out with this issue. Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Regards Surbhi Aggarwal On Mon, Sep 14, 2020 at 7:31 PM Ray Daniel Zimmerman <r...@cornell.edu<mailto:r...@cornell.edu>> wrote: You would have to confirm the cause of the non-convergence before drawing any conclusions. You may find FAQ #5<https://matpower.org/doc/faq/#pfconvergence> helpful. Ray On Sep 12, 2020, at 6:43 AM, Surbhi Aggarwal <surbhiaggarwa...@gmail.com<mailto:surbhiaggarwa...@gmail.com>> wrote: Hello If the system does not converge within 10 iterations, can we say that it is one of the most severely affected lines in contingency studies? Or any other explanation for this non-convergence? Regards Surbhi Aggarwal On Sat, Sep 12, 2020 at 10:40 AM Surbhi Aggarwal <surbhiaggarwa...@gmail.com<mailto:surbhiaggarwa...@gmail.com>> wrote: Hello Thank you all for your help and I'll definitely look out for the solutions. Thank you so much Regards Surbhi Aggarwal On Fri, Sep 11, 2020 at 7:16 PM Ray Daniel Zimmerman <r...@cornell.edu<mailto:r...@cornell.edu>> wrote: Hi Surbhi, It looks like you were using the default Newton power flow. To increase the number of iterations, you can increase the pf.nr.max_it option. However, this is rarely necessary. If the Newton method is going to converge, it usually does so in less than the default 10 iterations. If it is not converging within 10 iterations, it is likely that the problem has no solution or that you need a better starting point. Ray On Sep 11, 2020, at 5:48 AM, Mirish Thakur <mirishtha...@gmail.com<mailto:mirishtha...@gmail.com>> wrote: Hi Surbhi, First of all why do you want to increase the number of steps/ iterations in the calculation? Is there any specific need? Actually the number of iterations depends on the power flow optimization tolerance. Maybe you can focus on the mpoption function and read it carefully. You can update something like here: mpopt = mpoption('pf.alg', 'FDXB', 'pf.tol', 1e-4) I'm not 100 % sure about it but it may help you. Thanks. Kind regards Mirish On Fri, Sep 11, 2020 at 8:41 AM Surbhi Aggarwal <surbhiaggarwa...@gmail.com<mailto:surbhiaggarwa...@gmail.com>> wrote: Hello Everyone, I tried your suggestions and it somehow helped me. Now my doubt is " If I want to increase my number of iterations, what steps I have to do"? Please help me out Regards Surbhi Aggarwal On Wed, Sep 9, 2020 at 7:50 PM Ray Daniel Zimmerman <r...@cornell.edu<mailto:r...@cornell.edu>> wrote: I also suspect that the network becomes disconnected by removing certain lines. You can check the connectivity using case_info()<https://matpower.org/docs/ref/matpower7.0/lib/case_info.html> or find_islands()<https://matpower.org/docs/ref/matpower7.0/lib/find_islands.html>. If there are islands, then each island needs to have its own reference bus, and the local generation needs to be able to meet the local load in each island. Ray On Sep 9, 2020, at 4:57 AM, Mirish Thakur <mirishtha...@gmail.com<mailto:mirishtha...@gmail.com>> wrote: Hi Surbhi, Could you please check the continuity of the system before running opf? I mean to say if you run opf for (n-1) criteria the system should be still well connected. There shouldn't be any island within the system itself. Perhaps you can test your model by removing any one parallel line (which you can observe later for analysis) and run the model under (n-1) criteria. This could be the easiest way to test the model. Thank you. Kind regards Mirish On Wed, Sep 9, 2020 at 9:11 AM Surbhi Aggarwal <surbhiaggarwa...@gmail.com<mailto:surbhiaggarwa...@gmail.com>> wrote: Dear All, I am working on optimal power flow to find the contingency (N-1) situation of the system. For some of the line outage, certain output pops out MATPOWER Version 7.0, 20-Jun-2019 -- AC Power Flow (Newton) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Warning: Matrix is singular to working precision. > In mplinsolve (line 75) In newtonpf (line 110) In runpf (line 260) Newton's method power flow (power balance, polar) did not converge in 10 iterations. >>>>> Did NOT converge (7.38 seconds) <<<<< What to do to correct this error. It would be of great help. Regards Surbhi Aggarwal