Hola, no borren este post, me lo vienen pidiendo desde hace mucho la gente de youtube. ¿Quieren colaborar con el aprendizaje o me van a borrar el post?
¿Donde caxo está mi post de MATLAB? Todos los links eran míos.
%Tiro oblicuo con y sin resistencia del viento:
%m, Vi, theta y k, se ingresan desde Command Window poniendo: mruv(dato1,dato2,dato3,dato4)
% m: masa
% Vi: velocidad de arranque
% theta: angulo respecto de la horizontal
% k: constante de proporcionalidad. f=kv (fuerza de rozamiento del aire)
function [] = mruv(m,Vi,theta,k)
clc
g=9.8;
%Este segmento no considera la resistencia del viento
tp=2*Vi*sin(theta)/9.8;
ym=(tp/2)*Vi*sin(theta)-(1/2)*g*(tp/2).^2;
xm=Vi*cos(theta)*tp;
fprintf('velocidad inicial: %.2fnn',Vi);
fprintf('tiempo total en movimiento sin rozamiento: %.2fn',tp);
fprintf('x maximo sin rozamiento: %.2fn', xm);
fprintf('y maximo sin rozamiento: %.2fnn',ym);
axis([0 xm 0 ym])
for c=0:.3:xm
x=linspace(c,c+.5,25);
y=tan(theta)*x-.5*g*(x.^2)/((Vi*cos(theta)).^2);
if y<0
break
end
hold on
pause(.01)
plot(x,y)
end
%Este segmento considera la resistencia del viento:
n=(Vi*sin(theta)+m*g/k);
tm=-(m/k)*log(m*g/(n*k));
xmax=(m/k)*Vi*cos(theta)*(1-exp(-(k/m)*(2*tm)));
ymax=((m/k)*g+Vi*sin(theta))*(m/k)*(1-exp(-(k/m)*tm))-(m/k)*g*tm;
fprintf('tiempo maximo con rozamiento: %.2fn',2*tm);
fprintf('x maximo con rozamiento: %.2fn',xmax);
fprintf('y maximo con rozamiento: %.2fn',ymax);
for c=0:.05:2*tm
t=linspace(c,c+.3,25);
x=(m/k)*Vi*cos(theta)*(1-exp(-(k/m)*t));
y=((m/k)*g+Vi*sin(theta))*(m/k)*(1-exp(-(k/m)*t))-(m/k)*g*t;
hold on
pause(.03)
plot(x,y,'r')
end
end
Yapa: Es una papa finalizar este programa. La idea es que la cámara tome imágenes secuencialmente, encuentre diferencias en ella y te alerte.
function intRectangle()
clc
img=imread('img5.jpg');
imgb=im2bw(img,.35);
[fil col]=size(imgb);
[ dir ]=bwlabel(imgb); acá va todo junto el dir entre corchetes
propied=regionprops(dir);
for n=1:length(propied)
xcen=propied(n).Centroid(1);
ycen=propied(n).Centroid(2);
y=linspace(1,ycen,15);
i=1;
for x=linspace(1,xcen,15)
pause(.01)
imshow(imgb)
rectangle('Position',[x,y(i),col*(1+(1-x)/xcen),fil*(1+(1-y(i))/ycen)],'EdgeColor','r','lineWidth',2)
i=i+1;
end
end
imshow(imgb)
end
Voy volver a subir más cosas de estas.
¿Donde caxo está mi post de MATLAB? Todos los links eran míos.
%Tiro oblicuo con y sin resistencia del viento:
%m, Vi, theta y k, se ingresan desde Command Window poniendo: mruv(dato1,dato2,dato3,dato4)
% m: masa
% Vi: velocidad de arranque
% theta: angulo respecto de la horizontal
% k: constante de proporcionalidad. f=kv (fuerza de rozamiento del aire)
function [] = mruv(m,Vi,theta,k)
clc
g=9.8;
%Este segmento no considera la resistencia del viento
tp=2*Vi*sin(theta)/9.8;
ym=(tp/2)*Vi*sin(theta)-(1/2)*g*(tp/2).^2;
xm=Vi*cos(theta)*tp;
fprintf('velocidad inicial: %.2fnn',Vi);
fprintf('tiempo total en movimiento sin rozamiento: %.2fn',tp);
fprintf('x maximo sin rozamiento: %.2fn', xm);
fprintf('y maximo sin rozamiento: %.2fnn',ym);
axis([0 xm 0 ym])
for c=0:.3:xm
x=linspace(c,c+.5,25);
y=tan(theta)*x-.5*g*(x.^2)/((Vi*cos(theta)).^2);
if y<0
break
end
hold on
pause(.01)
plot(x,y)
end
%Este segmento considera la resistencia del viento:
n=(Vi*sin(theta)+m*g/k);
tm=-(m/k)*log(m*g/(n*k));
xmax=(m/k)*Vi*cos(theta)*(1-exp(-(k/m)*(2*tm)));
ymax=((m/k)*g+Vi*sin(theta))*(m/k)*(1-exp(-(k/m)*tm))-(m/k)*g*tm;
fprintf('tiempo maximo con rozamiento: %.2fn',2*tm);
fprintf('x maximo con rozamiento: %.2fn',xmax);
fprintf('y maximo con rozamiento: %.2fn',ymax);
for c=0:.05:2*tm
t=linspace(c,c+.3,25);
x=(m/k)*Vi*cos(theta)*(1-exp(-(k/m)*t));
y=((m/k)*g+Vi*sin(theta))*(m/k)*(1-exp(-(k/m)*t))-(m/k)*g*t;
hold on
pause(.03)
plot(x,y,'r')
end
end
Yapa: Es una papa finalizar este programa. La idea es que la cámara tome imágenes secuencialmente, encuentre diferencias en ella y te alerte.
function intRectangle()
clc
img=imread('img5.jpg');
imgb=im2bw(img,.35);
[fil col]=size(imgb);
[ dir ]=bwlabel(imgb); acá va todo junto el dir entre corchetes
propied=regionprops(dir);
for n=1:length(propied)
xcen=propied(n).Centroid(1);
ycen=propied(n).Centroid(2);
y=linspace(1,ycen,15);
i=1;
for x=linspace(1,xcen,15)
pause(.01)
imshow(imgb)
rectangle('Position',[x,y(i),col*(1+(1-x)/xcen),fil*(1+(1-y(i))/ycen)],'EdgeColor','r','lineWidth',2)
i=i+1;
end
end
imshow(imgb)
end
Voy volver a subir más cosas de estas.