Matematické Fórum

Dobrý den, zde přikládám jeden skriptík a mám problém s vykreslením grafu. Stále se mi nechce vykreslit, a tak se obracím na někoho z vás, kdo by se mi na to mohl mrknout a kdyžtak mi to poopravit. Jsem začátečník v MatLabu, a tak se s tím ještě seznamuji. A pokud mo6no pokud n2kdo ovl8d8 kinematiku kulisov0ho mechanismu kzvn0ho, tak jestli by se mi někdo mohl podívat na průběh výpočtu.
http://www.sdilej.eu/pics/a6b017e87ada2 … 194425.jpg
Předem moc děkuji.


function kulmech3_GUI()
figure('Units','Pixels','Position',[200 90 800 650], 'Color', [1 1 1],'Resize', 'off');
uicontrol('Style','text','Position',[0 550 300 25], 'Backgroundcolor',[1 1 1],'String','Rychlost, zrychlení a poloha kulisového mechanismu kyvného','FontSize',15,'FontWeight','demi','Foregroundcolor','black');

uicontrol('Style','text','Position',[25 500 120 30], 'String', 'Úhlové zrychlení: ', 'FontSize', 12, 'FontWeight', 'demi', 'Backgroundcolor',[1 1 1],'Foregroundcolor','black');
Omega = uicontrol('Style','edit','Position',[145 505 25 25], 'String', 0 ,'FontSize', 12, 'FontWeight', 'demi', 'Backgroundcolor','white','Foregroundcolor','black');

axes('Units','Pixels','Position',[75 190 700 250]);

uicontrol('Style','Pushbutton','Position',[40 450 130 30], 'String', 'Zobraz','FontSize', 12, 'FontWeight', 'demi', 'Backgroundcolor',[0 0.5 1], 'Foregroundcolor','black','Callback', @kulmech3);

function kulmech3 (handle, evendata)
omega = str2double(get(Omega,'String'));
H = 2; % H = a - ze zadání
r = 0.5;
b = 3;
c = 2;
beta1 = 0;%výchozí úhel
pole = zeros(360, 3);

for i = 0 : 1 : 360;
    alfa = omega * i;
    L = sqrt((r * r + H * H) - 2 * r * H * cosd(alfa));
    beta = asind ((r / L) * sind (90 + alfa));
    R = b / cosd(beta);
   
    %rychlost - obvodová rychlost rotační páky s kamenem
    v0 = omega * r;
    va = v0 * sind(alfa + beta);
    vb = va * (R / L);
    v = vb * cosd(beta);
    pole (alfa + 1, 1) = v;

    %zrychleni - úhlová rychlost kyvné páky s kulisou
    omega1 = vb / R;
    anl = omega1 * omega1 * R;
    a = anl * sind (beta);
    pole (alfa + 1, 2) = a;
   
    %poloha táhla vůči výchozí pozici
    h = R * sind(beta1) - R * sind(beta); %výchozí úhel je nulový
    %h = (-1) * R *sind(beta);
    pole (alfa + 1, 3) = h;
end
plot(pole)
grid on
xlabel('osa x')
ylabel('osa y')
%legend('rychlost [m * s^-1]','zrychlení [m * s^-2]','poloha [m]')
end
end