- •Contents
- •1 LAboratory work # 1
- •Mathematical model
- •Stages of a program elaboration
- •Call Desktop matlab
- •Script-files and Function-files
- •Enter of input data by awarding method. Comments
- •Organization of enter of the input data by a dialogue mode
- •Creation of Function-file
- •Graphical output
- •2 LAboratory work # 2
- •Debugging and verification of programs
- •Search of syntactic mistakes
- •Debugging with the help of Editor/Debugger
- •Verification of results of calculation
- •3 LAboratory work # 3
- •The task for fulfillment
- •Individual tasks
- •4 LAboratory work # 4
- •Mathematical model
- •The block-diagram of algorithm of calculation according to mathematical model
- •The task for fulfillment
- •5 LAboratory work # 5
- •The task for fulfillment
- •Individual tasks
- •6 LAboratory work # 6
- •Mathematical model
- •Determination of zero approximation
- •Program of calculation in matlab environment
- •Results of calculation
- •Individual tasks
- •The task for fulfillment
- •7 LAboratory work # 7
- •Mathematical model
- •Program of calculation at matlab environment
- •Results of calculation
- •Individual tasks
- •The task for fulfillment
- •8 LAboratory work # 8
- •Mathematical model
- •Results of calculation
- •Improvement of convergence of the Newton method
- •The task for fulfillment
- •9 LAboratory work # 9
- •Mathematical model
- •The program of calculation in matlab environment
- •Results of calculation
- •The task for fulfillment
- •10 LAboratory work # 10
- •The task for fulfillment
- •Individual tasks
- •LIst of literature
Organization of enter of the input data by a dialogue mode
Enter the input data with the help of the operator input:
R1=input (' input value R1 = ');
R2=input (' input value R2 = ');
R3=input (' input value R3 = ');
L2=input (' input value L2 = ');
C3=input (' input value C3 = ');
ph1=input (' input value ph1 = ');
ph2=input (' input value ph2 = ');
For this purpose:
Click M - file Circuit 1 and replace input operators by above-mentioned ones;
Transfer output results of calculation of arguments phI1, phI2, phI3 of currents I1, I2, I3 from radians into degrees:
disp ([' I1M = ', num2str (I1M), ' phI1 = ', num2str (phI1/pi*180)]);
disp ([' I2M = ', num2str (I2M), ' phI2 = ', num2str (phI2/pi*180)]);
disp ([' I3M = ', num2str (I3M), ' phI3 = ', num2str (phI3/pi*180)]);
Replace appropriate output operators;
Save the program as Circuit 2 into the same directory where Circuit1 is saved;
From Command Window click Circuit 2 for performance (or click Run icon on panel of tools of window of Editor/Debugger)(see1.9);
Enter numerical values of the input data by means of dialogue mode. After each entered numerical value, press ENTER.
Results of calculation according program Circuit 2 completely should coincide with results of calculation according program Circuit 1, except values of arguments phI1, phI2, phI3 of complex numbers I1, I2, I3 for which output values are transformed from radians into degrees.
> Circuit 2
input value R1=1
input value R2=2
input value R3=3
input value L2=1e-3
input value C3=1e-4
input value ph1=30
input value ph2=60
I1=1.31391+30.2501i
I2 =-0.904861+28.5767i
I3=2.2188+1.6734i
error calculation eps=3.5527e-015+7.1054e-015i
I1M=30.2786 phI1=87.5129
I2M=28.591 phI2=91.8136
I3M=2.7791 phI3=37.0242
Creation of Function-file
For the further output of time dependences of branch currents to the graphic, let’s calculate arrays of these functions with the help of MATLAB subroutine - function.
Mathematical dependence of a branch current from time is expressed
where I - the module of a current complex (effective value of the branch current),
- argument of the current complex.
Let’s realize subroutine - function of calculation of i (t) as current (IM, f, phI, t), where values of formal parameters are:
IM - the module of a current complex;
f - frequency;
phI - argument of a current complex ;
t - current time.
Open a new file in the text editor and type the following program - function:
function im=current (IM, f, phI, t)
im=sqrt (2) *IM*sin (2*pi*f*t+phI);
Save the created subroutine into the same personal directory with name “current”.
Graphical output
Click M - file Circuit 2;
Write down at the end of the program the fragment realizing filling of arrays of instantaneous values of branch currents and their graphical view:
% GRAPH
h=0.04/200;
t (1)= 0;
i1 (1) =current (I1M, f, phI1, t (1));
i2 (1) =current (I2M, f, phI2, t (1));
i3 (1) =current (I3M, f, phI3, t (1));
for k=2:201
t (k) =t (k-1) +h;
i1 (k) =current (I1M, f, phI1, t (k));
i2 (k) =current (I2M, f, phI2, t (k));
i3 (k) =current (I3M, f, phI3, t (k));
end
plot (t, i1, t, i2, t, i3);
Save the program as Circuit 3;
Click Circuit 3 for performance (see .1.9).
The calculated graphical branch current dependences of the time are shown in Fig. 1-2.
Draw up diagrams (colors, line styles, line width, grid), as shown in Fig.1.2;
To change colors, line styles and line width click right button of mouse on the image of edited graphical curve. From the fallen out menu choose necessary option.
For appearance grid press Main menuEditAxes properties…Grid show (for x-axis)Grid show (for y-axis).
Figure 1.2 – Branch current dependences of the time
Write down values of the input data into the field of the graphics.
Start the program once again having changed the input data:
e1 (t) =200sin (t +/2) V; e2 (t) =80sin (t +/3) V;
R1=0.5 Ohm; R2=1 Ohm; R3=4 Ohm;
C3=2*10-4 F; L2=4*10-3 H.
Look, as it will be reflected in calculation results.