интеграции между кривыми

Вот участок с 2-х тенденций на нем:интеграции между кривыми

enter image description here

Я хочу найти площадь под черной кривой (черные линии являются частью одного тренда), но выше синяя кривая.

%w = wavelength array, Tf = filter Transmission, Ts = spectrum Transmission 
figure 
plot(w,Tf,'b') 
hold on 
plot(w,Ts,'k')

источник

2014-10-31 user4204005

Я бы сказал спонтанно: 'дТ = макс (Tf, Ts) -Tf; 'Для quadratur вы можете использовать, например, trapez then: 'area = sum (diff (w). * (dT (1: end-1) + dT (2: end)))/2;' – matheburg

@matheburg Мой мозг жарится сегодня. У вас было решение, но 'trapz' делает вторую половину. :) – chappjc

@chappjc в 'trapz.m' вы, вероятно, найдете линию, очень похожую на мою: P – matheburg

Выберите часть черной линии (Ts), которая находится выше синей линии (Tf) и интеграции (trapz):

d = max(Ts-Tf,0); % or max(Tf,Ts)-Tf; % thanks, matheburg! 
trapz(w,d)

источник

2014-10-31 21:36:31 chappjc

На самом деле это не место, оно также должно зависеть от 'w' (см. Мой комментарий выше). – matheburg

@matheburg True! Ответили слишком поспешно. Как насчет использования 'trapz'? Это работает? – chappjc

Должно быть одинаково ...;) – matheburg

ответ

Смежные вопросы