apocan07
Centipat
- Katılım
- 11 Mayıs 2020
- Mesajlar
- 3
Merhaba öncelikle konuyu nerede açacağımı bilmediğim için yazılımda açtım eğer uygun bir yer varsa taşırsanız sevinirim. Kusura bakmayın.
Bir octave kodum var şu şekilde;
Function f=turbfriction(x)
%X(1)=diameter.
%X(2)=velocity.
%X(3)=reynolds number.
%X(4)=friction.
F =[x(2)-0.35/(3.14*x(1)*x(1)/4);
X(3)-x(2)*x(1)/(1.655*10^(-5));
1/sqrt(x(4))+2.0*log10(2.51/(x(3)*sqrt(x(4))));
x(4)-20*x(1)*2*9.81/(150*x(2)*x(2))];
Options=optimset('tolx',1E-18,'maxfunevals',4000000,'tolfun',1E-18,'maxiter',200000)
[x,fval] = fsolve('turbfriction',[0.1 0.1 100000 0.1 ],options)
Function f=turbfriction(x)
%X(1)=diameter.
%X(2)=velocity.
%X(3)=reynolds number.
%X(4)=friction.
F=[x(2)-0.35/(3.14*x(1)*x(1)/4);
X(3)-x(2)*x(1)/(1.655*10^(-5));
1/sqrt(x(4))+2.0*log10(2.51/(x(3)*sqrt(x(4))));
x(4)-20*x(1)*2*9.81/(150*x(2)*x(2))];
Ama şöyle bir hata alıyorum:
Error: X(2): Out of bound 1 (dimensions are 1x1)
Error: Called from.
Turbfriction at line 8 column 2.
Bir octave kodum var şu şekilde;
Function f=turbfriction(x)
%X(1)=diameter.
%X(2)=velocity.
%X(3)=reynolds number.
%X(4)=friction.
F =[x(2)-0.35/(3.14*x(1)*x(1)/4);
X(3)-x(2)*x(1)/(1.655*10^(-5));
1/sqrt(x(4))+2.0*log10(2.51/(x(3)*sqrt(x(4))));
x(4)-20*x(1)*2*9.81/(150*x(2)*x(2))];
Options=optimset('tolx',1E-18,'maxfunevals',4000000,'tolfun',1E-18,'maxiter',200000)
[x,fval] = fsolve('turbfriction',[0.1 0.1 100000 0.1 ],options)
Function f=turbfriction(x)
%X(1)=diameter.
%X(2)=velocity.
%X(3)=reynolds number.
%X(4)=friction.
F=[x(2)-0.35/(3.14*x(1)*x(1)/4);
X(3)-x(2)*x(1)/(1.655*10^(-5));
1/sqrt(x(4))+2.0*log10(2.51/(x(3)*sqrt(x(4))));
x(4)-20*x(1)*2*9.81/(150*x(2)*x(2))];
Ama şöyle bir hata alıyorum:
Error: X(2): Out of bound 1 (dimensions are 1x1)
Error: Called from.
Turbfriction at line 8 column 2.
Son düzenleyen: Moderatör:
