Здравейте !
Нужна ми е помощ при решаване на интеграл с Матлаб . Случая е следния
Организирал съм следните масиви
Xi се изменя по следния начин Xi=0: L/20:L , където L е число което задавам . Причил съм снимка , където са показани стойностите на X,Y,Z .
DC(Xi,Xj) - дадена ми е формула , от която намирам DC за всяко J и i
DC1=[..] -
DC2=[..]
............
DC20=[..]
DC21=[..]
Z=[..] - задават се стойностите от таблицата
Z приема по-големи стойности от DCj при Z>10.40
Y i
Y1=[..] - задават се стойностите от таблицата .
Y2=[..]
.........
Y20=[..]
Y21=[..]
т.е. за величината DC за всяка стойност на J ще имаме стойност на DC при различно положение Xi .
Пример за J=0
Xi DC
0.00 10.4000
5.69 10.2966
11.38 9.9967
17.06 9.5296
22.75 8.9410
28.44 8.2887
34.13 7.6365
39.81 7.0483
45.50 6.5817
51.19 6.2825
56.88 6.1800
62.56 6.2842
68.25 6.5849
73.94 7.0526
79.63 7.6415
85.31 8.2940
91.00 8.9461
96.69 9.5339
102.38 9.9998
108.06 10.2983
113.75 10.4000
http://picbg.net/img...9efd0a86aff.jpg
http://picbg.net/img...cd664e8278f.jpg
Нужно е да се реши интеграла Omega като се използва метода на трапеците ( ако някой от другите методи за числено интегриране е по-лесен за записване може и с него ) . Все още съм начинаещ с Матлаб и не ми е ясно как да "обвържа" горната граница на интеграла Dc с Z . Нужно ми е да направя алгоритъм при който
да се интегрира У при съответното Xi от 0 до DCj за конкретния случай така че DCj<Z .
Ако имате някакви идеи , примери подобни на моя случай или конкретна литература ( eng/bg ) ще бъда благодарен !
пп - смятам , че ще е нещо от вида
ако за улеснение в примера използвам
Z=[0,1,2,3,4,5,6,7,8,9,10]
DCJ=5.85
Y=[0,1,2,3,4,5,6,7,8,9,10]
то тогава би трябвало за да се "отсеят" нужните У за интегрирането да се използва формула от вида
if(Z(1)<DCJ)
Y(1)
elseif(Z(2)<DCJ)
Y(2)
elseif(Z(3)<DCJ)
Y(3)
elseif(Z(4)<DCJ)
Y(4)
elseif(Z(5)<DCJ)
Y(5)
elseif(Z(6)<DCJ)
Y(6)
elseif(Z(7)<DCJ)
Y(7)
elseif(Z(8)<DCJ)
Y(8)
.. повтаря се до изчерпване на броя Z
end
но дава грешка , неправилно използване на командата if
пп2 - ако някой има идея как процеса да се осъществи в Maple също би ми било доста полезно .