jagomart
digital resources
picture1_Processing Pdf 176835 | 4354 Item Download 2023-01-28 20-28-18


 136x       Filetype PDF       File size 0.70 MB       Source: download.atlantis-press.com


File: Processing Pdf 176835 | 4354 Item Download 2023-01-28 20-28-18
proceedings of the 2012 2nd international conference on computer and information application iccia 2012 an application of matlab in practical problems xie huiyang bi qiuxiang college of sciences fixed income ...

icon picture PDF Filetype PDF | Posted on 28 Jan 2023 | 2 years ago
Partial capture of text on file.
                              Proceedings of the 2012 2nd International Conference on Computer and Information Application (ICCIA 2012)
                                      An Application of MATLAB in practical problems 
                                                    
                                   XIE Huiyang                                                        Bi Qiuxiang 
                                 College of Sciences                                            Fixed-Income Department 
                             Beijing Forestry University                                         GF securities CO., LTD 
                                   Beijing, China                                                  Guangzhou, China 
                                xhyang@bjfu.edu.cn                                                  gfbqx@tom.com
                                                                                                  12ŠŠ1                  101
              Abstract—This article gives an example of using Matlab to                         
              solve practical problems                                                          
                                                                                           AB==012, 022
                                                                                     Let                                             ,  
                                                                                                
                 Keywords-Linear Algebra; Matlab                                                
                                                                                                  Š364                  351
                                                                                                
                                 I.  INTRODUCTION                                    calculate  AB ,  AŠ1B. 
                 Matlab(Matrix Laboratory) is a commercial                          
              mathematical software developed by MathWorks for                     >> A = [1 2 -1; 0 1 2; -3 6 4]; 
              numerical computation. One of the nice features of Matlab            >> B = [-1 0 1; 0 2 2; 3 5 1]; 
              is its ease of computations with vectors and matrices.               >> A * B 
              Matrix operations, numerical analysis, signal processing,            Result of  AB  returned: 
              graphics, and user interface creation can be easily                  ans = 
              manipulated in Matlab. It also provides many                             -4    -1     4 
              implementations of algorithms and interface to other                      6    12     4 
              programming language. With help of additional toolboxes                  15    32    13                 
              and package, Matlab even supports symbolic computation,              >> inv(A) * B 
              graphics manipulation and model-base designing for                              AŠ1B
              dynamic and embedded systems.                                        Result of         returned: 
                 Matlab is widely used in all areas of science and                 ans = 
              technology in education and research at universities. Many               -1.0000    0.1304    1.3478 
              courses of mathematics and engineering are using Matlab                   0        0.3478    0.2609 
              for education on campus.                                                  0        0.8261    0.8696 
                 Linear algebra is a branch of mathematics, which is               (The numbers could be expressed as fractions either.)   
              central to modern mathematics and incredibly useful in the            
              modern world. It plays important roles in the researches of          Sample 2. 
              mathematics, mechanics, physics and other subjects, and it           Find the solution of equations 
              is widely applied to solve the real world problems. For                +++=
                                                                                   xxxx
                                                                                       22 0
              most non-linear problems, the usual method is to                   1234
              approximate them by linear problems. People are only                    +ŠŠ=
                                                                                    xxxx
                                                                                2220
              really good at solving linear problems by now due to the           1234
                                                                                      ŠŠŠ=
                                                                                   xxxx
              development and widespread use of the computer in areas                        430
              that apply mathematics. By using Matlab, the complicated           1234. 
              calculations could be avoided. With Matlab, we can utilize            
              linear algebra into real life to solve practical problems. In        >> A = [1  2  2  1; 2  1  -2  -2; 1  -1  -4  -3];  
              the past two decades, the applications of linear algebra to           
              real world problems have mushroomed. In this paper,                  >> format rat  // Let the solution outputs in the 
              some examples on solving linear problems with Matlab are          rational form. 
              illuminated and how to use Matlab in practice is discussed.          >> B = null(A, 'r')  // Get the rational basis of space 
                                                                                A. 
                      APPLICATION IN LINEAR ALGEBRA OPERATIONS                     B = 
                 II.                                                                  2            5/3 
                   Linear algebra operations such as Matrix                           -2           -4/3 
              multiplication, finding inverse of a matrix, solving a                  1             0 
              system of linear equations and diagonalizing quadratic                  0             1 
              forms would lead to boring and complicated calculations.             The column vectors in B form the basic set of 
              However, these operations can be performed very simply               solutions. Apply the following commands to get the 
              using matlab.                                                     general solution. 
                   Sample 1.                                                        
                                                                                   >> syms k1 k2    // define the symbols k1 and k2 
                                                                                   >> X = k1 * B(: , 1) + k2 * B(: , 2)  // write the 
                                                                                general solution 
                                                          Published by Atlantis Press, Paris, France. 
                                                                      © the authors 
                                                                          1631
                              Proceedings of the 2012 2nd International Conference on Computer and Information Application (ICCIA 2012)
                 >> pretty(X)                       // make X                      >> zeroelement(1, 3) 
              more pretty                                                          >> zeroelement(1, 4) 
                                                                                   >> wamparray(2, 3) 
                 We usually reduce coefficient matrix into row echelon             >> zeroelement(2, 3) 
              form to use general solution to solve linear equations.              >> zeroelement(2, 4) 
              Command “rref” would provide row echelon form matrix                 >> zeroelement(3, 4) 
              in Matlab.                                                            
                                                                                   Each step implemented, the result will be displayed to 
                 >> B = rref(A)                                                help users observe. You can attach a semicolon at the end 
                 B =                                                           of the row to omit the intermediate output. 
                      -5/3    1         0       -2                                III.  APPLICATION IN SOLVING PRACTICAL PROBLEMS 
                       4/3    0        1        2                
                       0      0        0        0                                  Sample 1. There’re five cities with its own airport and 
                                                                               some one-way airlines between the cities. As figure i 
                 To observe intermediate data, sometimes, the                  shows, a directed edge from node A to node B represents 
              user-defined functions are needed, showing the                   an airline from city A to city B. 
              information expected. In sample 3, reducing matrix to its 
              row echelon form is implemented by user-defined M 
              function. 
                  
                 Sample 3. 
                            14115
                                                                                                                           
                                                                                                      Figure i 
                             2867
                      A= The cities and the airlines could be represented by an 
                           
                             32128                                             adjacency matrix            . If there’s an airline from city 
                                                                                                 A=()a
                                                                                                      ij
                            715660
                            i to city j, let a =1, otherwise a =0. Figure I can be 
                 Let                           , get the row echelon                           ij                 ij
              form of A.                                                       represented as matrix A. 
                                                                                           01101
                 We will program two M-functions for the intermediate                                       
              data. Given the matrix A and row index i and column                                           
                                                                                           10011
              index j, the M-function named zeroelement.m is going to                                       
                                                                                     A=                     
                                       a =0               a ≠0                             10010
              reduce A by making  ij           ,  when  ii         and                                      
              i < j                                                                        01100
                    .                                                                                       
                                                                                                            
                 function    [A]= zeroelement(A, i, j)                                     01000
                 if  i>j | i==j                                                                               
                    error('i is not less than j')                                  In  A2 , aij is the number of route from city i to city j 
                 else                                                          though exact one transfer. Enter the following commands 
                       A(j, :)=A(j, :)-A(i, :)*(A(j, i)/A(i, i));              in Matlab. 
                 end                                                               A=[0,1,1,0,1;1,0,0,1,1;1,0,0,1,0;0,1,1,0,0;0,1,0,0,0]
                                                                                A1*= AA
                 The other function named swaparray.m will swap the                A1 =  
              i-th row and j-th row in A. 
                                                                                        2     1     0     2     1
                 function    [A] = swaparray(A, i, j)   
                 if  i>n|j>n                                                            0     3     2     0     1
                    error('out of range') 
                 else                                                                   0     2     2     0     1
                    temp1=A(i, :); 
                    temp2=A(j, :);                                                      2     0     0     2     1
                    A(i, :)=temp2; 
                    A(j, :)=temp1;                                                      1     0     0     1     1 
                 end                                                               Let’s take an element in A1 for verification. A1[2][3] = 
                                                                               2 means there’re two different strategies to reach city 3 
                 Invoking these functions repeatedly, we’ll finally get        from city 2 though a transfer. Enumerating all the paths, 
              the row echelon form of matrix A.                                we can find exactly two ways which are 2→4→3, 
                                                                               2→1→3. 
                 >> A=[1 4 11 5; 2 8 6 7; 3 2 12 8; 1 15 66 0];                    The matrix keeps the number of routes between each 
                 >> zeroelement(1, 2) 
                                                          Published by Atlantis Press, Paris, France. 
                                                                      © the authors 
                                                                         1632
                              Proceedings of the 2012 2nd International Conference on Computer and Information Application (ICCIA 2012)
              two cities though at most one transfer could be simply                                         1
                                          2                                                            
                                                                                                         10
              described as B =+AA. Enter the command                                                   
                                                                                            a                2
                                                                                            
              BA=+A1                                                                          n
                            in Matlab.                                               ()                    1
                                                                                      n                
                                                                                                 ,01
                                                                                   XbM
                 B =                                                                     ==
                                                                                              n
                                                                                            
                       2     2     1     2     2                                                           2
                                                                                            c          
                                                                                              n
                                                                                            000
                       1     3     2     1     2                                                       
                                                                                                       
                                                                                                       
                       1     2     2     1     1                                                                      
                                                                                     ()nn(Š1)               ()nn(0)
                       2     1     1     2     1                                   X     =MX             X     =MX 
                                                                                  There’re three eigenvalues to matrix M, so M can be 
                       1     1     0     1     1                               diagonalized which means there is a diagonal matrix D and 
                 Rest can be deduced by analogy.                               a invertible matrix P satisfy M = PDPŠ1 .So 
                 Sample 2.                                                        ()nnŠ1 (0)
                                                                                XP= DPX
                 In autosomal genetic, individual inherits one gene for                              should be satisfied. The diagonal 
              each gene pair from its parents to form its own unique           elements in D are the eigenvalues of M and P is the 
              gene pairs. The inherited gene is selected randomly. If a        eigenvector correspondingly. 
              parent has gene pair Aa, the descendant will have the same          First of all, we’ll calculate the eigenvalues and 
              probabilities to inherit gene A or gene a. For a parent with     eigenvector of M. 
              gene pair aa and another with gene pair Aa, the descendant          >> A = [1 0.5 0;0 0.5 1;0 0 0]; 
              can be with gene pair Aa or aa in the same probability.             >> [P, D] = eig(A);  // eigenvalues and eigenvector of 
              Enumerating all the parents gene pairs, the probabilities of     matrix A 
                                                                                        1 Š0.7071    0.4082         1   0   0
              each case are listed below.                                              
                                                                                       
                                                                                   PD=Š0    0.7071    0.8165 ,   =0 0.5 0
                                                                                       
                                      Parents’ gene pairs                              
                                                                                        0     0      0.4082         0   0   0
                                                                                       
                                AAAŠ a          AaAŠ a           aaŠaa                                                         ,  
                        AAAŠ A          AAaŠ a           AaaŠ a                   Combined with the results returned, the distributions at 
                                                                               each generation can be calculated by 
             s  DescenAA  1 1/2 0 1/4 0 0 
              gene                                                                                                           a
                                                                                                                            
                                                                                                                              0
                    Aa                                                                                                 (0)  
                                                                                                                     Xb
                d         0 1/2 1 1/2 1/2 0                                                                               =   0
             p  a                                                                                                           
             ai nt’                                                               ()nnŠ1 (0)                                
             r      aa                                                                                                       c
                                                                                                                              0
                                                                                XP= DPX                                     
                          0 0 0 1/4 1/2 1                                                                , with                     and 
                                                                               ab++c=1
                 There’re some crops consist of three                            000. The multiplication of matrix in Matlab is 
              genotypes              based on a certain distribution. All      helpful here.   
                        AAA,,aaa
              the crops are pollinated by the crops with  AA genotype.                                 UMMARY 
                                                                                                      S
              Let ab,,c  be the portion of crops with                             Linear algebra is an efficient tool to help do scientific 
                    nnn
              genotype              in n-th generation.               .        calculations on immense amounts of data, especially when 
                       AAA,,aaa                          (0n =  ,1,2)          it is combined with computers. This article gives some 
              ab,,c are the original distribution and                          examples to demonstrate the advantages of using Matlab 
               000                                                             in linear algebra. Teaching linear algebra with Matlab not 
              satisfying              . The n-th generation distribution 
                       ab++c=1
                        000                                                    only enhances students’ learning interest, but also 
              satisfies the following equations.                               improves their ability of applying theories to solve 
                             1                                                practical problems.   
                   aa=+b
                  nnŠŠ11n
                              2                                                                    CKNOWLEDGMENT 
                            1                                                                   A
                 
                   bc=+b                                                          Funded projects: (863) "large number DNA 
                  nnŠŠ11n
                            2                                                 decomposition computer model research project (Code: 
                 c =0                                                         2009AA01Z413) 
                  n
                                                                                                      EFERENCES 
                 Let                                                                                 R
                                                                               [1] Chen Huaichen, Gong Jiemin, and MATLAB introductory linear 
                                                                               algebra practice (Second Edition), publishing house of electronics 
                                                                               industry, 2009 
                                                                               [2] Zhou Jianxing, will Tsu Take, MATLAB Xingming, from entry to 
                                                                               the master, the people post and Telecommunications Press, 2008 
                                                                               Ḥ
                                                          Published by Atlantis Press, Paris, France. 
                                                                     © the authors 
                                                                         1633
The words contained in this file might help you see if this file matches what you are looking for:

...Proceedings of the nd international conference on computer and information application iccia an matlab in practical problems xie huiyang bi qiuxiang college sciences fixed income department beijing forestry university gf securities co ltd china guangzhou xhyang bjfu edu cn gfbqx tom com abstract this article gives example using to solve ab let keywords linear algebra i introduction calculate matrix laboratory is a commercial mathematical software developed by mathworks for numerical computation one nice features b its ease computations with vectors matrices operations analysis signal processing result returned graphics user interface creation can be easily ans manipulated it also provides many implementations algorithms other programming language help additional toolboxes package even supports symbolic inv manipulation model base designing dynamic embedded systems widely used all areas science technology education research at universities courses mathematics engineering are campus bran...

no reviews yet
Please Login to review.