· 27 3 0 sin" exp 3 R+$ C eCE; S log10 R :UTS logD+$ sqrt round-(L H18 B GK min QA:OMNPC 4min(A)...

28

Transcript of  · 27 3 0 sin" exp 3 R+$ C eCE; S log10 R :UTS logD+$ sqrt round-(L H18 B GK min QA:OMNPC 4min(A)...

���� �������

MATLAB84/5ZMATLAB6:-

“Matrix Laboratory” 9�^� �#% 92;9,�!NXBTPYB$'Z��

• )�0��J[F�

• OAKV\_��)�]-"�\`��)�]

• �&�CDIQ

• (�7+�T@MTU6BTL?GAH[V

• �9�!>�32LU[ELK9Scilab-Octave0.<

Z��1=<�7�*��% -CPRW[CSY-���-BTL?GAD�

1

��•������• �&,"��,"�• �����)$!

•���%�!*�• �(#� �� plot• $+�(',�• for�• if�• while�• switch�

• �����-������.

2https://jp.mathworks.com/help/matlab/getting-started-with-matlab.html↓こちらのテキストも参考に

�� ���

3

u $+0'! 0'!�)/#-,������.&10�

>> a=1, b=2, c=a+b �a =

1. b =

2. c =

3.

>> a=1+2+3 ...�+4+5 �

a =15

u ���2��(."'���������

2��(."'���%*1%���

RYIVTZI

4

[KHWRN8?ARYIVTZIMATLAB9��B�<4&)1*5RYIVUB�A��!��.A6(: JSZOFDZOF8 QCEX�B�.A!

[��B��.A+68?ARYIVTZI�-#��B��.A+68?34 �=�=��6�8��.A+6'5(A!

u M-file8RYIVUB� -4�� \JZPEX>WZH:��7#]

�]���6QCEX�:��,/A��'"A!06%; ���'func17@ QCEX�:func1.m9?$8.A!

>> edit filename �6JSZOFDZOF8�.A6 GMDL'�2�'A

��������

5

a=1;b=2;c = a + b

>> edit test1.m� �a�b����� ������

function c = test1(a,b)c = a + bend

�!�� �a�b����� ���

�����������>> test1c=

3

�����������>> c=test1(1,2)c=

3

>> edit test1.m

ORGR �5����

1. EJBI?�'" >> edit test22. �5��?��1�'S�MRH1��T"

3. ��+>SLACP�6!test2.m5<$4+>T"4. FNQKDBQKD1 �?�/09>">> a = 5; b = 2; c = 1; [r,s] = test2(a,b,c)

��TORGR �5�%=!�9�9 �;�5ORGR �? 7�+(2:1&8+"

6

function [x,y] = test2(a,b,c)x = a + b;y = a – c;

�TFNQKDBQKD1�/.��?-588 �5�1�$(261&8,@"S�*�5��4)0:!#5:543>T

AC?C��4��

7

function [x,y] = test2(a,b,c)x = a + b;y = a – c;

��: test2D@=>B1��*-;E

function5��4���3��2��

�E��5"6',%#1�87,!

F����3��<�(!

�,;��D���1�%E5 '2/)�:-;!

�,;��D���1�%E9 '2/)�:-;!./+ ��4�0�&�7;��&$;!

8

u -14/$#

n % -14/,��6%�����(!-14/%$*7�

n 8 ��&���

&�+)%.0-34�8�, �*

�7 >> a=1, b=3;�

n ��"��

MATLAB!'�C��"� %��"��,���*

��252

��+��� �'&*)(%���#

u a�b%����"�a + b, a – b, a � b, a / b ��%������# test3%�����!������������$#'&*)( �����!��

9

10

>> x=[1; 2; 3]x =

1. 2. 3.

�� ��#! ��(,3(4"�!%1-5�� ��#! 6 �).'01�

���linspace$��������� ������$���+&*/$�>> linspace(1,9,5)ans =1 3 5 7 9

>> A=[1 2 3; 4 5 6; 7 8 10]A =

1. 2. 3. 4. 5. 6. 7. 8. 10.

+&*/2�� ���

11

� A�i�j���� A(i,j)���A =

1 2 3 4 5 67 8 10

>> A(2,2)ans =

5� A�m,:)�����A(:,n)��������>> A(2,:)ans =

4 5 6>> A(:,1)ans =

14 7

�������

����7�� (��8

��!� : A% B /�#>> A + B

��� : A�, B /�>> A – B

�� � : A.*B) A% B(����(�/�+.>> A .* B

��-� : ��� A(i,j)/B(i,j)(��/�+.>> A ./B

�A(*��: A.^B)��� A(i,j)( B(i,j)�$�.��/�+.>> A.^2>> A.^B

��(�'5614/� .9203��"$&�%��$�&�($��!

12

����

13

����� : �� A& B(����(��A(��) B(�� �#%'"+*��%!'�(%� >> A *B

���(���- xA = B (���� A& B(��)�#%'"+*��%!'�(%� >> A /B

���(��- Ax = B (�%$��� A& B(��)�#%'"+*��%!'�(%� >> A¥B

��� (�����(���������,>> A’

14

A =1. 2. 3. 4.

����>> R=inv(A)R =- 2. 1. 1.5 - 0.5

� m�n������ �1���>> d=ones(2,3)d =

1. 1. 1.1. 1. 1.

� m�n������ �%���>> d=zeros(2,1)d =

0.0.

� m�n������>> E=eye(2)E =

1. 0. 0. 1.

������"��

��!���$����"�������� �#�

�:�0!(%����;

15

$483652)1#." ,-1��

�� �� ���

for ����/�0!( for 794��<794/���

end

while ����/�0!( while�����

end

if, elseif, else ���� if ����

end

switch, case ���. *'%case�+��2�0�&

switch ��case ���

otherwise�

end

16

�* $ ��

�'���

U @A�1.

S @A�4.

UT ��

ST ��

TT �5.

~T �52:.P��<$=HLGQ

%����

&& /8

|| >7= Pshift + ¥FODR �Q

+���D��6B);�.B-EJKMIP���&Q0"�<��,��=# 39;(�4C,!�?"�9:B-

�'���N%����

‘for�

17

for�1�5�)��8����5�(0��7(2)"

���.� 9<: X7Y%4Z2-1*,�'+.&4��7�5�)

for X=Y:Z9��:

end

���.� 9=: X7Y%4Z2-W*,�'+.&4��7�5�)for X=Y:W:Z9��:

end

9�: i71%4102-1*,�'+.&4i0��7�5�)"for i=1:10disp(i); % disp(x)1�x7��/��)6��

end

9�:• for i=1:2:10 91%4102-2�3-�'+!107�$.#�1�5�):

• for i=10:-1:1 910%412--1�3-��'+;1/.62-�5�):

for•������������������

• !��#�����

18

for i=1:ndisp(i);

end

for i=n:-1:1disp(i);

end

for2.m

function s = sum1(a)n = size(a,1);s = 0;for i=1:n

s = s + a(i);end

sum1.m>> x = [1 4 5 2 7];>> t = sum1(x)

� � >> n = 10;>> for1

>> n = 10;>> for2

for1.m$��" ���% $��" ���%

$����%

‘while

19

Eps = 1; while 1+Eps > 1

Eps = Eps/2;endEps = Eps*2

�while ,�����3��&�����+�1�&3� -'�9�:eps1.m

Eps"1/1�# �,�Eps/23��&�%2

91.0 !0($*�# �����-)+��3�!.254876:

���>> eps1Eps =

2.2204e-16

‘if�

20

$���6%��3���5 9�?/3-1"�,>8.&A���6p.9��BA�B:-�� x *10��5�%k5�@C�,/=

if x >= 10k = k + 1;

end

$if�@,<4�)+�-0(��%else ; elseif@��.=

:-��C*�2'>7% �C@�!%��C*�2'>7%��D@#1%�2'>7% �D@�!%,<4��D:�2'>7% �E@�!

if ��C �C

elseif ��D �D

else �E

end

if�

if1.m ���������

if2.m ������

21

if a > 4disp(b);

end

function c = if2(a,b)if a > b

c = a + b;else

c = a - b;end

>> a = 3; b = 5;>> c = if2(a,b)

���

>> a = 3; b = 5;>> if1

���

⇓� �������������

⇓� ��������

‘switch�

22

function evenodd(n)

switch mod(n,2)case 0

disp(’even number');case 1

disp(‘odd number');otherwise

disp(’error');end

��'$� case �% �/�,��*

3�4evenodd.m ��&��/��+-��

��(��mod(x,y) )�x mod y/��*

2%�#"���5(��even3��4���6(��odd3��4

!.��)012& -

23

>> x = linspace(0, 2*pi);>> y = sin(x); >> plot(x,y,'-o')

�plot&�=��+387/���387,Figure<'� 10961*��$.-��plot+�*��"$.-�

�Figure</�Save as�&�:��&�;�jpg� )(&�%-#'!���

��+��'�+��,�axis��/� �>> axis([0 2*pi -1.5 1.5])

=��3870524

3��������

24

>> t=linspace(-pi,pi,10) %-π� π��10��>> [X,Y]=meshgrid(t) %����t� ��XY���>>Z=sin(X)+cos(Y);>> surf(Z) >> mesh(Z) >> contourf(Z)

25

4�� 1

zeros .*:?G@OSE-�M��

ones .*:?G@TE-�M��

eye ��-�M��

diag :QNPRE��F�0-�M�� :-�E��F6�0.*M�

magic 1 9J n^2 HAE��M�'>@6-��B���E�:)><CK n - n �E-�M��

rand A = rand(n) F62�3 (0,1) �E"$�!��9J�;=LK(����M�I n - n �E-�M��

randn A = randn(n) F6"$#/��9J�;=LK(����M�I n - n �E-�M��

linspace y = linspace(a,b) F6a B b. E3A6+�D)35C 100 %E-QNPR y M��y = linspace(a,b,n) F6a B b E3A6+�D)35C n %E-QNPR y M��

7�,8-�&�4�

26

%�� $�

¥(mldivide)

X = A ¥ B 4&����� AX = B @ X 72*3"-' X = inv(A) * B 5;<�/'R¥JQA9&Windows��49HEACMEBKT4��.?>S

/(mrdivide) X = A / B4&����� XB = A @ X 72*3"-' X = A * inv(B) 5;<�/'

det ���@#�

rank ��8MPA8#�

lu [L,U,P] = lu(A) 9&P*A = L*U 56>��!�� L&��!�� U&���� P @#�

qr [Q,R] = qr(A) 9&A ,m � n �50>5& A = Q*R 56>=+6&m � n �8��!�� R 5m �m �8����R �8�9LGDN��S Q @#�

eig d = eig(A) 9&�� A 80:38 ��@#�'[X,D] = eig(A) 9& �IAFO@�:1�� X 5 ��@�:1�!�� D @#�'

(�)��%�

27

3� 0�

sin "�

exp ��3�R+$��C�eCE;�Slog10 ��R�:UTS�logD+$��sqrt ��

round -(L�H18��B�GK

min���QA:OMNPC��4min(A) D A C��-(L2=QA :,�C��4min(A) D4A C�LOMNP?<>J�84��C��-(L�F,OMNPL2=

max���QA :OMNPC��4max(A) D A C��-(L2=QA :,�C��4max(A) D4A C�LOMNP?<>J�84��C��-(L�F,OMNPL2=5

mean���QA :OMNPC��4mean(A) D A C���L2=QA :,�C��4mean(A) D4A C�LOMNP?<>J�84���9IAK,OMNPL2=

std!#��Qs = std(X) D X :OMNPC?;4�/C (1) L�%<>!#��L.'QX:,�C��4std(X) D X C��C-(C!#��L�F,OMNPL2=

sum)�QA :OMNPC��4sum(A) D-(C�L.'QA :,�C��4sum(A) D A C�LOMNP?<>�84��C�L,OMNP?<>2=

6*7�&3�A@