mat lab slides

8/13/2019 Mat Lab Slides

1/63

C LF P. A


2/63

2

? G


3/63

3

? G


4/63

4

( ( ( ( ))))

2

( )

Matlabinterpretor

Signalprocessing

M-files

Control System

system

identification

Toolboxes Simulinkgraphic

interface

DSP Power system

Blocksets

Graphicwindows(outputs)

Commandwindow


5/63

5

? G


6/63

6

( )

C


7/63

7

( ) ( ) ( ) ( )

Command window

Start menus: Simulink

Command history

Variables being used


8/63

8

( ). ( ). ( ). ( ).


9/63

9

? G


10/63

10

( ) ( ) ( ) ( )

: : 1. 2. * /3. +

E>>8+3*5

= 23>>(8+3)*5

= 55


11/63

11

( ) ( ) ( ) ( )

A :=

>> = 3 >> = +2 2 ,

C

>>6 = >> +2=20


12/63

12

( ) ( ) ( ) ( )

Complex numbers

; suppresses screen printing


13/63

13

1. A =

.A =15 =8 . 20%?

2. G = 5+9 =6 2 , + =1+7 , = 12+64 / = 1.2+1.1 . C M ( + )

.

3. :

4. C R :

63% = 10 ,

R=10 C=10 F? C RC.


14/63

14

? G

A F


15/63

15

/ ( )/ ( )/ ( )/ ( )

or name =[start_val:increment:end_val]

Define a vector name =[start_val:end_val]

Names are case sensitive v2 V2

Multiplication element wise: .

X

=

Not

X

=


16/63

16

C = 1 2 4 5 7 3 (1,2,4,5,7,3) = (1, ) (0,0,0 .,0)

= (1, ) (1,1,1 .,1) =( : : ) ( , + , +2 , , + )

+ < +( +1) ( : )

= 1 2 1 2 = 1+ 2 1 2 = 1.* 2 1 2


17/63

17

/ / / /

E : >> = 1 0 2 5

+ ( )

Use function root to ge the roots of apolynomialsExample

>>roots(p)ans=

-2.09461.0473+1.1359i

1.0473-1.1359i

+


18/63

18

; C

(:,N) N S ( ) (

) D ( )

I ( )


19/63

19

( ) ( ) ( ) ( )

, = ( ), .

E :2 0.1 2

. 2+2

:/ / ,

B : ! I 2

2 2, 1.99, 1.98 1.99, 2.


20/63

20

? G


21/63

21

( ) ( ) ( ) ( )

, /

( , )

3 ,

( 10)

( 10)

Sampled! Interpolated!


22/63

22

( ) ( ) ( ) ( ) A ( )

. M

NOT the same

X-axis vector

same

Y-axis vector


23/63

23

( ) ( ) ( ) ( )

: !!!

10 samples/period 40 samples/period


24/63

24

( ) ( ) ( ) ( )

: >> ( , , , , , ) C :>> ( ) A : >> ( ) A : >> ( )

A : >> ( , ) C :>> ( ) : >> ( , , ) A : >>


25/63

25

( ) ( ) ( ) ( )

.

,


26/63

26

? G


27/63

27

I .

( ) ( DOS (. ) NI ,

) , ( ,

!) NI ,


28/63

28

, : >>

Figure options

2 curves on same plot


29/63

29

A , : C

A , !

function y = myfunc(param1, param2, )

myfunc returns y and accepts param1 andparam2 as input parameters

name the file after the functions name


30/63

30

F L F =1: ,

< >

I I ( )

< 1>

< 2>

A


31/63

31

A

N , .

I

( !!!!) , . .

M .

D .

2222


32/63

32

2 2 2 21.

2 :A = 1 2 3 4 5 6 7 8 9 10

2. , +1

.3.

( 1 13).

4. B L P F (B : B

). :

G ( B) : .

()

2

10)(10)( f H Log f G =

3333


33/63

33

3 3 3 3 2 . 2 . 2 . 2 .

Figure 1 Figure 2

3333 ( )( )( )( )


34/63

34

3333 ( )( )( )( ) A

) : = 100 H A = 0.3 = 30 . !

) A :

) R 0.1 0.1 ( )

) D .

) C ( )

,

4 ( )4 ( )4 ( )4 ( )


35/63

35

4 ( ) 4 ( ) 4 ( ) 4 ( )

.

. 2

.


36/63

36

? G


37/63

37

/ E O O

E


38/63

38

C .

.

D :

( )

>> .

=

x0x1x2...

xn-1xn

y0y1y2...

yn-1yn


39/63

39

C . . D

>> .

, . ,

>> . 1 2


40/63

40

S , . . H : =

( F ) .

F :

F ( )

D

time freq

T

1/T

FT


41/63

41

F ,

D F

A 1/ H : N = 1/(2)

( = N/2). A > :

s(k) : k th time sampleS(n) : nth frequency sample

Ts : sampling periodf d : DFT frequency step 1/(N.T s )

( )( )( ) ( )


42/63

42

( ) ( ) ( )( )

A signal is defined as followed (4 time samples given)if k =0 ,1

otherwise

Compute the 4 frequency samples using:

0)(

1)(

=

=

k s

k s

0 4321

1

k

( ) ( ) ( ) ( )


43/63

43

( )( )( )( )


44/63

44

:

( : !)

. ,

. =2

!

where


45/63

45Frequency not @ 50Hz -> plotted versus theindex of the samples, not versus frequency

2 peaks!!! And symetrical withrespect to vectors length/2

Magnitude of fft


46/63

46

,

, 1 2 2

, . . 2 . 0 ( /2) /

>> = 0:N/2 1 * /NI FF (

2 1.5 6 0.2 0.3 0.6 21.560.20.30.6

1 NN/2+1

2 1.5 6 0.2 0.3 0.6

4 3 12 0.4 0.6 1.2

+=

Fold

vectorafter FFT

sample index

vectorto plot


47/63

47

ZOOM

Frequency ok

Magnitude not ok

(x1024 too large) Divide by N to getCorrect magnitude

5555


48/63

48

6. 1, . ( ),

. .


49/63

49

C E E D ( D)


50/63

50


51/63

51

B : >> = (A)*

B G ( ): >> = (A, )( A )

B >> >>

>> >>

/ / / /


52/63

52

6 6 6 6


53/63

53


54/63

54

G , .

E D : L , ,

. C

C .


55/63

55

Click on Start

Simulink librairies


56/63

56

Block diagram of firstmodel


57/63

57

, . , , .

,

. ,

.

/ / / /


58/63

58

, . , ,

.

. , .

, .

, .

. , . , ,

. , ,

, . ,

.


59/63

59

, . A

. S . , E .

, .

, . ,

F =>F . I

OK. , .

,

. , . . I

, .


60/63

60

D .

/ . A :

N E E .

.

7777


61/63

61

. 1, 100 .

. .

1, 1 , 0.3 100 .

/ 16 (4 ). . 1 .

.


62/63

62

comparators

cos(2 f0t)

sin(2 f0t)

comparators

I

Q

x(t)

b3

b2

b0

b1


63/63

63

10 samples/periode 40 samples/periode

mat lab slides

Documents