Difference between revisions of "Scilab"
(→Advanced Level Scilab) |
(→Advanced Level Scilab) |
||
Line 102: | Line 102: | ||
==Advanced Level Scilab== | ==Advanced Level Scilab== | ||
# File Handling- Scilab File handling | # File Handling- Scilab File handling | ||
+ | #* Writing to a file using write() | ||
+ | #* Reading from a file using read() | ||
+ | #* Opening an existing file using mopen() | ||
+ | #* Closing an already opened file using mclose() | ||
# File Handling- User Defined Input and Output in Scilab | # File Handling- User Defined Input and Output in Scilab | ||
+ | #* Input Function | ||
+ | #* mprintf() | ||
+ | #* save() and load() | ||
+ | #* Used to quit scilab midway through calculation and continue at later stage | ||
# Numerical methods- Integration | # Numerical methods- Integration | ||
# Numerical methods- Solving Non- linear Equations | # Numerical methods- Solving Non- linear Equations | ||
Line 111: | Line 119: | ||
# Numerical methods- ODE- Applications | # Numerical methods- ODE- Applications | ||
# Optimization Using Karmarkar Functions | # Optimization Using Karmarkar Functions | ||
+ | #* About Optimization | ||
+ | #* Use of Scilab function Karmarkar in Optimization | ||
# Optimization of Non-linear Functions | # Optimization of Non-linear Functions | ||
# Digital Signal Processing | # Digital Signal Processing | ||
Line 120: | Line 130: | ||
# Discrete systems | # Discrete systems | ||
# Calling User Defined Functions in XCOS | # Calling User Defined Functions in XCOS | ||
− | + | #* Write a squaring function | |
− | #* | + | #* Use of scifunc block in XCOS |
− | #* | + | #* Use of MUX block |
− | #* | + | #* Call functions having multiple inputs and output |
− | #* | + | |
− | + | ||
==Xcos== | ==Xcos== |
Revision as of 11:45, 3 May 2013
Template:Spoken Tutorial Banner
Scilab (http://scilab.org, http://scilab.in, http://scilab.cn) is an open source scientific software package for numerical computations. It has an extremely reliable and efficiently coded numerical library. It is also a high productivity tool: Through its interpreted language, one can quickly develop the code required to solve problems. Typically, if it takes ten lines of C code for some calculation, Scilab would require only one for the same purpose.
Scilab runs on all popular operating systems. Since 1994, it has been distributed freely along with the source code via the Internet. It is useful for students at school, college and roupas da moda research level. It is also useful for researchers and professionals.
IIT Bombay is leading the effort to popularise Scilab in India. This is part of the Free and Open source Software for Science and Engineering Education (FOSSEE) project, supported by the National Mission on Education through ICT of MHRD (http://spoken-tutorial.org/NMEICT-Intro).
IIT Bombay is using Spoken Tutorials (http://spoken-tutorial.org) to create learning material for FOSS. This is the main page for the organisation of the scripts required for Scilab spoken tutorials. We invite the Scilab user community to participate in this activity.
Contents
Note: Each numbered topic corresponds to a single spoken tutorial. Each bulleted point corresponds to a command or topic that must be covered organic nuts in the given spoken tutorial.
General Information
This topic will include all functionality in Scilab that iphone photography is not domain specific- functionality that is required across several domains.
Why Scilab
Basic Level Introduction to Scilab
This level will include a set of tutorials that are required to be known in order to qualify as "Scilab Literate". The tutorials here will teach programming weight loss pills fundamentals and the most commonly used Scilab functionality- Matrix operations and Plotting.
- Installing
- Show where to download from and how to decide which version to choose (OS and 32/64bit) (www.scilab.org/download)
- Windows installation (Internet Connection is necessary)
- Linux installation (using package manager- show only Debian/Ubuntu as example (sudo apt-get install scilab) as well as generic binary)
- Mac
- Compilation from source can come as a part of a more advanced tutorial
- Getting Started: Basic Housekeeping
- Expressions: Show mathematical expressions with numbers (42+4^2-64/4).
- Variables: Assign scalars to variables- a=1, b=2, c=3.
- Diary:CD, diary open and close and show transcript of some calculations given below
- Define symbolic constants. Use %i as example. Verify showing %i*%i to be -1. Also show sqrt(-1). Define a complex number (5*%i) and show a complex number expression (10+5*%i)*(2*%i). Show %pi and %e.
- Basic functions Show sin(%pi/2), cos(%pi/2), exp(1) versus %e, exp(2) versus %e^2. Define a matrix.
- suppressing output(;)
- help,clc
- --------------------------------------------------------------------------------------------------------------------------------------------
- Matrix Operations
- Vector Operations
- Define p = [1 2 3], q = [2 3 4]
- size(p), length(p)
- p+q, q-p, p*q, p', p'*q, p*q'
- Define P = [1 2 3; 4 5 6]
- size(P), P'
- Matrix Operations
- Square matrices Q = [1 5; 2 4], S = [2 0; 4 1]
- det(Q), diag(Q), T=Q+S, Q-S, Q*S, E = 2*Q + 3*S
- Matrix generation: zeros(3,4), ones(2,5), eye(4,4), diag([1 2 3]), rand(2,3). Also show how one can generate these matrices (zeros, ones, eye, rand) of the same size as another matrix by passing that other matrix as an argument- zeros(A), rand(A).
- Ranges: 1:4, 2:2:8, linspace(1, 9, 5) (linspace is to be explained in Plotting 2D graphs tutorial)
- Referring: P(2), P(2,1), P(6), P(3,2), P(1,$), P(2,:), P(:,2), P([2 3],[1 3]), P(1,2)=8.5
- Elementary row operations:P(2,:) = P(2,:) - 4*P(1,:)
- Solving equations: P = [1 2 -1; -2 -6 4 ; -1 -3 3], q = [1; -2; 1], x=P\q
- --------------------------------------------------------------------------------------------------------------------------------------------
- Conditional Branching
- Explain booleans: %T and %F- they form the basis for the if statement.
- First explain 'if' and 'then' with the example: n = 42, if (n == 42) then disp("The number is forty two"), end
- Now explain the use of the 'else' keyword: n = 54, if (n == 42) then disp("The number is forty two"), else disp("The number is not forty two"), end
- Now explain the use of the 'elseif' keyword: n = 54, if (n == 42) then disp("The number is forty two"), elseif (n == 54) then, disp("The number is fifty four"), else disp("The number is neither forty two nor fifty four"), end
- Say that if there are several branches, it may be clearer to use the 'select' keyword.
- Give example for select: n = 100, select n, case 42, disp("The number is forty two"), case 54, disp("The number is fifty four"), else disp("The number is neither forty two nor fifty four"), end
- Branching can be based on satisfaction of combination of multiple conditions as well: a=1, b=2; if (a>0 & b>0), disp('product of a and b would be positive'), elseif (a==0 | b==0) disp('product of a and b would be zero'), end
- Iteration
- Explain syntax of 'for' statement- tell that the variable iterates over a list/vector/matrix (or an expression that evaluates to any of these).
- Give example: for i = 1:5, disp (i), end
- Then explain break condition. Use example: for i = 1:10, disp(i), if (i==5), break, end, end
- Then explain continue condition. Use example: for i = 1:10, if (i<=5) then continue, else disp(i), end, end
- Explain while condition.
- Give example: i = 0; while(i <=5), i = i + 1; disp(i), end
- Then explain break, continue condition with example in while loop.
- Plotting 2D graphs
- About linspace: linspace is a linearly spaced vector.
- Plot a simple graph: x=linspace(12,34,10), y=linspace(-.1,2,10), plot(x,y)
- plot2d(x,y,style=1), style can have different values, style=1 gives black coloured graph line, style=2 gives blue coloured graph line and so on
- Show that if one tries to plot another graph without closing the first one, the plot will appear over the first graph. This may or may not be desirable- one may want to plot the second graph in a fresh window. Then explain how scf() can be used to alternate between graph windows.
- Using clf() clear the graphic window.
- Configure the title for the plot using title("My title"), Configure the axis separately using xlabel("X"), ylabel("Y")
- Configure a legend to know what curve is associated with y1 = x^2 and what curve is associated with y2 = 2x^2.
- Divide a graphic window into a matrix of sub-windows using subplot(mnp)
- Scripts and Functions
- Explain that one often repeats a set of commands- in which case it is helpful to save that set of commands for future or repeated use. The commands can be saved as scripts or functions.
- Change directory to the desktop. Open the scilab editor and type the commands (each on a new line): disp("Hello world"); a=1; b=2; c=a+b; disp(c); disp("Goodbye, world");
- Save the above file to the desktop. Then open the same file using a regular text editor such as notepad to show that it is indeed a text file. Now load the file into scilab using the scilab editor's execute menu option.
- Change the value of a to 5 in the editor, save and close it. Now execute the script directly from the scilab interpreter using exec.
- Functions: Show the syntax of functions, explain the function keyword, input arguments and the structure of output arguments when there is more than one argument.
- Show the following function in the editor: function [degrees] = radians2degrees(radians) degrees = radians*(180/%pi); endfunction. Save the function file and exec it from scilab. Then show radians2degrees(%pi/2), radians2degrees(%pi/4).
- Inline functions: show how to define degrees2radians, i.e. inverse of the above function, as an inline function.
- .sce versus .sci: These are just conventions- sci is used for scripts containing only function definitions- sce is used for all other scripts.
- Input and Output
- Console input and output:
- input function: age = input("Enter your age: "). Then try entering name: name = input("Enter your name"). Enter name without quotes- will give an error. Then enter name with qutoes- that'll work. Explain that if string input is expected, then you can use: name = input("Enter your name", "string")
- Output function: Show disp(age) and disp(name). Formatted output: mprintf("Your name is %s and your age is %d", name, age)
- Reading and writing in files: Do a = rand(20, 1)
- Text files:
- Output- write(): write('my-random-numbers.txt', a). Open and show the text file.
- Input- read(): b = read('my-random-numbers.txt', 20, 1). Now do b - a.
- Binary files: Back up a thus: a_backup = a. Now clear a.
- Output- save(): save('my-random-numbers.dat', a). Try opening in notepad- show it's really binary.
- Input- load(): load('my-random-numbers.dat', 'a'). Now do a_backup - a.
- Say that you can load only some variables if needed: save('personaldetails.dat', name, age), clear, disp(name) //error, disp(age) //error, load('personaldetails.dat', 'name'), disp(name) //no error, disp(age) //error
- Text files:
- Console input and output:
Advanced Level Scilab
- File Handling- Scilab File handling
- Writing to a file using write()
- Reading from a file using read()
- Opening an existing file using mopen()
- Closing an already opened file using mclose()
- File Handling- User Defined Input and Output in Scilab
- Input Function
- mprintf()
- save() and load()
- Used to quit scilab midway through calculation and continue at later stage
- Numerical methods- Integration
- Numerical methods- Solving Non- linear Equations
- Numerical methods- Gauss Elimination Method and Gauss Jordan Method
- Numerical methods- Iterative Methods
- Numerical methods- Interpolation
- Numerical methods- ODE- Euler methods
- Numerical methods- ODE- Applications
- Optimization Using Karmarkar Functions
- About Optimization
- Use of Scilab function Karmarkar in Optimization
- Optimization of Non-linear Functions
- Digital Signal Processing
- Filter Design- Windowing Techniques
- Filter Design- Window Based FIR Filter
- Filter Design- IIR Digital Filter
- Filter Design- Applying Digital Filter
- Control systems
- Discrete systems
- Calling User Defined Functions in XCOS
- Write a squaring function
- Use of scifunc block in XCOS
- Use of MUX block
- Call functions having multiple inputs and output
Xcos
Signals and Systems
Control Systems
Outline for a set of spoken tutorials on Control Systems
Representing control systems and their properties
- First order systems (1)
- Defining a system
- Simulating the system
- Step test
- Sine test (at two different frequencies)
- Root locus
- Bode plot
- Nyquist plot (1)
- Second order systems (2)
- Underdamped system
- Defining the system
- Roots of denominator
- Step test
- Underdamped system
- Other system representations (3)
- Discrete time systems (3)
- Defining a system
- Simulating the system (step test)
- Bode plot
- Root locus
- Discretizing continuous time systems (3)
- Discrete time systems (3)
- State space representation (4)
- Defining the system and initial state
- Simulating the system
- Plotting output state and internal state
- Root locus
- Conversion to transfer function
- Compare roots of denominator of answer with eigenvalues of A matrix.
- Multiple subsystems (5)
- Modeling a system (6)
- Obtain ODE representation of a system
- Laplace transform representation of the system
- Comparison with a second order model (obtain natural freq. and damping factor)
- Analyzing the system (7)
- Open loop step response
- Step response with unity feedback
- Steady state error: open loop system as well as unity feedback system
- Root locus and finding the gain from a point on the root locus.
- Controlling the system (8)
- Proportional control
- Deciding the parameters to be achieved (overshoot/rise time, etc)
- Determining corresponding system parameters (natural freq., damping factor)
- Obtaining location of corresponding point on root locus and gain there
- Building and simulating the proportional controller
- Other Controllers (9)
- PI Control
- PD Control
- PID Control
- Comparing all controllers
Robotics
- rtss Toolbox
Signal Processing
- Signal Processing Basics
- Audio processing
- Image Processing
- Fourier Transforms
Other Topics
Ordinary Differential Equations
(Spoken Tutorial-1 Scripts)- Introduction to ODEs
Use <math>\frac{dx}{dt} = sin(2t)</math> as the first order example and solve using ode.
- System:
- 2nd order differntial equation, Lotka Volterra equation example,
- Using champ to show arrows.
- With input (nonlinear input/state system). van der Pol oscillator?
- 3rd order, show example and plot.
- One example for discrete time system (difference equation) (of order 2)
- Show points evolving in plane.
- Linear system: syslin
- Solution:
- Using ODEOPTIONS
- Implicit Differential Equations
Boundary Value Problems
Difference Equations
Differential Algebraic Equations
Partial Differential Equations
Linear Algebra
Revise linsolve.
spec, det, generalized eigenvalue (spec)
inv and pinv (pseudo inverse), svd, companion-form
trace, sum of eigenvalues almost same as trace (similarly determinant and product of eigenvalues)
Optimization
Probability and Statistics
Scilab Toolboxes
Symbolic toolbox Scripts Scilab Wavelet toolbox Scripts
Special Functions
- Bessel functions
- Beta function
- Definition
- beta(m,n)=beta(n,m)
- beta(m,n)=gamma(m)*gamma(n)/gamma(m+n)
- Gamma function
- Definition
- gamma(1), gamma(0)
- gamma(n+1)=n!
- gamma(11/2), gamma(1/2)
- Error function
- Legendre functions
Least squares fit, splines
Polynomial interpolation
Close relation between interpolation and van der Monde matrix
Subset of columns of van der Monde matrix will correspond to lower degree fit.
Solvability of over-determined equations perhaps not possible.
Pseude inverse corresponds to "least squares fit" of specified degree.
Data Handling
- Databases
- XMLLab
==Web Interface for Scilab (W3 Scilab)==
Contributors and Content Editors
Lavitha Pereira, Minal, Nancyvarkey, PoojaMoolya, Pratik kamble