Matematické Fórum

PanTau · 20. 03. 2016 13:01 — Editoval PanTau (20. 03. 2016 13:02)

Ahoj, pomůže mi prosím někdo přepsat tento kód do javy?

MATLAB:

Skrytý text:

Code:

% simulace drahy naboje ve vetru
% souradny system: xy = podlaha, z = vyska nad podlahou
%
% pocitaji se rychlost v a polohovy vektor x
% soucasne se pro kontrolu pocitaji vcista a xcista (bez pritomnosti vzduchu)

% doba simulace [s]
T = 30;

% casovy krok simulace [s]
deltaT = 0.01;

% rychlost vetru [m/s]
vw = [0, 10, 0];

% gravitacni zrychleni [m/s^2]
g = [0, 0, -10];

% pocatecni rychlost [m/s]
v0 = [100, 0, 100];

% odpor naboje ve vzduchu [1/s]
c = 0.05;

% pocatecni poloha [m]
x0 = [0, 0, 0];

%%%%%%%%%%%%
% simulace
t = (0:deltaT:T);
v = zeros(length(t), 3);
x = zeros(length(t), 3);
vcista = zeros(length(t), 3);
xcista = zeros(length(t), 3);

% definice pocatecni rychlosti a polohy; indexy do pole zacinaji na hodnote 1
v(1, :) = v0;
x(1, :) = x0;
vcista(1, :) = v(1, :);
xcista(1, :) = x(1, :);

for ti = 2:length(t)
  v(ti, :) = v(ti-1, :) + g*deltaT + (vw - v(ti-1, :))*c*deltaT;
  x(ti, :) = x(ti-1, :) + v(ti-1, :)*deltaT;
  vcista(ti, :) = vcista(ti-1, :) + g*deltaT;
  xcista(ti, :) = xcista(ti-1, :) + vcista(ti-1, :)*deltaT;
end

%%%%%%%%%%%%%
% vizualizace
figure(1);
plot(t, x(:, 3), 'r', t, xcista(:, 3), 'b', t, 0*t, 'k');
legend('vyska (ve vzduchu)', 'vyska (ve vakuu)', 'teren');
xlabel('t [s]');
ylabel('z [m]');
title('vyska nad terenem v case');

figure(2);
plot(x(:, 1), x(:, 3), 'r', xcista(:, 1), xcista(:, 3), 'b', xcista(:,1), 0*xcista(:,1), 'k');
legend('trajektorie (ve vzduchu)', 'trajektorie (ve vakuu)', 'teren');
xlabel('x [m]');
ylabel('z [m]');
title('bokorys trajektorie');
axis('equal');

figure(3);
plot(x(:, 1), x(:, 2), 'r', xcista(:, 1), xcista(:, 2), 'b');
legend('trajektorie (ve vzduchu)', 'trajektorie (ve vakuu)');
xlabel('x [m]');
ylabel('y [m]');
title('pudorys trajektorie');
axis('equal');

figure(4);
plot3(x(:, 1), x(:, 2), x(:,3), 'r', xcista(:, 1), xcista(:, 2), xcista(:,3), 'b');
legend('trajektorie (ve vzduchu)', 'trajektorie (ve vakuu)');
xlabel('x [m]');
ylabel('y [m]');
title('trajektorie 3D');
axis('equal');
grid('on');
box on;
zlim([0 max(max(xcista(:,3)))]);
xlim([min(min(xcista(:,1))) max(max(xcista(:,1)))]);
ylim([min(min(xcista(:,2))) max(max(xcista(:,2)))]);

Zvládl jsem tohle, ale dál si nevím rady:

Code:

package test;

public class testovaciTrida {

    public static void main(String[] args) {
        new testovaciTrida();
    }

    public testovaciTrida() {
        // doba simulace [s]
        double T = 30;

        // casovy krok simulace [s]
        double deltaT = 0.01;

        //rychlost vetru [m/s]
        Vector3 vw = new Vector3(0, 10, 0);

        //gravitacni zrychleni [m/s^2]
        Vector3 g = new Vector3(0, 0, -10);

        //pocatecni rychlost [m/s]
        Vector3 v0 = new Vector3(100, 0, 100);

        //odpor naboje ve vzduchu [1/s]
        double c = 0.05;

        //pocatecni poloha [m]
        Vector3 x0 = new Vector3(0, 0, 0);

        //////////////////////////////////////
        // simulace
        Vector3 t = new Vector3(0, deltaT, T);
        Vector3 v = new Vector3(0, 0, 0);
        Vector3 x = new Vector3(0, 0, 0);
        Vector3 vcista = new Vector3(0, 0, 0);
        Vector3 xcista = new Vector3(0, 0, 0);

    }
}

class Vector3 {

    public final double x, y, z;

    public Vector3(double x, double y, double z) {
        this.x = x;
        this.y = y;
        this.z = z;
    }
    
}

Nepotřebuji vykreslit grafy, stačí vypsat výsledky do konzole.

Děkuji za pomoc

PanTau · 21. 03. 2016 11:09

Problém jsem již vyřešil, děkuji za pomoc kaja.marik.

http://forum.matweb.cz/viewtopic.ph … 69#p510469

Matematické Fórum

#1 20. 03. 2016 13:01 — Editoval PanTau (20. 03. 2016 13:02)

Pomoc při psání jednoduchého matlab kódu do javy

Code:

Code:

#2 21. 03. 2016 11:09

Re: Pomoc při psání jednoduchého matlab kódu do javy

Zápatí