Численные методы - Octave Использование командной строки

Материал из Xgu.ru

Начнем с элементарного примера -- построения эпюры изгибающих моментов консольной балки, нагруженной распределенной нагрузкой (см. рис.).

В командной строке Octave последовательно набираем команды (для тех кто впервые видит консольное приложение -- в конце каждой строчки после символа ";" следует нажимать Enter):

> L = 100;
> q = 2000;
> x = 0:0.1:L;
> y = -(q/2)*(x.^2);
> plot(x,y)

В результате должно получится нечто похожее на график.

И так, в первой строке переменной L (длина балки) присваивается значение 100, во второй -- переменной q (нагрузка), значение 2000. Третьей строкой создается вектор (массив) узлов от 0 до L с шагом 0.1, в которых будет рассчитываться момент.

Допускается точку с запятой ";" в конце каждой строки не ставить, но в таком случае будет выводится результат вычислений. Для сравнения:

Вариант 1				Вариант 2
> q = 2000				> q = 2000;
q = 2000				>
>

или

Вариант 1				Вариант 2
> x = 0:10				> x = 0:10;
x = 					>
	0 1 2 3 4 5 6 7 8 9 10
>

или

Вариант 1				Вариант 2
> 2+2					> 2+2;
ans = 4				>
>

Значения изгибающих моментов y (4-я строка) рассчитываются по известной формуле: (поскольку приложенная к каждой рассматриваемой отсеченной части нагрузка приводит к сжатию нижних волокон балки, то согласно правилу знаков, момент является отрицательным). Выражение x.^2 приводит к возведению в квадрат каждого i-го элемента массива x.

И наконец, последней строкой будет выведен на экран график распределения моментов по длине балки.

KorenkovVladimir 10:36, 19 октября 2007 (EEST)