질문이 좀 깁니다. 훌륭하신 분들의 많은 조언을 구합니다.ㅎ

안녕하세요 고수님들.

원래는 C를 사용하다가 여러가지 유용성을 확인한 후에 Matlab으로 갈아타기를 시도 중인 학생입니다.

현재 Matlab으로 ODE 방정식을 푸는 프로그래밍을 해보고 있는데요.

제가 풀고자 하는 내용이 너무 복잡해서 전체에 관해서 물어보는 것보다 제가 간단히 설명드리고 질문을 하는 것이 나을듯 합니다.

간단히 말씀드리면 세 개의 미분 방정식을 푸는 건데요, 식에서 독립변수 앞의 계수가 상수가 아니고 time step에 따라 바뀌는 구조가 되겠습니다.

식을 간단히 써보면,,

(1) dC/dt=-k*C                                   where,   k=const

      I=k*N*C(t)                                     where,   N=const

(2) dx/dt=beta11*x*y

(3) dy/dt=(beta12*x*y)+I

where, beta11=dc1+dc2

            beta12=dc2+dc2

            dc1=x/y

            dc2=y/x

Initial condition은 다음과 같습니다.

(1) C(0)=0.1

(2) x(0)=0

(3) y(0)=0

구하고자 하는 시간 구간은 0~10 sec이고요.

이렇게 되어 있습니다.

위에 써놓은 방정식은 제가 질문의 요점만 넣어서 간단히 만든식이기 때문에 물리적 혹은 수학적으로는 말도 안되는 구조가 될수도 있겠네요.

요점만 말씀드리면 미분방정식에 들어가 있는 계수가 시간에 따라 변화하는 독립변수에 따라 변화하는 구조가 되겠습니다.

제가 간단히 생각한 m-file은 다음과 같습니다.

기본적으로 m-file을 2개를 만들어야 하는데, (3)번 식에서 (1)에서 구한 I가 필요하기 때문에 m-file을 3개를 만들었습니다.

=========================================================================================

<monomer.m>

function dCdt=monomer(t,C)
global k N I

dCdt=-(k.*C);

<f.m>

function dydt=f(t,y)
global k N I

dc1=y(1)/y(2);

dc2=y(2)/y(1);

beta11=dc1+dc2;

beta12=dc2+dc2;

dydt=[beta11.*y(1).*y(2); beta12.*y(1).*y(2)+I];

<main.m>

clc

global k N I

k=(8.26*10.^4)*exp(-88800/(8.314*1100));      % reation const. [/s]
N=6.023*10^23;                                               % Avogadro Number [#/mole]

tspan=[0 10];

Czero=0.1;

yzero=[0;0];

[t,C]=ode15s(@monomer,tspan,Czero);

I=k.*N.*C;

figure;

plot(t,y);

[t,y]=ode15s(@f,tspan,yzero);

figure;

plot(t,dc1);

hold on;

plot(t,dc2);

=========================================================================================

위와 같이 프로그래밍 해보았는데요, 알 수 없는 error가 뜨네요..

어디가 문제인지 완전!!! 고수님들이 좀 도와주세요....ㅠ(완전 고수 아니셔도 할 수 있는 건가요..?ㅜㅠ)

더불어 한가지 더 궁금한 것은 ode15s를 사용하면 tspan내에서 자동적으로 time step이 결정되어서 계산이 되게 되자나요,,,

그럼 식(1)을 계산하여 나온 I 값이 <f.m> 파일에서 달라진 time step에 대해서 유연하게 들어갈 수 있을까요?

언뜻 생각하기엔 불가능할 거 같은데 말이죠,,,,

이상  matrix기반의 matlab의 개념이 아직 덜 잡힌 왕초보의 질문이었습니다.