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