Metody Obliczeniowe
Działania macierzowe w systemie MATLAB
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski
Operacje na macierzach
• Definiowanie
A=[2 2 2 1; 1 2 3 1];
A=[min:krok:max]
A=eye(3)
A =
1 0 0
0 1 0
0 0 1
A=ones(2,3)
A =
1 1 1
1 1 1
A=zeros(3,2)
A =
0 0
0 0
0 0
Met
Działania macierzowe w systemie MATLAB
ody Obliczeniowe
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski
• Dostęp do elementów macierzy
» A(i,j)
- odwołanie do elementu w wierszu i i kolumnie j;
• Działania na macierzach i działania tablicowe
Suma, różnica, iloczyn
A+B, A-B, A*B
- obowiązują reguły z matematyki
UWAGA: CZĘSTY BŁĄD
2
y
x
=
,
[0 1 2 3 4 5 6 7 8 9 10]
x
=
1 10
1 10
x
x
x
x
⋅
ŹLE !!!
x=[0:1:10];
y=x.^2=;
y=x.*x;
Metody Obliczeniowe
Działania macierzowe w systemie MATLAB
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski
Funkcje
Function
[pole,obwod] = prostokat(a,b)
% funkcja obliczajaca pole prostokata
%-------------------------------------------------
% WEJSCIE
% a - długość jednego z bokow prostokata
% t - długość drugiego z bokow prostokata
%-------------------------------------------------
% WYJSCIE:
% pole - pole prostokata
% obwod – obwod prostokata
pole = a*b;
obwod=2*(a+b);
Funkcję zapisujemy w pliku prostokat.m
Funkcję wywołujemy w sposób następujący:
[pole,obwod] = prostokat(2,3)
>> pole = 6
>> obwod = 10
Metody Obliczeniowe
Działania macierzowe w systemie MATLAB
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski
function
[Ke] = ke_beam(EJ,L)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% funkcja generuje lokalna macierz sztywnosci elementu belkowego
%-----------------------------------------------------------------
% WEJSCIE:
% EJ = sztywnosc gietna EJ
% L = dlugosc elementu
%-----------------------------------------------------------------
% WYJSCIE:
% Ke = macierz sztywnosci 4x4 wzgledem przemieszczen:
% v_a,fi_a,v_b,fi_b
%-----------------------------------------------------------------
% AUTOR:
% M. HIRSZ KATEDRA MECHANIKI BUDOWLI I MOSTOW, POLITECHNIKA
GDANSKA,PAZDZIERNIK 2006
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%******************************
% V_i Fi_i V_k Fi_k
%******************************
Ke=[ 12 6*L -12 6*L ;
6*L 4*L^2 -6*L 2*L^2;
-12 -6*L 12 -6*L ;
6*L 2*L^2 -6*L 4*L^2]*EJ/L^3;
Metody Obliczeniowe
Działania macierzowe w systemie MATLAB
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski
Format plików
•
.m - zawiera ciąg instrukcji (skrypt, kod ASCII)
•
.mat - zawiera wyniki (kod binarny)
Zasady nadawania nazw plików .m
•
Nazwa pliku nie może zaczynać się od cyfry
•
W nazwie pliku nie mogą znajdować się
•
spacje,
•
polskie znaki,
•
oraz znaki typu + - : , . ; / \ itd.
•
Nie wolno nadawać nazw takich, jakie mają wbudowane funkcje matlabowe, np. sin.m
Metody Obliczeniowe
Działania macierzowe w systemie MATLAB
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski
okno komend
skrypt
(command window)
(m-file)
program
uruchamiamy
klawiszem
tu oglądamy wyniki
tu piszemy program
F5
oba okna musza być widoczne w czasie pracy na ekranie !!!!
Metody Obliczeniowe
Działania macierzowe w systemie MATLAB
Politechnika Gdańska
Wydział Inżynierii Lądowej i Środowiska
Katedra Mechaniki Budowli i Mostów, studia dzienne, sem. 5
Magdalena Rucka, Wojciech Witkowski