|
Trapezoidal rule에서와 같이 two point를 연결하는 직선으로, 함수를 근사하여 적분을 하지만 Gauss-Legendre formula는 적분구간이 two point에 closed 돼있지 않고 범위를 확장하여 positive, negative error를 통해 더 정확한 적분 값을 얻을 수 있다.
Gauss-Legendre form
|
|
|
trapezoidal rule]을 이용
It1=0;
It2=0;
k=1.5/0.025; % 사출 지연 시간을 제외한 추력발생 시간/time interval
for i=1:k
It1=It1+(0.5*h)*(Y1(1,i)+Y1(1,i+1)); %total impulse 계산
It2=It2+(0.5*h)*(Y2(1,i)+Y2(1,i+1));
end
T1a=It1/1; %평균 추력을 구함.
T2a=It2/1; %평균 추력을 구함.
mp=0.02; %
|
|
메뉴펼치기
