Приближение с помощью кривых

Приближение кривой, имеющей шум (случайные отклонения), с помощью модели асимметричного пика с итеративным процессом (алгоритм Гаусса — Ньютона с переменным коэффициентом затухания α).
Сверху: исходные данные и модель.
Снизу: изменение нормализованной суммы квадратов отклонений.

Приближение с помощью кривых [1][2] — это процесс построения кривой или математической функции, которая наилучшим образом приближается к заданным точкам [3] с возможными ограничениями на кривую [4]. Для построения такого приближения может использоваться либо интерполяция [5], где требуется точное прохождение кривой через точки, либо сглаживание[англ.][6][7], когда «сглаживающая» функция проходит через точки приближённо. Связанный раздел — регрессионный анализ[8][9], который фокусируется, главным образом, на вопросах статистического вывода, таких как, какая неопределённость заключена в кривой, которая приближает данные с некоторыми случайными ошибками. Построенные кривые могут быть использованы для визуализации данных [10][11], для вычисления значений функции в точках, в которых значение не задано [12] и для определения связи между двумя и более переменными[13]. Экстраполяция означает использование полученной кривой за пределами данных, полученных из наблюдения [14], и порождает некоторую неопределённость[15], поскольку может зависеть от метода построения кривой.

Различные типы приближения с помощью кривых

Приближение заданных точек функциями

Чаще всего ищется приближение в виде y=f(x).

Приближение заданных точек линейными и полиномиальными функциями

Приближение синуса полиномиальными функциями
Приближение точек, сгенерированных функцией синуса, полиномиальными функциями.
Красная прямая — многочлен первой степени, зелёная линия — второй степени, оранжевая линия — третьей степени , а синяя — четвёртой степени

Начнём приближения с многочлена первой степени:

Это прямая с угловым коэффициентом a. Через любые две точки можно провести прямую, так что многочлен первой степени проходит через любые две точки с различными абсциссами.

Если порядок уравнения увеличить до многочленов второй степени, получим:

Эта функция описывает параболу. Параболу можно провести через любые три точки.

Если увеличить порядок многочлена до третьей степени, получим:

Такую кривую можно построить для любых четырёх заданных точек. Более того, как правило, можно построить такую кривую, если заданы в точности четыре ограничения. Каждое ограничение может быть точкой, углом, или кривизной (которая является обратной величиной от радиуса соприкасающейся окружности). Ограничения в виде угла и кривизны чаще всего добавляется к концам кривой, и такие ограничения часто называются граничными условиями. Одинаковые граничные условия часто используются для обеспечения гладкого перехода между полиномиальными кривыми внутри сплайна. Ограничения более высокого порядка, такие как скорость изменения кривизны, также могут быть заданы. Это, например, может использоваться в строительстве автострадных развязок для расчёта скорости изменения сил, действующих на автомобиль (см. рывок) во время прохождения развязки и расчёта максимальной допустимой скорости.

Многочлен первой степени можно также построить, если заданы одна точка и угол, в то время как многочлен третьей степени может быть построен для двух точек, заданного углового коэффициента и заданной кривизны. Возможны другие комбинации ограничений для этих и более высоких степеней многочлена.

Если задано более чем n + 1 условий (где n — степень многочлена), можно всё равно попытаться построить полиномиальную кривую, удовлетворяющую этим условиям. Однако в общем случае такую кривую построить не удастся (например, многочлен первой степени по трём точкам можно построить, только если эти точки коллинеарны). Поэтому нужны некоторые методы осуществления приближения. Метод наименьших квадратов является одним из них.

Имеется несколько причин для получения приближённого решения, когда просто увеличение степени многочлена дало бы точное прохождение через точки:

  • Даже если точное решение существует, из этого не следует, что его легко найти. В некоторых алгоритмах можем получить расходящуюся последовательность и точное решение может оказаться невычислимым, в других случаях может потребоваться слишком много компьютерного времени для нахождения точного решения. В этих ситуациях может оказаться более приемлемым приближённое решение.
  • Эффект усреднения ненадёжных данных в выборке, может быть предпочтителен точному следованию точкам выборки, при котором может возникнуть выгибание кривой.
  • Феномен Рунге: при интерполяции полиномами высоких степеней может возникнуть эффект нежелательных осцилляций. Если кривая проходит через точки A и B, от кривой ожидается, что она пройдёт где-то близко к середине отрезка AB. Это может оказаться неверным в случае многочленов высокой степени — отклонение может быть очень велико. Для многочленов малой степени кривая, скорее всего, пройдёт рядом с серединой отрезка (а в случае многочлена первой степени обязательно пройдёт через середину).
  • Многочлены низкой степени обычно «гладкие», а многочлены высокой степени, как правило, «волнистые». Говоря точнее, максимальное число точек перегиба полиномиальной кривой равно n-2, где n — порядок многочлена. Точка перегиба — это точка, где кривизна кривой меняет знак. Заметим, что многочлены высокой степени не обязательно «волнистые», они также могут быть и «гладкими», но нет никаких гарантий «гладкости», в отличие от многочленов малого порядка. Многочлен десятой степени может иметь до восьми точек перегиба, но может иметь и меньше, или вообще не иметь.

Степень многочлена, более высокая, чем необходимо для точного прохождения кривой через точки, нежелательна по всем перечисленным выше причинам, но, кроме того, приводит к бесконечному числу решений. Например, многочлен первой степени (прямая) с ограничением одной точкой вместо обычных двух приводит к бесконечному числу решений. Отсюда вырастает проблема — как сравнивать и выбрать только одно решение, и это может стать проблемой как для программ, так и для людей. По этой причине лучший выбор — как можно меньшая степень для точного выполнения всех ограничений, и, может быть, даже меньшая степень, если приближённое решение допустимо.

Приближение заданных точек другими функциями

Могут быть также использованы в некоторых случаях и другие типы кривых, такие как тригонометрические функции (например, синус и косинус).

В спектроскопии данные могут быть приближены нормальным распределением, распределением Коши, контуром Фойгта[англ.] и связанными функциями.

Алгебраическое приближение и геометрическое приближение кривыми

Для алгебраического анализа данных, «приближение» обычно означает поиск кривой, минимизирующей вертикальное (по оси y) отклонение точки от кривой (к примеру, метод наименьших квадратов). Для графических приложений и при работе с изображениями геометрическое приближение ищет наилучшее визуальное приближение, что, обычно, означает попытку минимизировать расстояние до кривой (к примеру, метод наименьших полных квадратов) или минимизируются отклонения по обеим координатам. Геометрическое приближение непопулярно, поскольку, обычно, вовлекает нелинейные и/или рекуррентные вычисления, хотя оно и даёт эстетически более приемлемый и геометрически более аккуратный результат [16][17][18].

Приближение заданных точек плоскими кривыми

Если функция не может быть задана в виде , можно попытаться приблизить с помощью плоской кривой.

В некоторых случаях могут быть использованы и другие типы кривых, такие как конические сечения (дуги окружности, эллипса, параболы и гиперболы) или тригонометрические функции (такие как синус и косинус). Например, траектории объектов под воздействием силы тяжести являются параболами (если не принимать во внимание сопротивление воздуха). Таким образом, привязка (экспериментальных) точек траектории к параболической кривой будет иметь смысл. Приливы соответствуют синусоидальной модели, так что данные о приливах следует сопоставлять с синусом или суммой двух синусов различных периодов, если принимать во внимание влияние как луны, так и солнца.

В случае параметрической кривой эффективно рассматривать каждую координату как отдельную функцию от длины кривой. Если исходные данные можно упорядочить, можно использовать длину хорды[19].

Геометрическое приближение заданных точек окружностью

Приближение окружностью по методу Куупа (Coope). Точки соответствуют дуге окружности с центром в (1, 1) и радиусом 4.
Различные модели приближения эллипсом
Приближение эллипсом с минимизацией алгебраического расстояния (метод Фитцгиббона).

Куп[20] попытался решить задачу поиска наилучшего визуального приближения окружностью точек на плоскости. Метод элегантно преобразует нелинейную задачу в линейную, которую уже можно решить, не прибегая к рекурретным методам, а потому результат получается быстрее, чем при использовании предыдущих подходов.

Геометрическое приближение эллипсом

Вышеупомянутая техника была распространена на эллипсы общего вида[21] путём добавления нелинейного шага, в результате чего получен быстрый метод, находящий, однако, визуально привлекательные эллипсы произвольной ориентации и расположения.

Приложение к поверхностям

Заметим, что хотя до сих пор обсуждение относилось к плоским кривым, большинство результатов распространяется на поверхности в трёхмерном пространстве, каждый кусочек которых определяется сеткой кривых в двух параметрических направлениях. Поверхность может состоять из одного и более кусков в обоих направлениях.

Программы

Множество пакетов обработки статистических данных[англ.], таких как R, и пакетов программ численного анализа[англ.], таких как GNU Scientific Library, MLAB, DataMelt, Maple, MATLAB, SciPy и OpenOpt, включают средства приближения кривыми в различных сценариях. Существуют также программы, специально написанные для приближения кривыми. Их можно найти в статьях «Пакеты обработки статистических данных[англ.]» и «Пакеты программ численного анализа[англ.]».

См. также

Примечания

  1. Arlinghaus, 1994.
  2. Kolb, 1984.
  3. Halli, Rao, 1992, с. 165.
  4. Silver, 2012.
  5. Kiusalaas, 2005, с. 21.
  6. Guest, 2012, с. 349.
  7. См. также: Сглаживающий оператор
  8. Пакет Prism[англ.]кампании. Документация: «Fitting Models to Biological Data Using Linear and Nonlinear Regression» (Harvey Motulsky, Arthur Christopoulos).
  9. Freund, Wilson, Sa, 2006, с. 269.
  10. Daud, Sagayan, Yahya, Najwati, 2009, с. 689.
  11. Hauser, 2009, с. 227.
  12. Williams, 1976, с. 150.
  13. Salkind, 2010, с. 266.
  14. Klosterman, 1990, с. 1.
  15. Yoe, 1996, с. 69.
  16. Ahn, 2008.
  17. Chernov, Ma, 2011, с. 285–302.
  18. Liu, Wang, 2008, с. 384–397.
  19. Ahlberg, Nilson, Walsh, 1967, с. 51.
  20. Coope, 1993, с. 381.
  21. Sheer, 1997.

Литература

  • Sandra Lach Arlinghaus. Practical Handbook of Curve Fitting. — CRC Press, 1994. — ISBN 0849301434.
  • William M. Kolb. Curve Fitting for Programmable Calculators. — 3. — Syntec, Incorporated, 1984. — ISBN 0943494028.
  • John R. Hauser. Numerical Methods for Nonlinear Engineering Models. — Springer, 2009. — ISBN 978-1-4020-9919-9.
  • I.D. Coope. Circle fitting by linear and nonlinear least squares // Journal of Optimization Theory and Applications. — 1993. — Т. 76, вып. 2. — С. 381. — doi:10.1007/BF00939613.
  • Encyclopedia of Research Design / Neil J. Salkind. — SAGE Publications, 2010. — Т. 1. — ISBN 978-1-4129-6127-1.
  • Rudolf J. Freund, William J. Wilson, Ping Sa. Regression Analysis / Statistical Modeling of a Response Variable. — 2. — Elsevier, 2006. — ISBN 0-12-088597-2.
  • Jaan Kiusalaas. Numerical Methods in Engineering with MATLAB®. — Cambridge University Press, 2005. — ISBN 0-521-85288-9.
  • Richard E. Klosterman. Community Analysis and Planning Techniques. — Rowman & Littlefield Pub Inc, 1990. — ISBN 084767651X.
  • Hanita Daud, Vijanth Sagayan, Noorhana Yahya, Wan Najwati. Visual Informatics: Bridging Research and Practice (IVIC 2009) / Halimah Badioze Zaman, Peter Robinson, Maria Petrou, Patrick Olivier, Heiko Schröder, Timothy K. Shih. — Berlin, Heidelberg, NewYork: Sprintger, 2009. — Т. 5857. — (Lecture Notes in Computer Science). — ISBN 3-642-05035-2.
  • An Introduction to Risk and Uncertainty in the Evaluation of Environmental Investments / Charles E. Yoe. — West Chester, Pensylvania: The Greeley-Polhemus Group, Inc., 1996.
  • P. G. Guest. Numerical Methods of Curve Fitting. — Cambridge Academ, 2012. — ISBN 978-1-107646-5-7.
  • S.S. Halli, K.V. Rao. Advanced Techniques of Population Analysis. — 1992. — С. 165. — ISBN 0306439972.
  • Sung-Joon Ahn. Geometric Fitting of Parametric Curves and Surfaces // Journal of Information Processing Systems. — 2008. — Декабрь (vol. 4 (вып. 4). — doi:10.3745/JIPS.2008.4.4.153. Архивировано 13 марта 2014 года.
  • J. H. Ahlberg, E. N. Nilson, J. L. Walsh. The theory of splines and their applications. — New York, London: Academic Press, 1967.
  • N. Chernov, H. Ma. Computer Vision / Sota R. Yoshida. — Nova Science Publishers, 2011. — С. 285–302. — ISBN 9781612093994.
  • Nate Silver. The Signal and the Noise: Why So Many Predictions Fail-but Some Don't. — Penguin Group, 2012. — ISBN 978-1-59-420411-1.
  • Dudley Williams. Spectroscopy / Claire Marton.. — Academic Press, 1976. — Т. 13, Part 1. — (Methods of Experimental Physics). — ISBN 0124759130.
  • Yang Liu, Wenping Wang. Advances in Geometric Modeling and Processing / F. Chen, B. Juttler. — 2008. — Т. 4975. — С. 384–397. — ISBN 978-3-540-79245-1. — doi:10.1007/978-3-540-79246-8_29.
  • P. Sheer. A Software Assistantfor Manual Stereo Photometrology. — University of the Witwater-srand, 1997.

Литература для дальнейшего чтения

  • N. Chernov (2010), Circular and linear regression: Fitting circles and lines by least squares, Chapman & Hall/CRC, Monographs on Statistics and Applied Probability, Volume 117 (256 pp.). [1] Архивная копия от 22 октября 2020 на Wayback Machine

Read other articles:

العلاقات المارشالية الدومينيكية جزر مارشال دومينيكا   جزر مارشال   دومينيكا تعديل مصدري - تعديل   العلاقات المارشالية الدومينيكية هي العلاقات الثنائية التي تجمع بين جزر مارشال ودومينيكا.[1][2][3][4][5] مقارنة بين البلدين هذه مقارنة عامة ومرجع�...

 

1994 World JuniorChampionships in AthleticsTrack events100 mmenwomen200 mmenwomen400 mmenwomen800 mmenwomen1500 mmenwomen3000 mwomen5000 mmen10,000 mmenwomen100 m hurdleswomen110 m hurdlesmen400 m hurdlesmenwomen3000 msteeplechasemen4 × 100 m relaymenwomen4 × 400 m relaymenwomen5000 m walkwomen10,000 m walkmenRoad events20 km walkmenField eventsHigh jumpmenwomenPole vaultmenLong jumpmenwomenTriple jumpmenShot putmenwomenDiscus throwmenwomenHammer throwmenJavelin throwmenwomenCombined event...

 

Artikel ini mengenai Provinsi Batman di Turki. Untuk kegunaan lain silakan melihat: Batman (disambiguasi). Lokasi Provinsi Batman Batman adalah sebuah provinsi Turki yang terletak di sisi tenggara Anatolia. Sungai Batman mengalir melalui wilayah ini. Provinsi ini penting karena muatan reserve minyak buminya. Kilang minyak pertama didirikan pada tahun 1955. Dari Batman ada pipa minyak ke Iskenderun sepanjang 494 kilometer. Katun merupakan komoditas dagang utama pertanian di daerah ini. Sebuah ...

American actor Tony RobertsRoberts in 2019BornDavid Anthony Roberts (1939-10-22) October 22, 1939 (age 84)Manhattan, New York City, U.S.OccupationActorYears active1962–present From the original Broadway cast of Play It Again, Sam. Third from left is Tony Roberts; fourth from left is Woody Allen. Diane Keaton is on the far right. (1969) Tony Roberts (born David Anthony Roberts; October 22, 1939)[1] is an American actor. He is known for his roles in six Woody Allen movies—...

 

Yearly period of high rainfall Rainy season redirects here. For the film, see Wet Season (film). For other uses, see Rainy Season (disambiguation). Part of a series onWeather Temperate and polar seasons Winter Spring Summer Autumn Tropical seasons Dry season Harmattan Wet season Storms Cloud Cumulonimbus cloud Arcus cloud Downburst Microburst Heat burst Derecho Lightning Volcanic lightning Thunderstorm Air-mass thunderstorm Thundersnow Dry thunderstorm Mesocyclone Supercell Tornado Anticyclon...

 

Artikel ini sebatang kara, artinya tidak ada artikel lain yang memiliki pranala balik ke halaman ini.Bantulah menambah pranala ke artikel ini dari artikel yang berhubungan atau coba peralatan pencari pranala.Tag ini diberikan pada Januari 2023. Kepanembahan Cirebon1679–1775Ibu kotaKota CirebonBahasa yang umum digunakanBahasa Cirebon 1679-1773, Belanda 1679-1773Agama IslamPemerintahanMonarki (kepanembahan)Gusti Panembahan • 1679 (Deklarasi Pakungwati - didirikannya Peguron C...

Postseason college football bowl game College football game2021 Goodyear Cotton Bowl ClassicCollege Football Playoff Semifinal86th Cotton Bowl Classic Cincinnati Bearcats Alabama Crimson Tide (13–0) (12–1) AAC SEC 6 27 Head coach: Luke Fickell Head coach: Nick Saban APCoachesCFP 444 APCoachesCFP 111 1234 Total Cincinnati 3030 6 Alabama 710010 27 DateDecember 31, 2021Season2021StadiumAT&T StadiumLocationArlington, TexasMVPBrian Robinson Jr. (RB, Alabama)Will Anderso...

 

Olacaceae Olax imbricata Klasifikasi ilmiah Kerajaan: Plantae (tanpa takson): Tracheophyta (tanpa takson): Angiospermae (tanpa takson): Eudikotil (tanpa takson): Eudikotil inti Ordo: Santalales Famili: OlacaceaeR.Br.[1] Genera Lihat teks Olacaceae adalah salah satu suku anggota tumbuhan berbunga. Menurut Sistem klasifikasi APG II suku ini dimasukkan ke dalam bangsa Santalales, klad dikotil sejati inti (core Eudikotil) namun tidak termasuk ke dalam dua kelompok besar, Rosidae dan aste...

 

Romanian academic and politician (1883–1972) Onisifor GhibuBorn(1883-05-31)May 31, 1883Szelistye, Austria-HungaryDiedOctober 3, 1972(1972-10-03) (aged 89)Sibiu, Socialist Republic of RomaniaResting placeMunicipal Cemetery, SibiuAlma materUniversity of Bucharest Eötvös Loránd UniversityUniversity of JenaOccupation(s)Academic, educatorEmployerBabeș-Bolyai UniversitySpouse Veturia Ghibu [ro] ​ ​(m. 1911; died 1959)​Webs...

Junior role in the British government For the unpaid aide role, see Parliamentary Private Secretary to the Cabinet Office. United Kingdom Parliamentary Secretary for the Cabinet OfficeRoyal Arms of His Majesty's GovernmentIncumbentAlex Burghartsince 27 October 2022Cabinet OfficeStyleMinisterReports toChancellor of the Duchy of Lancaster Secretary of State in the Cabinet OfficeNominatorPrime Minister of the United KingdomAppointerThe Monarchon advice of the Prime MinisterTerm lengthAt His...

 

Libertarian socialist political and militant group in southern Mexico Zapatista Army of National LiberationEjército Zapatista de Liberación Nacional (EZLN)Flag of the EZLNAlso known asZapatistasLeaders Subcomandante Insurgente Galeano Comandanta Ramona Subcomandante Elisa FoundationNovember 17, 1983 (1983-11-17)Dates of operation1994–presentCountryMexicoHeadquartersSan Cristóbal de las Casas, ChiapasActive regionsChiapasIdeology Neozapatismo Indigenismo Anti-imperiali...

 

Legislation or statute defining public toilet access based on gender identity The examples and perspective in this article deal primarily with North America and do not represent a worldwide view of the subject. You may improve this article, discuss the issue on the talk page, or create a new article, as appropriate. (May 2018) (Learn how and when to remove this message) Part of a series onTransgender topics      OutlineHistoryTimeline Gender identities Androgyne Bissu...

The Fifth EstateSutradaraBill CondonProduser Steve Golin Brad Dorros Michael Sugar SkenarioJosh SingerBerdasarkanInside WikiLeaksoleh Daniel Domscheit-BergWikiLeaksoleh David Leigh dan Luke HardingPemeran Benedict Cumberbatch Daniel Brühl Anthony Mackie David Thewlis Alicia Vikander Stanley Tucci Laura Linney Penata musikCarter BurwellSinematograferTobias A. SchliesslerPenyuntingVirginia KatzPerusahaanproduksi DreamWorks Pictures Reliance Entertainment Participant Media Anonymous Conte...

 

74th United States Secretary of the Treasury (born 1946) Henry PaulsonOfficial portrait, 200674th United States Secretary of the TreasuryIn officeJuly 10, 2006 – January 20, 2009PresidentGeorge W. BushDeputyRobert M. KimmittPreceded byJohn W. SnowSucceeded byTim Geithner Personal detailsBornHenry Merritt Paulson Jr. (1946-03-28) March 28, 1946 (age 78)Palm Beach, Florida, U.S.Political partyRepublicanSpouseWendy JudgeChildren2, including MerrittEducationDartmouth College (BA)H...

 

Artikel ini tidak memiliki referensi atau sumber tepercaya sehingga isinya tidak bisa dipastikan. Tolong bantu perbaiki artikel ini dengan menambahkan referensi yang layak. Tulisan tanpa sumber dapat dipertanyakan dan dihapus sewaktu-waktu.Cari sumber: Universitas Dharmawangsa – berita · surat kabar · buku · cendekiawan · JSTOR Universitas DharmawangsaNama sebelumnyaSekolah Tinggi Ilmu Sosial & Ilmu Politik Dharmawangsa Sekolah Tinggi Hukum Dharmaw...

Este artículo trata sobre el sistema de transporte. Para la unidad de medida, véase Historia del metro (unidad de medida). Estación de King's Cross St. Pancras en el Metro de Londres Vista del Metro de Moscú Estación de Chamartín del Metro de Madrid En 1843, el inglés Charles Pearson propuso, como parte de un plan de mejora para la ciudad de Londres, abrir túneles subterráneos con vías férreas. En 1853, tras diez años de debates, el parlamento inglés autorizó la propuesta y en 1...

 

Brisbane BaylandsThis picture, facing west, shows the city of Brisbane at the south end (left side), which is sheltered by the prominent San Bruno Mountain rising in the background. The proposed Brisbane Baylands includes the land between Bayshore Boulevard and U.S. 101 (closer to the camera) north of the triangular Brisbane Lagoon.ProjectDeveloperUniversal Paragon CorporationWebsitebrisbanebaylands.comLocationPlaceCoordinates: 37°42′25″N 122°24′14″W / 37.707°N 122....

 

This article includes a list of references, related reading, or external links, but its sources remain unclear because it lacks inline citations. Please help improve this article by introducing more precise citations. (January 2011) (Learn how and when to remove this message) Record label Duophonic Ultra High Frequency Disks LimitedFounded1991FounderStereolabGenreVariousCountry of originUnited KingdomLocationLondon Duophonic Ultra High Frequency Disks Limited (also known as Duophonic Records...

Giordano Bruno ZilliInformazioni personaliArbitro di Calcio Federazione Italia SezioneReggio Emilia Attività nazionale AnniCampionatoRuolo 1940-19501942-1950Serie BSerie AArbitroArbitro Giordano Bruno Zilli (Finale Marina, 11 aprile 1909[1] – Reggio Emilia, 10 novembre 1981[2]) è stato un arbitro di calcio italiano. Carriera Nato a Finale Marina, oggi Finale Ligure, si trasferisce a Reggio Emilia dove inizia nel 1928[3] anche la carriera di arbitro nel locale ...

 

Anna d'AlverniaContessa d'Alvernia e di LauragaisStemma In carica1501 - 1524 PredecessoreGiovanni II SuccessoreCaterina de' Medici Contessa di LauragaisIn carica1501 - 1524 PredecessoreGiovanni II SuccessoreGiovanni Stewart Nome completoAnna de La Tour Altri titolisignora de La Tour Nascita1496 MorteSaint-Saturnin (Puy-de-Dôme), giugno 1524 Luogo di sepolturaabbazia di Bouschet-Vauluisant Yronde-et-Buron DinastiaCasato di La Tour d'Auvergne PadreGiovanni III de La Tour d'Auvergne Madre...