Scribd
Upload
123 views

Examples

The document provides an introduction to MATLAB programming for engineers. It covers MATLAB basics, vectors and matrices, loops, plots, and examples of using MATLAB for electrical circuits problems. Examples are provided to demonstrate solving circuits problems and plotting impedance versus frequency for RLC circuits using MATLAB code.

Uploaded by

venki249
Copyright
© © All Rights Reserved
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
123 views

Examples

The document provides an introduction to MATLAB programming for engineers. It covers MATLAB basics, vectors and matrices, loops, plots, and examples of using MATLAB for electrical circuits problems. Examples are provided to demonstrate solving circuits problems and plotting impedance versus frequency for RLC circuits using MATLAB code.

Uploaded by

venki249
Copyright
© © All Rights Reserved
Available Formats
Download as PPTX, PDF, TXT or read online on Scribd
You are on page 1/ 32

Matlab Programming for Engineers

Introduction to Matlab
Matlab Basics
Vectors and Matrices
Loops
Plots
MATLAB examples

Program :
>> r=10;
>> i=0:2:10;
>> v=i.*r;
>> p=(i.^2)*r;
>> sol=[p v i]
Output :
sol =
Columns 1 through 6
0

40

160

360

640

60

80

1000

Columns 7 through 12
0

20

40

100

Columns 13 through 18
0

7/31/15 09:48:50 AM

10

MATLAB PRESENTATION

Use MATLAB in Electrical


Circuits :

7/31/15 09:48:50 AM

MATLAB PRESENTATION

7/31/15 09:48:50 AM

MATLAB PRESENTATION

7/31/15 09:48:50 AM

MATLAB PRESENTATION

Solution usng MATLAB :


Z = [40 -10 -30;
-10 30 -5;
-30 -5 65];
V = [10 0 0]';
% solve for the loop currents
I = inv(Z)*V;
% current through RB is calculated
IRB = I(3) - I(2);
fprintf('the current through R is %f Amps \n',IRB)
% the power supplied by source is calculated
PS = I(1)*10;
fprintf('the power supplied by 10V source is %f watts
\n',PS)
7/31/15 09:48:50 AM

MATLAB PRESENTATION

MATLAB answers are :


the current through R is 0.037037 Amps
the power supplied by 10V source is
4.753086 watts

7/31/15 09:48:50 AM

MATLAB PRESENTATION

Vectors and Matrices

Arithmetic operations Matrices


Example:
-j5

290o

10

j10

1.50
o

Solve for V1 and V2


7/31/15 09:48:50 AM

MATLAB PRESENTATION

Vectors and Matrices

Arithmetic operations Matrices


Example (cont)

(0.1 + j0.2)V1 j0.2V2 = -j2


- j0.2V1

+ j0.1V2 = 1.5

0.1 j0.2 j0.2


j0.2

j0.1

V1
V
2

j2
=

1
.
5

Vectors and Matrices

Arithmetic operations Matrices


Example (cont)
>>>A=[(0.1+0.2j)0.2j;0.2j0.1j]
A=
0.1000+0.2000i00.2000i
00.2000i0+0.1000i
>>>y=[2j;1.5]
y=
02.0000i
1.5000
*
>>>x=A\y
A\B is the matrix division of A in
x=
to B, which is roughly the same as
14.0000+8.0000i
INV(A)*B
28.0000+1.0000i
*
>>>

Vectors and Matrices

Arithmetic operations Matrices


Example (cont)
>>>V1=abs(x(1,:))
V1=
16.1245
>>>V1ang=angle(x(1,:))
V1ang=
0.5191
V1 = 16.1229.7o V

7/31/15 09:48:51 AM

MATLAB PRESENTATION

12

HOW TO IMPORT LOAD DATA

SNO

se

re

1
2
3
4
5

1
2
3
4
5

2
3
4
5
6

7/31/15 09:48:51 AM

Resistance

0.0005
0.0005
0.0015
0.0251
0.3660

MATLAB PRESENTATION

Reactance

0.0012
0.0012
0.0036
0.0294
0.1864

13

ldata=[
1 1 2 0.0005 0.0012
2 2 3 0.0005 0.0012
3 3 4 0.0015 0.0036
4 4 5 0.0251 0.0294
5 5 6 0.3660 0.1864
];

7/31/15 09:48:51 AM

MATLAB PRESENTATION

14

clear all
ldata=[
1 1 2 0.0005 0.0012
2 2 3 0.0005 0.0012
3 3 4 0.0015 0.0036
4 4 5 0.0251 0.0294
5 5 6 0.3660 0.1864
];
nb=6;
nl=nb-1;
se=ldata(:,2);
re=ldata(:,3);
r=ldata(:,4);
x=ldata(:,5)
for i=1:nl
z(i)=r(i)+j*x(i);
end
for i=1:5
disp(z(i))
end
7/31/15 09:48:51 AM

MATLAB PRESENTATION

15

If R = 10 Ohms and the current is increased from 0 to 10


A with increments of 2A, write a MATLAB program to
generate a table of current, voltage and power
dissipation.

7/31/15 09:48:51 AM

MATLAB PRESENTATION

16

M-files : script and function files (script)

Example RLC
circuit

R = 10

eg4.m
eg5_exercise1.m

+
V

Exercise 1:
Write an mfile to plot Z, Xc and XLversus
frequency for R =10, C = 100 uF, L = 0.01 H.
7/31/15 09:48:51 AM

MATLAB PRESENTATION

17

M-files : script and function files (script)

Example RLC
circuit
Total impedance is given
by:

Z R j L

When

X C XL

Z R

1
o
LC
7/31/15 09:48:51 AM

MATLAB PRESENTATION

18

M-files : script and function files (script)

eg4.m
eg5_exercise1.m

Example RLC
circuit
120

Z
Xc
Xl

100

80

60

40

20

7/31/15 09:48:51 AM

200

400

600

800

1000

1200

1400

MATLAB PRESENTATION

1600

1800

2000
19

eg6.m

M-files : script and function files (script)

Example RLC
circuit
R = 10

+
V

For a given values of C and L, plot the following versus the frequency
a)

the total impedance ,

b)

Xc and XL

c)

phase angle of the total impedance


for 100 < < 2000

7/31/15 09:48:51 AM

MATLAB PRESENTATION

20

M-files : script and function files (script)

Example RLC
circuit100

eg6.m

Magnitude
Mag imp
Xc
Xl

80
60
40
20
0

200

400

600

800

1000

1200

1400

1600

1800

2000

1200

1400

1600

1800

2000

Phase

100
50
0
-50
-100

7/31/15 09:48:51 AM

200

400

600

800

1000

MATLAB PRESENTATION

21

Simulink

Model simplified representation of a system e.g. using


mathematical equation
We simulate a model to study the behavior of a system
need to verify that our model is correct expect results
Knowing how to use Simulink or MATLAB does not
mean that you know how to model a system

7/31/15 09:48:51 AM

MATLAB PRESENTATION

22

Simulink

Problem: We need to simulate the resonant circuit


and display the current waveform as we change the
frequency dynamically.
10

i
Varies from 0
to 2000 rad/s

100 uF

+
v(t) = 5 sin t

0.01 H

Observe the current. What do we expect ?

The amplitude of the current waveform will become


maximum at resonant frequency, i.e. at = 1000 rad/s
7/31/15 09:48:52 AM

MATLAB PRESENTATION

23

Simulink

How to model our resonant circuit ?


i

10

100 uF

+
v(t) = 5 sin t

0.01 H

Writing KVL around the loop,

7/31/15 09:48:52 AM

di 1
v iR L
idt
dt C

MATLAB PRESENTATION

24

Simulink

Differentiate wrt time and re-arrange:


2

1 dv di R d i
i

2
L dt dt L dt LC
Taking Laplace transform:

sV R
I
2
sI s I
L L
LC
sV
1
2 R
I s s
L
L
LC

7/31/15 09:48:52 AM

MATLAB PRESENTATION

25

Simulink

Thus the current can be obtained from the voltage:

I V

7/31/15 09:48:52 AM

s(1/ L)

R
1
2
s s
L
LC

s(1/ L )
R
1
2
s s
L
LC
MATLAB PRESENTATION

26

Simulink

Start Simulink by typing simulink at Matlab prompt


Simulink library and untitled windows appear

It is where we obtain
the blocks to
construct our model

7/31/15 09:48:52 AM

It is here where we
construct our model.

MATLAB PRESENTATION

27

Simulink

Constructing the model using Simulink:


Drag and drop block from the Simulink library
window to the untitled window
1
Sine Wave

7/31/15 09:48:52 AM

s+1
Transfer Fcn

MATLAB PRESENTATION

simout
To Workspace

28

Simulink

Constructing the model using Simulink:

s(1/ L)
R
1
2
s s
L
LC

s(100)
2
6
s 1000s 1 10
100s
s2+1000s+1e6

Sine Wave

Transfer Fcn

i
To Workspace

v
To Workspace1

7/31/15 09:48:52 AM

MATLAB PRESENTATION

29

Simulink

eg8_sim.mdl

We need to vary the frequency and observe the current


5
Amplitude

Ramp

v
To Workspace3

w
1
1000
Constant

To Workspace2

s
Integrator
Dot Product3

100s
s2+1000s+1e6

sin
Elementary
Math

Dot Product2

Transfer Fcn1

i
T o Workspace

From initial problem definition, the input is 5sin(t).


You should be able to decipher why the input works, but you
do not need to create your own input subsystems of this
form.
7/31/15 09:48:52 AM

MATLAB PRESENTATION

30

Simulink

1
0.5
0
-0.5
-1

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

-5

7/31/15 09:48:52 AM

MATLAB PRESENTATION

31

Simulink

eg9_sim.mdl

The waveform can be displayed using scope similar to


the scope in the lab

5
Constant1

2000
Constant

0.802
Slider
Gain

7/31/15 09:48:52 AM

100s
sin

s
Dot Product2
Integrator Elementary
Math

MATLAB PRESENTATION

s2+1000s+1e6
Transfer Fcn

Scope

32

You might also like