Cyfroteka.pl

klikaj i czytaj online

Cyfro
Czytomierz
00365 008448 11062970 na godz. na dobę w sumie
Excel w nauce i technice. Receptury - książka
Excel w nauce i technice. Receptury - książka
Autor: Liczba stron: 432
Wydawca: Helion Język publikacji: polski
ISBN: 83-246-0477-4 Data wydania:
Lektor:
Kategoria: ebooki >> komputery i informatyka >> podręczniki szkolne >> książki okołoszkolne
Porównaj ceny (książka, ebook, audiobook).

Excel to nie tylko arkusz kalkulacyjny
-- to potężne narzędzie obliczeniowe

Większość użytkowników komputerów kojarzy Excel z programem wykorzystywanym w biurach i urzędach do tworzenia zestawień, tabelek, wykresów i przeprowadzania obliczeń matematycznych. Tymczasem jego możliwości są o wiele większe. Excel, dzięki wbudowanym funkcjom, dodatkowym pakietom i językowi VBA może służyć jako narzędzie do wykonywania złożonych operacji obliczeniowych, przydatnych naukowcom i inżynierom. Dziś do przeprowadzenia zaawansowanych obliczeń i symulacji nie trzeba już wyspecjalizowanych stacji roboczych i trudnych w obsłudze aplikacji -- wystarczy komputer i Excel.

Książka 'Excel w nauce i technice. Receptury' to zbiór sposobów, które ułatwią odkrycie wszystkich możliwości obliczeniowych tego programu. Opisano w niej metody rozwiązywania różnych, nawet bardzo skomplikowanych zadań matematycznych, związanych z porządkowaniem danych i tworzeniem wykresów. Przedstawiono sposoby przeprowadzania analiz statystycznych i analiz ciągów czasowych, aproksymowania wykresów funkcji, obliczeń macierzowych, rozwiązywania równań liniowych, nieliniowych i różniczkowych oraz numerycznego różniczkowania i całkowania.

Dzięki tej książce przekonasz się, że Excel jest narzędziem przydatnym w każdym biurze projektowym i pracowni naukowej.

Znajdź podobne książki Ostatnio czytane w tej kategorii

Darmowy fragment publikacji:

IDZ DO IDZ DO PRZYK£ADOWY ROZDZIA£ PRZYK£ADOWY ROZDZIA£ SPIS TREœCI SPIS TREœCI KATALOG KSI¥¯EK KATALOG KSI¥¯EK KATALOG ONLINE KATALOG ONLINE ZAMÓW DRUKOWANY KATALOG ZAMÓW DRUKOWANY KATALOG TWÓJ KOSZYK TWÓJ KOSZYK DODAJ DO KOSZYKA DODAJ DO KOSZYKA CENNIK I INFORMACJE CENNIK I INFORMACJE ZAMÓW INFORMACJE ZAMÓW INFORMACJE O NOWOœCIACH O NOWOœCIACH ZAMÓW CENNIK ZAMÓW CENNIK CZYTELNIA CZYTELNIA FRAGMENTY KSI¥¯EK ONLINE FRAGMENTY KSI¥¯EK ONLINE Wydawnictwo Helion ul. Koœciuszki 1c 44-100 Gliwice tel. 032 230 98 63 e-mail: helion@helion.pl Excel w nauce i technice. Receptury Autor: David Bourg T³umaczenie: Zbigniew Waœko ISBN: 83-246-0477-4 Tytu³ orygina³u: Excel Scientific and Engineering Cookbook Format: B5, stron: 432 Excel to nie tylko arkusz kalkulacyjny — to potê¿ne narzêdzie obliczeniowe Wiêkszoœæ u¿ytkowników komputerów kojarzy Excel z programem wykorzystywanym w biurach i urzêdach do tworzenia zestawieñ, tabelek, wykresów i przeprowadzania obliczeñ matematycznych. Tymczasem jego mo¿liwoœci s¹ o wiele wiêksze. Excel, dziêki wbudowanym funkcjom, dodatkowym pakietom i jêzykowi VBA mo¿e s³u¿yæ jako narzêdzie do wykonywania z³o¿onych operacji obliczeniowych, przydatnych naukowcom i in¿ynierom. Dziœ do przeprowadzenia zaawansowanych obliczeñ i symulacji nie trzeba ju¿ wyspecjalizowanych stacji roboczych i trudnych w obs³udze aplikacji — wystarczy komputer i Excel. Ksi¹¿ka „Excel w nauce i technice. Receptury” to zbiór sposobów, które u³atwi¹ odkrycie wszystkich mo¿liwoœci obliczeniowych tego programu. Opisano w niej metody rozwi¹zywania ró¿nych, nawet bardzo skomplikowanych zadañ matematycznych, zwi¹zanych z porz¹dkowaniem danych i tworzeniem wykresów. Przedstawiono sposoby przeprowadzania analiz statystycznych i analiz ci¹gów czasowych, aproksymowania wykresów funkcji, obliczeñ macierzowych, rozwi¹zywania równañ liniowych, nieliniowych i ró¿niczkowych oraz numerycznego ró¿niczkowania i ca³kowania. (cid:129) Podstawy obs³ugi Excela (cid:129) Programowanie w jêzyku VBA (cid:129) Importowanie danych do arkuszy (cid:129) Sortowanie i filtrowanie (cid:129) Obliczanie przedzia³ów ufnoœci (cid:129) Analiza i predykcja ci¹gów czasowych (cid:129) Regresja liniowa (cid:129) Aproksymacja funkcji (cid:129) Rozwi¹zywanie równañ i uk³adów równañ (cid:129) Ca³kowanie i ró¿niczkowanie (cid:129) Optymalizacja (cid:129) Obliczenia finansowe Dziêki tej ksi¹¿ce przekonasz siê, ¿e Excel jest narzêdziem przydatnym w ka¿dym biurze projektowym i pracowni naukowej. Wstęp .............................................................................................................................. 9 1. Podstawy pracy z programem Excel .............................................................................15 15 15 21 23 26 29 32 34 35 35 36 39 42 44 46 47 48 50 53 1.0. Wprowadzenie 1.1. Poznawanie interfejsu 1.2. Wprowadzanie danych 1.3. Ustawianie typu danych dla komórki 1.4. Zaznaczanie wielu komórek 1.5. Wprowadzanie formuł 1.6. Styl odwołań W1K1 1.7. Odwoływanie się do więcej niż jednej komórki 1.8. Zrozumienie priorytetu operatorów 1.9. Stosowanie potęg w formułach 1.10. Funkcje wbudowane 1.11. Formatowanie arkuszy 1.12. Definiowanie własnych stylów formatowania 1.13. Korzystanie z poleceń Kopiuj, Wytnij, Wklej i Wklej specjalnie 1.14. Używanie nazw komórek (jak zmiennych w programowaniu) 1.15. Kontrolowanie poprawności danych 1.16. Stosowanie makr 1.17. Wstawianie komentarzy i równań 1.18. Uzyskiwanie pomocy 2. Poznawanie języka Visual Basic for Applications (VBA) ........................................... 55 55 56 59 63 64 2.0. Wprowadzenie 2.1. Nawigowanie po edytorze VBA 2.2. Pisanie funkcji i podprogramów 2.3. Typy danych 2.4. Definiowanie zmiennych 3 2.5. Definiowanie stałych 2.6. Używanie tablic 2.7. Komentowanie kodu 2.8. Wpisywanie długich instrukcji w kilku liniach 2.9. Używanie instrukcji warunkowych 2.10. Wykorzystanie pętli 2.11. Uruchamianie programów VBA 2.12. Poznawanie funkcji wbudowanych 2.13. Poznawanie obiektów Excela 2.14. Tworzenie własnych obiektów w VBA 2.15. Korzystanie z pomocy VBA 65 66 67 68 69 70 72 75 76 81 84 3. Gromadzenie i porządkowanie danych ......................................................................85 85 85 90 91 92 94 97 97 100 102 105 109 116 3.0. Wprowadzenie 3.1. Importowanie danych z plików tekstowych 3.2. Importowanie danych z plików tekstowych delimitowanych 3.3. Importowanie danych metodą „przeciągnij i upuść” 3.4. Importowanie danych z baz danych Accessa 3.5. Importowanie danych ze stron internetowych 3.6. Konwersja tekstu na kolumny 3.7. Usuwanie dziwnych znaków z zaimportowanego tekstu 3.8. Zamiana jednostek 3.9. Sortowanie danych 3.10. Filtrowanie danych 3.11. Wyszukiwanie wartości w tabelach 3.12. Pobieranie danych z plików XML 4. Tworzenie wykresów .................................................................................................. 119 119 119 126 128 129 132 134 138 140 140 145 148 151 4.0. Wprowadzenie 4.1. Tworzenie prostych wykresów 4.2. Typy wykresów — krótki przegląd 4.3. Formatowanie wykresów 4.4. Modyfikowanie osi wykresu 4.5. Ustawianie skali logarytmicznej lub półlogarytmicznej 4.6. Tworzenie wykresów o większej liczbie osi 4.7. Zmienianie typu wykresu 4.8. Łączenie wykresów różnych typów 4.9. Tworzenie wykresów typu Powierzchniowy 3-W 4.10. Tworzenie wykresów konturowych 4.11. Opisywanie wykresów 4.12. Zapisywanie własnych typów wykresów 4 | Spis treści 4.13. Kopiowanie wykresów do Worda 4.14. Wyświetlanie słupków błędów 151 152 5. Analiza statystyczna .................................................................................................. 153 153 154 158 161 162 166 168 172 174 175 5.0. Wprowadzenie 5.1. Obliczanie statystyk podsumowujących 5.2. Wykreślanie rozkładu częstości 5.3. Obliczanie przedziałów ufności 5.4. Korelowanie danych 5.5. Wyznaczanie rang i percentyli 5.6. Przeprowadzanie testów statystycznych 5.7. Przeprowadzanie analizy ANOVA 5.8. Generowanie liczb losowych 5.9. Pobieranie próbek 6. Analiza szeregów czasowych .....................................................................................177 177 177 178 180 186 191 194 197 200 202 204 6.0. Wprowadzenie 6.1. Wykreślanie szeregów czasowych 6.2. Dodawanie linii trendu 6.3. Obliczanie średnich ruchomych 6.4. Wygładzanie danych za pomocą średnich ważonych 6.5. Centrowanie danych 6.6. Usuwanie trendu z szeregów czasowych 6.7. Szacowanie wskaźników wahań sezonowych 6.8. Usuwanie wahań sezonowych 6.9. Prognozowanie 6.10. Zastosowanie dyskretnej transformaty Fouriera 7. Funkcje matematyczne ............................................................................................... 215 215 215 216 217 219 221 222 223 223 224 225 227 7.0. Wprowadzenie 7.1. Korzystanie z funkcji sumujących 7.2. Dzielenie 7.3. Mnożenie 7.4. Przegląd funkcji wykładniczych i logarytmicznych 7.5. Używanie funkcji trygonometrycznych 7.6. Kontrolowanie znaków 7.7. Pierwiastkowanie 7.8. Zaokrąglanie i obcinanie liczb 7.9. Zamiana systemów liczbowych 7.10. Manipulowanie macierzami 7.11. Wykonywanie działań na wektorach Spis treści | 5 7.12. Wykorzystywanie funkcji arkuszowych w kodzie VBA 7.13. Wykonywanie działań na liczbach zespolonych 230 231 8. Dopasowywanie krzywej i regresja .......................................................................... 233 8.0. Wprowadzenie 233 8.1. Przeprowadzanie liniowego dopasowywania krzywej za pomocą wykresów 233 8.2. Przeprowadzanie dopasowania liniowego za pomocą funkcji arkusza 237 240 8.3. Liniowe dopasowywanie krzywej przy użyciu jednej funkcji arkuszowej 8.4. Przeprowadzanie wielokrotnej regresji liniowej 243 8.5. Generowanie nieliniowego dopasowania krzywej przy użyciu wykresów Excela 8.6. Dopasowywanie krzywych przy użyciu dodatku Solver 8.7. Ocenianie jakości dopasowania 8.8. Wyznaczanie przedziałów ufności 246 248 252 258 9. Rozwiązywanie równań ............................................................................................ 261 261 269 271 274 281 286 287 9.0. Wprowadzenie 9.1. Rozwiązywanie równań metodą graficzną 9.2. Iteracyjne rozwiązywanie równań nieliniowych 9.3. Automatyzowanie żmudnych zadań za pomocą VBA 9.4. Rozwiązywanie układów równań liniowych 9.5. Rozwiązywanie układów równań nieliniowych 9.6. Rozwiązywanie równań metodami klasycznymi 10. Numeryczne całkowanie i różniczkowanie .............................................................. 293 293 294 298 10.0. Wprowadzenie 10.1. Obliczanie całek oznaczonych 10.2. Implementacja metody trapezów w VBA 10.3. Zastosowanie całkowania numerycznego do wyznaczania środka ciężkości obszaru 10.4. Obliczanie momentu drugiego rzędu dla danego obszaru 10.5. Obliczanie całek podwójnych 10.6. Różniczkowanie numeryczne 300 303 304 307 11. Rozwiązywanie równań różniczkowych zwyczajnych ........................................... 315 315 315 11.0. Wprowadzenie 11.1. Rozwiązywanie zagadnień początkowych pierwszego rzędu 11.2. Zastosowanie metody Rungego-Kutty do rozwiązywania zagadnień początkowych drugiego rzędu 11.3. Rozwiązywanie układów równań sprzężonych 11.4. Rozwiązywanie zadań brzegowych metodą strzałów 321 326 332 6 | Spis treści 12. Rozwiązywanie równań różniczkowych cząstkowych ........................................... 337 337 339 341 345 12.0. Wprowadzenie 12.1. Rozwiązywanie równań różnicowych za pomocą Excela 12.2. Iteracyjne rozwiązywanie równań różnicowych za pomocą dodatku Solver 12.3. Rozwiązywanie zagadnień początkowych 12.4. Wykorzystanie Excela do rozwiązywania zagadnień sformułowanych przy użyciu metody elementów skończonych 349 13. Przeprowadzanie analizy optymalizacyjnej w Excelu ............................................. 353 353 354 357 13.0. Wprowadzenie 13.1. Tradycyjne programowanie liniowe z wykorzystaniem Excela 13.2. Analizowanie zagadnień z zakresu optymalizacji alokacji zasobów 13.3. Uzyskiwanie bardziej realistycznych wyników przy użyciu ograniczeń całkowitoliczbowych 13.4. Rozwiązywanie zadań kłopotliwych 13.5. Optymalizowanie projektów technicznych 13.6. Korzystanie z raportów Solvera 13.7. Zastosowanie algorytmu genetycznego do optymalizacji 362 364 370 373 377 14. Wprowadzenie do obliczeń finansowych .................................................................393 393 394 394 395 396 397 397 399 400 402 403 404 14.0. Wprowadzenie 14.1. Obliczanie wartości bieżącej 14.2. Obliczenie wartości przyszłej 14.3. Określanie wymaganej stopu zwrotu 14.4. Podwajanie zasobów pieniężnych 14.5. Ustalanie miesięcznych płatności 14.6. Analiza przepływów finansowych 14.7. Uzyskiwanie zakładanej wartości przyszłej 14.8. Wyznaczanie wartości bieżącej netto 14.9. Szacowanie stopy zwrotu 14.10. Rozwiązywanie zagadnień odwrotnych 14.11. Wyznaczanie progu rentowności Skorowidz ................................................................................................................... 407 Spis treści | 7 ROZDZIAŁ 9. 9.0. Wprowadzenie W tym rozdziale skoncentrujemy się na rozwiązywaniu równań za pomocą Excela. Rozwiązy- wanie równań w sensie ogólnym może polegać na znajdowaniu pierwiastków pojedynczych równań, wyznaczaniu wartości zmiennych niezależnych, przy których zmienna zależna przyj- muje zadaną wartość, lub rozwiązywaniu układów równań nieliniowych. Istnieje wiele tra- dycyjnych, ręcznych i komputerowych, metod rozwiązywania równań. Metoda Newtona ite- racyjnego znajdowania pierwiastków równań nieliniowych czy metoda eliminacji Gaussa z podstawianiem wstecznym dla rozwiązywania układów równań liniowych to przykłady ta- kich klasycznych metod. W swojej, zaliczanej już do klasyki, książce Introduction to Numerical Analysis Hildebrand przedstawił kilka klasycznych metod rozwiązywania równań. W książ- kach z serii Numerical Recipes możemy znaleźć algorytmy rozwiązywania równań zapisane w różnych językach programowania1. Algorytmy prezentowane w tych i innych książkach poświęconych metodom numerycznym są skuteczne i dobrze służą naukowcom oraz inży- nierom. W tym rozdziale chciałbym jednak pokazać, jak łatwo można, wykorzystując możli- wości Excela, rozwiązywać równania przy małej ilości programowania, a w niektórych przy- padkach bez żadnego programowania poza tworzeniem formuł arkuszowych. Po tym, jak zobaczymy, co Excel ma nam do zaoferowania, jeśli Czytelnik nadal będzie chciał napisać program dla jakiejś klasycznej metody, pokażę implementację kilku takich metod przy uży- ciu języka VBA i Excela. Niewiele można zyskać przez zastosowanie Excela lub innego tego typu narzędzia do roz- wiązywania równań, które można z łatwością rozwiązać ręcznie, stosując proste przekształ- cenia algebraiczne. Dlatego w prezentowanych tu przykładach skoncentrujemy się niemal wy- łącznie na równaniach nieliniowych. Należy również dodać, że opisywane tutaj techniki mogą być stosowane nie tylko do równań mających postać wyrażeń matematycznych, mimo iż większość przykładowych równań jest tutaj prezentowana w takiej właśnie postaci. Możemy na przykład mieć do czynienia z „rów- naniem” złożonym z kilku formuł realizujących wyszukiwanie w tabelach lub podobnych, które można łatwo utworzyć za pomocą arkusza kalkulacyjnego. Taki arkusz może potem 1 Patrz: F.B. Hildebrand, Introduction to Numerical Analysis, Dover Publications, 1974 oraz książki Numerical Re- cipes wydane przez Cambridge University Press, których autorami są: Press, Teukolsky, Vetterling i Flannery. 261 reprezentować równanie nieliniowe, które będziemy chcieli rozwiązać. Możemy mieć rów- nież do czynienia z obliczeniami rozłożonymi na wiele arkuszy. Ostatecznie, niezależnie od tego, czy mamy do czynienia z czysto matematycznym wyrażeniem, czy też ze skompliko- wanym arkuszem, zawsze możemy wskazać zmienne niezależne i zmienne zależne. Rozwią- zywanie tego typu równań polega zwykle na znajdowaniu takich wartości zmiennych nieza- leżnych, przy których zmienne zależne przyjmują wartości zadane. W pewnych przypadkach zmienne zależne mogą być tylko pośrednio związane z rozwiązywanym równaniem, a wte- dy możemy próbować znaleźć pierwiastek równania dla pewnej miary efektu (z taką sytu- acją mamy do czynienia na przykład przy dopasowywaniu krzywej metodą najmniejszych kwadratów). Ten typ obliczeń występuje również w zagadnieniach (lub częściowo pokrywa się z nimi) związanych z analizą „co-jeśli” oraz optymalizacją (patrz rozdział 13.) W obu tych przypad- kach często zachodzi potrzeba stosowania tych samych technik iteracyjnych, jakich używa się do rozwiązywania równań nieliniowych. Ponieważ w niniejszym rozdziale bardzo często będziemy korzystać z takich narzędzi Excela, jak polecenie Szukaj wyniku i dodatek Solver, chciałbym je pokrótce przedstawić, pokazując ich możliwości oraz występujące między nimi różnice. Zarówno polecenie Szukaj wyniku, jak i dodatek Solver są użytecznymi narzędziami umożli- wiającymi przeprowadzanie obliczeń iteracyjnych bez konieczności pisania odpowiednich programów. Ja wykorzystuję je bardzo często do rozwiązywania różnych problemów, począw- szy od rozwiązywania układów równań nieliniowych, przez analizę optymalizacyjną, po pro- gnozowanie. Dodatek Solver wykorzystałem nawet do rozwiązania równań ruchu zbudowa- nego i testowanego przeze mnie wodolotu. Z pozoru może wydawać się, że obydwa narzędzia robią to samo — pozwalają znaleźć w spo- sób iteracyjny wartość w komórce docelowej przez zmianę wartości w innej, powiązanej z nią komórce. Istnieją jednak między tymi narzędziami pewne różnice, o których należy pamię- tać, decydując się na wybór jednego z nich w celu rozwiązania określonego problemu. Dalsza część niniejszego wprowadzenia zawiera opis tych różnic. Polecenie Szukaj wyniku Polecenie Szukaj wyniku jest bardzo łatwe w użyciu. Wystarczy po prostu podać docelową wartość dla komórki docelowej i określić niezależną komórkę, której zawartość będzie zmie- niana, aby osiągnięta została zadana wartość docelowa. Na rysunku 9.1 pokazany jest inter- fejs tego narzędzia. Aby uzyskać do niego dostęp, należy z głównego menu wybrać polecenie Narzędzia/Szukaj wyniku. Rysunek 9.1. Okno dialogowe polecenia Szukaj wyniku 262 | Rozdział 9. Rozwiązywanie równań Pole Ustaw komórkę zawiera odwołanie do komórki docelowej, której zawartość ma być osta- tecznie równa wartości podanej w polu Wartość. Odwołanie do komórki, której zawartość ma być zmieniana, podajemy w polu Zmieniając komórkę. Komórka docelowa musi zawierać for- mułę, która z kolei musi bezpośrednio lub pośrednio odwoływać się do komórki podanej w polu Zmieniając komórkę. Polecenie Szukaj wyniku może być wywołane również z podprogramu napisanego w języku VBA. Nie musimy więc robić tego ręcznie, co ma szczególne znaczenie w sy- tuacji, gdy polecenie to musi być wywoływane wielokrotnie. Przykład takiej sytuacji zawiera receptura 9.3. Posługując się poleceniem Szukaj wyniku, możemy zmieniać zawartość tylko jednej komórki, co oznacza, że możemy go stosować tylko w przypadkach z jedną zmienną niezależną. Po- nadto musimy podać wartość docelową w sposób jawny, tzn. w postaci konkretnej liczby. Po- lecenie Szukaj wyniku nie pozwala nam również nakładać żadnych ograniczeń na zawartość komórki niezależnej (tzn. tej, którą podajemy w polu Zmieniając komórkę). Zgodnie z tym, co można przeczytać w krótkim artykule na stronie pomocy technicznej Mi- crosoftu (artykuł nr 100782), polecenie Szukaj wyniku wykorzystuje algorytm przeszukiwania liniowego, przyjmując wartości początkowe zbliżone do wartości podanej w komórce nieza- leżnej. Oznacza to, że wartość znajdująca się w komórce niezależnej w chwili uruchomienia polecenia Szukaj wyniku służy jako wartość początkowa mająca charakter prognozy wyniku. Należy o tym pamiętać szczególnie podczas rozwiązywania problemów, które mają więcej niż jedno rozwiązanie. Na przykład, jeśli próbujemy znaleźć pierwiastki równania trzeciego stopnia, podana przez nas wartość początkowa może prowadzić do jednego rozwiązania, podczas gdy inna wartość początkowa może dać inne rozwiązanie. Wynika stąd również, że jeśli polecenie Szukaj wyniku nie znajduje żadnego rozwiązania, wówczas należy spróbować zastosować inną wartość początkową. Polecenie Szukaj wyniku kończy obliczenia, gdy maksymalna wielkość zmian między iteracja- mi spada poniżej określonego progu. Przy ustawieniach domyślnych próg ten wynosi 0,001. Obliczenia zostaną przerwane również wtedy, gdy liczba iteracji przekroczy ustaloną war- tość (domyślnie 100). Kryteria te możemy zmienić, wybierając z głównego menu polecenie Narzędzia/Opcje, które otwiera okno dialogowe Opcje pokazane na rysunku 9.2. W oknie tym należy otworzyć zakładkę Przeliczanie. W jej środkowej części zobaczymy pole opcji Iteracja oraz dwa pola edycyjne Maksymalna liczba iteracji i Maksymalna zmiana. W polach tych należy wpisać odpowiednie wartości i zaznaczyć opcję Iteracja, aby nowe kryteria zaczę- ły obowiązywać. W przypadkach prostych problemów iteracyjnych polecenie Szukaj wyniku sprawdza się bar- dzo dobrze i jest łatwe w użyciu. Jednak gdy mamy do czynienia z bardziej złożonym pro- blemem lub potrzebujemy większej kontroli nad procesem iteracyjnym, należy sięgnąć po bardziej wyrafinowane narzędzie, jakim jest dodatek Solver. Dodatek Solver Dodatek Solver jest podobny w swoim działaniu do polecenia Szukaj wyniku — również po- zwala znaleźć wartość docelową w komórce docelowej przez iteracyjną zmianę wartości w ko- mórce niezależnej. Jednak jego możliwości są o wiele większe, co znajduje odzwierciedlenie w bardziej rozbudowanym interfejsie (rysunek 9.3). Dodatek Solver jest dostępny w menu 9.0. Wprowadzenie | 263 Rysunek 9.2. Okno dialogowe Opcje Narzędzia. Jeśli nie ma go wśród opcji tego menu, należy wybrać polecenie Narzędzia/Dodatki i na liście dostępnych dodatków zaznaczyć pozycję Dodatek Solver. Po wykonaniu tych czyn- ności można go uruchomić, wybierając polecenie Narzędzia/Solver. Rysunek 9.3. Okno dialogowe Solver - Parametry Dodatek Solver możemy wykorzystać do rozwiązywania problemów o wielu zmiennych. W takich przypadkach wartość docelowa w komórce docelowej będzie uzyskiwana przez zmianę wartości w kilku komórkach niezależnych. Komórki, których zawartość będzie zmie- niana, podajemy w polu Komórki zmieniane. Możemy przy tym zastosować dowolny styl od- wołań oraz odwołania do wielu komórek oddzielone średnikami. Wersja Solvera dołączana do Excela pozwala na podanie maksymalnie 200 komórek niezależnych. Dodatek Solver nie wymaga nawet podawania wartości docelowej. Zamiast tego możemy wybrać opcję minimalizowania lub maksymalizowania tej wartości. Odwołanie do komórki docelowej umieszczamy w polu Komórka celu (musi to być pojedyncza komórka zawierająca formułę). Przyciski wyboru opcji Równa pozwalają zdecydować, czy zawartość komórki do- celowej ma być minimalizowana, maksymalizowana lub ustalona przez przypisanie jej kon- kretnej wartości liczbowej. 264 | Rozdział 9. Rozwiązywanie równań Klikając przycisk Odgadnij, możemy pozwolić Solverowi „odgadnąć”, które komórki powinien uznać za niezależne dla danego problemu. Ja jednak unikam korzystania z tej możliwości. Sprowadza się ona do wybrania wszystkich komórek, do których od- wołuje się formuła w komórce docelowej i które same nie zawierają formuł. W przy- padku prostej formuły z niewielką liczbą zmiennych niezależnych, z których wszyst- kie chcemy zmieniać, takie rozwiązanie może być przydatne. Stwierdziłem jednak, że w większości przypadków, z jakimi miałem do czynienia, wybranych zostało zbyt dużo komórek — nawet te, których nie zamierzałem w ogóle zmieniać (bo zawierały np. stałe lub inne wartości, które z pewnych względów powinny pozostać niezmien- ne). Z drugiej strony, jeżeli mamy do czynienia z dużą liczbą zmiennych, ułatwieniem może być zlecenie Solverowi wybrania ich wszystkich, a następnie ręczne usunięcie z pola Komórki zmieniane tych, których nie chcemy zmieniać. To jeszcze nie wszystkie możliwości Solvera. Na komórki będące częścią rozwiązywanego problemu możemy nakładać określone ograniczenia. Mogą one mieć charakter równości lub granic górnych i dolnych. Jest to przydatne, gdy interesuje nas rozwiązanie z określonego ob- szaru lub gdy przeprowadzamy proces optymalizacji z ograniczeniami. Receptura 9.4 zawie- ra przykład wykorzystania tych ograniczeń przy rozwiązywaniu układu równań za pomocą Solvera. Z kolei rozdział 13. zawiera receptury z przykładami zastosowania Solvera do prze- prowadzania optymalizacji z ograniczeniami. Dodatek Solver (w wersji dołączanej do Excela) pozwala utworzyć maksymalnie 100 ograniczeń dla zagadnień nieliniowych i 200 dla zagad- nień liniowych. Solver korzysta również z większej liczby algorytmów znajdowania rozwiązań niż polecenie Szukaj wyniku. Do rozwiązywania problemów liniowych wykorzystuje metodę sympleks. Solver może także rozwiązywać problemy nieliniowe, wykorzystując metodę uogólnionego gradientu zredukowanego. W przypadku, gdy rozwiązywany problem wymaga zmiennych o wartościach całkowitych, Solver korzysta także z algorytmu podziałów i ograniczeń (ang. branch and bound). Niewielką kontrolę nad tymi algorytmami umożliwia okno dialogowe Solver - Opcje. Aby otworzyć to okno (pokazane na rysunku 9.4), należy kliknąć przycisk Opcje w oknie dialogo- wym Solver - Parametry. Rysunek 9.4. Okno dialogowe Solver - Opcje 9.0. Wprowadzenie | 265 Jak widać, jest tu znacznie więcej opcji w porównaniu z tylko dwiema dostępnymi dla pole- cenia Szukaj wyniku. Opcje Solvera to: Maksymalny czas Opcja ta reprezentuje maksymalny czas, jaki Solver może przeznaczyć na poszukiwanie wyniku. Po upływie tego czasu proces zostanie przerwany, a Solver wyświetli okno z in- formacją, że nie mógł znaleźć rozwiązania w wyznaczonym czasie. Wartość tę możemy ustalić maksymalnie na 32 767 sekund (nieco ponad 9 godzin), ale chyba nigdy nie zacho- dzi potrzeba ustawiania aż tak długiego czasu. Liczba iteracji Pozwala określić maksymalną liczbę iteracji, jakie Solver może wykonać w poszukiwaniu wyniku. Jeżeli ta liczba zostanie osiągnięta, zanim Solver znajdzie rozwiązanie, proces zo- stanie przerwany, a my zostaniemy poinformowani, że rozwiązanie nie zostało znalezio- ne przy zadanej liczbie iteracji. Wartość tę możemy ustalić maksymalnie na 32 767, chociaż w mojej praktyce nigdy nie musiałem stosować wartości większych niż 1000, a najczęściej pozostawiam wartość domyślną, czyli 100. Dokładność Jest to liczba dziesiętna z przedziału od 0 do 1, a jej wartość domyślna wynosi 1,0×10-6. Solver wykorzystuje ją do oceny, czy nałożone ograniczenia są spełnione. Im mniejsza jest ta liczba, tym większa precyzja. W praktyce bardzo rzadko zmieniam tę wartość. Tolerancja Tolerancja jest używana do określenia, czy ograniczenie wykorzystujące wartości całko- wite jest spełnione. Wyrażana jest w procentach, a jej wartość domyślna wynosi 5 . Jeśli nie stosujemy ograniczeń z wartościami całkowitymi, tolerancja nie jest używana. Zbieżność Ten parametr jest używany do określenia, czy Solver znalazł rozwiązanie. Jeśli dla pięciu ostatnich iteracji zmiana wartości w komórce docelowej nie przekracza wartości podanej w polu Zbieżność, Solver uznaje, że rozwiązanie zostało znalezione, i wyświetla okno ze sto- sowną informacją. Obliczenia będą kontynuowane do momentu osiągnięcia maksymalnej liczby iteracji lub maksymalnego czasu. Jeżeli zostanie napotkane którekolwiek z zadanych ograniczeń, Solver wyświetli komunikat z informacją, że nie może odnaleźć rozwiązania. Przyjmij model liniowy Jeżeli wiemy, że model rozwiązywanego problemu ma charakter liniowy, możemy zalecić Solverowi zastosowanie metody sympleks zamiast algorytmu uogólnionego gradientu zre- dukowanego. Gdy zaznaczymy tę opcję, Solver przeprowadzi kilka testów sprawdzających, czy według jego kryteriów nasz model jest rzeczywiście liniowy. W przypadku negatyw- nego wyniku tych testów zostanie wyświetlony komunikat z ostrzeżeniem, że model nie spełnia kryteriów liniowości. Przyjmij nieujemne Zaznaczenie tej opcji oznacza automatyczne nałożenie ograniczeń w postaci dolnej grani- cy wartości dla wszystkich komórek zmienianych. Jeśli wiemy, że wszystkie zmienne po- winny być zawsze nieujemne, i chcemy, aby rzeczywiście nie przyjmowały wartości ujem- nych, wówczas powinniśmy zaznaczyć tę opcję. W ten sposób możemy uniknąć ręcznego nakładania ograniczeń na każdą zmienną. 266 | Rozdział 9. Rozwiązywanie równań Automatyczne skalowanie W niektórych zagadnieniach wartości zmiennych niezależnych mogą znacznie różnić się pod względem wielkości od wartości zmiennych zależnych. W takich przypadkach do- brze jest tak przeskalować model, aby wartości wejściowe i wyjściowe były tego samego rzędu. Jeżeli wcześniej nie wykonaliśmy takiego skalowania, możemy zaznaczyć tę opcję, aby Solver zrobił to za nas. Wartości wejściowe i wyjściowe zostaną wówczas przeskalo- wane w wyniku podzielenia ich przez wartości początkowe podane w komórkach zmie- nianych i docelowych. W takich sytuacjach szczególnej wagi nabiera dobór odpowiednich wartości początkowych. Zawsze najlepiej jest ustalać te wartości zgodnie ze zdrowym roz- sądkiem i realiami rozważanego problemu. Pokaż wyniki iteracji Jeśli chcemy na bieżąco śledzić pracę Solvera podczas rozwiązywania danego problemu, możemy tę opcję zaznaczyć. Pojawi się wówczas okno dialogowe informujące nas o tym, że Solver przerwał pracę i że aktualne wyniki iteracji są wyświetlane w aktywnym arku- szu. Opcja ta przydaje się, gdy chcemy prześledzić wyniki poszczególnych etapów itera- cji, na przykład wtedy, gdy chcemy ustalić, dlaczego Solver nie może znaleźć rozwiąza- nia. W zwykłych warunkach nie zaznaczam tej opcji, ponieważ taki tryb pracy wymaga za każdym razem kliknięcia odpowiedniego przycisku nakazującego Solverowi przejście do następnego etapu. W sytuacji, gdy rozwiązanie danego problemu wymaga dużej licz- by iteracji, może to być bardzo żmudne i czasochłonne. Estymaty Tutaj możemy określić sposób, w jaki Solver będzie szacować wartości początkowe zmien- nych niezależnych. Opcja Styczna oznacza wykorzystanie do tego celu interpolacji liniowej, a opcja Kwadratowa — interpolacji kwadratowej, dającej lepsze rezultaty (szybsze znajdo- wanie rozwiązania) w przypadku problemów nieliniowych. Szczerze mówiąc, przy szyb- kości współczesnych procesorów trudno jest zauważyć istotną różnicę między tymi opcja- mi, jeśli chodzi o czas znajdowania rozwiązania typowych problemów. Pochodne Do obliczania gradientów Solver wykorzystuje metodę różnic skończonych, dając nam możliwość wyboru jednego z dwóch schematów tej metody. Opcja W przód pozwala wy- brać schemat różnicy przedniej (progresywnej), a opcja Centralne — schemat różnicy cen- tralnej. Różniczkowanie w oparciu o schemat różnicy centralnej wymaga większej ilości obliczeń, ale jest bardziej dokładne. I znów, jeśli chodzi o czas trwania obliczeń, różnica między tymi opcjami jest praktycznie niezauważalna na współczesnych komputerach. Ja najczęściej wybieram różniczkowanie centralne. (Analizę różnic między różniczką przed- nią a różniczką centralną zawiera receptura 10.6). Szukanie Tutaj możemy nakazać Solverowi rozwiązywanie problemu metodą Newtona lub metodą sprzężonego gradientu. Solver używa tych metod do określania kierunku poszukiwań podczas każdej iteracji. Metoda Newtona wymaga mniejszej liczby obliczeń niż metoda sprzężonego gradientu, ale za to angażuje więcej zasobów pamięciowych. Jeśli pamięć jest kwestią istotną — na przykład, gdy rozwiązujemy rozbudowany problem z dużą liczbą zmiennych i ograniczeń lub gdy zasoby pamięciowe naszego komputera są ograniczone — wówczas możemy wybrać metodę sprzężonego gradientu. 9.0. Wprowadzenie | 267 Istnieją jeszcze dwie inne opcje Solvera, które, mimo że nie wpływają bezpośrednio na proces obliczeniowy, są bardzo użyteczne. Opcjami tymi są: Załaduj model i Zapisz model. Ujmując rzecz krótko: Solver pozwala nam zapisywać modele i ponownie je wykorzystywać. To może być przydatne podczas sprawdzania różnych kombinacji zmiennych, ograniczeń i opcji Solvera pod kątem najlepszego dopasowania modelu do danego zagadnienia. Daje to również możli- wość odtworzenia modelu w dowolnym momencie (bez konieczności pamiętania wszystkich ustawień), jeśli zdecydujemy się na ponowne wykonanie tych samych obliczeń. Aby zapisać model, należy kliknąć przycisk Zapisz model. Pojawi się wówczas małe okno dia- logowe wzywające do określenia zakresu komórek, w których dane modelu mają być zapi- sane. Zakres ten powinien być położony z dala od innych danych i formuł znajdujących się w arkuszu, aby nie zostały one zastąpione danymi zapisywanego modelu. Do określenia tego zakresu wystarczy podanie jednej komórki — Solver umieści dane modelu w kolumnie, roz- poczynając od podanej komórki i zajmując trzy komórki plus jeszcze tyle, ile jest ograniczeń w zapisywanym modelu. Aby później odtworzyć zapisany model, wystarczy otworzyć okno Solver - Opcje i kliknąć przycisk Załaduj model. Zobaczymy wówczas małe okno dialogowe wzywające nas do zazna- czenia zakresu komórek, w których model został zapisany. Tym razem musimy zaznaczyć ten zakres w całości — zaznaczenie tylko pierwszej komórki nie wystarczy. Jeszcze jednym przydatnym elementem Solvera — jakiego nie ma w przypadku polecenia Szukaj wyniku — jest okno dialogowe Solver - Wyniki. Po znalezieniu rozwiązania Solver in- formuje nas o tym, otwierając okno dialogowe, takie jak to pokazane na rysunku 9.5. Rysunek 9.5. Okno dialogowe Solver - Wyniki Mamy tutaj do wyboru zachowanie wyniku uzyskanego przez Solver lub przywrócenie ory- ginalnych wartości początkowych. Ponadto możemy wydać Solverowi polecenie wygenero- wania jednego lub kilku spośród trzech raportów: Wyników, Wrażliwości i Granic. W tym celu należy kliknąć odpowiednią nazwę raportu. Możemy wybrać jeden, wszystkie lub dowolną kombinację dwóch raportów. Po kliknięciu przycisku OK Solver utworzy nowy arkusz dla każdego zaznaczonego raportu. Raporty te przydają się podczas interpretowania uzyskanych wyników, szczególnie w przypadku optymalizacji z ograniczeniami. Zagadnienia związane z tymi raportami i optymalizacją są szerzej opisane w rozdziale 13. Zawsze możemy skorzystać z pomocy, jaką dodatek Solver oferuje nam po kliknięciu przyci- sku Pomoc dostępnego we wszystkich jego oknach dialogowych. Czytelnikom pragnącym po- szerzyć swoją wiedzę na temat historii i wewnętrznych mechanizmów Solvera polecam arty- kuł pt. Design and Use of the Microsoft Excel Solver, autorstwa Daniela Fylstry, Leona Lasdona, Johna Watsona i Allana Warena, opublikowany na łamach pisma „Interfaces”, nr 5 (1998). 268 | Rozdział 9. Rozwiązywanie równań 9.1. Rozwiązywanie równań metodą graficzną Problem Chcemy graficznie wyznaczyć pierwiastki równania. Rozwiązanie Należy zapisać równanie w postaci y = g(x), obliczyć wartości zmiennej y dla zadanego zakre- su wartości zmiennej x i, wykorzystując możliwości Excela, utworzyć wykres tych wyników. Informacje na temat tworzenia wykresów w Excelu można znaleźć w rozdziale 4. Analiza Jeśli musimy znaleźć pierwiastki równania, którego nie możemy rozwiązać ręcznie, powinni- śmy rozpocząć od wykonania wykresu dla tego równania. Uzyskamy w ten sposób cenną in- formację na temat natury samego równania i położenia jego pierwiastków. Jeśli tylko możemy przedstawić równanie w postaci y = g(x), gdzie x oznacza zmienną niezależną, a y — zmien- ną zależną, wówczas z łatwością możemy zestawić kolumny z obliczeniami zmiennej y dla różnych wartości zmiennej x. To wszystko wydaje się bardzo proste i prawdę mówiąc, w kategoriach operacji arkuszowych, nie jest to trudne do wykonania. Jednak gdy mamy do czynienia z jakimś szczególnie zawi- kłanym równaniem, które na dodatek może mieć kilka pierwiastków, lepiej byłoby dyspono- wać wiedzą o położeniu tych pierwiastków. Wynika to stąd, że większość iteracyjnych metod rozwiązywania równań wymaga od nas podania przypuszczalnego położenia tych pierwiast- ków. Niektóre metody wymagają podania dwóch wartości początkowych, a inne — określenia przedziałów, w których pierwiastki mogą się znajdować. Sukces w wyznaczeniu pierwiast- ków za pomocą każdej z tych metod zależy więc od jakości naszych przewidywań i ustala- nych na tej podstawie wartości początkowych. (Konkretne przykłady można znaleźć w pozo- stałych recepturach z tego rozdziału). Rozważmy następujące równanie wielomianowe trzeciego stopnia: y a += bx 2 cx + 3 dx + To równanie ma oczywiście trzy pierwiastki, co oznacza, że istnieją trzy wartości zmiennej x, dla których zmienna y przyjmuje wartość równą zero. Bez trudu można to wykazać, wykre- ślając krzywą zależności y od x. Na rysunku 9.6 został przedstawiony wykres powyższego równania wykonany w Excelu. Jak widać, rzeczywiście krzywa przecina oś x w trzech punktach, a to dowodzi, że dla trzech wartości zmiennej x zmienna y przyjmuje wartość równą zero. Aby wykreślić tę krzywą, ze- stawiłem dwie kolumny: jedna zawiera wartości zmiennej x, a druga — obliczone wartości zmiennej y. Na rysunku 9.6 są to kolumny, odpowiednio: B i C. W tym przykładzie przy- jąłem dowolne wartości współczynników wielomianu i umieściłem je w kolumnach od C2 do C5. Na pasku formuły można zobaczyć formułę, jakiej użyłem do obliczania wartości y 9.1. Rozwiązywanie równań metodą graficzną | 269 Rysunek 9.6. Wykres wielomianu trzeciego stopnia odpowiadających poszczególnym wartościom zmiennej x. (Na zrzucie ekranu zaznaczona jest komórka C8 i na pasku formuły widoczna jest formuła zawarta w tej komórce). Formuła ta ma następującą postać: =$C$2+$C$3*B8+$C$4*B8^2+$C$5*B8^3. Teraz widzimy wyraźnie, jak zachowuje się nasze równanie i gdzie leżą jego pierwiastki. Wy- znaczenie dokładnych wartości zmiennej x, dla których y = 0 możemy zrealizować na kilka sposobów. Jeden z nich polega na aproksymowaniu wartości zmiennej x przez interpolację wartości w kolumnach B i C. W tym celu należy odszukać takie wartości zmiennej x, dla któ- rych zmienna y zmienia swój znak, a następnie należy dokonać interpolacji między tymi war- tościami. Inny sposób wyznaczenia miejsc zerowych może polegać na wykorzystaniu technik iteracyjnych. Przykłady takich technik zawarte są w recepturach 9.2 i 9.3. Z kolei receptura 9.6 pokazuje, jak wyznaczyć pierwiastki opisywanego tu równania trzeciego stopnia, stosując metodę Newtona i metodę siecznych. Jednak co robić w sytuacji, gdy mamy do czynienia z równaniem nie dającym się zapisać w postaci y = g(x)? Podczas próby wykonania wykresu takiego równania w opisany wyżej sposób napotkamy oczywiście problemy. Na szczęście Excel oferuje narzędzia, które możemy wykorzystać, aby ominąć te trudności. Wyjaśnienie tego, co mam na myśli, znajduje się w re- cepturze 9.2. 270 | Rozdział 9. Rozwiązywanie równań 9.2. Iteracyjne rozwiązywanie równań nieliniowych Problem Użytkownik chciałby rozwiązać w Excelu równanie nieliniowe metodą iteracyjną, ale nie wie, jak się do tego zabrać. Rozwiązanie Należy wykorzystać polecenie Szukaj wyniku lub dodatek Solver. Informacje na temat różnic w korzystaniu z tych narzędzi oraz zalet i wad każdego z nich zawarte zostały we wprowa- dzeniu do niniejszego rozdziału. Analiza W tej recepturze chcę pokazać na konkretnym przykładzie, jak rozwiązać równanie nielinio- we za pomocą takich narzędzi jak polecenie Szukaj wyniku i dodatek Solver. Nasze rozważa- nia skoncentrujemy na następującym równaniu: , 0 242 C f = log( CR N f ) Równanie to jest wykorzystywane w obliczeniach oporów ruchu statków do szacowania war- tości współczynnika oporu tarcia Cf jako funkcji liczby Reynoldsa RN. Istnieją jeszcze inne równania służące do tego samego celu, ale wybrałem ten klasyczny przypadek, ponieważ nie można go zapisać w postaci y = g(x). Musimy więc sięgnąć po metody iteracyjne, aby obliczyć wartość Cf odpowiadającą danej wartości RN. Pokażę teraz, jak wykorzystać polecenie Szukaj wyniku i dodatek Solver do rozwiązania tego problemu. W obu przypadkach musimy najpierw przekształcić to równanie do następującej postaci: 0 = , 0 242 C f - log( CR N f ) Teraz możemy zastosować jedno z wymienionych wyżej narzędzi do iteracyjnego poszuki- wania takiej wartości Cf, która przy zadanej wartości RN sprawi, że prawa strona tego równa- nia będzie równa 0. Wcześniej jednak musimy prawą stronę tego równania umieścić w arku- szu w postaci formuły, tak jak zostało to pokazane na rysunku 9.7. Formuła reprezentująca prawą stronę naszego równania znajduje się w komórce C5. Postać tej formuły jest następująca: =0,242/PIERWIASTEK(C4)-LOG(C3*C4;10), co zresztą widać na pasku formuły na rysunku 9.7. Komórka C3 zawiera zadaną liczbę Reynoldsa, a komórka C4 — współczynnik tarcia wyznaczony za pomocą Solvera lub polecenia Szukaj wyniku. Aby sko- rzystać z tych narzędzi, musimy w komórce C4 umieścić wartość początkową, którą powinna być przewidywana przez nas wartość współczynnika tarcia. Jak już wskazywałem we wstępie do niniejszego rozdziału, taka wartość początkowa powinna być dobierana bardzo starannie. 9.2. Iteracyjne rozwiązywanie równań nieliniowych | 271 Rysunek 9.7. Przykład równania nieliniowego Wyznaczanie Cf za pomocą polecenia Szukaj wyniku Aby znaleźć rozwiązanie za pomocą tego narzędzia, należy z głównego menu wybrać po- lecenie Narzędzia/Szukaj wyniku. Zostanie okno dialogowe Szukanie wyniku pokazane na ry- sunku 9.8. Rysunek 9.8. Okno dialogowe Szukanie wyniku dla przykładowego równania nieliniowego W polu Ustaw komórkę należy wpisać C5 lub kliknąć małą ikonę na prawo od pola edycyjne- go, aby tymczasowo przejść do arkusza i zaznaczyć tam komórkę C5 (po jej zaznaczeniu na- leży wcisnąć klawisz Enter, aby powrócić do okna Szukanie wyniku). Przypominam, że komór- ka C5 zawiera formułę reprezentującą prawą stronę naszego równania i chcemy, aby wynik tej formuły był równy zero. Dlatego w polu Wartość należy wpisać 0. Jest to nasza wartość docelowa. Teraz w polu Zmieniając komórkę należy umieścić odwołanie do komórki C4, która zawiera wartość współczynnika tarcia. Po kliknięciu przycisku OK zawartość komórki C4 będzie dotąd zmieniana, aż formuła w komórce C5 osiągnie wartość dostatecznie bliską 0. Wynik tych obliczeń jest pokazany na rysunku 9.7. Widać tu, że reszta w komórce C5 wynosi 9,6×10–4. 272 | Rozdział 9. Rozwiązywanie równań Powinienem dodać, że przyjęta przeze mnie wartość początkowa dla komórki C4 wynosiła 1×10–7. Mogłem przyjąć wartość 0, ponieważ wiem, że współczynnik tarcia powinien być bar- dzo małą liczbą, ale tego nie zrobiłem, bo to wywołałoby błąd dzielenia przez zero. Wystą- pienie takiego błędu podczas pracy polecenia Szukaj wyniku powoduje awaryjne zakończenie jego działania i wyświetlenie komunikatu o następującej treści: Formuła w komórce musi dawać w wyniku liczbę. Reszta wynosząca 9,6×10–4 nie oznacza złego wyniku, ale możemy uczynić go jeszcze lepszym. Jak wspominałem we wprowadzeniu, możemy zmienić ustawienia zbieżności dla polecenia Szukaj wyniku w oknie dialogowym Opcje (rysunek 9.2). Aby go otworzyć, należy z głównego menu wybrać polecenie Narzędzia/Opcje. Na zakładce Przeliczanie zaznaczyłem opcję Iteracja i ustawiłem parametr Maksymalna liczba iteracji na 5 000, a parametr Maksymalna zmiana na 1×10–8. Przy takich ustawieniach polecenie Szukaj wyniku znalazło rozwiązanie ze znacznie mniejszą resztą wynoszącą 3,4×10–9. Współczynnik tarcia dla zadanej liczby Reynoldsa wy- niósł przy tym 0,0124. Wyznaczanie Cf za pomocą Solvera Polecenie Szukaj wyniku sprawdza się dobrze w tym prostym przykładzie i w praktyce nie musielibyśmy tutaj stosować bardziej zaawansowanego narzędzia, jakim jest dodatek Solver. Jednak w celach ilustracyjnych chcę pokazać, jak można go wykorzystać do rozwiązania tego samego problemu. Arkusz należy przygotować tak samo jak poprzednio, natomiast z główne- go menu należy wybrać polecenie Solver zamiast Szukaj wyniku. Okno dialogowe Solver - Pa- rametry dla tego konkretnego przykładu zostało pokazane na rysunku 9.9. Rysunek 9.9. Okno dialogowe Solver - Parametry dla przykładu z równaniem nieliniowym W polu Komórka celu umieściłem odwołanie do komórki C5. Spośród opcji Równa wybrałem Wartość i w polu obok wpisałem 0. Tak jak poprzednio, komórka C5 zawiera prawą stronę naszego równania, która ma przyjąć wartość równą 0. Następnie w polu Komórki zmieniane umieściłem odwołanie do komórki C4, która zawiera wartość początkową współczynnika tarcia. Kliknięcie przycisku Rozwiąż rozpoczyna proces rozwiązywania równania. Solver rzeczywiście znajduje rozwiązanie, dając wartość współczynnika tarcia równą 0,0124 przy reszcie wynoszącej –5,45×10–7. Są to wartości porównywalne z uzyskanymi za pomocą polecenia Szukaj wyniku. 9.2. Iteracyjne rozwiązywanie równań nieliniowych | 273 Spróbujmy teraz, trochę dla zabawy, rozwiązać ten sam problem, ale z większą wartością po- czątkową współczynnika tarcia, np. 1 lub 5 albo jeszcze więcej. W takim przypadku Solver nie potrafi znaleźć rozwiązania! Przy takich wartościach początkowych Solver przeskakuje rze- czywiste rozwiązanie i w efekcie próbuje przyjąć dla dalszych iteracji wartość współczynnika tarcia równą 0, co oczywiście oznacza w tym przypadku błąd dzielenia przez zero i przed- wczesne zakończenie działania Solvera. Ten prosty przykład świadczy dobitnie o tym, jak wielkie znaczenie ma właściwy dobór wartości początkowych. Ponadto pokazuje, że jeżeli pierwsza próba jest nieudana, należy zmienić wartość początkową i spróbować ponownie. 9.3. Automatyzowanie żmudnych zadań za pomocą VBA Problem Wiemy już, jak stosować polecenie Szukaj wyniku i dodatek Solver do rozwiązywania równań nieliniowych, ale kłopot polega na tym, że potrzebujemy rozwiązać takie równanie dla okre- ślonego przedziału wartości, a ręczne wykonanie tego zadania staje się niezwykle żmudne, bo wymaga wielokrotnego stosowania tych narzędzi. Rozwiązanie Można zautomatyzować ten proces, wykorzystując język VBA, który umożliwia programowe wywoływanie Solvera lub polecenia Szukaj wyniku. Analiza Dodatek Solver i polecenie Szukaj wyniku są wygodnymi narzędziami do rozwiązywania roz- maitych problemów w Excelu. Ale jeszcze lepsze jest to, że możemy używać tych narzędzi bezpośrednio w kodzie VBA, nawet bez ręcznego otwierania ich okien dialogowych, zazna- czania komórek itp. Oznacza to, że możemy zautomatyzować wiele zadań wymagających użycia Solvera lub polecenia Szukaj wyniku, oszczędzając w ten sposób wiele cennego czasu. Rozważmy ponownie przykładowe równanie z receptury 9.2. Dla wygody przytaczam je tu- taj jeszcze raz: , 0 242 C f = log( CR N f ) Powstaje pytanie: jak utworzyć wykres tego równania? Gdybyśmy chcieli wykreślić krzywą zależności Cf od RN, musielibyśmy dla każdego punktu wyznaczającego tę krzywą użyć Solvera lub polecenia Szukaj wyniku. Nawet gdybyśmy nie zamierzali tworzyć wykresu, lecz tylko samą tabelę z wartościami Cf dla różnych wartości RN, wówczas również czekałby nas trud ręczne- go, wielokrotnego stosowania tych narzędzi dla każdej wartości RN. A można sobie wyobra- zić jeszcze bardziej złożony przypadek, gdy równanie zawiera jeszcze jeden parametr, który również chcemy systematycznie zmieniać. To musiałoby być naprawdę nudne zajęcie! 274 | Rozdział 9. Rozwiązywanie równań Na szczęście twórcy Solvera i polecenia Szukaj wyniku umożliwili wykorzystanie możliwo- ści tych narzędzi z poziomu języka VBA, dzięki czemu możemy tego typu zadania zautoma- tyzować. Automatyzowanie polecenia Szukaj wyniku Aby pokazać zautomatyzowany sposób wykorzystania polecenia Szukaj wyniku, rozpocznie- my od wprowadzenia pewnych modyfikacji do arkusza pokazanego na rysunku 9.7. Naszym celem będzie utworzenie podprogramu VBA wywołującego polecenie Szukaj wyniku, które z kolei będzie obliczać wartości Cf dla określonego zakresu wartości RN. Aby móc uruchomić podprogram VBA, dodałem do arkusza element sterujący w postaci przycisku. Kliknięcie te- go przycisku będzie rozpoczynało proces rozwiązywania. Wstawiłem również tabelę, która będzie automatycznie wypełniana obliczonymi wartościami Cf i RN. Na koniec dodałem jesz- cze wykres krzywej reprezentującej uzyskane wyniki. Zmodyfikowany w ten sposób arkusz został przedstawiony na rysunku 9.10. Rysunek 9.10. Przykład automatyzacji polecenia Szukaj wyniku Jak widać, jest tutaj kilka nowych elementów. Każdy z nich postaram się objaśnić. Elementy pochodzące z oryginalnego arkusza zostały wyróżnione czcionką pochyłą (patrz rysunek 9.10), natomiast cała reszta jest tutaj nowa. 9.3. Automatyzowanie żmudnych zadań za pomocą VBA | 275 Pierwszą rzeczą, jaką dodałem do tego arkusza, było wpisanie w komórkach C1 i C2 wartości określających początek i koniec rozważanego przedziału wartości RN. Na lewo od tych komó- rek umieściłem odpowiednie etykiety objaśniające. Następnie w komórkach od A9 do C30 skonstruowałem tabelę przeznaczona do przechowy- wania wartości RN i odpowiadających im wartości Cf . Początkowo tabela ta zawierała tylko nagłówki kolumn w wierszu 9 i numery wierszy w kolumnie A. Kolumny Cf i RN pozostały na razie puste — zostaną wypełnione po uruchomieniu podprogramu VBA, który wkrótce zaprezentuję. Po zestawieniu tabeli utworzyłem wykres, który będzie prezentował zależność Cf od RN. Jako że tabela danych na razie jest pusta, na wykresie nie zobaczymy żadnej krzywej. Jednak przy- gotowanie wykresu już na tym etapie oznacza, że po zakończeniu obliczeń wykonywanych przez podprogram VBA zostanie on automatycznie zaktualizowany i oczekiwana krzywa po- jawi się. (Jeśli Czytelnik potrzebuje dodatkowych informacji na temat tworzenia wykresów, proponuję zajrzeć do rozdziału 4.). Ostatnim elementem, jaki musimy jeszcze dodać do arkusza, jest przycisk, który na rysunku 9.10 nosi nazwę Oblicz Cf. Excel pozwala nam dodawać do arkuszy różne elementy sterujące, umożliwiając w ten sposób konstruowanie nawet bardzo rozbudowanych interfejsów gra- ficznych (ang. GUI — Graphical User Interface) podobnych do tych, które znamy ze standar- dowych programów działających w systemie Windows. Dla celów niniejszego przykładu po- każę, jak wstawić do arkusza prosty przycisk i powiązać go z kodem VBA, tak aby kliknięcie tego przycisku wywoływało jakąś akcję. W tym przypadku kod VBA obliczy wartości Cf dla całego zakresu wartości RN i wypełni tabelę danych. Zanim przystąpimy do umieszczenia przycisku w arkuszu, musimy upewnić się, że pasek narzędzi Przybornik formantów jest widoczny w Excelu. Pasek ten został pokazany na rysun- ku 9.11. Rysunek 9.11. Pasek narzędzi Przybornik formantów Jeśli jest niewidoczny, należy z głównego menu Excela wybrać polecenie Widok/Paski narzędzi/ Przybornik formantów. Aby dodać przycisk do arkusza, wystarczy kliknąć ikonę Przycisk polecenia (rysunek 9.11), a następnie kliknąć w obrębie arkusza, tam, gdzie chcemy przycisk umieścić. Po wykonaniu tych czynności na naszym arkuszu powinien pojawić się nowy element w postaci przycisku. Początkowo będzie on otoczony uchwytami (kółeczkami), które możemy przeciągać, zmienia- jąc w ten sposób rozmiary przycisku. Możemy również przeciągać sam przycisk, aby zmienić jego położenie. W momencie dodawania nowego elementu sterującego do arkusza Excel przechodzi w tryb projektowania, co jest sygnalizowane wciśnięciem przycisku Tryb projektowania na pasku Przy- bornik formantów. W tym trybie możemy zmieniać rozmiary elementów kontrolnych i prze- mieszczać je. Po wyjściu z trybu projektowania kliknięcie elementu kontrolnego spowoduje 276 | Rozdział 9. Rozwiązywanie równań jego uaktywnienie. Tryb projektowania można na przemian włączać i wyłączać, klikając wspo- mniany już przycisk Tryb projektowania na pasku narzędzi. Elementy sterujące dostępne na pasku Przybornik formantów są kontrolkami ActiveX. Jeśli w trybie projektowania klikniemy taki element prawym przyciskiem myszy, wówczas z menu podręcznego możemy wybrać polecenie Właściwości, które otwiera okno udostępniające wszystkie właściwości tego elementu. W tym oknie możemy dokonać edycji wielu charakterystyk wybranego elementu kontrolnego. Na przykład zmiana właściwości Caption oznacza zmianę tekstu wyświetlanego na samym elemen- cie kontrolnym. W naszym przykładzie ten tekst został zmieniony na Oblicz Cf, co widać na rysunku 9.10. Po dodaniu przycisku musimy powiązać go jakimś kodem VBA, jeśli kliknięcie tego przyci- sku ma powodować podjęcie jakiejkolwiek akcji. A zatem następny etap będzie polegał na napisaniu odpowiedniego podprogramu i powiązaniu go z nowym przyciskiem. Przykład 9.1 przedstawia procedurę, jaką przygotowałem na potrzeby naszego przykładu. (Czytelników, którzy poszukują informacji na temat dołączania takich jak ta procedur VBA do arkuszy Excela, odsyłam do rozdziału 2.). Przykład 9.1. Procedura ComputeCf Public Sub ComputeCf() Dim inc As Double With Worksheets( Równanie nieliniowe ) inc = (.Range( Rn_2 ) - .Range( Rn_1 )) / 20 For i = 0 To 20 .Range( Rn ) = .Range( Rn_1 ) + (inc * i) .Range( Fx ).GoalSeek goal:=0, ChangingCell:=.Range( Cf ) .Cells(10 + i, 2) = .Range( Rn ) .Cells(10 + i, 3) = .Range( Cf ) Next i End With End Sub Procedura, której nadałem nazwę ComputeCf, pobiera granice przedziału zmienności RN za- pisane w komórkach C1 i C2, a następnie, wykorzystując polecenie Szukaj wyniku, oblicza 20 wartości Cf . Po obliczeniu wszystkich wartości Cf wyniki są zapisywane w tabeli obejmującej komórki arkusza od A9 do C30. Pierwsza linia procedury zawiera deklarację zmiennej lokalnej o nazwie inc, która będzie przechowywać przyrost zmiennej RN (przyrost ten jest obliczany w dalszej części procedury), czyli różnicę między kolejnymi wartościami tej zmiennej. W następnej linii występuje instrukcja With. W językach zorientowanych obiektowo, takich jak VBA, gdzie dostęp do składników obiektu realizowany jest poprzez składnię „z kropką” (widoczną w powyższym przykładzie), instrukcje With pełnią bardzo użyteczną rolę. Pozwa- lają one określić, który obiekt ma być przyjęty jako domyślny przy odwoływaniu się do po- szczególnych składników. Dzięki temu możemy zaoszczędzić pisania, opuszczając w tych odwołaniach nazwę obiektu-rodzica. Na przykład w tym podprogramie użyłem instrukcji With Worksheets( Równanie nieliniowe ), informując w ten sposób VBA, że chcę uzyskać dostęp do obiektu typu Worksheet o nazwie Równanie nieliniowe. Jest to nazwa arkusza 9.3. Automatyzowanie żmudnych zadań za pomocą VBA | 277 wykorzystywanego w naszym przykładzie i można ją zobaczyć na zakładce arkusza w dolnej części rysunku 9.10. Teraz za każdym razem, gdy używam instrukcji takiej jak .Range( Rn_2 ), odwołując się w tym przypadku do zakresu komórek, VBA przyjmuje domyślnie, że jest to odwołanie do składnika typu Range obiektu Worksheet o nazwie Równanie nieliniowe. Bez użycia instrukcji With musielibyśmy napisać: Worksheets( Równanie nieliniowe ).Range ( Rn_2 ). Instrukcji With musi towarzyszyć instrukcja zamykająca End With, tak jak w przy- kładzie 9.1. W ten sposób działanie instrukcji With zostaje ograniczone do tych in- strukcji, które są zawarte między With a End With. Następna instrukcja po With oblicza wartość zmiennej inc, odwołując się do komórek zawie- rających granice przedziału zmiennej RN za pomocą obiektów typu Range, tak jak to zostało wcześniej opisane. Jako odwołania do tych komórek zastosowałem ich nazwy, Rn_2 i Rn_1. Uważam, że nazwy są bardziej czytelne niż standardowe odwołania. O nadawaniu nazw ko- mórkom można przeczytać w recepturze 1.14. Po obliczeniu wartości inc podprogram wchodzi w pętlę For przebiegającą przez wszystkie wartości RN. Dla każdej wartości RN obliczana jest iteracyjnie wartość Cf za pomocą polecenia Szukaj wyniku. Oto objaśnienie każdej linii tej pętli: .Range( Rn ) = .Range( Rn_1 ) + (inc * i) W tej linii obliczana jest kolejna wartość RN, która zapisywana jest w komórce o nazwie Rn (komórka C3). .Range( Fx ).GoalSeak goal:=0, ChangingCell:=.Range( Cf ) Tutaj następuje wywołanie polecenia Szukaj wyniku (ang. Goal Seak) dla komórki o nazwie Fx (komórka C4) jako komórki docelowej, której zawartość ma osiągnąć wartość 0. War- tość docelowa jest tutaj ustalana przez goal:=0. Parametr Zmieniając komórkę jest określa- ny przez ChangingCell:=.Range( Cf ). Wszystko to stanowi informację dla polecenia Szukaj wyniku, że zawartość komórki Cf ma być tak zmieniana, aby w komórce Fx pojawi- ła się wartość równa zero. .Cells(10 + i, 2) = .Range( Rn ) W tej linii zawartość komórki Rn jest zapisywana w tabeli danych, począwszy od komórki B10. Numer wiersza jest zwiększany o 1 przy każdym przejściu pętli. Dla wygody zasto- sowałem tutaj styl odwołań W1K1. (Więcej informacji na temat tego stylu zawiera recep- tura 1.6). .Cells(10 + i, 3) = .Range( Cf ) Wartość Cf wyznaczona przez polecenie Szukaj wyniku jest zapisywana do tabeli danych. Tutaj również zastosowałem styl W1K1. Teraz możemy powiązać napisaną właśnie procedurę ComputeCf z przyciskiem dodanym wcześniej do arkusza. W tym celu należy najpierw włączyć tryb projektowania (w sposób opisany wcześniej), a następnie kliknąć ten przycisk prawym przyciskiem myszy. Z otwarte- go w ten sposób menu podręcznego należy wybrać polecenie Wyświetl kod, które przeniesie nas do edytora VBA, gdzie powinniśmy zobaczyć nową procedurę przygotowaną dla nasze- go przycisku. Procedura ta została pokazana w przykładzie 9.2. 278 | Rozdział 9. Rozwiązywanie równań Przykład 9.2. Procedura przycisku Private Sub CommandButton1_Click() End Sub Jest to procedura obsługująca zdarzenie polegające na kliknięciu przycisku. Ponieważ nasz przycisk nosi domyślną nazwę CommandButton1, stąd nazwą tej procedury jest CommandBut- ton1_Click. Nazwę elementu sterującego typu ActiveX można zmienić przez zmianę właściwości Name w oknie właściwości tego elementu. Aby dokonać takiej zmiany, należy włączyć tryb projektowania, kliknąć wybrany element prawym przyciskiem myszy i z menu podręcznego wybrać polecenie Właściwości. W otwartym w ten sposób oknie Properties można zmodyfikować parametry wybranego elementu. Opisywana procedura jest na razie pusta. Umieścimy w niej wywołanie napisanej wcześniej procedury ComputeCf, tak jak to zostało pokazane w przykładzie 9.3. Przykład 9.3. Wywołanie procedury ComputeCf Private Sub CommandButton1_Click() ComputeCf End Sub Jeśli teraz klikniemy przycisk w Excelu (po wyłączeniu trybu projektowania), spowoduje to wykonanie procedury obliczającej wartości Cf dla zadanego zakresu wartości RN. Wyniki tych obliczeń zobaczymy w tabeli oraz na wykresie. Jeżeli zechcemy wyznaczyć wartości Cf dla innego zakresu wartości RN, jedyne, co musimy zrobić, to wpisać granice nowego zakresu (w komórkach C1 i C2) i ponownie kliknąć przycisk. Excel zawiera jeszcze jeden zestaw elementów sterujących dostępnych na pasku For- mularze (widocznym po wybraniu polecenia Widok/Paski narzędzi/Formularze). Elemen- ty te wyglądają tak jak stosowane w tym przykładzie elementy ActiveX, ale dają nam większą kontrolę i różnią się od nich sposobem modyfikowania właściwości oraz przypisywania funkcji. Ogólnie, elementy typu ActiveX są nowsze i lepiej oprogra- mowane, a elementy formularzowe są starsze i utrzymywane dla zachowania zgod- ności ze starszymi wersjami programów. Automatyzowanie Solvera Wykorzystanie Solvera możemy zautomatyzować w sposób analogiczny do tego, jaki zasto- sowaliśmy w celu zautomatyzowania użycia polecenia Szukaj wyniku. W rzeczywistości, aby w rozważanym przykładzie zamiast polecenia Szukaj wyniku zastosować dodatek Solver, wy- starczą niewielkie modyfikacje kodu VBA. Wszystkie elementy arkusza mogą pozostać takie same jak poprzednio, a zmodyfikować musimy tylko procedurę ComputeCf. Nową wersję tej procedury przedstawia przykład 9.4. Przykład 9.4. Procedura ComputeCf z użyciem Solvera Public Sub ComputeCf() Dim inc As Double With Worksheets( Równanie nieliniowe ) inc = (.Range( Rn_2 ) - .Range( Rn_1 )) / 20 9.3. Automatyzowanie żmudnych zadań za pomocą VBA | 279 For i = 0 To 20 .Range( Rn ) = .Range( Rn_1 ) + (inc * i) Tutaj rozpoczyna się nowy kod: .Range( Cf ) = 0.001 SolverOK SetCell:=Range( Fx ), MaxMinVal:=3, ValueOf:=0, ByChange:=Range( Cf ) SolverSolve UserFinish:=True SolverFinish KeepFinal:=1 .Cells(10 + i, 2) = .Range( Rn ) .Cells(10 + i, 3) = .Range( Cf ) Next i End With End Sub Zmiana w stosunku do poprzedniej wersji obejmuje cztery linie kodu następujące po komen- tarzu Tutaj rozpoczyna się nowy kod:. Te nowe linie to: .Range( Cf ) = 0.001 W tej linii ustawiana jest wartość początkowa zmiennej Cf. W tym przypadku będzie ona wynosić 0,001 i będzie służyć jako wartość wstępna, ustawiana przed wywołaniem Solvera. SolverOK SetCell:=Range( Fx ), MaxMinVal:=3, ValueOf:=0, ByChange:=Range( Cf ) W tej linii następuje wywołanie wewnętrznej procedury Solvera o nazwie SolverOK w celu inicjalizacji modelu Solvera. Komórka docelowa (jest nią komórka o nazwie Fx) jest okre- ślana przez SetCell:=Range( Fx ). Ponieważ chcemy jej przypisać konkretną wartość, ustalamy MaxMinVal:=3. Gdyby chodziło nam o maksymalizowanie wartości docelowej, należałoby napisać MaxMinVal:=1, a w przypadku minimalizowania — MaxMinVal:=2. Ko- lejny zapis ValueOf:=0 jest informacją dla Solvera, że wartość docelowa wynosi 0. Z kolei zapis ByChange:=Range( Cf ) informuje Solver, że ma poszukiwać rozwiązania przez zmianę zawartości komórki o nazwie Cf. SolverSolve UserFinish:=True Wywołanie procedury SolverSolve uruchamia Solvera. Przypisanie UserFinish:=True ma na celu poinformowanie Solvera, aby po odnalezieniu rozwiązania nie wyświetlał okna dialogowego Solver - Wyniki. Wyświetlanie tego okna może być przydatne np. podczas debugowania tego typu procedur, a wówczas należy ustawić ten parametr na False. SolverFinish KeepFinal:=1 Ta linia nakazuje Solverowi zachować znalezione wyniki. Po dokonaniu tych zmian w kodzie procedury i uruchomieniu jej przez kliknięcie przycisku na przykładowym arkuszu, powinniśmy ujrzeć wyniki bardzo podobne do tych, jakie uzyska- liśmy, stosując polecenie Szukaj wyniku. Jeśli podczas uruchamiania nowej procedury wystąpi błąd sygnalizowany komuni- katem o treści Sub or Function not Defined, należy sprawdzić, czy Solver jest zareje- strowany w środowisku VBA. W tym celu należy z głównego menu edytora VBA wybrać polecenie Tools/References i w otwartym w ten sposób oknie dialogowym Re- ferences upewnić się, że zaznaczona jest pozycja SOLVER. Na koniec należy kliknąć przycisk OK. 280 | Rozdział 9. Rozwiązywanie równań Informacje dodatkowe Za pomocą kodu VBA można kontrolować znacznie więcej aspektów Solvera, niż tutaj poka- załem. Więcej informacji można znaleźć w pomocy Excela oraz na stronie internetowej po- mocy technicznej Microsoftu http://support.microsoft.com, gdzie należy odszukać artykuły na temat Solvera (wśród wielu tematów dotyczących innych produktów tej firmy). Dla pragną- cych zgłębić tajniki stosowania Solvera najbardziej pomocnym może okazać się artykuł pt. „How to create Visual Basic macros by using Excel Solver in Excel 97” (artykuł nr 843304)2. 9.4. Rozwiązywanie układów równań liniowych Problem Chcielibyśmy rozwiązać układ równań liniowych za pomocą Excela. Rozwiązanie Układy równań liniowych można rozwiązywać w Excelu, stosując różne podejścia. Na przy- kład możemy napisać program w języku VBA będący implementacją jednego ze standardo- wych algorytmów rozwiązywania układów równań. Takie podejście jest niewątpliwie skutecz- ne, ale czasochłonne — wymaga napisania programu i jego przetestowania w celu usunięcia ewentualnych błędów. Jako że Excel udostępnia nam narzędzia do wykonywania działań na macierzach, możemy układ równań przedstawić w postaci macierzowej i znaleźć rozwiązanie przez odwrócenie odpowiedniej macierzy. Jeszcze łatwiejszy sposób polega na wykorzysta- niu dodatku Solver. W poniższe
Pobierz darmowy fragment (pdf)

Gdzie kupić całą publikację:

Excel w nauce i technice. Receptury
Autor:

Opinie na temat publikacji:


Inne popularne pozycje z tej kategorii:


Czytaj również:


Prowadzisz stronę lub blog? Wstaw link do fragmentu tej książki i współpracuj z Cyfroteką: