Лабораторная работа №1

Дисциплина «Компьютерные технологии и программирование»

Выполнил:

студент гр. И – 2013 – МА

Чернобай И.В.

Проверил:

ст. преподаватель

Койфман А.А.

Мариуполь, 2014

Лабораторная работа № 2

ТЕМА: «ПРОСТЫЕ ВЫЧИСЛЕНИЯ В SCILAB»

ЦЕЛЬ РАБОТЫ:

По указанной литературе изучить:

• системное меню редактора SciLab,

• основные системные команды,

• правила ввода команд и данных,

• ранжированные переменные,

• правила вывода результатов,

• правила вывода результатов в виде двумерных графиков,

• правила отладки программ.

Разработать алгоритмы решения задач из варианта задания.

Составить программы решения задач.

Задача 1

• Ввести текст в виде комментария, как заглавие программы.

• Ввести исходные данные.

• Задать изменение аргумента.

• Вычислить значения функций 1 и 2 для аргумента в заданном интервале.

• Вывести графики функций одновременно на одном графике в декартовых координатах. Для разных графиков использовать разный тип линий.

Задача 2

• Пункты 1...4 задачи 1.

• Вывести графики функций в двух подокнах на одном графике. Графики сделать в столбиковом формате.

Задача 3

Провести обработку данных с целью определения коэффициентов линейной модели. Для этого:

• Ввести прямоугольную матрицу X размерностью 10х3 и Y размерностью 10х1. Дополнить матрицу Х столбцом, заполненным единицами.

• Выполнить вычисление коэффцициентов линейной модели, используя матричную формулировку системы нормальных уравнений:

Вектор b (параметры модели) найти при помощи встроенного алгоритма решения СЛАУ

b = (X’*X)\(X’*Y)

Вычислить расчетные значения отклика модели Yr = b*X, построить график зависимости Y(представить точками) и Yr (представить линией) от X (в качестве независимой переменной выбрать один из столбцов) и друг от друга.

Задание

Функция 1 y = |sin(x)|^exp(x/2)

Функция 2 z = 5x - x^1.5+sin(x)

Начальное значение аргумента a = 0

Конечное значение аргумента b = 5

Шаг изменения аргумента h = 0,5

Код Scilab

Задача 1

// Задача 1

// Диапазон и шаг

a=0;

b=5;

h=0.5;

// Задание аргумента

X=a:h:b;

// Расчет функций

Y= abs(sin(X)).*exp(X /2);

Z= 5.*X -X ^1.5+sin(X);

// Вывод графиков с одинаковым типом линии в окно 1

scf(1);

plot(X,Y,X,Z);

//Включим координатную сетку

xgrid(1)

//Вывод графиков с разными типами линии в окно 2

scf(2);

plot(X,Y,'-',X,Z,':');

//Включим координатную сетку

xgrid(2)

Экспорт графика из Scilab

Задача 2

// Задача 2

// Диапазон и шаг

a=0;

b=5;

h=0.5;

// Задание аргумента

X=a:h:b;

//Расчет функций

Y= abs(sin(X)).*exp(X/2);

Z= 5.*X-X^1.5+sin(X);

// Вывод графика 1 в виде столбиков в подокно 1

subplot(2,1,1),plot2d3(X,Y);

xgrid(1)

// Вывод графика 2 в виде столбиков в подокно 2

subplot(2,1,2),plot2d3(X,Z);

xgrid(1)

Экспорт графика из Scilab

Задача 3

clc;

clear;

//X=[1:3;1:3;1:3;1:3;1:3;1:3;1:3;1:3;1:3;1:3]

X=rand(10,3)

X1=[1 1 1 1 1 1 1 1 1 1]

Y=rand(10,1)

X=[X X1']

b=(X'*X)\(X'*Y)

Yr=X*b

Y

Министерство образования и науки Украины

Государственное высшее учебное заведение

Приазовский государственный технический университет

Кафедра автоматизации технологических процессов и производств