algoritme of Simpson
1 次查看(过去 30 天)
显示 更早的评论
hej,
I got a question about making the algoritme myself, because the funtion i am going to calculate is too big for matlab (or it takes several days to complete).
My problem: I am trying to create a loop in matlab
%Constant
Cpl=1.006;
Cpd=1.86;
Cpw=4.183;
Ro=2500;
Pa=1.01325;
%IMS
Tnb=22;
Ti=34;
Tu=31;
Mw=5;
Ml=2;
%methode van Simpson
n=10; %interval
h=(Ti-Tu)/(2*n);
%hli-formules to make hli
hli=Cpl*Tnb+wTnb*(Ro+Cpd*Tnb);
PsTnb=(6.1121*exp(((18.678-(Tnb/234.5))*Tnb)/(257.14+Tnb)))/1000;
wTnb=0.622*(PsTnb/(Pa-PsTnb));
%formules needed for simpson
Ps=(6.1121*exp(((18.678-(Tw/234.5))*Tw)/(257.14+Tw)))/1000;
w=0.622*(Ps/(Pa-Ps));
hs=Cpl*Tw+w*(Ro+Cpd*Tw);
hl=hli+((Mw*Cpw)/Ml)*(Tw-Tu);
F=(1/(hs-hl));
I wil try to explain the loob: I need "F" for every temperature "Tw" starting with temperature "Ti" and ending with temperature "Tu" and with steps of "h".
For this example:(this is what the loob needs to be doing)
F0=(1/(hs(Ti)-hl(Ti))) --> store in matrix first row
F1=(1/(hs(Ti+h)-hl(Ti+h))) --> store in matrix second row
...
F2n=(1/(hs(Ti+2nh)-hl(Ti+2nh)))=1/(hs(Tu)-hl(Tu)) --> store in matrix last row
So if i run it, it will give me a matrix with all the values of "F"
with these matrix i can calculate the integral by using
y=(h/3)*((F0+F2n)+4*(F1+F3+F5+...+F2n-1)+2*(F2+F4+F6+...+F2n-2))
thanks in advance
Stijn
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!