- •Численные методы моделирования
- •Численные методы моделирования
- •Первичная обработка информации
- •1.Для 8 значений X из диапазона 5..8 заполнить и распечатать массив y[I]
- •Численное дифференцирование
- •Численное интегрирование
- •Численное решение дифференциальных уравнений
- •Основные команды мatlab
- •1.Требования к идентификатору
- •2.Команды ввода исходных данных
- •3. Арифметические операторы
- •4. Операторы отношения
- •5.Логические операторы
- •9..Графика
- •2.Пример 2.
- •3.Функции
Основные команды мatlab
1.Требования к идентификатору
Состоит из латинских букв и цифр
Прописные и строчные буквы различаются
2.Команды ввода исходных данных
а=2 %
M=[1 3 5 1]
M2=[23 3
2]
x=1:5
y=2:.1:3
a=(1:.1:2)%строка 11 элементов
b=(1:.1:5)’%столбец 41 элемент
3. Арифметические операторы
приоритеты
уровень 1
(^)% возведение в степень
(.^)% поэлементное возведение в степень
(‘)% транспонирование
уровень 2
(*) % умножение чисел (матриц согласованного размера)
(.*)% поэлементное умножение матриц одинакового размера
(/) % деление чисел
(./) % правое поэлементное деление матриц
(.\) % левое поэлементное деление матриц
Задание
|
Программа
|
Результат
|
1. |
a=2 b=4 a/b a\b c=a*b |
a = 2 b = 4 ans = 0.5000 ans = 2 c= 8 |
2. |
A=[1 2 -1 0] B=[2 1 -1 1] C=A*B D=A.*B M=A/B N=B\A K=A./B |
A = 1 2 -1 0 B = 2 1 -1 1 C = 0 3 -2 -1 D = 2 2 1 0 M = 1.0000 1.0000 -0.3333 0.3333 N = 0.6667 0.6667 - 0.3333 0.6667 K = 0.5000 2.0000 1.0000 0 |
уровень 3
(+) %сложение чисел (матриц согласованного размера)
(-)%вычитание чисел (матриц согласованного размера)
уровень 4
(:) %оператор формирования массивов
4. Операторы отношения
(<)
(<=)
(>)
(>=)
(= =)
(~=)
5.Логические операторы
(&)И
(|)ИЛИ
(~)НЕf
Оператор условия
If…else elseif end
if (усл)
P1
end
if (усл)
P1
else
P2
end
if (усл)
P1
else
P2
end
if (усл)
P1
elseif (усл1)
P2
else
p3
end
Оператор цикла с неопределенным числом операций
while..end
while <усл>
p
end
Оператор цикла с определенным числом операций
for..end
for<перем.цикла>=<н.знач>:<приращение>:<к.знач
p1
end
Операторы
abs(x)
acos(x)
sin(x)
asin(x)
cos(x)
acos(x)
tan(x)
atan(x)
cot(x)
acot(x)
log(x)
log2(x)
log10(x)
sqrt(x)
exp(a)поэлементное вычислениеexp
sqrt(a)% поэлементное вычисление корня квадратного
a^p% возведение в степен
a.^p% поэлементное возведение в степень
round(x)% округление round(10.6) = 11
floor(x)%округление до меньшего целогоfloor(10.6) = 10
mod%oстаток от деления с учетом знака mod(10,3) = 1
sum(a)% суммирование элементов вектора(матрица по столбцам)
cumsum(a)% куммулятивное суммирование элементов вектора
max(а)
min(а)
mean(а)
rand(3,4)% генерация матрицы 3*4
a=[]%стирает матрицу
zeros(3,4) нулевая матрица
diag(a) диагональная матрица
tril(a)-треугольная левая
triu(a)-треугольная правая
inv(a)-обращение матриц
a(1,:) %первая строка
b(:,3)% третий столбец
с(1:2,1:10)% блок матрицы
9..Графика
plot(x,F)
plot(x,F,x1,F1)
plot(x,F,x1,F1,’+c’) % с=[y,m,c,r,g,b,w,k] –цвет линии% .+=[ o + * - :v<>]-маркер
bar(x,F)
loglog
semilogx
smilogy
mesh
figure
subplot(3,3,1)
hold on
hold of
title(‘заголовок’)
xlabel
ylabel
zlabel
text(x,y,’text’)
2.Пример 2.
Задание
|
Программа
|
ответ
|
Построить график функции y=x*x (-1<x<1) |
x=(-1:.1:1) for i=1:21 y(i)=x(i)^2 end plot(x,y) xlabel('x') ylabel('y') title('y=x^2') |