Numerical solution of fuzzy differential equations by Milne s predictor-corrector method and the dependency problem

Similar documents
Application of Bernoulli wavelet method for numerical solution of fuzzy linear Volterra-Fredholm integral equations Abstract Keywords

Method For Solving Fuzzy Integro-Differential Equation By Using Fuzzy Laplace Transformation

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

On Control Problem Described by Infinite System of First-Order Differential Equations

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

Overview. Overview Page 1 of 8

AN EVOLUTIONARY APPROACH FOR SOLVING DIFFERENTIAL EQUATIONS

Fuzzy Laplace Transforms for Derivatives of Higher Orders

STUDY OF THE STRESS-STRENGTH RELIABILITY AMONG THE PARAMETERS OF GENERALIZED INVERSE WEIBULL DISTRIBUTION

An Automatic Door Sensor Using Image Processing

Fuzzy Hv-submodules in Γ-Hv-modules Arvind Kumar Sinha 1, Manoj Kumar Dewangan 2 Department of Mathematics NIT Raipur, Chhattisgarh, India

Int. J. Computers & Electrical Engineering, vol. 30, no. 1, pp , AN OBSERVER DESIGN PROCEDURE FOR A CLASS OF NONLINEAR TIME-DELAY SYSTEMS

CS537. Numerical Analysis

Comparison between the Discrete and Continuous Time Models

Heat Conduction Problem in a Thick Circular Plate and its Thermal Stresses due to Ramp Type Heating

Representing Knowledge. CS 188: Artificial Intelligence Fall Properties of BNs. Independence? Reachability (the Bayes Ball) Example

7 Wave Equation in Higher Dimensions

On two general nonlocal differential equations problems of fractional orders

An S-type singular value inclusion set for rectangular tensors

MEI Structured Mathematics. Module Summary Sheets. Numerical Methods (Version B reference to new book)

û s L u t 0 s a ; i.e., û s 0

Final Exam. Tuesday, December hours, 30 minutes

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

Numerical solution of the first order linear fuzzy differential equations using He0s variational iteration method

Additional Exercises for Chapter What is the slope-intercept form of the equation of the line given by 3x + 5y + 2 = 0?

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

The sudden release of a large amount of energy E into a background fluid of density

TMA4329 Intro til vitensk. beregn. V2017

COSC 3361 Numerical Analysis I Ordinary Differential Equations (I) - Introduction

Some algorthim for solving system of linear volterra integral equation of second kind by using MATLAB 7 ALAN JALAL ABD ALKADER

Higher Order Difference Schemes for Heat Equation

Numerical Solution of Linear Emden Fowler Boundary Value Problem in Fuzzy Environment

ln y t 2 t c where c is an arbitrary real constant

International Journal of Mathematical Archive-5(6), 2014, Available online through ISSN

On The Estimation of Two Missing Values in Randomized Complete Block Designs

Extremal problems for t-partite and t-colorable hypergraphs

Chapter Finite Difference Method for Ordinary Differential Equations

Online Completion of Ill-conditioned Low-Rank Matrices

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

Fig. 1S. The antenna construction: (a) main geometrical parameters, (b) the wire support pillar and (c) the console link between wire and coaxial

Numerical Solution of Fuzzy Fractional Differential Equations by Predictor-Corrector Method

Discretization of Fractional Order Differentiator and Integrator with Different Fractional Orders

CONSIDERATIONS REGARDING THE OPTIMUM DESIGN OF PRESTRESSED ELEMENTS

Computer Propagation Analysis Tools

Probabilistic Models. CS 188: Artificial Intelligence Fall Independence. Example: Independence. Example: Independence? Conditional Independence

Degree of Approximation of a Class of Function by (C, 1) (E, q) Means of Fourier Series

In many engineering and other applications, the. variable) will often depend on several other quantities (independent variables).

Lecture 20: Riccati Equations and Least Squares Feedback Control

Research on the Algorithm of Evaluating and Analyzing Stationary Operational Availability Based on Mission Requirement

EFFECT OF PERMISSIBLE DELAY ON TWO-WAREHOUSE INVENTORY MODEL FOR DETERIORATING ITEMS WITH SHORTAGES

Department of Chemical Engineering University of Tennessee Prof. David Keffer. Course Lecture Notes SIXTEEN

02. MOTION. Questions and Answers

The Production of Polarization

A STOCHASTIC MODELING FOR THE UNSTABLE FINANCIAL MARKETS

BMOA estimates and radial growth of B φ functions

The Combined Effect of Chemical reaction, Radiation, MHD on Mixed Convection Heat and Mass Transfer Along a Vertical Moving Surface

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

NUMERICAL SIMULATION FOR NONLINEAR STATIC & DYNAMIC STRUCTURAL ANALYSIS

Stress Analysis of Infinite Plate with Elliptical Hole

Servomechanism Design

DonnishJournals

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

Ordinary dierential equations

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997)

Numerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method

Functions Defined on Fuzzy Real Numbers According to Zadeh s Extension

An Application of Fuzzy Linear System of Equations in Economic Sciences

Tutorial on Strehl ratio, wavefront power series expansion, Zernike polynomials expansion in small aberrated optical systems By Sheng Yuan

Faα-Irresolute Mappings

Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

Mathematics Paper- II

, the tangent line is an approximation of the curve (and easier to deal with than the curve).

Low-complexity Algorithms for MIMO Multiplexing Systems

Exponential and Logarithmic Equations and Properties of Logarithms. Properties. Properties. log. Exponential. Logarithmic.

AN EFFICIENT INTEGRAL METHOD FOR THE COMPUTATION OF THE BODIES MOTION IN ELECTROMAGNETIC FIELD

An introduction to the theory of SDDP algorithm

ON DIFFERENTIABILITY OF ABSOLUTELY MONOTONE SET-VALUED FUNCTIONS

Review - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y

Hydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods

336 ERIDANI kfk Lp = sup jf(y) ; f () jj j p p whee he supemum is aken ove all open balls = (a ) inr n, jj is he Lebesgue measue of in R n, () =(), f

International Journal of Pure and Applied Sciences and Technology

On Energy-Efficient Node Deployment in Wireless Sesnor Networks

Estimation of Lag Time Between Onset of and Death from an Occult. Tumor

Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1

Dynamic Estimation of OD Matrices for Freeways and Arterials

Linear Dynamic Models

2-d Motion: Constant Acceleration

New and Faster Filters for. Multiple Approximate String Matching. University of Chile. Blanco Encalada Santiago - Chile

The Global Trade and Environment Model: GTEM

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra

A Novel Iterative Numerical Algorithm for the Solutions of Systems of Fuzzy Initial Value Problems

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2 " P 1 = " #P L L,

Monochromatic Wave over One and Two Bars

Problem Set 5: Universal Law of Gravitation; Circular Planetary Orbits

A HYBRID DEA/AHP MODEL OF R&D PROJECT SELECTION CONSIDERING ORDINAL EVALUATION FACTORS

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Bayes Nets. CS 188: Artificial Intelligence Spring Example: Alarm Network. Building the (Entire) Joint

International Journal of Mathematical Archive-3(12), 2012, Available online through ISSN

( ) exp i ω b ( ) [ III-1 ] exp( i ω ab. exp( i ω ba

Transcription:

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 Numeical soluion o uzz dieenial equaions b Milne s pedico-coeco meod and e dependenc poblem Kanagaajan K Indakuma S Muukuma S Depamen o MaemaicsSi Ramakisna mission Vidalaa College o As & Science Coimbaoe 6400 ABSTRACT Te sud o is pape suggess on dependenc poblem in uzz compuaional meod b using e numeical soluion o Fuzz dieenial equaions(fdes) in Milne s pedico-coeco meod. Tis meod is adoped o solve e dependenc poblem in uzz compuaion. We solve some uzz iniial value poblems o illusae e eo. KEYWORDS Fuzz iniial value poblem Dependenc poblem in uzz compuaion Milnes pedico-coeco meod.. INTRODUCTION Fuzz Dieenial Equaions (FDEs) ae used in modeling poblems in science and engineeing. Mos o e poblems in science and engineeing equie e soluions o FDEs wic ae saisied b uzz iniial condiions eeoe a Fuzz Iniial Value Poblem(FIVP) occus and sould be solved. Fuzz se was is inoduced b Zade[]. Since en e eo as been developed and i is now emeged as an independen banc o Applied Maemaics. Te elemena uzz calculus based on e exension pinciple was sudied b Dubois and Pade [4]. Seikkala[] and Kaleva[6] ave discussed FIVP. Buckle and Feuing[] compaed e soluions o FIVPs wic wee obained using dieen deivaives. Te numeical soluions o FIVP b Eule's meod was sudied b Ma e al.[8]. Abbasband and Allvianloo [ ] poposed e Talo meod and e ou ode Runge-Kua meod o solving FIVPs. Palligkinis e al.[0] applied e Runge-Kua meod o moe geneal poblems and poved e convegence o n-sage Runge-Kua meod. Allavianloo e. al.[8] and Banabas Bed [0] o solve e numeical soluion o FDEs b pedico-coeco meod. Te dependenc poblem in uzz compuaion was discussed b Amad and Hasan[4] and e used Eule's meod based on Zade's exension pinciple o inding e numeical soluion o FIVPs. Oma and Hasan[7] adoped e same compuaion meod o deive e ou ode Runge-Kua meod o FIVP. Lael Amad and Hasan[4] invesigae e dependenc poblem in uzz compuaion based on Zade exension pinciple. In is pape we sud e dependenc poblem in uzz compuaions b using Milne's pedico-coeco meod. 65

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04. PRELIMINARY CONCEPTS In is secion we give some basic deiniions. Deiniion. Subse à o a univesal se Y is said o be a uzz se i a membesip uncion µ à () akes eac objec in Y ono e ineval [0]. Te uncion µ à () is e possibili degees o wic eac objec is compaible wi e popeies a caaceized e goup. A uzz se can also be pesened as a se o odeed pais () Te suppo e coe and e eig o A ae especivel :> () : () sup. (4) Deiniion. A uzz numbe is a convex uzz subse A o R o wic e ollowing condiions ae saisied: (i) is nomalized. i.e. ; (ii) ae uppe semiconinuous; (iii) : ae compac ses o 0< and (iv) : ae also compac ses o 0<. Deiniion. I is e se o all uzz numbes and we can caaceize b is α-levels b e ollowing closed-bounded inevals: [] : [ ] 0< (5) [] : [ ] 0< (6) Opeaions on uzz numbes can be descibed as ollows: I en o 0<. [ ];. [] [ ];. [ ]; 4. ee 0 [] ; 5. [] [] wee s is scala and 6. o 0<. Deiniion.4 A uzz pocess is a mapping : wee I is a eal ineval [7]. Tis pocess can be denoed as: [] [ ] 0<. (7) 66

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 Te uzz deivaive o a uzz pocess x() is deined b [] 0<. (8) Deiniion.5 Tiangula uzz numbe ae ose uzz ses in F(R) in wic ae caaceized l c l c l b an odeed iple ( ) R wi U 0 and c [ U ] [ ] o an en α R. FUZZY INITIAL VALUE PROBLEM Te FIVP can be consideed as ollows suc a [ ] [ ] [ ] α c c l c c [ U ] ( α )( ) ( α )( ) d( ) ~ ( ( ) ) (0) Y 0 (0) d Wee : is a coninuous mapping and wi -level ineval [ ] [ ] 0<. () Wen is a uzz numbe e exension pinciple o Zade leads o e ollowing deiniion: I ollows a () [] [ ] 0< () Wee [ ] 0< (4) [ ] 0<. (5) Teoem. Le sais 0 (6) Wee : is a coninuous mapping suc a is non deceasing e IVP 0 (7) (9) 67

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 68 Has a soluion on o >0 and a 0 is e onl soluion o equaion (7) o 0. Ten e FIVP (0) as a unique uzz soluion. Poo.See [7] In e uzz compuaion e dependenc poblem aises wen we appl e saigowad uzz ineval aimeic and Zade's exension pinciple b compuing e ineval sepaael. Fo e dependenc poblem we ee [7]. 4. THE MILNE S PREDICTOR-CORRECTOR METHOD IN DEPENDENCY PROBLEM We conside e IVP in equaion (0) bu wi cisp iniial condiion and [ ]. Te omula o Milne s pedico-coeco meod is ollows: [ ] [ ] ) ( ) ( ) ( ) ( 4 4 P C P (8) Wee. 0 0 N N T K We conside e ig-and side o equaion (8) we modi e Milne s pedico-coeco meod b using dependenc poblem in uzz compuaion as one uncion [ ] 4 ) ( P V (9) B e equivalen omula [ ]. 4 4 C (0) Now le e omula ; 0 () Can exend equaion (0) in e uzz seing. Le [ ] epesen e -level o e uzz numbe deined in equaion (). We ewie equaion () using e -level as ollows:

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 [ ] () B appling equaion () in (8) we ge [ ] () Wee Teeoe [ ] (4) [ ]. (5) α α [ ( ( )) 4 ( ( )) ( ( ))] [ ] α min ( ) P (6) α α [ ( ( )) 4 ( ( )) ( ( ))] [ ] max ( ) P α (7) B using e compuaional meod poposed in [5] we compue e minimum and maximum in equaions (6) (7) as ollows (8) 9 5. NUMERICAL EXAMPLES In is secion we pesen some numeical examples including linea and nonlinea FIVPs. Example 5. Conside e ollowing FIVP. [0]; 0 <0.5; 4 0.5 0.5; (0) 0 >0.5; 69

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 Te exac soluion o equaion (0) is given b []. () Te absolue esuls o e numeical uzz Milne's pedico-coeco meod appoximaed soluions a 0. See Table and Figue and. TABLE Te eo o e obained esuls wi e exac soluion a. Figue : Te appoximaion o uzz soluion b Milne's pedico-coeco (0.) 70

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 Figue : Compaison beween e exac Milne's pedico-coeco pedico-coeco In is example e compaison o e absolue local eo beween Milne's pedico-coeco meod wi e uzz exac soluion is given in Table o vaious values o α -level α 00. K0.9 and ixed value o ( 0 ). Te esuls sows a Milne's pedicocoeco meod is moe accuae an pedico-coeco meod sows e gapical compaison o a uzz soluion beween exac Milne's pedico-coeco pedico-coeco a ixed ( 0 ). Te beaviou soluions o e end poins o e uzz inevals o a uzz exac soluion Milne's pedico-coeco and pedico-coeco uzz appoximaed soluions ae ploed and compaed in Figue a α 0. Figue cleal sow a Milne's pedicocoeco povides a moe accuae esuls an pedico-coeco meod. Example 5. Conside e ollowing FIVP 4 [0]; 0 <0.5; 4 0.5 0.5; 0 >0.5; () Te exac soluion o equaion () is given b []. () Te absolue esuls o e numeical uzz Milne's pedico-coeco meod appoximaed soluions a 0. See Table and Figue and 4. 7

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 TABLE Te eo o e obained esuls wi e exac soluion a. Figue : Te appoximaion o uzz soluion b Milne's pedico-coeco (0.) 7

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 Figue 4: Compaison beween e exac Milne's pedico-coeco pedico-coeco In is example we compae e soluion obained b Milne's pedico-coeco meod wi e exac soluion and pedico-coeco. We ave given e numeical values in Table ixed value o 0. and o dieen values o α. 6. CONCLUSION In is pape we used e Milne's pedico-coeco meod o solving FIVP b consideing e dependenc poblem in uzz compuaion. We compaed e soluions obained in wo numeical examples. ACKNOWLEDGEMENTS Tis wok as been suppoed b Tamilnadu Sae Council o Science and Tecnolog Tamilnadu India. REFERENCES [] S. Abbasband T. Allavianloo Numeical soluions o uzz dieenial equaions b Talo Meod Compuaional Meods in Applied Maemaics (00) -4. [] S. Abbasband T. Allavianloo Numeical soluion o Fuzz dieenial equaion b Runge- Kua meod Nonlinea Sudies (004) 7-9. [] M. Aamad M. Hasen A new appoac o incopoae unceaini ino Eule meod Applied Maemaical Sciences 4(5) (00) 509-50. [4] M. Aamed M. Hasan A new uzz vesion o Eule's meod o solving dienial equaions wi uzz iniial values Sians Malasiana 40 (0) 65-657. [5] M. Amad M. Hasan Incopoaing opimizaion ecnique ino Zade's exension pinciple o compuing non-monoone uncions wi uzz vaiable Sains Malasiana 40 (0) 64-650. [6] N. Z. Amad H. K. Hasan B. De Baes A new meod o compuing coninuous uncion wi uzz vaiable Jounal o Applied Sciences (7) (0) 4-49. [7] A. H. Alsonosi Oma Y. Abu Hasan Numeical soluion o uzz dieenial equaions and e dependenc poblem Applied Maemaics and Compuaion 9 (0) 6-7. [8] T. Allavianloo N. Amad E. Amad Numeical soluions o uzz dieenial equaions b pedico-coeco meod Inomaion Sciences 77(7) (007) 6-647. 7

Applied Maemaics and Sciences: An Inenaional Jounal (MaSJ ) Vol. No. Augus 04 [9] T. Allavianloo N. Amad E. Amad Eaum o Numeical soluions o uzz dieenial equaions b pedico-coeco meod Inomaion Sciences 77(7) (007) 6-647 Inomaion Sciences 78 (008) 780-78. [0] Banabas Bede Noe on Numeical soluions o uzz dieenial equaions b pedico-coeco meod Inomaion Sciences 78 (008) 97-9. [] A. Bonaini G. Bonempi A Qualiaive simulaion appoac o uzz dnamical modelsacm Tans. Model. Compu. Simula 4 (994) 85-. [] R. en Algoims o Minimizaion wiou Deivaives Dove Pubns 00. [] J. J. Buckle T. Feuing Fuzz dieenial equaions Fuzz Ses and Ssems0 (000) 4-54. [4] D. Dubois H. Pade Towads uzz dieenial calculus pa : dieeniaion Fuzz Ses and Ssems 8 (98) 5-. [5] E. Isaacson H. B. Kelle Analsis o Numeical Meods Wile New Yok 966. [6] Kaleva Osmo Fuzz dieenial equaions Fuzz Ses and Ssems 4 (987) 0-7. [7] O. Kaleva Inepolaion o uzz daa Fuzz Ses and Ssems 60 (994) 6-70. [8] M. Ma M. Fiedman A. Kandel Numeical soluions o uzz dieenial equaions Fuzz Ses and Ssems 05 (999) -8. [9] R. E. Mooe Ineval Analsis Paenice-Hall Englewood clis N. J 966. [0] S. Palligkinis G. Papageogiou I. Famelis Runge-Kua meods o uzz dieenial equaions Applied Maemaics and Compuaion 09 (009) 97-05. [] S. Seikkala On e uzz iniial value poblem Fuzz Ses and Ssems 4() (987) 9-0. [] L. A. Zade Fuzz Ses Inomaion and Conol 8(965) 8-5. 74