Have you looked at the documentation
<http://www.pserc.cornell.edu//matpower/docs/ref/> for the case_info()
<http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.0/case_info.html>,
find_islands()
<http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.0/find_islands.html>,
and extract_islands()
<http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.0/extract_islands.html>
functions? Btw, when using extract_islands()
<http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.0/extract_islands.html>
you will need to make sure that the island you extract has a reference bus.
Ray
> On Feb 19, 2015, at 7:01 PM, Bijay Hughes <[email protected]> wrote:
>
> Dear All,
>
> I am modeling blackout in matpower. I have some success as long as no
> isolated buses are formed. However when I have isolated buses, I receive the
> following command error:
>
>
> Index exceeds matrix dimensions.
>
> Error in dcpf (line 46)
> Va([pv; pq]) = B([pv; pq], [pv; pq]) \ (Pbus([pv; pq]) - B([pv; pq], ref) *
> Va0(ref));
>
> Error in runpf (line 161)
> Va = dcpf(B, Pbus, Va0, ref, pv, pq);
>
> Error in rundcpf (line 93)
> [varargout{1:nargout}] = runpf(casedata, mpopt, fname, solvedcase);
>
> Error in untitled3 (line 70)
> results=rundcpf(mpc);
>
> I model blackout by removing the attacked transmission lines in each loop.
> After some loops islands are formed. Then I cannot go further. I tried using
> the followoing command: extract_islands(mpc,1) but still gives the same
> command error. How do I solve this problem? How do I find out isolated buses
> and set its index to 4 (meaning non-used)? Any help is appreciated.
>
> best,
>
> Bijay
>
