Лабораторная работа Вывод графиков функций в Delphi

Для построения графиков функций используем метод «как в школе учили». Строить график будем так же , как на бумаге. Тем самым оставим те же недостатки «бумажного» метода построения, и даже усугубим их.

Рассмотрим пример программы приведенный ниже.

unit Graf;

interface

uses

  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms, Dialogs;

type

  TForm1 = class(TForm)

  procedure FormPaint(Sender: TObject);

  procedure FormResize(Sender: TObject);

  private

  { Private declarations }

  public

  { Public declarations }

  end;

var

  Form1: TForm1;

implementation

{$R *.dfm}

 Function f(x:real):real;

begin

f:=2*Sin(x)*exp(x/5); // наша функция

end;

// строит график функции

procedure GrOfFunc;

var

x1,x2:real; // границы изменения аргумента функции

y1,y2:real; // границы изменения значения функции

x:real; //аргумент функции

y:real; // значение функции в точке х

dx:real; // приращение аргумента

l,b:integer; // левый нижний угол области вывода графика

w,h:integer; // ширина и высота области вывода графика

mx,my:real; // масштаб по осям X и Y 

x0,y0:integer; // точка - начало координат

begin // область вывода графика

l:=10; // X - координата левого верхнего угла

b:=Form1.ClientHeight-20; //У - координата левого верхнего угла

h:=Form1.ClientHeight-40; // высота

w:=Form1.Width-40; // ширина

x1:=0; // нижняя граница диапазона аргумента

x2:=25; // верхняя граница диапазона аргумента

dx:=0.01; // шаг аргумента

// найдем максимальное и минимальное значения функции на отрезке [x1,x2]

y1:=f(x1); // минимум

y2:=f(x1); //максимум

x:=x1;

repeat

y := f (x);

if y < y1 then y1:=y;

if y > y2 then y2:=y;

x:=x+dx; until (x >= x2);

// вычислим масштаб

my:=h/abs(y2-y1); // масштаб по оси Y

mx:=w/abs(x2-x1); // масштаб по оси X

x0:=1;

y0:=b-Abs(Round(y1*my)) ;

with form1.Canvas do

begin // оси

MoveTo(l,b);LineTo(l,b-h);

MoveTo(x0,y0);LineTo(x0+w,y0);

TextOut(l+5,b-h,FloatToStrF(y2,ffGeneral,6,3));

TextOut(l+5,b,FloatToStrF(y1,ffGeneral,6,3));

// построение графика

x:=x1; repeat

y:=f(x);

Pixels[x0+Round(x*mx),y0-Round(y*my)]:=clRed;

x:=x+dx;

until (x >= x2);

end;

end;

procedure TForm1.FormPaint(Sender: TObject); 

begin

GrOfFunc;

end;

procedure TForm1.FormResize(Sender: TObject); 

begin

// очистить форму

form1.Canvas.FillRect(Rect(0,0,ClientWidth,ClientHeight));

// построить график

GrOfFunc;

end;

end.

Программа хорошо выводит графики, когда функция имеет как положительные, так и отрицательные значения. Причем весь график помещается в указанном прямоугольнике.

Недостатки:

Теперь посмотрим, а сколько вычислений значений функции делает программа? В данном случае (25-0)/0.01=2500. Для любого прямоугольника вывода. Чем был обусловлен выбор шага dx? Скорее всего, непрерывностью линии графика. Который, кстати, так и остался прерывистым на некоторых участках, там, где функция меняется быстро. Для борьбы этим напрашивается решение уменьшение значения dx, например - в 100 раз, доводя dx до 0.0001. это приведет к возрастанию вычислений функции до 250000. Но график все равно получиться прерывистые. Причем большая часть результатов вычислений просто пропадут в результате округления Round . Отсюда вытекает, что функцию нужно считать в количестве точек равных размерам окна по горизонтали, а отрезки вертикальных прямых можно нарисовать карандашом. И dx нужно выбирать например (25-0)/600= 0,0416666. График получится самый качественный, какой только возможно получить.

Вторым недостатком является сам метод построения (вычисление значений функции с шагом dx) работает как фильтр, отсекая высокочастотные гармоники, т.е если к функции f(x) добавить что-то вроде g(x)*sin(2*pi/dx*x), то результат вывода будет плачевным. Этот элемент никак не изменит предыдущий график. Хотя он может являться основным носителем информации о функции. И уж конечно очень непросто вывести на экран график дискретной функции (имеется в ввиду универсальными программами общего пользования, подобными приведенной). Если взять f(x)=2*Sin(x)*exp(x/5)+ exp(x*x)*sin(2*pi/dx*x), то данная программа второе слагаемое не заметит, хотя будет тратить время на расчет f(x)=2*Sin(x)*exp(x/5)* exp(x*x)*sin(2*pi/dx*x). Приведенная программа, некорректно отобразит его.

А вот если взять TAB MathGrapher 1.0 (распространена в Интернете) и просто ввести 5* Sin(200*pi*x), то мы получим чистый ноль. Вместо 5, понятно, можно написать любую функцию, да и вместо Sin(200*pi*x) любую периодическую с кратной частотой, и программа выдаст неверный график.

Как с этими недостатками бороться? Очевидно, нужно, чтобы программа давала возможность рассматривать функцию на небольшом отрезке, с dx=(х2-х1)/w (ширина вывода в пикселях). Т.е. программа должна автоматически менять шаг с уменьшением отрезка(значения х2, х1 нужно вводить вручную).