-->x = 12 * (6 : 9);
-->y = 1 : 4;
-->z = x ./ y
z=
72. 42. 32. 27.
Ответ к упражнению 4.5 (Векторизованное возведение в степень ) Получим вектор
x21; x22; x23; x24 при x = (1; 2; 3; 4), выполнив поэлементное возведение в степень с использованием оператора ” .ˆ”:
-->x = 1 : 4; -->y = x .^ 2
y=
1. 4. 9. 16.
Ответ к упражнению 4.6 (Применение функции к вектору ) Получим вектор (sin(x1); sin(x2);
: : : ; sin(x10)) при x 2 [0; ], применив функцию sin к вектору, полученному в результате вызова функции linspace:
-->x |
= |
linspace (0 , |
%pi , 10); |
|
|
|
|
||
-->y = sin (x) |
|
|
|
|
|
|
|
||
y |
= |
|
|
|
|
|
|
|
|
|
|
column |
1 |
to |
6 |
|
|
|
|
|
0. |
0.3420201 |
0.6427876 |
0.8660254 |
0.9848078 |
||||
0.9848078 |
|
|
|
|
|
|
|
||
|
|
column |
7 |
to |
10 |
|
|
|
|
|
0.8660254 |
0.6427876 |
0.3420201 |
1.225 D -16 |
|||||
Ответ к упражнению 4.7 (Векторизованные функции ) |
Вычислим значения y = f(x) |
||||||||
функции f, заданной уравнением |
|
|
|
|
|
||||
|
|
|
f(x) = log10 (r=10x + 10x) |
|
(3) |
||||
при r = 2:220:10 16 и x 2 [ 16; 0]. Следующий пример демонстрирует, как это можно сделать, используя операторы поэлементного деления ” ./” и возведения в степень ” .ˆ”:
-->r = 2.220 D -16;
-->x = linspace ( -16 , 0, 100);
-->y = log10 (r ./ 10 .^ x + 10 .^ x );
Данная функция используется при вычислении оптимального шага численного дифференцирования, который равен h=sqrt(%eps).
Список литературы
[1]Atlas - automatically tuned linear algebra software. http://math-atlas. sourceforge.net.
[2]Cecill and free software. http://www.cecill.info.
[3]Stephen L. Campbell, Jean-Philippe Chancelier, and Ramine Nikoukhah.
Modeling and Simulation in Scilab/Scicos . Springer, 2006.
86
[4] J.-P. Chancelier, F. Delebecque, C. Gomez, M. Goursat, R. Nikoukhah, and
´
S. Steer. Introduction a` Scilab, Deuxi`eme Edition. Springer, 2007.
[5] The Scilab Consortium. Scilab. http://www.scilab.org.
[6] Intel. Intel math kernel library. http://software.intel.com/en-us/ intel-mkl/.
[7] Cleve Moler. Numerical computing with MATLAB.
[8] Flexdock project. Flexdock project home. https://flexdock.dev.java.net/.
87
Предметный указатель
библиотеки, 62 библиотеки функций, 62 булев тип, 24 целые числа, 26
элементарные функции, 23 функции, последовательность вызова,
58 функции, тело, 59
функции, заголовок, 59 имя переменной, 22 комментарии, 20, 23 комплексные числа, 25 комплексное сопряжение, 44 консоль, 12
левосторонние аргументы, 58 логический тип, 24 матрицы, 32 модуль, 62 оператор ”^”, 22
оператор ”:”, 37 оператор ”;”, 22 оператор ” ’ ”, 44 пакетная обработка, 19
поэлементные операции, 43 правосторонние аргументы, 58 приглашение, 12 продолжение строки, 23 расположение панелей, 15 размер матрицы, 32 строки, 30 точка, 23
транспонирование, 44 ^, 22
’, 44
.’, 44
.., 23 //, 23 ;, 22
SCIHOME, 64 contour, 71 disp, 12 exec, 61 feval, 74
function, 58, 59 genlib, 62 help, 8
lib, 62 linspace, 71 plot, 70, 71 size, 35 testmatrix, 38 title, 75 xtitle, 75
:, 37
SCIHOME, 18
%i, 25
%pi, 24 ans, 30
Intel, 7
Linux, 7
Mac OS, 8
MKL, 7
Windows, 7
88