MATLAB: How coding drag coefficient

dragfor loop

Hi all, I need kindly help to write a code to calculate drag coefficient as the attached figure.

thank you in advance

Best Answer

See if this does what you want:

tau = @(U10) 1.22.*[1.2E-3.*((4 < U10) & (U10 < 11)) + 1E-3*(0.49 + 0.065*U10).*((11 <= U10) & (U10 < 25))].*U10.^2;
x = linspace(0, 24.9);
figure(1)
plot(x, tau(x))
xlabel('\itU_{10}\rm')
ylabel('\bf\tau\rm')
grid

Related Solutions

MATLAB: There is the error that “insufficient number of inputs”

Try :
integral2(@Tuu,0,pi/2,0,pi/4) % calling if the function
function U=Tuu(x,y)
d=2e-9;
mu=8;
delta=10;
vh=16;
HBAR=1.05457266e-34;
ME=9.1093897e-31;
ELEC=1.60217733e-19;
Kh=2.116e10;
kc=sqrt(2.*ME.*ELEC./HBAR.^2);
ku=kc.*sqrt(mu+delta);
kd=kc.*sqrt(mu-delta);
puu1=sqrt(ku.^2-ku.^2.*sin(x).^2+kc.^2.*vh);
puu2=sqrt(ku.^2-ku.^2.*sin(x).^2-kc.^2.*vh);
quu1=sqrt(ku.^2-(Kh-ku.*sin(x).*cos(y)).^2-ku.^2.*sin(x).^2.*sin(y).^2+kc.^2.*vh);
quu2=sqrt(ku.^2-(Kh-ku.*sin(x).*cos(y)).^2-ku.^2.*sin(x).^2.*sin(y).^2-kc.^2.*vh);
U=1/4.*ku.*sin(x).*(real(puu1).*exp(-2.*imag(puu1).*d)+real(puu2).*exp(-2.*imag(puu2).*d)+real(quu1).*exp(-2.*imag(quu1).*d)+real(quu2).*exp(-2.*imag(quu2).*d)+((real(puu1)+real(puu2)).*cos((real(puu1)-real(puu2)).*d)-(imag(puu1)-imag(puu2)).*sin((real(puu1)-real(puu2)).*d)).*exp(-(imag(puu1)+imag(puu2)).*d)-((real(quu1)+real(quu2)).*cos((real(quu1)-real(quu2)).*d)-(real(quu1)-imag(quu2)).*sin((real(quu1)-real(quu2)).*d)).*exp(-(imag(quu1)+imag(quu2)).*d));
end

MATLAB: Solving an equation in Matlab

These have analytic solutions:

syms  u(t) v(t) f taux tauy Hm rho R u0 v0 
du = diff(u);
dv = diff(v);
DE = [du - f*v == taux/(Hm*rho) - R*u;  dv - f*u == taux/(Hm*rho) - R*v]
Soln = dsolve(DE, u(0)==u0, v(0)==v0)
u_Soln = Soln.u;
v_Soln = Soln.v;
u_Soln = simplify(u_Soln, 'Steps',500)
v_Soln = simplify(v_Soln, 'Steps',500)
u_fcn = matlabFunction(u_Soln)
v_fcn = matlabFunction(v_Soln)

Note that ‘u_fcn’ and ‘v_fcn’ are anonymous functions you can use outside of the Symbolic Math Toolbox.

Best Answer

Related Solutions

MATLAB: There is the error that “insufficient number of inputs”

MATLAB: Solving an equation in Matlab

Related Question