목차
1.개요
2.항공기 비행조건
3.프로그램 설명
4.결과 예측 및 특성 파악
5.프로그램 코드[MATLAB]
2.항공기 비행조건
3.프로그램 설명
4.결과 예측 및 특성 파악
5.프로그램 코드[MATLAB]
본문내용
j-50)+12000;
elseif j>=75 && j<150
Th(j,1)=17000;
elseif j>=150 && j<175
Th(j,1)=-300*(j-150)+17000;
elseif j>=175 && j<250
Th(j,1)=100/3*(j-175)+9500;
elseif j>=250 && j<302
Th(j,1)=12000;
end
L=n*w;
cl=(n*w)/(0.5*rho*v(j)^2*s);
cd=0.024+0.12*cl^2;
D=(0.5*rho*v(j)^2*s)*cd;
T=Th(j,1);
[t, y]=ode45(\'newton_eqn\', [j-1:0.5:j], y_0(j,:));
y_0(j+1,:)=y(1,:);
v(j+1)=y(1,1);
tf(j+1,1)=t(3,1);
yf(j+1,:)=y(3,:);
h(j+1)=v(j)*sin(y(3,3))+h(j);
end
figure(3)
subplot(311)
plot(tf,Th); xlabel(\'Time [sec]\');title(\'Time reponse\'); ylabel(\'T\'); grid;
subplot(312)
plot(tf,yf(:,1));xlabel(\'Time [sec]\'); ylabel(\'Velocity [m/s]\');grid;
subplot(313)
plot(tf,h);xlabel(\'Time [sec]\');ylabel(\'height [m]\');grid;
%추력, 롤각도, 하중배수 3가지를 모두 한꺼번에 바뀔 경우
for j=1:1:301
n=N(j,1);
L=n*w;
cl=(n*w)/(0.5*rho*v(j)^2*s);
cd=0.024+0.12*cl^2;
D=(0.5*rho*v(j)^2*s)*cd;
Phi=phi(j)*pi/180;
T=Th(j,1);
[t, y]=ode45(\'newton_eqn\', [j-1:0.5:j], y_0(j,:));
y_0(j+1,:)=y(1,:);
v(j+1)=y(3,1);
tf(j+1,1)=t(3,1);
yf(j+1,:)=y(3,:);
h(j+1)=v(j)*sin(y(3,3))+h(j);
X(j+1)=v(j)*cos(y(3,3))*cos(y(3,2))+X(j); % X방향 거리를 구하기 위한 배열
Y(j+1)=v(j)*cos(y(3,3))*sin(y(3,2))+Y(j); % Y방향 거리를 구하기 위한 배열
end
figure(4)
subplot(311)
plot(tf,Th,\'y-^\',tf,N*10000,\'b-+\',tf,phi*500,\'r--\'); xlabel(\'Time [sec]\');title(\'Time reponse\'); grid;
legend(\'Thrust\', \'n * 10000\', \'phi *500\');
subplot(312)
plot(tf,yf(:,1));xlabel(\'Time [sec]\'); ylabel(\'Velocity [m/s]\');grid;
subplot(313)
plot(tf,h);xlabel(\'Time [sec]\');ylabel(\'height [m]\');grid;
figure(5)
plot3(X,Y,h);xlabel(\'X\'); ylabel(\'Y\');zlabel(\'Z\'); grid;
%%%%%% newton_eqn.m %%%
function dy=newton_eqn(t,y)
global g
global w
global Phi;
global L;
global D;
global T;
v_t=y(1);
psi=y(2);
gamma=y(3);
dy(1)=g*(((T-D)/w)-sin(gamma));
dy(2)=g/v_t*L/w*sin(Phi)/cos(gamma);
dy(3)=g/v_t*(L*cos(Phi)/w-cos(gamma));
dy=dy\';
elseif j>=75 && j<150
Th(j,1)=17000;
elseif j>=150 && j<175
Th(j,1)=-300*(j-150)+17000;
elseif j>=175 && j<250
Th(j,1)=100/3*(j-175)+9500;
elseif j>=250 && j<302
Th(j,1)=12000;
end
L=n*w;
cl=(n*w)/(0.5*rho*v(j)^2*s);
cd=0.024+0.12*cl^2;
D=(0.5*rho*v(j)^2*s)*cd;
T=Th(j,1);
[t, y]=ode45(\'newton_eqn\', [j-1:0.5:j], y_0(j,:));
y_0(j+1,:)=y(1,:);
v(j+1)=y(1,1);
tf(j+1,1)=t(3,1);
yf(j+1,:)=y(3,:);
h(j+1)=v(j)*sin(y(3,3))+h(j);
end
figure(3)
subplot(311)
plot(tf,Th); xlabel(\'Time [sec]\');title(\'Time reponse\'); ylabel(\'T\'); grid;
subplot(312)
plot(tf,yf(:,1));xlabel(\'Time [sec]\'); ylabel(\'Velocity [m/s]\');grid;
subplot(313)
plot(tf,h);xlabel(\'Time [sec]\');ylabel(\'height [m]\');grid;
%추력, 롤각도, 하중배수 3가지를 모두 한꺼번에 바뀔 경우
for j=1:1:301
n=N(j,1);
L=n*w;
cl=(n*w)/(0.5*rho*v(j)^2*s);
cd=0.024+0.12*cl^2;
D=(0.5*rho*v(j)^2*s)*cd;
Phi=phi(j)*pi/180;
T=Th(j,1);
[t, y]=ode45(\'newton_eqn\', [j-1:0.5:j], y_0(j,:));
y_0(j+1,:)=y(1,:);
v(j+1)=y(3,1);
tf(j+1,1)=t(3,1);
yf(j+1,:)=y(3,:);
h(j+1)=v(j)*sin(y(3,3))+h(j);
X(j+1)=v(j)*cos(y(3,3))*cos(y(3,2))+X(j); % X방향 거리를 구하기 위한 배열
Y(j+1)=v(j)*cos(y(3,3))*sin(y(3,2))+Y(j); % Y방향 거리를 구하기 위한 배열
end
figure(4)
subplot(311)
plot(tf,Th,\'y-^\',tf,N*10000,\'b-+\',tf,phi*500,\'r--\'); xlabel(\'Time [sec]\');title(\'Time reponse\'); grid;
legend(\'Thrust\', \'n * 10000\', \'phi *500\');
subplot(312)
plot(tf,yf(:,1));xlabel(\'Time [sec]\'); ylabel(\'Velocity [m/s]\');grid;
subplot(313)
plot(tf,h);xlabel(\'Time [sec]\');ylabel(\'height [m]\');grid;
figure(5)
plot3(X,Y,h);xlabel(\'X\'); ylabel(\'Y\');zlabel(\'Z\'); grid;
%%%%%% newton_eqn.m %%%
function dy=newton_eqn(t,y)
global g
global w
global Phi;
global L;
global D;
global T;
v_t=y(1);
psi=y(2);
gamma=y(3);
dy(1)=g*(((T-D)/w)-sin(gamma));
dy(2)=g/v_t*L/w*sin(Phi)/cos(gamma);
dy(3)=g/v_t*(L*cos(Phi)/w-cos(gamma));
dy=dy\';
추천자료
6시그마의 서비스업에 대한 적용에 관한 연구
청소년 약물중독 문제- 영화 <바스켓볼다이어리>를 보고
24~36개월 영아의 발달적 특징과 교육적 적용
[독점자본주의][국가독점자본주의][역사적 자본주의][지력자본주의]독점자본주의, 국가독점자...
되치기본
테니스(tennis) 포핸드(Forehand)
[인터넷방송][방송정책]인터넷방송의 정의, 인터넷방송의 특징, 인터넷방송의 종류, 인터넷방...
[경영학]대한항공 기업분석
[스쿼시] 스쿼시의 경기방식과 규칙 및 기술과 전략
경영ㆍ회계의 유비쿼터스 적용
학교생활지도의 의의, 학교생활지도의 목표, 학교생활지도의 계획, 학교생활지도의 기초생활...
대한항공 기업분석
소개글