매트랩 수치해석 질문드려요

게시물 수정 , 삭제는 로그인 필요

안녕하세요? 

실제 데이터로 다음과 같이 그래프를 얻었습니다.

여기서 이구간에 대해서 적분하고하는데 수치해석으로 하니깐 오차가 심한거 같아서요

실제 사용한 코드는 다음과 같이 사용했습니다.

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

X = xlsread('711.5nF_10KHz_coating_electrode.csv');

Y = xlsread('711.5nF_10KHz_monitoring capacitor.csv');

X1 = xlsread('711.5nF_10KHz_coating_electrode.csv','A3:A1282');

X2 = xlsread('711.5nF_10KHz_coating_electrode.csv','B3:B1282');

Y1 = xlsread('711.5nF_10KHz_monitoring capacitor.csv', 'A3 : A1282');

Y2 = xlsread('711.5nF_10KHz_monitoring capacitor.csv', 'B3 : B1282');

abs_X1 = abs(X1)

time = abs_X1 + X1

find(X2==0) % 0이되는 값찾기

first = 0이 되는 첫번째값

third = 0이 되는 세번째값

range = third - first +1

t1 = time(first)

t2 = time(third)

T = t2-t1

C = 711.5*10^-9

x= linspace(t1,t2,range)

Z = X2.*Y2

Z(first:third,:)

y=Z(first:third, :)*C/T

trapz(x,y)

figure

plot(x,y,'-*')

grid on

title('power')

xlabel('Time(s)')

ylabel('power(w)')

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

실제 1대1 대응되는 데이터 값만 있는데 이런경우는 어떻게 적분을 해야할까요?