МИНИСТЕРСВО ОБРАЗОВАНИЯ И НАУКИ
РОССИЙСКОЙ ФЕДЕРАЦИИ
БЛАГОВЕЩЕНСКИЙ ГОСУДАРСТВЕННЫЙ
ПЕДАГОГИЧЕСКИЙ УНИВЕРСИТЕТ
Физико-математический факультет
Кафедра информатики
РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ, ИСПОЛЬЗУЮЩИХ МЕТОДЫ ИНТЕГРИРОВАНИЯ, В СРЕДЕ MATLAB
Курсовая работа
Выполнил: студент курса
Научный руководитель:
кандидат физико-
математических наук, доцент
Благовещенск 2008
СОДЕРЖАНИЕ
ВВЕДЕНИЕ.. 3
1. ИНТЕГРИРОВАНИЕ ФУНКЦИЙ В MATLAB.. 5
1.1 Численный метод. 9
1.2 Символьный метод. 11
2. MATLAB – СРЕДА МОДЕЛИРОВАНИЯ.. 15
3. РЕАЛИЗАЦИЯ ЭКОНОМИЧЕСКОЙ МОДЕЛИ ВЗАИМОРАСЧЁТОВ ПРЕДПРИЯТИЙ В СРЕДЕ MATLAB.. 16
ЗАКЛЮЧЕНИЕ.. 19
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ... 20
ПРИЛОЖЕНИЯ.. 21
ВВЕДЕНИЕ
Невозможно представить себе современную науку без широкого применения математического моделирования. Сущность этого метода состоит в замене реального объекта его «образом» - математической моделью. Этот метод позволяет быстро и «безболезненно» изменить объект, изучить его свойства и поведение в различных средах и т.д. Неудивительно, что математическое моделирование бурно развивается и проникает во все сферы знаний.
![](data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEAYABgAAD/2wBDAAoHBwkHBgoJCAkLCwoMDxkQDw4ODx4WFxIZJCAmJSMgIyIoLTkwKCo2KyIjMkQyNjs9QEBAJjBGS0U+Sjk/QD3/wAALCAD7ATgBAREA/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/9oACAEBAAA/APZqKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKQnAyelZs/iPSbc7ZNQtywOCqNvIP0GacuuWsn+pS5l/wBy3c/0pf7XiHW2vQPX7M/+FRN4l0uNlWW58ot08yNlH5kVetr22vULWlxDOoOCYnDAflU9FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFQ3V3BZQ+bcyrGmcAnufQDufYVnm61PUOLOBbKDP+uuVy7D/ZjHT6sfwoHh61lO7UHmv3/6eHyn4IML+laMFtBbLtghjiX0RQv8qlopCARgjI96oXGhabdPvks4hJ2kQbHH0ZcH9ag/s/UbHnT73z488wXh3cf7Ljkfjuqe21eN5kt7yJ7O5c4WOUjDn/ZYcN/P2rQoooooooooooooooooooooooooooqtf3i2Fo0xRpG4VI0xukY9FGe5qnp+lSGcX+qss18fuqOY7cH+FB/Nup+nFatFFFFFFFQXlnb39s1vdRLLE3VW/mPQ+9Z1lLdaZeiwvpGnt5M/Zbpz8x/6Zv6sOx7geo52KKKKKKKKKKKKKKKKKKKKKKKKKKKyYlF/4jnkcZj09VjjB6eY43M34KVA+rVrUVm61rUWh28M9zFK8MkqxM8YB8vccbjk9B3xRea1DZ6xY6aYpZJrzdtK42ptGTkk+gPTNaOR+VGR6ijI9RRketLRVXUrFNSsJbaQ43D5WHVGHKsPcHBpmkXjX+lW9xIMSMuJAOzjhh+YNXaKKKKKKKKKKKKKKKKKKKKKKKKKytJ/d6pq8Dfe89ZR7qyDB/NWH4Vq0Vm6xaNfxw2ptxLbyMyzEsBtUowzg9etY8Oj6r9u0K5ulikmtjIbqRXwOYyi4B69ifqaRtN1e1hhKKJ5I4JIZcMP9LkbAV2zwAMEnOeuKpWnhzUo5J4TG8MZmt5vllGJo0jCPGT2y2Wx0Nalzo91LZNbR2aeXHBIYWkn3vvYMNuSOAMggg+3aqNvpGuWXkS2cEKXEenw2+WkUq0qn5i3Gce45NdiM4G7Ge+KWkJwMngCsvw3zokcuMCeSWdf915GcfowrVooooooooooooooooooooooooorK1SN7S4j1W3DEwrsuI1UsZIuvAHdTyPxHetC3uIru3jnt5FkikUMjqcgipaKKKKKKKKyNUnN/cf2PayYkdQ10458qI9s9mboPbJ7VqRxpDEscahUQBVUdABT6KKKKKKKKKKKKKKKKKKKKKKKKKKyW0yfTZprjSSu2Ql5LNzhGbuVP8BP4g+neprTWrW5lEEha2uu9vONr/h2Ye4zWhRRRRRUVzdQWcJluZo4Yx1Z2Cj9azTe3eqgppqtb25BzdzJgn/cQ9fqePY1d0/ToNMtzFbhiWYvJI53PIx6sx7mrVFFFFFFFFFFFFFFFFFFFFFFFFFFFFQXdlbX8PlXcEcyZyA65wfUeh9xWbNp0umwSTWurS20EalmW6/fRqB7khgAP9qs7wt4m1PxHpP9ox2Fu1s0rpEwlKNIqnG/BBxkjpmtk31+n39Lds/885lP88VVsvET6ibgWumzubeZreT94g2uuMg8+4q0brU2GF06JSe7XHA/IVnWF1rGrXF5BcXEGnvay+W8cMfmMykZVwzcYI/2fUdq0LfQrSGZZ5hJdXC8iW5bzCv+7nhfwArSooooooooooooooooooooooooooooooorzD4vW3iPWWtdG0REltpY2mlhif8AevtI5YHjbyMc8n6V33h7S00Xw9YafEMLbwKnPrjk/nmrlzOLa1lnYErGhcgDJOBmvMvg/eap/aGtw6raXEH2yY3kZljKgsThx9R8vFepVga9v0jUINdj3eRGvk36j/niTkSY/wBg8/QtW8rK6hlIZWGQQcgilooooooooooooooooooooooooooooooqK5uI7S2luJ2CRRIXdj2AGTWR4ctJn+0axfxlL2/IYRt1ghH+rj/I5PuxrcoqhpP+on/6+pv/AEM1fpskaSxtHIodHBVlYZBB6g1g+GxJpU1xoM5LLajzLN2bJe3J4H1Q/L9NvrXQUUUUUUUUUUUUUUUUUUUUUUUUUUUUUUVz2rE63rUOjRk/Zrfbc37DuAcxxf8AAiMn2X3roaKKoaR/qLj/AK+pv/QzV+isTxLZTtDBqenoXv8ATmMkaD/lqh/1kf8AwIDj3ArUsryHUbKG7tX3wzIHRvUGp6KKKKKKKKKKKKKKKKKKKKKKKKKKKKpatqA0vTJ7rZvdFxHHnG9zwq/iSBUOg6W2l6cFnZZLyZjNdSgffkbk/gOg9gK06KKoaR/qLj/r6m/9DNX6KKwLCP8AsTxDNYgqtnqBa4tl6bJRzKg9jw4/4FW/RRRRRRRRRRRRRRRRRRRRRRRRRRRRXOxn/hIPErSfe07SnKoO0tzjk/RAcf7xPpXRUUUVQ0j/AFFx/wBfU3/oZq/RRWbrulf2tpzRxv5V1ERLbTY5ilH3T9Ox9QSKXQ9V/tfTEnaPyp1JiuIc8xSrwy/n+hBrRoooooooooooooooooooooooooorK8QX0lrZJb2kipfXri3tiRnDHq2P9lct+FW9M0+DStOgsrUERQrgFjksepJPckkkn1NWqKKKoaR/qLj/AK+pv/QzV+iiisC7CaH4iivtxS11JlguB/Cs2MRv7ZA2H/gNb9FFFFFFFFFFFFFFFFFFFFFFFFFISFBJIAHUmuf0Mf23qL6/Kv7na0FgpHSLPzSfVyP++QPWuhooooqhpP8AqJ/+vqb/ANDNX6KKKrajYQapp89ldLuhnQo3qPce46is/wAPahNMk+nX7Z1CwYRynp5qn7kg9mH6gjtWzRRRRRRRRRRRRRRRRRRRRRRRRWF4hY6jLDoMEpR7sF7kr1W3Bw30LZCj6n0rajjSGJI4lCRoAqqowAB0Ap9ULzWbW0nFvl57k9IIF3v+IH3R7nAqN5tXuAPs9tbWoPUzuXYD/dXjP40waTfSgfa9auTznbbxpEv8if1oHhy1P+tuL+b/AK6Xkhx+tInhjTowfLFymSW+W5kHJ5z97rS/2CU/1GqanEc5yZ/M/wDQw1At9atmzHe212mfuzxeWwH+8vH6U6TWWsiP7Ss5reM/8tk/exj6kcr9SMe9aMM0dxEssMiyRsMqyHIP40+sHxAn9m3Nvr8e4fZAY7pVGd8DHnj/AGThvoG9a3EdZEV0YMjDIYHIIp1FFFFFFFFFFFFFFFFFFFFFFV76+g02ylu7uQRwxKWZj6D+ZrN8O2U6rcanqCFL6/YOyH/ljGPuR/gOT7k1sO6xozyMFRRlmY4AHqaxhcXevf8AHozWumEf8fHSWcf7H91f9o8nt61p2djbWEPlWkKRJ1O0ck+pPUn3NWKKKKKKKyrnSHgkNzo7x2twTl4yP3U3+8B0P+0OfrU2naol88kEkbW95CB5tu5yVz0IP8SnsR/OrrKroVcBlYYIPQisHQGbSLyXQJ2JSIGWwdj9+DP3M+qE4+m2ugoooooooooooooooooooooornLi3t/E3iHy5laWx0lw2M/JJc8EfXYO3q3tXR1k6gDqWox6aMG3RRNdgjIZc/In4kEn2X3rVAwMDpS0hIAJJwB3oBDAFSCDyCKWiiiiiszWLaTy1vrNR9stcsox/rE/ijP1HT3Aq9bXEd3bRTwtujlUOp9QRkVm+IrKSa0jvbSIyX1g/nwKpwXxwyf8CXI/L0q/YX0Gp2EN5avvhmUMp/ofQjoRViiiiiiiiiiiiiiiiiiiiisnxDqcthZJDZANqF4/kWqH++erH/ZUZY/T3qzpOmRaPpkNlASyxjl2+87HlmPuSST9au1k6P8ANqGsO33zdBT9BGgH6VrUVha07tr2jwT5FhI0hc5wrSgDy1b2+8cHqQKydRvZbC+1WCxMq2bvaxlo2IEUsjlX2n+HC7ScdM9jUtxrk+j6PfvaQea9jdLDKZ5pJN5YoBtYkk/f6exFX7nxG73VvaabAJ55DKJC/wAixtGBlTnBydwOPTJqoniq6fW302SC3t3KQGNnYsgd1LMjMOM/KQMdavHVNS/4SI6cIbXZ5H2gHc27bvC49M45rcoorK8N8aOFH3UnnRP90SuF/QCtWudt/wDin/EjWp+XT9VdpYfSK46un0cfMPcN610VFFFFFFFFFFFFFFFFFFFFYGiSLrep3Gt7AbdQbaxYn7yA/PIPZmAA9lHrW/RWTIf7N8QCQj9xqAWMt2WZQcZ/3l4/4CPWtaimuiyKVdQynqCMikWGNIvLWNBHjG0Dj8qjeztpIRC9vE0Q5CFAVH4Uv2S38kReRF5QOdmwYz9Ka+n2kjMz2sDM+NxMYJOOmfpTja25uBcGCPzgMeZtG78+tTUVR1e+axsS0K77mUiKBP70h6fh3PsDUunWa6fp1vaqdwiQLu/vHufxNWao6xpw1TS5rYMEkIDROf4JAcq34ECmaFqh1bTFllj8q5jYxXMX/POVeGH07j1BFaNFFFFFFFFFFFFFFFFFFc/4lmlv5rfQLOQxzXoL3EinDRW4Pzke7cKPqT2rcggitbeOCCNY4o1CIijAUDgAVJRUN1bRXls8E67o3GCM4P4EdDWXBfXGjuLXVnaS3zthvj0I7LL/AHW/2uh9jxWyCCARyDS0UUUUUVT1HU4tORQyvLPJkRQRDLyH0A/qeBVfTrC5a5a/1NwbphiOFDmO3X+6PVj3bv2wK1KKK57UwdC1uLVoxizu2W3vlHRWJxHL+BIU+xHpXQ0UUUUUUUUUUUUUUUUVDeXcNhZzXV1II4IULu56ADrWX4ct55IpdWvkKXd+Q+xhgwxD/Vx/UA5PuTW1RRWRqPiKC1u/sFnG19qRGRbQnlB6u3RB7n8AarpoNzqibvEdytwjEN9igytuvs3eT/gXHtVv+x3tIwNJuntQDkQsPMi+gU8qPoRim/2lqFpxfaY8ij/lrZt5gP8AwE4Yfhn60qeJ9JJxJeJbtnGLgGLn/gQFXYr60n/1NzBJ3+SQGpDPEASZUAHcsKpza9pVvjzdRtFJ6L5ykn6DOTUA8QxTkiwtLy8PqkJRc/7z4FOEer3rqZpIrCHqUhPmSN7FiMD8AfrUcvhizeT7RDLcwXwBAu0lJkOexzkEdPlIx7VU/tnUvD5ZfEMaz2Q+7qNtGcKP+msYyV/3hkfSt+3uIbuBJ7aVJYnGVdGDKw9iKlorN8RW11eeHr+3sBCbmWFljEwypJHQ1wPwl8X6tq2parpXiGeR72Da6LKoVlA+VhjH+6fxNeoUUUUUUUUUUUVmza9ZpI0VuZLydTgx2qeYQfcjhfxIpi3GtXDArZWtrH386Yu35KMfrUrW2puBnUIUx/zztuv5sahOl6iTn+3roewgh/8AiKemn6lH/wAxiST/AK6W8f8A7KBVLWdJ1bUo7eP7TZSwxzLLJC0TRiXacqpOW4zg9OcVa/tme2ONS06eFf8AnrD++T8dvzD8qv2l7bX8Xm2k8cyZwSjZwfQ+lcf4/wDiXaeCxHbQxLeahICfJD4Ea9i349vrVfwXr2r+P/D3m3VyunrG3lzG1UiWTvkEjCAg9snjqK7PTdKs9It/IsbdIUJyxHLOfVieWPuauUUU1kVxhlBHuKpy6Jpc/wDrtNs5Mc/PAp/pUY8O6MCCNI08EdxbJ/hVyKzt4M+TbxR5OfkQCpqKKKwrjw+llPJfaNO1hMcvLEib4Zv96Md+Oq4P1rz7w58ZLn/hJZLDxLapbW08pEU20p5HYBgRyvueRmvXA6NGJFZShG4MDwR65rNfxBasxSxSa/kBwRbJuUH3f7o/OsKTwzf3fjC08RQw2mm3EcbRTAsZTOpGBuAAGR1Bz2xXQNaanJ11KOP/AK52w/8AZiai/srUf+g9d/8AfiH/AOIp6WOpR4/4mzSY/wCetunP124okGtRYMTWFx6hg8R/DlqZ/bjW3/IUsbi0H/PRR5sf/fS9PxArRt7mG7hEttLHLG3R0YMD+IqWiiiquoahFp0AkkDO7sEjiQZeRj0VR6//AK6pf2feamxbVJfKtyOLOFuD/vv1b6DA+tacFvDaxCK3iSKNeiooAH4CpKKKKKKzr3RYbpnmt3ks7xgQLiDAb8R0b8Qa5XT/AAvoPhm6V9a0+Ca5uHx/ac+ZRK5P8W7Plsfy967lEWNQqKFUdABgCnUUUUUUUUUUUUVieJLfSbuCO21LT4dQmnOyG3ZAXc9Tgn7oHUntVXRfCK2Fqbe5nlayDZh09ZWaGEf3cn5mHsTj0FdFHGkMaxxIqIowqqMAD2FPooooorMuNFjMjT2EjWN0TkyRAbXP+2nRv5+4otdTljuo7LU40huZAfLdGzHNjrtzyD32n9a06KjnnjtreSeZgscalmY9gOTWbplm9zcf2tfIRcyAiCNv+XeM9AB/ePVj+HQVrUUUUUUUUUyWKOeJ4pkV43BVlYZDA9QRWXpqyaXeHTJGd7YrvtJGOSFHWMn24x7fSteiiiiiiiiiiiq99eJYWctzKCVjH3V6segA9ycCqmlac8LPfXuH1C4A8xs5Ea9o19FH6nmtOiiiiiiiiquoafBqdo1vcqSpIZWU4ZGHRlPYj1qDSbqZxLZ3pBu7UhWYDAkU/df8QOfcGtGsrxD89hDbn7lxcwxP/ulwSPxAx+Natc5qsKTeNtHR95RrW5LKHYAkGPGQDjjJ/Os7RdTk0rS4hHEZjdaxPaDzJW+T94+05OeAFxite28SC8ulsoIcXe6ZXDt8g8plDYbHOdwx+OelNtvEU9/qAgsrZHRrOO8Us5UlWYrt6dflPtTYfE7ESSz2oSCPUDYEq+WLbgoYDHTJ6dR1pY/FAlvJ7WOJDILeS4hcMdkiq+372Oc5B4yPypkHimS5uNMtkswk928sUwd/9Q8YJKnA56H9PWoU8YyOltGbEpdSXotJUZ/ljDZKvnHIYYI/H0q7Z+JYrnWlsGRcTCVoJUYsriMgHnAGfm6DOMHNblY3iS8OnQ2N2kMkzJdovlxDLMGypAHrzXDeJvi5qekq62vhW+gIz+9v0KqPfA/xrs9N8RvPpdrNNY38skkKOzxW3ysSAcjnpVr+3R/0DNT/APAf/wCvR/bo/wCgZqf/AID/AP16P7dH/QM1P/wH/wDr0f26P+gZqf8A4D//AF6P7dH/AEDNT/8AAf8A+vR/bo/6Bmp/+A//ANej+3R/0DNT/wDAf/69H9uj/oGan/4D/wD16P7dH/QM1P8A8B//AK9H9uj/AKBmp/8AgP8A/XqrdaoL26sLb7FeQh7hSWmh2r8oLY/St+iuW1nVJrTxDY3SG4+xwTC0uFCHyz5g++T0+Vtg9smruu6/LpENxNFFDPHAqsy7yGBLBcdOPvA02LX7gy6kZbaMQabLsnZGJYjyw5YDHPDdKr33i/7HaCX7MskjQRXKIH4McjhBzjggsPrzT73xS1nOYGgTzY7yK1lJc7R5gyrA4/Mf402bxPdw6FfakbFNtm0odDIQSEwPTqc/hirN5r0+mXFst9bItvK+ySdHyI8nCMcjoSQD6E1PoesPrAvma38lba5a3ALZLYAOT6Z3DitWuc1i7ex8aaAsY/4//Ot5f91UMg/UfrXR1Q1q1kvNKmSAAzoVliB7ujBlH5gCrFleRX9nFcwHMci7h6j1B9weKSSwtZbyO7kt42uYgVSUr8yg9QDWbfeF7O6jtIoFS1ht7sXbRxxjEjjPX65NaB02zMEcP2WHyozlE2DC/SlTT7SO8N2ltEtwUEfmBQG2joufT2qjpfh22024uZ2K3Es1y9yrvGAYy/UA1M/h/SZH3vp1qWCsoPlDgN1H0NSHRtOZo2Nlb7oiSh8sZUkYJH4cU3+w9M2RJ9gt9sLBox5Y+UjoR9MmlXRdOS6juVsbdZ42LJIIxlSepB7Zq9WZeOl1rNnaL8xgzcyD+6MFVz7knj/dNaLIrqVdQynggjINCIsaBEUKqjAAGAKdRRRRRRRRRWdrYZLD7QiF2tnWbaOpCnn9M1eilSeJJYmDxuoZWU5BB6Gn1SbR9PayktGs4DbSPveIoNrNnOSPXPNEuj6dPHLHLZW7pMVaVWjGHI4BPrjAqS20+0s5JntreKJ5jukZFALnGMn1qjq3h221LSprGLbapMyF2jjGSFYMB9MirLaNp8kUcctlbuqSCUAxjAf+99felbRtOa1ltmsoDBMxeSMoNrsepI706TSrGUoZLSF9ieWu5AcLnOPpkCn2tja2RlNrBHCZnMkmxcb2Pc+pqxWNFEmp+JHuuGi09TDGf+mp+/8AkMD8T6Vs0VjXUM2kXj39orSWkhzdWyrkg/8APRAO/qO/Xr107W7gvrdJ7WVJYnGVdDkGpqKKKKKKKyr/AFabzzZaTCt1eD/WEtiO3B6M5/8AZRyfarOm6cNPgYNI09xId807j5pG9fYdgOwq5RRRRRRRRRRRRWHIJfDtw8yAyaRIS0iAZNoe7KB1Q9x26jjONmKWOeJZYXWSNhlWU5BHsafRRRRRRRWLeX02qTvp+kzbdp23N2vIhH91T0Ln/wAd6nsK1LW1israO3t0CxxjAH9T6n3qaiis640dGla4spXs7ljkvH91z/tp0b+fvUI1HUbIY1GwMqj/AJbWXzjHuh+Yfhuqe213TLtgsN7AXPGxm2t/3ycGtCiiobi7t7RN9zPFCn96Rwo/WqB8QW8xK6fFPfuOP3CfID7ucL+tILTUtRJN/MLSDPEFq53MP9qTAP4Lj61oWtpBZQiK2iSKMc4UdT6n1NTUUUUUUUUUUUUUUVlyaMIXM2lzNZy9SijMLn/aTp+IwaaNVvLRcanp8gA6zWv72P64+8PyP1qza6zp16QttewSOeiBwG/LrV2iiq9zfWtkAbq5hhB6eY4XP51R/t+O4yNMtbi+PTci7Iwfd2wPyyaF0+91DDarOI485FtasQp/334LfQYH1rSggitoVigjWONBhUQYAqSiiiiioLmxtbwYuraGYD/nogb+dUB4Z0tRiK3kgHpBPJEPyVhVc+GJMnbrWpqueB52cD0qWHwzAFP2u7v7ps8F7qRcD0wpAP41Yg8P6VbMXj0+33k5LsgZvzPNaAAAAAwB2FLRRRRRRRRRRRRRRRRRRVW602yvf+Pq0gm4xmSMN/Oqf/CM6aq4hjmg9BDcyRgfgGAqt/wjEv8A0G9T/wC/tTxeGrUJi6nvbp+7SXUgBHoVDAfpVm20LTLQ7oLC3Vs53lAWz9TzV+iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiv/Z)
Создание модели проходит в 3 этапа: модель – алгоритм – программа.
На первом этапе строится модель, наиболее полно отображающая свойства объекта. Модель исследуется теоретическими методами, что позволяет получить важные предварительные знания об объекте. Второй этап включает в себя разработку алгоритма, для реализации модели на компьютере. Модель представляется в форме, удобной для применения численных методов, определяется последовательность вычислительных и логических операций, которые необходимо провести для нахождения искомых величин с заданной точностью. На третьем этапе создаются программы, переводящие модель и алгоритм на доступный компьютеру язык. К ним предъявляются требования экономичности и адаптивности к особенностям решаемых задач и используемых компьютеров. Их можно назвать электронным эквивалентом изучаемого объекта, уже пригодным для непосредственного испытания на компьютере.
Целью данной курсовой работы является изучение приёмов численного и символьного интегрирования на базе математического пакета прикладных программ, а также реализация математической модели, основанной на методе интегрирования.
1. ИНТЕГРИРОВАНИЕ ФУНКЦИЙ В MATLAB
Возможны два различных подхода к определению определённого интеграла.
ОПРЕДЕЛЕНИЕ 1: приращение F(b)-F(a) любой из преобразованных функций F(x)+c при изменении аргумента от x=a до x=b называют определённым интегралом от a до b функции F и обозначается
.Причём функция F является первообразной для функции f на некотором промежутке D, а числа а и b принадлежат этому промежутку. Это можно записать следующим образом:
, это формула Ньютона-Лейбница.
ОПРЕДЕЛЕНИЕ 2: Если при любой последовательности разбиений отрезка [a;b] таких, что δ=maxΔxi→0 (n→∞) и при любом выборе точек
интегральная сумма σk=
f(εi) Δxi стремится к одному и тому же конечному пределу А, то это число А и есть определённый интеграл, т.е Δxi=A(2). Где Δхi=xi-xi-1 (i=1,2,…,n) ε=maxΔxi – начало разбиения
произвольная точка из отрезка [xi-1;xi]
сумма всех произведений f(εi)Δxi, (i=1,…,n). Простыми словами, определенный интеграл есть предел интегральной суммы, число членов которой неограниченно возрастает, а каждое слагаемое стремится к нулю. ГЕОМЕТРИЧЕСКИЙ СМЫСЛ:
Всякая непрерывная на отрезке [a,b] функция f интегрируема на отрезке [a,b], функция f неотрицательна, но определённый интеграл
численно равен S криволинейной трапеции, ограниченной графиком функции f, осью абсцисс и прямыми x=a и x=b,
.Рассмотрим основные методы интегрирования: метод трапеций, метод прямоугольников и метод Симпсона.
Формула прямоугольников
Теперь рассмотрим первый вид приближённого вычисления:
требуется вычислить определённый интеграл:
.Пусть на отрезке [a,b] задана непрерывная функция y=f(x). Разделим отрезок [a,b], аналогично как в формуле трапеций: точками a=x0, x1, x2,…, xn=b на n равных частей длины Δх, где Δх=(b-a)/n.
Обозначим через y0, y1 ,y2,…, yn-1, yn значение функции f(x) в точках x0, x1, x2…, xn, то есть, если записать в наглядной формуле:
Y0=f(x0), y1=f(x1), y2=f(x2)…yn,=f(xn).
В данном способе подынтегральную функцию заменяем функцией, которая имеет ступенчатый вид.
Составим суммы: y0Δx+ y1Δx1+ y2Δx2…+yn-1Δx; Y1Δx+ y2Δx+…+ynΔx.
В результате вычислений получаем конечную формулу прямоугольников:
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAQgAAAAwCAYAAADtqOSLAAAAAXNSR0ICQMB9xQAAAAlwSFlzAAAOxAAADsQBlSsOGwAAABl0RVh0U29mdHdhcmUATWljcm9zb2Z0IE9mZmljZX/tNXEAAAauSURBVHja7V09b9w2GGaXLk0Bj00RekunLj7eWk/BRQ6yejgfjA4NMgSBOmTwYKSyNwONf0AzNZsz9A80Gbo2/8BA8gfajC1QoIDboyTqg0eKpETqTrpneIBEvpOo9+X78P0ij5yfnxMAAAAVIAQAAEAQAACAIAAAAEEAAACCGCCOGHlFCPmPEPbhOEluQyYACAKo4emMPiN3orePLy5uQR4ACAIocHHx+FZ0h7yls6fPIA8ABAHUkCTHtxkhH9ji7CHkAYAggBrOFuwhzz/Mk/gu9yQIITc8JxGaMPLcxw1/9iKOHlBCPsKLAUAQG4Y0/0DpNSX0OoqTPRFyhEpaCo+FsKNX4v8Rpa/7ICUABAG0yD/IZJB5FfQjJ4yQ5FD3JlBFAQISRDapebnO/8Qebf4hjva4ay+v3KEIQkUEBUlJpAEA3giiWJlSgvDvqnJDYrP4kb0h0MshkJTIP8grd2rIirJnmTfI5KyExtCFjuQ8g+56Vx0A48WC0Rc6+9JO9FAJrnSVpdFrlx6BdFVcxtWbThKpwUsGbWuwvsgozYEYvJU2OgDGHRofUHp178npxIog+CQLQRCpsezef9lmYvLvTsjkzabG1br+B5330BUZEdQJQoQ4TfmHLjoARhwepwsZ+/X45ORLI0HwSd2WILKJm7nHKmPpEq7YEFeZO9G75+HyD1nlYvnsT3Q5Ar8EkXkK/Hnpey+JKJqSX/h7J/NoX+VFdNUBMF6o7MvaVbZ2sfPVUm43zjwA+q5LmJAZ4eSdyeB0ycI+BFzNKYTsQygrJvnzlvqqXVPoz4cOgBF7EcK+Kl6E0dC7uLwrf+/oatu2Mevic6wQ2CMCuNmXF4Iwrdg+S282YwsV949B+Sh/Alb2dXLyuZIgiolkYWByOVTnVpsy+TLBFMk2bWmw7h2UW6x5TL54UH2Wy71H7T5aVFNU/Rq+yLZNZSzEeEJW6ELKoFphCLkAFvaVhxkNsa29i25yXYVRqoQjJxVFsk13r2pyTtVmLAiLE4LrvccfX+r3Z5QkWyd5X4lWV+PQjWfByM9dxjMkgsjeNcsx+ZSBOVVQElJngrBxXU2TsxyYufJQEkTmKdi0GdveuzopXRuXhk4Qus9w/c7Y7Pk6PIgQ4xmaB5HL4M8QOumFIGwqBmIlt5mcpsqDeIHJ1+Q32zbjdVU1Ngk2OlDJj8sumif7voxDlIDb5k26jqfNOHwTxLpl0CtB2OwzsFm9bCsU2Qvs/LVDyN/yZ3VGgENc3HQgJmOXlcqqjdwwx3yMp+qmtx1HW/h4tk+drIUgbGJUK/d2Hu1TSq5NOYLsBT79VyalctyrZGV777YhhnkSrA8uOpATYLpmqz5de9/jGVqIEVInwQnCtnSm25Zcu0+0+M6mfyF7gc/+Wa1k0MvplF4pww7Le4/agzDoQJ6McTK/y9gi2RTj8DWeIRPEk9PDr3zqRG9fnqoYtrG9jkjSwSyvxTP2KK1KVNqVdTs4U2Ht7PxR3WXK78P/LTpAF4wlBwf0heu9xwxbMhcTZDqd/qTSv+gWbVOua2McuvFkG/jIG1fSdxmHKmS1vdaHDJz1byEvuYSqn0gW5UBT9+TKZ6V7iuqCEFqlVVi5xbzMJXz7XHxOfF/+ruu9twE2nZRCbioZpSTLji6Fx+ZCst2Mg9zo5sP9XfYyMEHcKMjAeK0vGbiSRJO8vHdSuuzZ6HMvBtBeB00xrvDU+nTVm8bThiAGqbcGGWQ2wX4/PD2cLEPI96YkdJO8VPbVmiDaGGvXbeQpU6NVuLMX0djRmnsIKohwTXw29IEzpvFsA0GYZHAWx1/MZrPvI8Z+nCfzPZ6j0JVRTfJS2ZeTKyqf2Ozq8nQ502FbVot+vIhSB6K7NI7vfWM6RKYPD8JlPGOdE3UZHFw1LtT5IcWiczgNv9NcQ6Wylttyk7zS8GJ3NcehJwjVduHiQJL2cbyJEW1WL8DfqlRsUbfZe5N/j69QoRK9LuMpJrx0yMkockWWMqgSNU/ON8miSV46+/KaTAkVKnQ5wAbwG66JBqhN0EfRiLTFO3eXpPADX/XzMvb7Jlno5MWv6/RZc1VEnTx0MwYAAMNAPQYRpUMkAntDTe4VZi+6OnGuBbAJBAFsQNxZOTau3K6O3yUBQBBbjeoP/laTfy6NaAAAghgrQeRt4HwvSbU6hKPzABAEoDzQJtQP7gAACGJgUG2ZR/4BAEEAWk8BR9QDIAig9otc4lp1azY6SAEQxBZD+RuboituQzoWARAEAAAACAIAABAEAAAgCAAAQuJ/TMSQzmm9RkcAAAAASUVORK5CYII=)
Формула трапеций
Возьмём определённый интеграл
, где
— непрерывная подынтегральная функция, которую мы для наглядности будем предполагать положительной. При вычислении интеграла с помощью формулы трапеций подынтегральная функция f заменяется функцией, график которой представляет собой ломанную линию звенья которой соединяют концы ординат yi-1 и yi (i=1,2,…,n).
Тогда площадь криволинейной трапеции, ограниченной линиями x=a, x=b, y=0, y=f(x), а значит (следуя из геометрического смысла), и значение нужного нам интеграла, приблизительно равна сумме площадей обычных трапеций с основаниями yi-1 и yi и высотой h=(b-a)/n, так как (если более привычно выражать для нас) h это Δx,a Δx=(b-a)/n при делении отрезка на n равных отрезков при помощи точек x0=a<x1<…<xn=b. Прямые x=xk разбивают криволинейную трапецию на n полосок. Принимая каждую из этих полосок за обыкновенную трапецию, получаем, что площадь криволинейной трапеции приблизительно равна сумме обыкновенных трапеций.
![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAASwAAADMCAMAAAD+kG5xAAAAAXNSR0ICQMB9xQAAAANQTFRFAAAAp3o92gAAAAF0Uk5TAEDm2GYAAAAJcEhZcwAADsQAAA7EAZUrDhsAAAAZdEVYdFNvZnR3YXJlAE1pY3Jvc29mdCBPZmZpY2V/7TVxAAAAUUlEQVR42u3BAQEAAACCIP+vbkhAAQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADwbu/cAAEWlzZdAAAAAElFTkSuQmCC)
Площадь крайней полоски слева равна произведению полусуммы основания на высоту
Итак, запишем сказанное выше в математическом виде:
– это и есть формула трапеций.Формула Симпсона (формула парабол).
Разделим отрезок [a;b] на чётное число равных частей n=2m. Площадь криволинейной трапеции, соответствующей первым двум отрезкам [x0,x1], [x1,x2] и ограниченной заданной кривой y=f(x), заменим площадью криволинейной трапеции, которая ограничена параболой второй степени, проходящей через три точки M0[x0,y0], M1[x1,y1], M2[x2,y2] и имеющей ось, параллельную оси Oy (рис). Такую криволинейную трапецию будем называть параболической трапецией.
Уравнение па
раболы с осью, параллельной оси Oy, имеет вид:
. Коэффициенты A, B и C однозначно определяются из условия, что парабола проходит через три заданные точки. Аналогичные параболы строятся и для других пар отрезков. Сумма параболических трапеций и даст приближённое значение интеграла. Сначала вычислим площадь одной параболической трапеции. И продолжая вычисления, получаем формулу Симпсона: ![](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAe8AAAApCAYAAAAcT6HQAAAAAXNSR0ICQMB9xQAAAAlwSFlzAAAOxAAADsQBlSsOGwAAABl0RVh0U29mdHdhcmUATWljcm9zb2Z0IE9mZmljZX/tNXEAAAvKSURBVHja7V09k9vGGYZc26kj5XCdkzbHPbWqPDwwltOdYxyGTaS5gqNhYXmGhSxB6jTjU2+psZxKzh/I2Kd0yYylP+BkdP4Dkf9CRGMBLLhY7iexAAHeUzwFcbxd7Lsv99n3Y98NHj16FAAAAAAAMBxACAAAAAAA8gYAAAAAAOQ9RMEGwRXIAQAAAAB5DwCPH5++H4XBeUbey3B85y5kAgAAAIC8B4ITEj6J5ukBZAEAAABslbwfJuRmsBe9PH38+H0IT400nV4dBaPzaZpehTwAAACArZI3JaXxeP4p4rkWmxxy8i2VU1uyarNtAAAAYMDkLRJEGkc3Dg/DF9mzdyR5eBMClOOEBN+S6fQ2CYILKivf3oo74/Bu3m4Qvh2PyX14QwAAAEDeK+sxCJYZQfzCYrc0jns9mX+WzqODkJw8gQDl3omMtH8OAnJB3ea5rILgF9lmh22OZJC1nSfD7QUvqVXP94XEOAAAAJA3R0DUuiusbOYyZ8QOwtBsekpy5WUpkje1zpl8peDaqP0P95yRObwgAAAAIG8pOadxfDBdLK4tFtNrSMZSI3eZc2RaWN7hf5tmnheeEPIzL/ei7dFrzAUAAADIO4+prpF34f59S+OsOAL1SOmxEDc2orXcZFPAt7NYnH7Au9ABAACAS07elCh8u8W3mR3dVb9ilnmRWFbEvn2RN223dMXTZLglTYyLY2ymfOknMvgBABg0efu06LgY+lrVMZcz0Zucn64S7zrKyK7Fsj3KsBxHGRMnF3EaHxQErt8c4My5tX5elPJdkjj9RP098gNkedl0I5vzxeJaq5v+MDo/XSw+2Eb/wNb16/vjxeJ3umdu5O2J7AryCF9TV7tY5EUWxzW2t0EMWRYGuAwo5XsBslGDJv0d7ZNnTEblMbylmAQIWV7iBTZbc/aD/f98NLs3am3jSNdFgbwrnQNx7/Y6fTf6Y6ZfP2X6RZjnT/asc/JWueCbWIQuiVpFrfHwu8sWp+/K4hbJr8/ELJVTHN3gdUM8jte192Kb8uyi7z7ri3nNIa9MRJrrz/7hU1vCTQh5MJmEZyJ5d2Fxu77r0Prd1vg2JfBMv37krW3ZMy15V4uXB/LWWda6uLos9ih+zi0kjVuafT//EewfPaNjsWl3V2CS77YW4y7lvSlRiGEjUw7Irsizj+Tdp9/nPAo/Z3knPgiDro/Uw5OvZQJ5JyR4TnVOaXV5kIsrufmai76S97Z1bZbp13u5fv39Pd0zM3k7urPlbdAiLxR1dyPvSpcTfl49rCoOUz3jNhQq65tOQBV3zr4/T65/Rn9wtu2KGwjbAir9s7rl8i1ls1wLX2y4WXNZjOmCJMq7zXr5m5CReHZeJ8s2xrRNefaNvGvjK4lNFx/ug/VtSxi5UUGSB3Q9yTcF3BgLqzt8pXLTS+UyJR+7ysWF3OJVn/9s0mdfyTtbF7+5Uo7v+PT0N/z42OeOrO9XEuv7lcz6boW8VRaMjUUuWtYqC1tW/KRKjGMVyCTVzUztcgTnVEClVy4Yg3zFv9M5H5PxfZcFn21k6JG1fDHOfiA2mxsffbf1bpXOhNF3LnkZvZBnuUB12XcXcyLGffP1iURfbqMcMItPyxIa+fExwtCNj7rLmb6I5G0T686JRSYXCyJ1fVcffcr6pYTUhkHE2mTkbdsPGx8jStfxedKvNwfx/T/XPMfCs9bJW1VZrCJQQ/uMeFXuSrE0qCxWKdvlm9r1oTTbhI18xaQ/vnqeDaqd/wpL281N077bfLdi41a/ytVeV4cpz6by8tWHamEVk8V868sm5J2tHV/w7+s6PjqmKE5v1I6UZuTNLLyczA3kLZNLFM2PTQQVG981dpoLmz6ZZbtpv44G4zcexkdcx+eTvHn9kj1rn7w1SWU2C6Kp9KdI3jILSSxVatNuE7e51lrvCFbkLWys5uPxp7LwQxuuVlPfPnfiru9G5SbqhRV5exyTb3kO2W0us3bb1BdTWxx539XeP2Bw1SrIPicZOs6qPoSOvDm55K53iVx8u6/X+oyiY9anq473MebNiLLp+Gz0TE/eMwl5z7oj7yoJQ+LesnFFzsfkdn6DmWLnLy5aootedSGIqd3L4DbnZUflRMbz2+uEFbyz8U40IRuxb/68tQ/PiMu7sZLA7EfKrCMbXfU5pibkJptLl0uE+k7eBtl+0WRxtWlL5zZvQkxrbnPBPW0jl+vR/BZzFU/C4Ac+lNIGefN9lmvt/2znYQjkXRtf+bypnuW6n83NlWxuZDFscfOgetYqebM2dC5vVRIQfUG2kFbxx0ypKeny3xct+1oFsvz/Dv9xGB6+SObRn5iLyqbdXYApyaqan8nkjGXiyxZ70YXsYzHW9c2sB/FoVqMfi8W7rW7Q41Hov0mWm46pC3lWhGS52ewjeRdHPYPzanxc/JFZRuKxqk30RdXWujfRT8KajrxNCWsmudi+g+u7Vn1G0RO+z4fz+W+p3tH3PtyPnppixH0lb9X4BN34nm/PVdd071TqV+1omOxZq+RtU0hFdvymuqe6XGy4rPV30u+uW9pvqyzzNP6QFT9IojC1bXdXoDvexM3xWtER3nVsc3tcpYyOZGMqeMI2DszypZswulk7vjf7ve1NajbvtrobXe1dYcd2TD962zF5l6ei74SQNJebg+Wt67vqq/SoaT8rFvBG49O6osMzTl8u4jvjj/MsXaov2f/rrGVVW1JdIclzX0fFauQtZG2bjorp5NI2eavmgnmrfJBq5UHYO8p1yfRZ08/XLuMrvRYXqrPVVDfYpqr0kLz5y2x8s8gIn/2B/r/MSrYZe3EUMX4uHhUTn5nJu8HxHdMZbBvr0Py/KPnZRL5JkkykCwnniZC5YX1A1XftO1F0i+kf3XxRT8kR2X+2t7f37yiJbrXxXnp3mt4Sch1T2/JkZ4jvzT4aUfIecj2DDfTlQRhOXqj0xZS/krX1V3Geui5T2kTn2rRsVX1WiXCjk78NWdem02mke/9cN/hNFqdrYRj+axyPb1NdY22IuqaaG5lVbyqTqrTaXMmbEbZL9bNNSp3qstiHhC7OjLPsSSrfKlSQnl7VWa6ulre1J8aib5F42HxPwslZnM4/ZP8rSyxrfTEVM235MVHLzmFMXuUp6XtlIXEJUQO6fa6JbEt9+SrXl/J/mb5I5PJO9BLI5skm1t2dzk1emOTik7yrPhen15gHo0myXj91rRifyWpmusHHo6mu8d4dpmv5kcjciud17ehlPneCjFxj3d7Je1UExa3+c9sXk/SRtFeu2vZrZTOLoXA7m690ZTHvotiN2c1r/QOxvE6Wv5mN30Cw96pipltYIPiLSZqMqUt5uiSs9WpBLcdnqieu15fkjD9XbrWGice6tnwxCN//KjFMLxfv5G3RJzNGRMu07yiKoJjHV+gGdV8Hlfuaus3DaP554d0qdI2dKzfNDfvOJha3lrxlmeK4KtHXD5HertbvCy5Ytn3X+QBi/JnuYGnclsqK7Xpdk7C2jeLMbn1MnRNhtrDs4u9pTbaZlcKKn8j0RbeGydoaqlxiSQW2Tvor5bZrXDErdOP/4hhjqmsZ6VIS53TtTXGWfHlFsbaWldyOXjat3KZeRCXJYDSzEbHmhsQ9IPclAAAA0E+suT6Yi4BZXbyr2iYRDdBZs8XZS1W8W3xu+gwAAABccvLmLUP+mkQ++9hXAtNltbrzM5llIQixiAJXJGLJin0Ucb/yvLHwGXIFAAAAedfOPtfPuJJ0lSzhlhkOFGBFQGgGLSVe1Vn66nhXEt2i36Vn1anXI0824z6DvAEAAEDeWlSlQsPR61EYvG7r+sbdJ29J3fXSyq49y+RMSPQVlXH+eTQ6Z5Wz2rw+EwAAANgh8gbaIW/ZefVio7S6b7yKk3MlYId+vh0AAAAAeQ8C0rvHhTCEeNuZ+D+oKgcAAACAvDtGdbVkmW0uvQWNVpUqXeKii5wVkNjFi1QAAAAAkHdvUbtqVDhyR8mdz+QXXeTbKpwCAAAAgLwBAAAAAAB5AwAAAADIGwAAAACAHuNXOVuT9sWt0cYAAAAASUVORK5CYII=)
Теперь рассмотрим методы решения интегралов с помощью программы Matlab.
1.1 Численный метод
Вычисление определенных интегралов.
Рассмотрим пример:
.