40 likes | 196 Views
실습문제 1. 구간 0 ≤ t≤3 에서 다음 미분방정식을 나타내는 함수를 작성하여 y(t) 에 대한 그래프생성. function dy =differentialequation1( t,y ) dy =-y-sin(3*t);. >>[t y]=ode45(@differentialequation1,[0 3],1). 실습문제 2. >> dy =@( time,y )-time*y+4*time/y; >>[t y]=ode45( dy ,[0,2],1). 위의 식에 대한 해를 구하여 그래프로 표현하시오. 실습문제 3.
E N D
실습문제1 구간 0≤t≤3에서다음 미분방정식을 나타내는 함수를 작성하여 y(t)에 대한 그래프생성 function dy=differentialequation1(t,y) dy=-y-sin(3*t); >>[t y]=ode45(@differentialequation1,[0 3],1)
실습문제2 >>dy=@(time,y)-time*y+4*time/y; >>[t y]=ode45(dy,[0,2],1) 위의식에 대한 해를 구하여 그래프로 표현하시오.
실습문제3 구간 0≤t≤50,(y1,y3)에대한 상평면 그래프 생성 function dy=differentialequation2(t,y) dy=[10*(y(2)-y(1)) 28*y(1)-y(2)-y(1)*y(3) y(1)*y(2)-8/3*y(3)]
실습문제4 다음 Roberson 상미분방정식 시스템은 3가지 화학물들 간의 반응을 나타낸다. 이 함수를 작성하고 초기조건 y(1)=1,y(2)=0,y(3)=0일때 구간 0 ≤t≤0에서 미분방정식을 풀어보시오. function dy=reac(t,y) dy=[-0.04*y(1)+1e4*y(2)*y(3) 0.04*y(1)-1e4*y(2)*y(3)-3e7*y(2)^2 3e7*y(2)^2] >> [T,Y] = ode45(@reac,[0 3],[1 0 0]'); >> plot(T,Y(:,2))