guilda/pss.m at main · moris512/guilda · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
classdef pss < handle
% モデル ：PSSの実装モデル
%         発電機モデルに付加するために実装されたクラス
%親クラス：handleクラス
%実行方法：obj = pss(parameter)
%　引数　：parameter : table型．「'Kpss','Tpss','TL1p','TL1','TL2p','TL2'」を列名として定義
%　出力　：pssクラスのインスタンス

    properties
        nx
        A
        B
        C
        D
        sys
    end

    methods
        function obj = pss(pss_in)
            if nargin < 1 || isempty(pss_in)
                pss_in = ss(0);
            end
            obj.set_pss(pss_in);
            sys = ss(obj.A, obj.B, obj.C, obj.D);
            sys.InputGroup.omega = 1;
            sys.OutputGroup.v_pss = 1;
            obj.sys = sys;
        end

        function nx = get_nx(obj)
           nx = obj.nx;
        end

        function name_tag = get_state_name(obj)
            name_tag = cell(1,obj.get_nx);
            if obj.get_nx ~= 0
                for i = 1:obj.get_nx
                    name_tag{i} = ['xi',num2str(i)];
                end
            end
        end

        function [dx, u] = get_u(obj, x_pss, omega)
            dx = obj.A*x_pss + obj.B*omega;
            u = obj.C*x_pss + obj.D*omega;
        end

        function x = initialize(obj)
           x = zeros(obj.get_nx, 1);
        end

        function sys = get_sys(obj)
            sys = obj.sys;
        end

        function set_pss(obj, pss)
            if istable(pss)
                Kpss = pss{:, 'Kpss'};
                Tpss = pss{:, 'Tpss'};
                TL1p = pss{:, 'TL1p'};
                TL1 = pss{:, 'TL1'};
                TL2p = pss{:, 'TL2p'};
                TL2 = pss{:, 'TL2'};
                obj.A = [
                    -1/Tpss, 0, 0;
                    -Kpss*(1-TL1p/TL1)/(Tpss*TL1), -1/TL1, 0;
                    -(Kpss*TL1p)*(1-TL2p/TL2)/(Tpss*TL1*TL2), (1-TL2p/TL2)/TL2, -1/TL2
                    ];
                obj.B = [
                    1/Tpss;
                    Kpss*(1-TL1p/TL1)/(Tpss*TL1);
                    (Kpss*TL1p)*(1-TL2p/TL2)/(Tpss*TL1*TL2)
                    ];
                obj.C = [-(Kpss*TL1p*TL2p)/(Tpss*TL1*TL2), TL2p/TL2, 1];
                obj.D = (Kpss*TL1p*TL2p)/(Tpss*TL1*TL2);
            else
                [obj.A, obj.B, obj.C, obj.D] = ssdata(pss);
            end
            obj.nx = size(obj.A, 1);
        end
    end
end