clear
>> x1=0;
x2=0;

t = [0:0.05:2.95]
k = [0: 1: 59]

A = [1 0.05: -0.65 0.85-0.005.*k];
B = [0 : 0.05*exp(-0.1.*(k+1))];
C = [13 0];
D = [0];

 

----------------

 

위처럼 변수를 지정해준다음에

 

x1(k+1)와 x2(k+1) 의 상태방정식에 관한 식인데요,

 

위 처럼 코드를 실행하면 요런 실행 결과가 나옵니다. (당연하겠지만요 ㅠ)

 

t =

  열 1 ~ 12

         0    0.0500    0.1000    0.1500    0.2000    0.2500    0.3000    0.3500    0.4000    0.4500    0.5000    0.5500

  열 13 ~ 24

    0.6000    0.6500    0.7000    0.7500    0.8000    0.8500    0.9000    0.9500    1.0000    1.0500    1.1000    1.1500

  열 25 ~ 36

    1.2000    1.2500    1.3000    1.3500    1.4000    1.4500    1.5000    1.5500    1.6000    1.6500    1.7000    1.7500

  열 37 ~ 48

    1.8000    1.8500    1.9000    1.9500    2.0000    2.0500    2.1000    2.1500    2.2000    2.2500    2.3000    2.3500

  열 49 ~ 60

    2.4000    2.4500    2.5000    2.5500    2.6000    2.6500    2.7000    2.7500    2.8000    2.8500    2.9000    2.9500

k =

  열 1 ~ 21

     0     1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16    17    18    19    20

  열 22 ~ 42

    21    22    23    24    25    26    27    28    29    30    31    32    33    34    35    36    37    38    39    40    41

  열 43 ~ 60

    42    43    44    45    46    47    48    49    50    51    52    53    54    55    56    57    58    59

 

 

제가 하고싶은건 

K = 0~59 일때의 x1(k+1)값, x2(k+1), y값이 출력되게 하는건데요....

예를 들면 다음과 같이..

 

               x1(k+1)          x2(k+1)      y(k+1)

k  = 0      0.09238098    0.21654         0.654612

k = 1            ~               ~                   ~ 

.

.

.

k= 59         ~                 ~                    ~

 

어떤 함수를 사용해야 위의 표처럼 값이 나올까요..?

3시간째 인터넷 뒤지고 있는데 원하는 것만 딱 얻기가 너무 힘들어서

염치불구하고 질문글 올립니다..

 

혹시나 드리는 말씀이지만 코드를 짜달라는 얘기가 아닙니다 ㅠㅠ

어떤 '함수'를 사용해야 되는지만 알려주시면 감사하겠습니다. 코드는 짜도 제가 짭니다!