Форум на Исходниках.RU https://forum.sources.ru Форум на Исходниках.RU Форум на Исходниках.RU https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1510688 Wed, 04 Apr 2007 14:49:24 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1510688 Arslan Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509928 Wed, 04 Apr 2007 07:05:22 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509928 Arslan, вот так, например (см. аттач)

Кроме этого (для удобства использования твоей программы) я бы тебе посоветовал сделать еще кое-что... Напиши маленькую функцию, которая будет запрашивать у пользователя число, и заносить его в переменную, но если он сразу нажмет Enter - то значение переменной не изменится... А сами tmp_1, tmp_2, ... опиши не как переменные, а как типизированные константы, и задай им некоторое значение по умолчанию...]]> volvo877 Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509727 Wed, 04 Apr 2007 02:45:59 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509727
    {$N+}
     
    Uses Graph,Crt;
     
    Type type_mayat=record {tip mayatnika}
     
          DlinaNiti:real;
     
          MassaNiti:real;
     
          MassaSharika:real;
     
          RadiusSharika:real;
     
          a_begin,a_new,a_old:real;
     
          J,w0:real;
     
         end;
     
    Var
     
       p_menu:pointer;
      
       Var T, curr_t: Integer;
     
       main_mayat:type_mayat; {Maytnik}
     
       dT,fT,g,b:real;{Shag po vremeni, Polnoe vremya, g, koef triniya}
     
       tmp_1:integer; {svobodnie peremennie}
     
       tmp_2,tmp_3,tmp_4,tmp_5:real;{svobodnie peremennie}
     
       Key:char;
     
       graphtmp2,graphtmp3:integer; {graficheskie peremenie}
     
     
     
    Procedure GraphOn; {vklyuchenie graphiki}
     
    Var
     
       Gd,Gm,ErrorCode:integer;
     
    Begin
     
         Gd:=detect;
     
         initgraph(Gd,Gm,'d:\pascal\bgi');
     
         ErrorCode:=GraphResult;
     
         if ErrorCode <> grOk then
     
         begin
     
          WriteLn('oshibka:');
     
          WriteLn(GraphErrorMsg(ErrorCode));
     
          Halt(1);
     
         end;
     
    End;
     
     
     
    Procedure GraphOff;{Otklu4enie graphiki}
     
    Begin
     
         CloseGraph;
     
    End;
     
     
     
    Procedure mayatnik(DlinaNiti_:integer;
     
      MassaNiti_,MassaSharika_,RadiusSharika_,a_begin_:real);
     
    Begin
     
         with main_mayat do
     
         begin
     
          DlinaNiti:=DlinaNiti_;
     
          MassaNiti:=MassaNiti_;
     
          MassaSharika:=MassaSharika_;
     
          RadiusSharika:=RadiusSharika_;
     
          a_begin:=a_begin_;
     
          a_new:=a_begin_;
     
          J:=(1/3)*MassaNiti*sqr(DlinaNiti)+MassaSharika*(2/5*sqr(RadiusSharika)+
     
      sqr(DlinaNiti));
     
          w0:=sqrt(g*(MassaSharika*DlinaNiti+MassaNiti*DlinaNiti/2)/J);
     
              {vi4islenie sobstvennoi 4astoti}
     
         end;
     
    End;
     
     
     
    Procedure dvizh_mayat; {phizika mayatnika}
     
    Var
     
       costmp:double;
     
    Begin
     
         with main_mayat do
     
         begin
     
          a_old:=a_new;
     
          a_new:=a_begin*cos(w0*fT);{novii ugol mayatnika}
     
          if a_begin>0 then
     
          begin
     
                costmp:=cos(a_begin)+(b*DlinaNiti*sqr(a_new-a_old))/(dT*MassaSharika*g);
     
                {vi4islenie novogo cos ugla amplitudi}
     
                if (costmp<>0) then
     
                  a_begin:=arctan(sqrt(1-sqr(costmp))/costmp);
     
                  {vi4islenie novoi amplitudi}
     
          end
     
          else a_begin:=0;
     
          fT:=fT+dT;
     
         end;
     
    End;
     
     
     
    Procedure ris_mayat; {risovanie mayatnika}
     
    Var
     
       xNow,yNow:real;
     
    Begin
     
         with main_mayat do
     
         begin
     
          xNow:=graphtmp2+DlinaNiti*cos(a_new-pi/2);
     
          yNow:=21-DlinaNiti*sin(a_new-pi/2);
     
          setfillstyle(0,0);
     
          bar(graphtmp2-round(DlinaNiti+RadiusSharika),21,
     
      graphtmp2+round(DlinaNiti+RadiusSharika),21+round(DlinaNiti+RadiusSharika));
     
          setcolor(7);
     
          setlinestyle(0,0,3);
     
          line(graphtmp2,21,round(xNow),round(yNow));
     
          setlinestyle(0,0,1);
     
          setcolor(1);
     
          setfillstyle(1,1);
     
          fillellipse(round(xNow),round(yNow),
     
      round(RadiusSharika),round(RadiusSharika));
     
         end;
     
    End;
     
     
     
    Begin
     
         GraphOFF;
     
         writeln('vvedite dlinu niti (recom. 250)');
     
         readln(tmp_1);
     
         writeln('vvedite massu niti (recom. 20)');
     
         readln(tmp_2);
     
         writeln('vvedite massu sharika (recom. 30)');
     
         readln(tmp_3);
     
         writeln('vvedite radis sharika (recom. 15)');
     
         readln(tmp_4);
     
         writeln('vvedite ugol otklonenia (recom. 0.5)');
     
         readln(tmp_5);
     
         writeln('vvedite shag po vremeni (recom. 0.1)');
     
         readln(dT);
     
         GraphOn;
     
         b:=0.1;
     
         fT:=0;
     
         Key:=#0;
     
         g:=9.8;
     
         GraphOn;
     
         graphtmp2:=round(getmaxx/2);
     
         graphtmp3:=round(getmaxx/3);
     
         cleardevice;
     
         mayatnik(tmp_1,tmp_2,tmp_3,tmp_4,tmp_5);
     
    setcolor(15);
     
          setfillstyle(3,15);
     
          bar(graphtmp3,0,2*graphtmp3,20);
     
          line(graphtmp3,20,2*graphtmp3,20);
     
        repeat
     
            if(main_mayat.a_begin<>0)then
     
            begin
     
              dvizh_mayat;
     
              ris_mayat;
     
              delay(500);
     
     
              inc(T); if T = 5 then begin
     
                T := 0; inc(curr_t);
     
                putpixel(curr_t,
     
                  (getmaxy - 240) + trunc(180 * cos(main_mayat.a_new)),
     
                  white );
     
              end;
     
     
            end;
     
             if keypressed then key:=readkey;
     
         until (key=#27);
     
        GraphOff;
     
    End.
]]> Arslan Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509605 Tue, 03 Apr 2007 18:57:46 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509605
Возможно ты что-то добавлял к тому что есть здесь, или убирал...]]> volvo877 Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509602 Tue, 03 Apr 2007 18:54:54 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=1509602 Arslan Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734584 Tue, 31 May 2005 20:04:23 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734584
    inc(T); if T = 10 then begin { только при каждом 10-м пересчете }
    T := 0; inc(curr_t);  { увеличиваем координату X }
    { и далее выводим точку: }
     
    { с новой координатой по оси времени (ось ОХ) }
    putpixel(curr_t,
     
    { по оси OY - значение cos угла маятника }
      (getmaxy - 240) + trunc(180 * cos(main_mayat.a_new)),
     
    { точка выводится красным цветом }
      red );


Коэффициент 180 при косинусе угла подбирался экспериментально (так же как и 240), чтобы график находился там, где он находится сейчас... Измени 240 на другое число - и график будет отображаться выше или ниже, 180 влияет на амплитуду графика...]]> volvo877 Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734562 Tue, 31 May 2005 19:25:19 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734562 volvo877 можешь пожалуйста обьяснить построение графика зависимости угла отклонения от времени :unsure:]]> Guest Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734237 Tue, 31 May 2005 12:58:59 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734237 Цитата Pourtous @
ИМХО теме пора в смехогрех.
Я так не думаю. Помочь можно.

Цитата Guest @
формулу расстояния я не смогла вывести...т.к. ничего похожего не нашла....
я хочу все таки сделать график зависимости периода от времени...
формулу периода я знаю... T=2*pi*sqrt(J/m*g*l), но это формула без учета радиуса шарика и массы нити....И я не знаю как её записать правильно....
J(- момент инерции) уже выражен через массу шарика, массу нити, длину нити и радиуса шарика...это есть в самом начале моей проги...

подскажите как правильно записать эту формулу?
Задай этот вопрос в разделе "Алгоритмы". Как только выяснишь формулу, построить график не проблема. ;)]]> Romtek Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734215 Tue, 31 May 2005 12:30:36 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734215 Pourtous Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734199 Tue, 31 May 2005 12:19:17 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734199 я хочу все таки сделать график зависимости периода от времени...
формулу периода я знаю... T=2*pi*sqrt(J/m*g*l), но это формула без учета радиуса шарика и массы нити....И я не знаю как её записать правильно....
J(- момент инерции) уже выражен через массу шарика, массу нити, длину нити и радиуса шарика...это есть в самом начале моей проги...

подскажите как правильно записать эту формулу?

пожалуйста помогите, просто о-о-очень срочно надо(завтра уже защита самого проекта :( )]]> Guest Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734084 Tue, 31 May 2005 10:28:46 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734084 Цитата Guest @
и как это сделать
Ну выведи формулу: как ты без неё собираешься график строить?]]> Romtek Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734072 Tue, 31 May 2005 10:19:49 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=734072 Цитата Romtek @
Цитата Guest @
а еще подскажите пожалуйста воообще какие еще можно графики зависимости сделать

Например, зависимость пройденного маятником расстояния от времени.

и как это сделать :blink: :wacko:]]> Guest Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733792 Tue, 31 May 2005 06:25:18 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733792 Цитата Guest @
а еще подскажите пожалуйста воообще какие еще можно графики зависимости сделать

Например, зависимость пройденного маятником расстояния от времени.]]> Romtek Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733677 Tue, 31 May 2005 01:15:25 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733677
Цитата
Можно и со вторым что-нибудь придумать...

я думала сделать график зависимости периода от времни, но по-моему так не получится.....
а еще подскажите пожалуйста воообще какие еще можно графики зависимости сделать :unsure:]]> Guest Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733126 Mon, 30 May 2005 11:47:43 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733126 Цитата Guest @
зависимость угла отклонения от времени

можно попробовать вывести вот так:
         { это для красоты... }
         setcolor(white);
         rectangle(10, getmaxy-100, getmaxx-10, getmaxy-10);
         T := 0; curr_t := 10;
     
        { измени свой Repeat ... Until на этот: }
        repeat
            if(main_mayat.a_begin<>0)then
            begin
              dvizh_mayat;
              ris_mayat;
              delay(500);
     
              inc(T); if T = 10 then begin
                T := 0; inc(curr_t);
                putpixel(curr_t,
                  (getmaxy - 240) + trunc(180 * cos(main_mayat.a_new)),
                  red );
              end;
     
            end;
             if keypressed then key:=readkey;
         until (key=#27);

Ну, и естественно, объяви
    Var T, curr_t: Integer;

Можно и со вторым что-нибудь придумать...

Добавлено:
Еще кое-что... В процедуре ris_mayat у тебя постоянно перерисовывается кусок потолка.
          setcolor(15);
          setfillstyle(3,15);
          bar(graphtmp3,0,2*graphtmp3,20);
          line(graphtmp3,20,2*graphtmp3,20);

По-моему, это - лишнее. Эту часть можно смело перенести в основную программу и поставить сразу после вызова процедуры mayatnik:
         cleardevice;
         mayatnik(tmp_1,tmp_2,tmp_3,tmp_4,tmp_5);
         { Вот сюда }
]]> volvo877 Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733048 Mon, 30 May 2005 10:39:14 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=733048 Цитата Romtek @
Цитата Guest @
но я не знаю как сделать так что бы внизу модели выводились разные графики зависимости (например: зависимость угла отклонения от времени,; я точно не знаю какие еще можно графики построить, т.к. мало в этом соображаю
Ты покажи формулы этих зависимостей, а мы тебе подскажем как строить графики для них. С ними, как раз, проблем не должно быть.

какие формулы зависимостей? на сколько я поняла графики строятся путем перевода значения физической величины значению экранной координаты....по общей формуле:
x=x0+(x1-x0)/(i1-i0)*(i-i0)
y=y0+(y1-y0)/(j1-j0)*(j-j0)
но я сама точно не знаю :unsure: :'(]]> Guest Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=732471 Sun, 29 May 2005 18:25:28 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=732471 Цитата Guest @
но я не знаю как сделать так что бы внизу модели выводились разные графики зависимости (например: зависимость угла отклонения от времени,; я точно не знаю какие еще можно графики построить, т.к. мало в этом соображаю
Ты покажи формулы этих зависимостей, а мы тебе подскажем как строить графики для них. С ними, как раз, проблем не должно быть.

Кстати,
         graphtmp2:=round(getmaxx/2);
         graphtmp3:=round(getmaxx/3);
можно заменить на
         graphtmp2 := GetMaxX div 2;
         graphtmp3 := GetMaxX div 3;
т.к. div выполняет целочисленное деление. ;)]]> Romtek Pascal https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=732296 Sun, 29 May 2005 12:19:51 +0000 построение графиков https://forum.sources.ru/index.php?showtopic=97151&view=findpost&p=732296
Вот сам маятник…
    {$N+}
    Uses Graph,Crt;
    Type type_mayat=record {tip mayatnika}
          DlinaNiti:real;
          MassaNiti:real;
          MassaSharika:real;
          RadiusSharika:real;
          a_begin,a_new,a_old:real;
          J,w0:real;
         end;
    Var
       p_menu:pointer;
       main_mayat:type_mayat; {Maytnik}
       dT,fT,g,b:real;{Shag po vremeni, Polnoe vremya, g, koef triniya}
       tmp_1:integer; {svobodnie peremennie}
       tmp_2,tmp_3,tmp_4,tmp_5:real;{svobodnie peremennie}
       Key:char;
       graphtmp2,graphtmp3:integer; {graficheskie peremenie}
     
     
    Procedure GraphOn; {vklyuchenie graphiki}
    Var
       Gd,Gm,ErrorCode:integer;
    Begin
         Gd:=detect;
         initgraph(Gd,Gm,'');
         ErrorCode:=GraphResult;
         if ErrorCode <> grOk then
         begin
          WriteLn('oshibka:');
          WriteLn(GraphErrorMsg(ErrorCode));
          Halt(1);
         end;
    End;
     
     
    Procedure GraphOff;{Otklu4enie graphiki}
    Begin
         CloseGraph;
    End;
     
     
    Procedure mayatnik(DlinaNiti_:integer;
      MassaNiti_,MassaSharika_,RadiusSharika_,a_begin_:real);
    Begin
         with main_mayat do
         begin
          DlinaNiti:=DlinaNiti_;
          MassaNiti:=MassaNiti_;
          MassaSharika:=MassaSharika_;
          RadiusSharika:=RadiusSharika_;
          a_begin:=a_begin_;
          a_new:=a_begin_;
          J:=(1/3)*MassaNiti*sqr(DlinaNiti)+MassaSharika*(2/5*sqr(RadiusSharika)+
      sqr(DlinaNiti));
          w0:=sqrt(g*(MassaSharika*DlinaNiti+MassaNiti*DlinaNiti/2)/J);
              {vi4islenie sobstvennoi 4astoti}
         end;
    End;
     
     
    Procedure dvizh_mayat; {phizika mayatnika}
    Var
       costmp:double;
    Begin
         with main_mayat do
         begin
          a_old:=a_new;
          a_new:=a_begin*cos(w0*fT);{novii ugol mayatnika}
          if a_begin>0 then
          begin
                costmp:=cos(a_begin)+(b*DlinaNiti*sqr(a_new-a_old))/(dT*MassaSharika*g);
                {vi4islenie novogo cos ugla amplitudi}
                if (costmp<>0) then
                  a_begin:=arctan(sqrt(1-sqr(costmp))/costmp);
                  {vi4islenie novoi amplitudi}
          end
          else a_begin:=0;
          fT:=fT+dT;
         end;
    End;
     
     
    Procedure ris_mayat; {risovanie mayatnika}
    Var
       xNow,yNow:real;
    Begin
         with main_mayat do
         begin
          xNow:=graphtmp2+DlinaNiti*cos(a_new-pi/2);
          yNow:=21-DlinaNiti*sin(a_new-pi/2);
          setcolor(15);
          setfillstyle(3,15);
          bar(graphtmp3,0,2*graphtmp3,20);
          line(graphtmp3,20,2*graphtmp3,20);
          setfillstyle(0,0);
          bar(graphtmp2-round(DlinaNiti+RadiusSharika),21,
      graphtmp2+round(DlinaNiti+RadiusSharika),21+round(DlinaNiti+RadiusSharika));
          setcolor(7);
          setlinestyle(0,0,3);
          line(graphtmp2,21,round(xNow),round(yNow));
          setlinestyle(0,0,1);
          setcolor(1);
          setfillstyle(1,1);
          fillellipse(round(xNow),round(yNow),
      round(RadiusSharika),round(RadiusSharika));
         end;
    End;
     
     
    Begin
         GraphOFF;
         writeln('vvedite dlinu niti (recom. 250)');
         readln(tmp_1);
         writeln('vvedite massu niti (recom. 20)');
         readln(tmp_2);
         writeln('vvedite massu sharika (recom. 30)');
         readln(tmp_3);
         writeln('vvedite radis sharika (recom. 15)');
         readln(tmp_4);
         writeln('vvedite ugol otklonenia (recom. 0.5)');
         readln(tmp_5);
         writeln('vvedite shag po vremeni (recom. 0.1)');
         readln(dT);
         GraphOn;
         b:=0.1;
         fT:=0;
         Key:=#0;
         g:=9.8;
         GraphOn;
         graphtmp2:=round(getmaxx/2);
         graphtmp3:=round(getmaxx/3);
         cleardevice;
         mayatnik(tmp_1,tmp_2,tmp_3,tmp_4,tmp_5);
        repeat
            if(main_mayat.a_begin<>0)then
            begin
              dvizh_mayat;
              ris_mayat;
              delay(500);
            end;
             if keypressed then key:=readkey;
         until (key=#27);
        GraphOff;
    End.


зарание благодарю...]]> Guest Pascal