กลับไปหน้าบทความ

อ่าน 49 นาที

ฟังก์ชันเบสเซล

ฟังก์ชันเบสเซลเป็นคลาสของฟังก์ชันพิเศษที่มักปรากฏในปัญหาที่เกี่ยวข้องกับการเคลื่อนที่ของคลื่นการนำความร้อนและปรากฏการณ์ทางกายภาพอื่นๆ ที่มี สมมาตร แบบวงกลมหรือทรงกระบอก

ฟังก์ชันเบสเซล

ฟังก์ชันเบสเซลอธิบายส่วนรัศมีของการสั่นสะเทือนของเยื่อวงกลม

ฟังก์ชันเบสเซลเป็นคลาสของฟังก์ชันพิเศษที่มักปรากฏในปัญหาที่เกี่ยวข้องกับการเคลื่อนที่ของคลื่นการนำความร้อนและปรากฏการณ์ทางกายภาพอื่นๆ ที่มี สมมาตร แบบวงกลมหรือทรงกระบอก ฟังก์ชันเหล่านี้ตั้งชื่อตามนักดาราศาสตร์และนักคณิตศาสตร์ ชาวเยอรมัน ฟรีดริช เบสเซลผู้ซึ่งศึกษาฟังก์ชันเหล่านี้อย่างเป็นระบบในปี ค.ศ. 1824 [ 1 ]

ฟังก์ชันเบสเซลเป็นคำตอบของสมการเชิงอนุพันธ์สามัญ ชนิดหนึ่ง โดย เฉพาะ โดยที่เป็นจำนวนที่กำหนดรูปร่างของคำตอบ จำนวนนี้เรียกว่าอันดับของฟังก์ชันเบสเซล และสามารถเป็นจำนวนเชิงซ้อนใดๆ ก็ได้ แม้ว่าสมการเดียวกันจะเกิดขึ้นสำหรับทั้งและแต่นักคณิตศาสตร์กำหนดฟังก์ชันเบสเซลแยกกันสำหรับแต่ละ เพื่อให้แน่ใจว่าฟังก์ชันมีพฤติกรรมที่ราบรื่นเมื่ออันดับเปลี่ยนไป

กรณีที่สำคัญที่สุดคือเมื่อเป็นจำนวนเต็มหรือครึ่งจำนวนเต็ม เมื่อเป็นจำนวนเต็ม ฟังก์ชันเบสเซลที่ได้มักเรียกว่าฟังก์ชันทรงกระบอกหรือฮาร์มอนิกทรงกระบอกเนื่องจากเกิดขึ้นเองตามธรรมชาติเมื่อแก้ปัญหา (เช่น สมการลาปลาส) ในพิกัดทรงกระบอกเมื่อเป็นครึ่งจำนวนเต็ม คำตอบจะเรียกว่าฟังก์ชันเบสเซลทรงกลมและใช้ในระบบทรงกลม เช่น ในการแก้สมการเฮล์มโฮลทซ์ในพิกัดทรงกลม

แอปพลิเคชัน

สมการของเบสเซลเกิดขึ้นเมื่อต้องการหาคำตอบที่แยกตัวแปรได้ของสมการลาปลาสและสมการเฮล์มโฮลทซ์ในระบบพิกัดทรงกระบอกหรือทรงกลม ฟังก์ชันเบสเซลจึงมีความสำคัญอย่างยิ่งสำหรับปัญหาการแพร่กระจายคลื่นและศักย์สถิตหลายๆ ปัญหา ในการแก้ปัญหาในระบบพิกัดทรงกระบอก จะได้ฟังก์ชันเบสเซลที่มีอันดับจำนวนเต็ม ( ); ในปัญหาในระบบพิกัดทรงกลม จะได้ฟังก์ชันเบสเซลที่มีอันดับครึ่งจำนวนเต็ม ( ) ตัวอย่างเช่น:

ฟังก์ชันเบสเซลยังปรากฏในสาขาอื่นๆ เช่นการประมวลผลสัญญาณ (เช่นการสังเคราะห์เสียง FM , หน้าต่างไคเซอร์หรือตัวกรองเบสเซล ) นอกจากนี้ยังปรากฏในคณิตศาสตร์บริสุทธิ์ในฐานะส่วนหนึ่งของการขยายฟูริเยร์ของรูปแบบมาสส์

คำจำกัดความ

เนื่องจากสมการนี้เป็นสมการเชิงอนุพันธ์เชิงเส้นจึงสามารถปรับขนาดคำตอบให้มีขนาดใดก็ได้ ขนาดที่เลือกใช้สำหรับฟังก์ชันนั้นมาจากงานวิจัยในยุคแรกๆ ที่ฟังก์ชันเหล่านี้ปรากฏเป็นคำตอบของอินทิกรัลจำกัด แทนที่จะเป็นคำตอบของสมการเชิงอนุพันธ์ เนื่องจากสมการเชิงอนุพันธ์เป็นอันดับสอง จึงต้องมี คำตอบ ที่เป็นอิสระเชิงเส้น สอง คำตอบ คือ คำตอบแบบที่หนึ่งและคำตอบแบบที่สอง อย่างไรก็ตาม ขึ้นอยู่กับสถานการณ์ การกำหนดรูปแบบต่างๆ ของคำตอบเหล่านี้ก็มีความสะดวก ตารางด้านล่างสรุปรูปแบบต่างๆ ไว้ และจะอธิบายในส่วนต่อไป โดยทั่วไปจะใช้ตัวห้อยnแทนเมื่อทราบว่า เป็นจำนวนเต็ม

พิมพ์ชนิดแรกชนิดที่สอง
ฟังก์ชันเบสเซล เจαวายα
ฟังก์ชันเบสเซลที่ดัดแปลงฉันαเคα
ฟังก์ชันของแฮงเคลชม(1) α= J α + iY αชม(2) α= J αiY α
ฟังก์ชันเบสเซลทรงกลมเจเอ็น
ฟังก์ชันเบสเซลทรงกลมดัดแปลง ในk n
ฟังก์ชัน Hankel ทรงกลมชม.(1) = j n + iy nชม.(2) = j niy n

ฟังก์ชันเบสเซลชนิดที่สองและฟังก์ชันเบสเซลทรงกลมชนิดที่สองบางครั้งจะถูกแทนด้วยN n และ n n ตามลำดับ แทนที่จะเป็นY nและy n [ 4 ] [ 5 ]

ฟังก์ชันเบสเซลชนิดแรก: J α

กราฟแสดงฟังก์ชันเบสเซลชนิดที่หนึ่งสำหรับลำดับจำนวนเต็ม
กราฟแสดงฟังก์ชันเบสเซลชนิดแรกในระนาบจากถึง

ฟังก์ชันเบสเซลชนิดแรก ซึ่งแสดงด้วยJ α ( x )คือคำตอบของสมการเชิงอนุพันธ์ของเบสเซล สำหรับ  α ที่เป็นจำนวนเต็มหรือบวก ฟังก์ชันเบสเซลชนิดแรกจะมีค่าจำกัดที่จุดกำเนิด ( x = 0 ) ในขณะที่สำหรับ αที่เป็นลบและไม่ใช่จำนวนเต็ม  ฟังก์ชันเบสเซลชนิดแรกจะลู่เข้าสู่ค่าอนันต์เมื่อxเข้าใกล้ศูนย์ เป็นไปได้ที่จะกำหนดฟังก์ชันโดยใช้ชุดอนุกรมแมคลาลิน คูณด้วย (โปรดทราบว่าαไม่จำเป็นต้องเป็นจำนวนเต็ม และไม่อนุญาตให้ใช้กำลังที่ไม่ใช่จำนวนเต็มในชุดอนุกรมเทย์เลอร์) ซึ่งสามารถหาได้โดยการใช้วิธี Frobeniusกับสมการของเบสเซล: [ 6 ] โดยที่Γ( z )คือฟังก์ชันแกมมาซึ่งเป็นการขยายทั่วไปของ ฟังก์ชันแฟ กทอเรียลไปยังค่าที่ไม่ใช่จำนวนเต็ม ผู้เขียนบางคนก่อนหน้านี้กำหนดฟังก์ชันเบสเซลชนิดแรกแตกต่างกัน โดยพื้นฐานแล้วไม่มีการหารด้วยใน; [ 7 ]คำจำกัดความนี้ไม่ได้ใช้ในบทความนี้ ฟังก์ชันเบสเซลชนิดแรกเป็นฟังก์ชันสมบูรณ์ (entire function)ถ้าαเป็นจำนวนเต็ม มิฉะนั้นจะเป็นฟังก์ชันหลายค่า (multivalued function)ที่มีจุดเอกฐานที่ศูนย์ กราฟของฟังก์ชันเบสเซลจะมีลักษณะคล้ายกับ ฟังก์ชัน ไซน์หรือโคไซน์ที่แกว่งไปมาและลดลงตามสัดส่วนของ(ดูรูปแบบเชิงเส้นกำกับด้านล่างด้วย) แม้ว่ารากของฟังก์ชันโดยทั่วไปจะไม่เป็นคาบ ยกเว้นในเชิงเชิงเส้นกำกับสำหรับx ที่มีค่ามาก (อนุกรมนี้แสดงว่าJ 1 ( x )คืออนุพันธ์ของJ 0 ( x )เช่นเดียวกับ−sin xคืออนุพันธ์ของcos xโดยทั่วไปแล้ว อนุพันธ์ของJ n ( x )สามารถแสดงในรูปของJ n ± 1 ( x ) ได้ ด้วยเอกลักษณ์ด้านล่าง )

สำหรับ αที่ไม่ใช่จำนวนเต็มฟังก์ชันJ α ( x )และJ α ( x )เป็นอิสระเชิงเส้น และดังนั้นจึงเป็นสองคำตอบของสมการเชิงอนุพันธ์ ในทางกลับกัน สำหรับลำดับจำนวนเต็มnความสัมพันธ์ต่อไปนี้ใช้ได้ (ฟังก์ชันแกมมามีขั้วเดี่ยวที่จำนวนเต็มที่ไม่เป็นบวกแต่ละตัว): [ 8 ]

นั่นหมายความว่าคำตอบทั้งสองไม่เป็นอิสระเชิงเส้นอีกต่อไป ในกรณีนี้ คำตอบที่เป็นอิสระเชิงเส้นตัวที่สองก็คือฟังก์ชันเบสเซลชนิดที่สอง ดังที่จะกล่าวถึงต่อไป

อินทิกรัลของเบสเซล

นิยามอื่นของฟังก์ชันเบสเซลสำหรับค่าจำนวนเต็มของnสามารถทำได้โดยใช้การแสดงอินทิกรัล: [ 9 ] ซึ่งเรียกอีกอย่างว่าสูตร Hansen-Bessel [ 10 ]

นี่คือแนวทางที่เบสเซลใช้[ 11 ]และจากคำจำกัดความนี้ เขาได้อนุมานคุณสมบัติหลายประการของฟังก์ชัน คำจำกัดความนี้สามารถขยายไปยังลำดับที่ไม่ใช่จำนวนเต็มโดยใช้อินทิกรัลของ Schläfli ตัวใดตัวหนึ่ง สำหรับRe( x ) > 0 : [ 9 ] [ 12 ] [ 13 ] [ 14 ] [ 15 ]

ความสัมพันธ์กับอนุกรมไฮเปอร์จีโอเมตริก

ฟังก์ชันเบสเซลสามารถแสดงได้ในรูปของอนุกรมไฮเปอร์จีโอเมตริกทั่วไปดังนี้[ 16 ]

นิพจน์นี้เกี่ยวข้องกับการพัฒนาฟังก์ชันเบสเซลในแง่ของฟังก์ชันเบสเซล-คลิฟฟอร์

ความสัมพันธ์กับพหุนาม Laguerre

ในแง่ของพหุนาม Laguerre L kและพารามิเตอร์t ที่เลือกโดยพล การ ฟังก์ชัน Bessel สามารถแสดงได้ดังนี้[ 17 ]

ฟังก์ชันเบสเซลชนิดที่สอง: Y α

กราฟแสดงฟังก์ชันเบสเซลชนิดที่สองสำหรับลำดับจำนวนเต็ม

ฟังก์ชันเบสเซลชนิดที่สอง ซึ่งแสดงด้วยY α ( x )บางครั้งอาจแสดงด้วยN α ( x ) แทน เป็นคำตอบของสมการเชิงอนุพันธ์เบสเซลที่มีจุดเอกฐานที่จุดกำเนิด ( x = 0 ) และมีหลายค่าบางครั้งเรียกว่าฟังก์ชันเวเบอร์ ตามที่ HM Weber  ( 1873 ) แนะนำและเรียกอีกอย่างว่าฟังก์ชันนอยมันน์ตามชื่อของCarl Neumann [ 18 ]

สำหรับαที่ไม่ใช่จำนวนเต็มY α ( x )มีความสัมพันธ์กับJ α ( x )โดย

ในกรณีที่อันดับn เป็นจำนวนเต็ม ฟังก์ชันจะถูกกำหนดโดยการหาลิมิตเมื่อα ซึ่งไม่ใช่จำนวนเต็ม เข้าใกล้n :

ถ้าnเป็นจำนวนเต็มที่ไม่เป็นลบ เราจะมีอนุกรม[ 19 ] โดยที่คือฟังก์ชันไดแกมมาซึ่งเป็นอนุพันธ์ลอการิทึมของฟังก์ชันแกมมา[ 5 ]

นอกจากนี้ยังมีสูตรอินทิกรัลที่สอดคล้องกัน (สำหรับRe( x ) > 0 ): [ 20 ]

ในกรณีที่n = 0 : (โดยที่ เป็นค่าคงที่ของออยเลอร์ )

กราฟแสดงฟังก์ชันเบสเซลชนิดที่สองในระนาบเชิงซ้อนตั้งแต่ถึง

Y α ( x )มีความจำเป็นในฐานะที่เป็นคำตอบอิสระเชิงเส้นลำดับที่สองของสมการเบสเซลเมื่อ αเป็นจำนวนเต็ม แต่ Y α ( x )มีความหมายมากกว่านั้น มันสามารถพิจารณาได้ว่าเป็นคู่ "ตามธรรมชาติ" ของ J α ( x )ดูเพิ่มเติมในหัวข้อย่อยเกี่ยวกับฟังก์ชันแฮงเคลด้านล่าง

ยิ่งไปกว่านั้น เมื่อαเป็นจำนวนเต็ม ความสัมพันธ์ต่อไปนี้ก็เป็นจริงเช่นเดียวกับกรณีของฟังก์ชันประเภทแรก:

ทั้งJ α ( x )และY α ( x )เป็นฟังก์ชันโฮโลมอร์ฟิกของxบนระนาบเชิงซ้อนที่ตัดตามแกนจริงลบ เมื่อαเป็นจำนวนเต็ม ฟังก์ชันเบสเซลJจะเป็นฟังก์ชันเอนไทร์ของxถ้าxถูกกำหนดให้มีค่าที่ไม่เป็นศูนย์ ฟังก์ชันเบสเซลจะเป็นฟังก์ชันเอนไทร์ของ α

ฟังก์ชันเบสเซลชนิดที่สอง เมื่อαเป็นจำนวนเต็ม เป็นตัวอย่างหนึ่งของคำตอบชนิดที่สองในทฤษฎีบทของฟุคส์

ฟังก์ชัน Hankel: H(1) α, ชม(2) α

กราฟแสดงฟังก์ชันแฮงเคลชนิดที่หนึ่งH(1) ( x )โดยที่n = −0.5ในระนาบเชิงซ้อนจาก−2 − 2 iถึง2 + 2 i
กราฟแสดงฟังก์ชันแฮงเคลชนิดที่สองH(2) ( x )โดยที่n = −0.5ในระนาบเชิงซ้อนจาก−2 − 2 iถึง2 + 2 i

อีกหนึ่งรูปแบบที่สำคัญของคำตอบเชิงเส้นสองคำตอบของสมการเบสเซลคือฟังก์ชันแฮงเคลชนิดที่หนึ่งและชนิดที่สอง H(1) α( x )และH(2) α( x )ถูกกำหนดเป็น[ 21 ] โดยที่iคือหน่วยจินตนาการการรวมเชิงเส้นเหล่านี้ยังเป็นที่รู้จักในชื่อฟังก์ชันเบสเซลชนิดที่สามซึ่งเป็นสองคำตอบที่เป็นอิสระเชิงเส้นของสมการเชิงอนุพันธ์ของเบสเซล ตั้งชื่อตามเฮอร์มันน์ ฮันเคิ

รูปแบบการรวมเชิงเส้น เหล่านี้ มีคุณสมบัติที่ดูเรียบง่ายมากมาย เช่น สูตรเชิงอะซิมโทติกหรือการแสดงแทนเชิงอินทิกรัล ในที่นี้ "เรียบง่าย" หมายถึงการปรากฏของตัวประกอบในรูปแบบe i f (x)สำหรับค่าจริงโดยที่ และเป็นค่าจริง ฟังก์ชันเบสเซลชนิดที่หนึ่งและชนิดที่สองคือส่วนจริงและส่วนจินตนาการตามลำดับของฟังก์ชันแฮงเคลชนิดแรก และส่วนจริงและส่วนจินตนาการลบของฟังก์ชันแฮงเคลชนิดที่สอง ดังนั้น สูตรข้างต้นจึงเป็นอนาล็อกของสูตรของออยเลอร์โดยแทนที่H(1) α( x ) , H(2) α( x )สำหรับและสำหรับ, , ดังที่แสดงไว้อย่างชัดเจนในการขยายอนุกรมเชิงอะซิมโทติก

ฟังก์ชัน Hankel ใช้ในการแสดงคำตอบของสมการคลื่นทรงกระบอกที่แผ่ออกไปด้านนอกและด้านในตามลำดับ (หรือในทางกลับกัน ขึ้นอยู่กับข้อกำหนดเรื่องเครื่องหมายสำหรับความถี่ )

โดยใช้ความสัมพันธ์ก่อนหน้านี้ สามารถแสดงออกมาได้ดังนี้

ถ้าαเป็นจำนวนเต็ม จะต้องคำนวณลิมิต ความสัมพันธ์ต่อไปนี้ใช้ได้ ไม่ว่าαจะเป็นจำนวนเต็มหรือไม่ก็ตาม: [ 22 ]

โดยเฉพาะอย่างยิ่ง ถ้าα = m + 1/2โดยที่ mเป็นจำนวนเต็มที่ไม่เป็นลบ ความสัมพันธ์ข้างต้นแสดงให้เห็นโดยตรงว่า

สิ่งเหล่านี้มีประโยชน์ในการพัฒนาฟังก์ชันเบสเซลทรงกลม (ดูด้านล่าง)

ฟังก์ชัน Hankel ยอมรับการแสดงอินทิกรัลต่อไปนี้สำหรับRe( x ) > 0 : [ 23 ] โดยที่ขอบเขตการอินทิเกรตบ่งชี้ถึงการอินทิเกรตตามเส้นโค้งที่สามารถเลือกได้ดังนี้: จาก−∞ถึง 0 ตามแกนจริงลบ จาก 0 ถึง± π iตามแกนจินตนาการ และจาก± π iถึง+∞ ± π iตามเส้นโค้งที่ขนานกับแกนจริง[ 20 ]

ฟังก์ชันเบสเซลที่ดัดแปลง: I α , K α

ฟังก์ชันเบสเซลยังคงใช้ได้แม้กับอาร์กิวเมนต์เชิงซ้อนxและกรณีพิเศษที่สำคัญคือกรณีของอาร์กิวเมนต์จินตภาพล้วน ในกรณีนี้ คำตอบของสมการเบสเซลเรียกว่าฟังก์ชันเบสเซลแบบดัดแปลง (หรือบางครั้ง เรียกว่า ฟังก์ชันเบสเซลไฮเปอร์โบลิก ) ชนิดที่หนึ่งและชนิดที่สองและกำหนดเป็น[ 24 ] เมื่อαไม่ใช่จำนวนเต็ม เมื่อαเป็นจำนวนเต็ม จะใช้ลิมิต ฟังก์ชันเหล่านี้ถูกเลือกให้เป็นค่าจริงสำหรับอาร์กิวเมนต์จริงและบวกx ดังนั้น การขยายอนุกรมสำหรับI α ( x )จึงคล้ายกับการขยายอนุกรมสำหรับJ α ( x )แต่ไม่มีตัวประกอบ สลับ (−1) m

สามารถแสดงได้ในรูปของฟังก์ชัน Hankel:

เมื่อใช้สูตรทั้งสองนี้ ผลลัพธ์ของ ซึ่งโดยทั่วไปเรียกว่า อินทิกรัลของนิโคลสัน หรือ สูตรของนิโคลสัน สามารถหาได้ดังนี้

โดยที่เงื่อนไขRe( x ) > 0เป็นไปตามที่กำหนด นอกจากนี้ยังสามารถแสดงได้ว่า เฉพาะเมื่อ| Re( α ) | < 1/2และ Re ( x ) ≥ 0แต่ไม่ใช่เมื่อ x = 0 [ 25 ]

เราสามารถแสดงฟังก์ชันเบสเซลอันดับแรกและอันดับสองในรูปของฟังก์ชันเบสเซลที่ดัดแปลงแล้วได้ (ซึ่งใช้ได้ก็ต่อเมื่อπ < arg z⁠)π/2) : [ 26 ]

I α ( x )และ K α ( x )เป็นคำตอบอิสระเชิงเส้นสองคำตอบของสมการเบสเซลที่แก้ไขแล้ว : [ 27 ]

แตกต่างจากฟังก์ชันเบสเซลทั่วไป ซึ่งเป็น ฟังก์ชันที่แกว่งไปมาตามจำนวนจริง ฟังก์ชันและเป็นฟังก์ชันที่เพิ่มขึ้นและลดลง แบบเอกซ์ponential ตามลำดับ เช่นเดียวกับฟังก์ชันเบสเซลทั่วไปฟังก์ชันจะมีค่าเป็นศูนย์ที่x = 0สำหรับα > 0และมีค่าจำกัดที่x = 0สำหรับα = 0ในทำนองเดียวกันจะ มีค่าลู่เข้าสู่ ค่าอนันต์ที่x = 0โดยมีจุดเอกฐานเป็นแบบลอการิทึมสำหรับ> 0และ1/2 Γ(| α |)(2/ x ) | α |มิฉะนั้น [ 28 ]

ฟังก์ชันเบสเซลแบบดัดแปลงชนิดแรกสำหรับ
ฟังก์ชันเบสเซลแบบดัดแปลงชนิดที่สองสำหรับ

สูตรอินทิกรัลสองสูตรสำหรับฟังก์ชันเบสเซลที่ดัดแปลงคือ (สำหรับRe( x ) > 0 ): [ 29 ]

ฟังก์ชันเบสเซลสามารถอธิบายได้ว่าเป็นการแปลงฟูริเยร์ของกำลังของฟังก์ชันกำลังสอง ตัวอย่างเช่น (สำหรับRe(ω) > 0 ):

สามารถพิสูจน์ได้โดยการแสดงความเท่ากับนิยามอินทิกรัลข้างต้นสำหรับK 0ซึ่งทำได้โดยการอินทิเกรตเส้นโค้งปิดในควอดแรนต์แรกของระนาบเชิงซ้อน

ฟังก์ชันเบสเซลที่ดัดแปลงชนิดที่สองอาจแสดงด้วยอินทิกรัลของ Bassett [ 30 ]

ฟังก์ชันเบสเซลที่ดัดแปลงK 1/3และK 2/3สามารถแสดงได้ในรูปของอินทิกรัลที่ลู่เข้าอย่างรวดเร็ว[ 31 ]

ฟังก์ชันเบสเซลที่ปรับปรุงแล้วมีประโยชน์ในการแสดงการแจกแจงลาปลาซเป็นส่วนผสมของการแจกแจงปกติในระดับเลขชี้กำลัง

ฟังก์ชันเบสเซลแบบดัดแปลงชนิดที่สองยังถูกเรียกด้วยชื่ออื่นๆ ดังต่อไปนี้ (ซึ่งปัจจุบันไม่ค่อยใช้แล้ว):

ฟังก์ชัน Bessel ทรงกลม: j n , y n

กราฟแสดงฟังก์ชันเบสเซลทรงกลมชนิดแรกj n ( z )โดยที่n = 0.5ในระนาบเชิงซ้อนจาก−2 − 2 iถึง2 + 2 i
กราฟแสดงฟังก์ชันเบสเซลทรงกลมชนิดที่สองy n ( z )โดยที่n = 0.5ในระนาบเชิงซ้อน จาก−2 − 2 iถึง2 + 2 i
ฟังก์ชันเบสเซลทรงกลมชนิดแรกสำหรับ
ฟังก์ชันเบสเซลทรงกลมชนิดที่สองสำหรับ

เมื่อแก้สมการเฮล์มโฮลทซ์ในพิกัดทรงกลมโดยใช้วิธีการแยกตัวแปร สมการเชิงรัศมีจะมีรูปแบบดังนี้

คำตอบสองคำตอบที่เป็นอิสระเชิงเส้นของสมการนี้เรียกว่าฟังก์ชันเบสเซลทรงกลมj nและy nและมีความสัมพันธ์กับฟังก์ชันเบสเซลธรรมดาJ nและY nโดย[ 33 ]

y nยังถูกแทนด้วยn nหรือη n ; ผู้เขียนบางคนเรียกฟังก์ชันเหล่านี้ว่าฟังก์ชัน นอยมันน์ทรงกลม

จากความสัมพันธ์กับฟังก์ชันเบสเซลทั่วไป จะเห็นได้โดยตรงว่า:

ฟังก์ชันเบสเซลทรงกลมสามารถเขียนได้อีกแบบหนึ่งว่า (สูตรของเรย์ลีย์ ) [ 34 ]

ฟังก์ชันเบสเซลทรงกลมลำดับที่ศูนย์j 0 ( x )เรียกอีกอย่างว่าฟังก์ชัน sinc (ที่ไม่ได้รับการทำให้เป็นมาตรฐาน) ฟังก์ชันเบสเซลทรงกลมแรกๆ ได้แก่: [ 35 ] และ[ 36 ]

รากที่ไม่เป็นศูนย์กลุ่มแรกๆ ของฟังก์ชันเบสเซลทรงกลมกลุ่มแรกๆ มีดังนี้:

รากที่ไม่เป็นศูนย์ของฟังก์ชันเบสเซลทรงกลม (ชนิดที่หนึ่ง)
คำสั่งรากที่ 1รากที่ 2รากที่ 3รูท 4รากที่ 5
3.1415936.2831859.42477812.56637115.707963
4.4934097.72525210.90412214.06619417.220755
5.7634599.09501112.32294115.51460318.689036
6.98793210.41711913.69802316.92362120.121806
8.18256111.70490715.03966518.30125621.525418
รากที่ไม่เป็นศูนย์ของฟังก์ชันเบสเซลทรงกลม (ชนิดที่สอง)
คำสั่งรากที่ 1รากที่ 2รากที่ 3รูท 4รากที่ 5
1.5707964.7123897.85398210.99557414.137167
2.7983866.1212509.31786612.48645415.644128
3.9595287.45161010.71564713.92168617.103359
5.0884988.73371012.06754415.31539018.525210
6.1978319.98246613.38528716.67662519.916796

ฟังก์ชันการสร้าง

ฟังก์ชันเบสเซลทรงกลมมีฟังก์ชันก่อกำเนิด[ 37 ]

การขยายอนุกรมจำกัด

ตรงกันข้ามกับฟังก์ชันเบสเซลจำนวนเต็มทั้งหมดJ n ( x ), Y n ( x )ฟังก์ชันเบสเซลทรงกลมj n ( x ), y n ( x )มีการแสดงออกเป็นอนุกรมจำกัด: [ 38 ]

ความสัมพันธ์เชิงอนุพันธ์

ต่อไปนี้f nคือj n , y n , h ใดๆ ก็ได้(1) , ชม(2) สำหรับn = 0, ±1, ±2, ... [ 39 ]

ฟังก์ชัน Hankel ทรงกลม: h(1) , ชม(2)

กราฟแสดงฟังก์ชันแฮงเคลทรงกลมชนิดที่หนึ่งh(1) ( x )โดยที่n = −0.5ในระนาบเชิงซ้อนจาก−2 − 2 iถึง2 + 2 i
กราฟแสดงฟังก์ชันแฮงเคลทรงกลมชนิดที่สองh(2) ( x )โดยที่n = −0.5ในระนาบเชิงซ้อนจาก−2 − 2 iถึง2 + 2 i

นอกจากนี้ยังมีฟังก์ชันแฮงเคล ในรูปแบบทรงกลมอีกด้วย :

มีสูตรสำเร็จรูปง่ายๆ สำหรับฟังก์ชันเบสเซล อันดับ ครึ่งจำนวนเต็มในรูปของฟังก์ชันตรีโกณมิติ มาตรฐาน และด้วยเหตุนี้จึงใช้ได้กับฟังก์ชันเบสเซลทรงกลมด้วย โดยเฉพาะอย่างยิ่งสำหรับจำนวนเต็มที่ไม่เป็นลบn : และh(2) คือค่าสังยุคเชิงซ้อนของสิ่งนี้ (สำหรับx ที่เป็นจำนวนจริง ) ดังนั้นจึงสรุปได้ว่าj 0 ( x ) = ไซน์x/xและ y 0 ( x ) = − cos x/xและอื่นๆ

ฟังก์ชันแฮงเคลทรงกลมปรากฏในปัญหาที่เกี่ยวข้องกับ การแพร่กระจาย คลื่นทรงกลมตัวอย่างเช่น ในการขยายแบบมัลติโพลของสนามแม่เหล็กไฟฟ้า

ฟังก์ชันริกคาติ–เบสเซล: S n , C n , ξ n , ζ n

ฟังก์ชัน Riccati -Bessel แตกต่างจากฟังก์ชัน Bessel ทรงกลมเพียงเล็กน้อย:

ฟังก์ชัน Riccati–Bessel Sn พล็อตเชิงซ้อนจาก -2-2i ถึง 2+2i
ฟังก์ชัน Riccati–Bessel Sn พล็อตเชิงซ้อนจาก −2 − 2 iถึง 2 + 2 i

พวกมันสอดคล้องกับสมการเชิงอนุพันธ์

ตัวอย่างเช่น สมการเชิงอนุพันธ์ประเภทนี้ปรากฏในกลศาสตร์ควอนตัมขณะแก้ส่วนประกอบรัศมีของสมการชโรดิงเกอร์ด้วยกำแพงศักย์อนันต์ทรงกระบอกสมมุติ[ 40 ]สมการเชิงอนุพันธ์นี้และคำตอบของ Riccati–Bessel ยังเกิดขึ้นในปัญหาการกระเจิงของคลื่นแม่เหล็กไฟฟ้าโดยทรงกลม ซึ่งรู้จักกันในชื่อการกระเจิงของ Mieตามคำตอบที่ตีพิมพ์ครั้งแรกโดย Mie (1908) ดูเช่น Du (2004) [ 41 ]สำหรับการพัฒนาล่าสุดและเอกสารอ้างอิง

ตามแนวคิดของDebye (1909) บางครั้งจะ ใช้ สัญลักษณ์ψ n , χ nแทนS n , C n

รูปแบบเชิงเส้นกำกับ

ฟังก์ชันเบสเซลมี รูปแบบ เชิงเส้นกำกับ ดังต่อไปนี้ สำหรับอาร์กิวเมนต์ขนาดเล็กจะได้เมื่อไม่ใช่จำนวนเต็มลบ: [ 6 ]

เมื่อαเป็นจำนวนเต็มลบ เราจะได้ว่า

สำหรับฟังก์ชันเบสเซลชนิดที่สอง เรามีสามกรณี โดยที่γคือค่าคงที่ออยเลอร์-มาสเชโรนี (0.5772...) สำหรับกรณีที่สอง (โดยที่เป็นจำนวนเต็มบวก) เทอมหนึ่งจะเด่นกว่าอีกเทอมหนึ่ง เว้นแต่ว่า จะเป็นจำนวนจินตนาการ

สำหรับอาร์กิวเมนต์จริงขนาดใหญ่z ≫ | α 21/4 |เราไม่สามารถเขียนรูปแบบเชิงเส้นกำกับที่แท้จริงสำหรับฟังก์ชันเบสเซลชนิดที่หนึ่งและชนิดที่สองได้ (เว้นแต่αจะเป็นจำนวนเต็มครึ่ง) เนื่องจากมีศูนย์ไปจนถึงอนันต์ ซึ่งจะต้องตรงกับการขยายเชิงเส้นกำกับใดๆ อย่างแม่นยำ อย่างไรก็ตาม สำหรับค่าarg z ที่กำหนด เราสามารถเขียนสมการที่มีพจน์อันดับ| z | −1 ได้ : [ 42 ]

(สำหรับα = 1/2(พจน์สุดท้ายในสูตรเหล่านี้จะหายไปโดยสิ้นเชิง ดูได้จากฟังก์ชันเบสเซลทรงกลมด้านบน)

รูปแบบเชิงอะซิมโทติกของฟังก์ชัน Hankel มีดังนี้:

สามารถขยายแนวคิดเหล่านี้ไปยังค่าarg z อื่นๆ ได้ โดยใช้สมการที่เชื่อมโยงH(1) α( ze im π )และH(2) α( ze im π )ถึงH(1) α( z )และH(2) α( z ) . [ 43 ]

เป็นที่น่าสนใจว่าถึงแม้ฟังก์ชันเบสเซลชนิดแรกจะเป็นค่าเฉลี่ยของฟังก์ชันแฮงเคลสองฟังก์ชัน แต่J α ( z )ไม่ได้เป็นค่าประมาณเชิงเส้นกำกับของค่าเฉลี่ยของรูปแบบเชิงเส้นกำกับทั้งสองนี้เมื่อzเป็นค่าลบ (เพราะค่าใดค่าหนึ่งจะไม่ถูกต้อง ขึ้นอยู่กับค่า arg zที่ใช้) แต่รูปแบบเชิงเส้นกำกับสำหรับฟังก์ชันแฮงเคลทำให้เราสามารถเขียนรูปแบบเชิงเส้นกำกับสำหรับฟังก์ชันเบสเซลชนิดแรกและชนิดที่สองสำหรับz ที่เป็นจำนวนเชิงซ้อน (ไม่ใช่จำนวนจริง) ตราบใดที่| z |เข้าสู่ค่าอนันต์ที่มุมเฟสคงที่arg z (โดยใช้รากที่สองที่มีส่วนจริงเป็นบวก):

สำหรับฟังก์ชันเบสเซลที่ดัดแปลงHankelได้พัฒนาการขยายอนุกรมอสิมโทติกด้วยเช่นกัน: [ 44 ] [ 45 ]

นอกจากนี้ยังมีรูปแบบเชิงอะซิมโทติก (สำหรับจำนวนจริงขนาดใหญ่) [ 46 ]

เมื่อα = 1/2เงื่อนไขทั้งหมด ยกเว้นเงื่อนไขแรก จะหายไป และเราจะได้

สำหรับข้อโต้แย้งเล็กๆ น้อยๆเรามี

คุณสมบัติ

สำหรับฟังก์ชันเบสเซลใดๆ ที่มีอันดับไม่เป็นจำนวนเต็มลบ อนุพันธ์ของฟังก์ชันสามารถกำหนดได้ดังนี้: [ 47 ]

หรือเทียบเท่ากัน

สูตรเหล่านี้สามารถใช้เพื่อกำหนดความสัมพันธ์เวียนเกิดสำหรับซึ่งมีรูปแบบทั่วไปมากกว่าที่ระบุไว้ด้านล่าง[ 47 ]

สำหรับลำดับจำนวนเต็มα = nนั้นJ nมักถูกกำหนดโดยใช้ชุดอนุกรมลอเรนต์สำหรับฟังก์ชันก่อกำเนิดซึ่ง เป็นแนวทางที่ PA Hansen ใช้ในปี 1843 (สามารถขยายไปสู่ลำดับที่ไม่ใช่จำนวนเต็มได้โดยใช้การอินทิเกรตตามเส้นโค้งหรือวิธีการอื่น ๆ)

อนุกรมอนันต์ของฟังก์ชันเบสเซลในรูปแบบที่เกิดขึ้นในระบบทางกายภาพหลายระบบและถูกกำหนดในรูปแบบปิดโดยอนุกรมซุง[ 48 ]ตัวอย่างเช่น เมื่อ N = 3: โดยทั่วไปแล้ว อนุกรมซุงและอนุกรมซุงสลับจะเขียนได้ดังนี้:

การขยายอนุกรมโดยใช้ฟังก์ชันเบสเซล ( อนุกรมแคปเทน ) คือ

ความสัมพันธ์ที่สำคัญอีกประการหนึ่งสำหรับลำดับจำนวนเต็มคือการขยายแบบ Jacobi–Anger : และ

วิธีหลังนี้เทียบเท่ากับ วิธีที่ใช้ในการขยายคลื่นระนาบให้เป็นผลรวมของคลื่นทรงกระบอกหรือใช้ในการหาอนุกรมฟูริเยร์ ของ สัญญาณ FMที่มีการปรับโทนเสียง

โดยทั่วไปแล้ว อนุกรมนี้ เรียกว่าการขยาย Neumann ของfสัมประสิทธิ์สำหรับν = 0มีรูปแบบที่ชัดเจน โดยที่O kคือพหุนาม Neumann [ 49 ]

ฟังก์ชันที่เลือกไว้จะยอมรับการแสดงแทนแบบพิเศษ เนื่องจาก ความสัมพันธ์เชิงตั้งฉาก

โดยทั่วไปแล้ว ถ้าfมีจุดแยกสาขาใกล้จุดกำเนิดในลักษณะที่ จากนั้น หรือ โดยที่คือการแปลงลาปลาของf [ 50 ]

อีกวิธีหนึ่งในการกำหนดฟังก์ชันเบสเซลคือสูตรการแสดงแทนแบบปัวซงและสูตรเมห์เลอร์-โซนีน: โดยที่ν > − 1/2และ zC [ 51 ] สูตร นี้ มีประโยชน์โดยเฉพาะเมื่อทำงานกับการแปลงฟูริเยร์

Because Bessel's equation becomes Hermitian (self-adjoint) if it is divided by x, the solutions must satisfy an orthogonality relationship for appropriate boundary conditions. In particular, it follows that: where α > −1, δm,n is the Kronecker delta, and uα,m is the mth zero of Jα(x). This orthogonality relation can then be used to extract the coefficients in the Fourier–Bessel series, where a function is expanded in the basis of the functions Jα(xuα,m) for fixed α and varying m.

An analogous relationship for the spherical Bessel functions follows immediately:

If one defines a boxcar function of x that depends on a small parameter ε as: (where rect is the rectangle function) then the Hankel transform of it (of any given order α > −1/2), gε(k), approaches Jα(k) as ε approaches zero, for any given k. Conversely, the Hankel transform (of the same order) of gε(k) is fε(x): which is zero everywhere except near 1. As ε approaches zero, the right-hand side approaches δ(x − 1), where δ is the Dirac delta function. This admits the limit (in the distributional sense):

A change of variables then yields the closure equation:[52] for α > −1/2. For the spherical Bessel functions the orthogonality relation is: for α > −1.

คุณสมบัติสำคัญอีกประการหนึ่งของสมการเบสเซล ซึ่งเป็นผลมาจากเอกลักษณ์ของอาเบลเกี่ยวข้องกับวรอนสเกียนของคำตอบ โดยที่A αและB αเป็นคำตอบสองคำตอบใดๆ ของสมการเบสเซล และC αเป็นค่าคงที่ที่ไม่ขึ้นกับx (ซึ่งขึ้นอยู่กับ α และฟังก์ชันเบสเซลเฉพาะที่พิจารณา) โดยเฉพาะอย่างยิ่ง สำหรับ α > −1

สำหรับα > −1ฟังก์ชันเอไรเซชันคู่ของจีนัส 1, x α J α ( x )จะมีเฉพาะศูนย์ที่เป็นจำนวนจริงเท่านั้น ให้ เป็นศูนย์บวกทั้งหมดของฟังก์ชันนี้ แล้ว

(ยังมีอินทิกรัลและเอกลักษณ์อื่นๆ อีกจำนวนมากที่เป็นที่รู้จัก แต่ไม่ได้แสดงไว้ในที่นี้ อย่างไรก็ตาม สามารถค้นหาได้จากเอกสารอ้างอิง)

ความสัมพันธ์เวียนเกิด

ฟังก์ชันJ α , Y α , H(1) αและH(2) αทั้งหมดเป็นไปตามความสัมพันธ์เวียนเกิด[ 53 ] และ โดยที่ZหมายถึงJ , Y , H (1)หรือH (2)เอกลักษณ์ทั้งสองนี้มักจะรวมกัน เช่น บวกหรือลบ เพื่อให้ได้ความสัมพันธ์อื่นๆ ที่หลากหลาย ด้วยวิธีนี้ ตัวอย่างเช่น เราสามารถคำนวณฟังก์ชันเบสเซลที่มีลำดับสูงกว่า (หรืออนุพันธ์ที่สูงกว่า) โดยกำหนดค่าที่ลำดับต่ำกว่า (หรืออนุพันธ์ที่ต่ำกว่า) โดยเฉพาะอย่างยิ่ง เป็นไปตาม[ 54 ]

โดยใช้ความสัมพันธ์ก่อนหน้านี้ เราสามารถหาความสัมพันธ์ที่คล้ายคลึงกันสำหรับ ฟังก์ชันเบสเซล ทรงกลมได้ :

และ

ฟังก์ชันเบสเซล ที่ดัดแปลงแล้ว มีความสัมพันธ์ที่คล้ายคลึงกันดังนี้: และ และ

ความสัมพันธ์เวียนเกิดมีดังนี้ โดยที่C αหมายถึงI αหรือe αi π K αความสัมพันธ์เวียนเกิดเหล่านี้มีประโยชน์สำหรับปัญหาการแพร่กระจายแบบไม่ต่อเนื่อง

การก้าวข้าม

ในปี ค.ศ. 1929 คาร์ล ลุดวิก ซีเกลพิสูจน์ว่าJ ν ( x ) , J ' ν ( x )และอนุพันธ์ลอการิทึมJ ' ν ( x )/J ν ( x )จำนวนอดิศัยคือจำนวนอดิศัยเมื่อ νเป็นจำนวนตรรกยะและ xเป็นจำนวนพีชคณิตและไม่เป็นศูนย์ [ 55 ]การพิสูจน์เดียวกันนี้ยังแสดงให้เห็นว่าเป็นจำนวนอดิศัยภายใต้สมมติฐานเดียวกัน [ 56 ]

ผลรวมที่มีฟังก์ชันเบสเซล

ผลคูณของฟังก์ชันเบสเซลสองฟังก์ชันยอมรับผลรวมต่อไปนี้: จากความเท่าเทียมกันเหล่านี้ สรุปได้ว่า และเป็นผลสืบเนื่องมาจาก

ผลรวมเหล่านี้สามารถขยายให้รวมตัวคูณพจน์ซึ่งเป็นฟังก์ชันพหุนามของดัชนีได้ ตัวอย่างเช่น

ทฤษฎีบทการคูณ

ฟังก์ชันเบสเซลปฏิบัติตามทฤษฎีบทการคูณ โดยที่λและνอาจถูกเลือกเป็นจำนวนเชิงซ้อนใดๆ ก็ได้[ 57 ] [ 58 ]สำหรับ| λ 2 − 1 | < 1 , [ 57 ]นิพจน์ข้างต้นยังคงเป็นจริงหากJถูกแทนที่ด้วยYเอกลักษณ์ที่คล้ายคลึงกันสำหรับฟังก์ชันเบสเซลที่ดัดแปลงและ| λ 2 − 1 | < 1คือ และ

จุดศูนย์ของฟังก์ชันเบสเซล

สมมติฐานของบูร์เกต์

เบสเซลเองเป็นผู้พิสูจน์แต่เดิมว่าสำหรับจำนวนเต็มที่ไม่เป็นลบnสมการJ n ( x ) = 0มีคำตอบอนันต์ในx [ 59 ] อย่างไรก็ตามเมื่อพล็อตฟังก์ชันJ n ( x )บนกราฟเดียวกัน จะไม่มีศูนย์ใดที่ตรงกันสำหรับค่าn ที่แตกต่างกัน ยกเว้นศูนย์ที่x = 0ปรากฏการณ์นี้เรียกว่าสมมติฐานของบูร์เกต์ ตามชื่อของ นักคณิตศาสตร์ชาวฝรั่งเศสในศตวรรษที่ 19 ผู้ศึกษาฟังก์ชันเบสเซล โดยเฉพาะอย่างยิ่งระบุว่าสำหรับจำนวนเต็มใดๆn ≥ 0และm ≥ 1ฟังก์ชันJ n ( x )และJ n + m ( x )ไม่มีศูนย์ร่วมกันอื่นใดนอกจากศูนย์ที่x = 0สมมติฐานนี้ได้รับการพิสูจน์โดยคาร์ล ลุดวิก ซีเกลในปี 1929 [ 60 ]

การก้าวข้าม

ในปี พ.ศ. 2462 ซีเกลพิสูจน์ว่าเมื่อνเป็นจำนวนตรรกยะ รากที่ไม่เป็นศูนย์ทั้งหมดของJ ν (x)และJ ' ν (x)เป็นจำนวนอดิศัย [ 61 ]เช่นเดียวกับรากทั้งหมดของK ν (x) [ 56 ] นอกจาก นี้ ยังเป็นที่ทราบ กันว่ารากทั้งหมดของอนุพันธ์อันดับสูงสำหรับn ≤ 18 เป็น จำนวนอดิศัย ยกเว้นค่าพิเศษและ[ 61 ]

วิธีการเชิงตัวเลข

สำหรับการศึกษาเชิงตัวเลขเกี่ยวกับค่าศูนย์ของฟังก์ชัน Bessel โปรดดูGil, Segura & Temme (2007) , Kravanja et al. (1998)และโมลเลอร์ (2004) .

ค่าตัวเลข

ศูนย์แรกใน J 0 (เช่น j 0,1 , j 0,2และ j 0,3 ) เกิดขึ้นที่อาร์กิวเมนต์ประมาณ 2.40483, 5.52008 และ 8.65373 ตามลำดับ[ 62 ]

ประวัติศาสตร์

ปัญหาเกี่ยวกับคลื่นและความยืดหยุ่น

ฟังก์ชันเบสเซลปรากฏครั้งแรกในงานของแดเนียล เบอร์นูลลีในปี ค.ศ. 1732 ขณะที่เขากำลังวิเคราะห์สายที่สั่นซึ่งเป็นปัญหาที่โยฮันน์ เบอร์นูลลีบิดา ของเขาเคยศึกษามาก่อน [ 1 ]แดเนียลพิจารณาโซ่ที่ยืดหยุ่นซึ่งแขวนจากจุดคงที่ด้านบนและเป็นอิสระที่ปลายด้านล่าง[ 1 ]การแก้สมการเชิงอนุพันธ์นำไปสู่การแนะนำฟังก์ชันที่ปัจจุบันถือว่าเป็น ฟังก์ชันเบสเซล เบอร์นูลลียังได้พัฒนาวิธีการหาค่าศูนย์ของฟังก์ชันอีกด้วย[ 1 ]

ในปี ค.ศ. 1736 เลออนฮาร์ด ออยเลอร์ได้ค้นพบความเชื่อมโยงระหว่างฟังก์ชันอื่นๆ (ซึ่งปัจจุบันรู้จักกันในชื่อพหุนามลากูร์ ) กับคำตอบของเบอร์นูลลี ออยเลอร์ยังได้แนะนำห่วงโซ่ที่ไม่สม่ำเสมอซึ่งนำไปสู่การแนะนำฟังก์ชันที่เกี่ยวข้องกับฟังก์ชันเบสเซลที่ดัดแปลง[ 1 ]

ในช่วงกลางศตวรรษที่สิบแปดฌอง เลอ รอนด์ ดาเลมแบร์ได้ค้นพบสูตรแก้สมการคลื่นภายในปี 1771 มีการโต้แย้งกันระหว่างเบอร์นูลลี ออยเลอร์ ดาเลมแบร์ และโจเซฟ-หลุยส์ ลากรองจ์เกี่ยวกับลักษณะของคำตอบของสายที่สั่น[ 1 ]

ในปี 1778 ออยเลอร์ทำงานเกี่ยวกับการโก่งงอโดยนำเสนอแนวคิดของภาระวิกฤตของออยเลอร์เพื่อแก้ปัญหานี้ เขาได้นำเสนออนุกรมสำหรับ[ 1 ] นอกจาก นี้ ออยเลอร์ยังได้หาคำตอบของเมมเบรน 2 มิติที่สั่นในพิกัดทรงกระบอกในปี 1780 เพื่อ แก้สมการเชิงอนุพันธ์ของเขา เขาได้นำเสนออนุกรมกำลังที่เกี่ยวข้องกับสำหรับจำนวนเต็มn [ 1 ]

ในช่วงปลายศตวรรษที่ 18 Lagrange, Pierre-Simon LaplaceและMarc-Antoine Parsevalก็ได้ค้นพบสิ่งที่เทียบเท่ากับฟังก์ชัน Bessel เช่นกัน[ 1 ]ตัวอย่างเช่น Parseval ได้ค้นพบการแสดงแทนแบบอินทิกรัลโดยใช้โคไซน์[ 1 ]

ในช่วงต้นทศวรรษ 1800 โจเซฟ ฟูริเยร์เคยแก้สมการความร้อนในปัญหาที่มีสมมาตรทรงกระบอก[ 1 ]ฟูริเยร์ได้รับรางวัลจากสถาบันวิทยาศาสตร์แห่งฝรั่งเศสสำหรับผลงานนี้ในปี 1811 [ 1 ]แต่รายละเอียดส่วนใหญ่ของงานของเขา รวมถึงการใช้ชุดอนุกรมฟูริเยร์ยังคงไม่ได้รับการตีพิมพ์จนกระทั่งปี 1822 [ 1 ]ปัวซง ซึ่งเป็นคู่แข่งกับฟูริเยร์ ได้ขยายงานของฟูริเยร์ในปี 1823 โดยแนะนำคุณสมบัติใหม่ของฟังก์ชันเบสเซล รวมถึงฟังก์ชันเบสเซลลำดับครึ่งจำนวนเต็ม (ปัจจุบันรู้จักกันในชื่อฟังก์ชันเบสเซลทรงกลม) [ 1 ]

ปัญหาทางดาราศาสตร์

ในปี ค.ศ. 1770 ลากรองจ์ได้นำเสนอการขยายอนุกรมของฟังก์ชันเบสเซลเพื่อแก้สมการของเคปเลอร์ซึ่งเป็นสมการอดิศัยในทางดาราศาสตร์ฟรีดริช วิลเฮล์ม เบสเซลได้เห็นวิธีแก้ปัญหาของลากรองจ์แล้ว แต่พบว่าจัดการได้ยาก ในปี ค.ศ. 1813 ในจดหมายถึงคาร์ล ฟรีดริช เกาส์เบสเซลได้ทำให้การคำนวณง่ายขึ้นโดยใช้ฟังก์ชันตรีโกณมิติ[ 1 ]เบสเซลตีพิมพ์ผลงานของเขาในปี ค.ศ. 1819 โดยนำเสนอวิธีการอนุกรมฟูริเยร์โดยอิสระโดยไม่ทราบถึงผลงานของฟูริเยร์ที่ตีพิมพ์ในภายหลัง[ 1 ] ในปี ค.ศ. 1824 เบสเซลได้ทำการตรวจสอบฟังก์ชันอย่างเป็นระบบ ซึ่งทำให้ฟังก์ชันเหล่านั้นได้รับชื่อของเขา[ 1 ]ในเอกสารเก่า ฟังก์ชันเหล่านี้เรียกว่าฟังก์ชันทรงกระบอกหรือแม้แต่ฟังก์ชันเบสเซล-ฟูริเยร์[ 1 ]

ดูเพิ่มเติม

หมายเหตุ

  1. ^ a b c d e f g h i j k l m n o p q r Dutka, Jacques (1995). "เกี่ยวกับประวัติศาสตร์ยุคแรกของฟังก์ชันเบสเซล". วารสารประวัติศาสตร์วิทยาศาสตร์ที่แม่นยำ 49 ( 2): 105– 134. doi : 10.1007/BF00376544 .
  2. ^ "การเลี้ยวเบนของเส้นใย"สถาบันวิจัยการแพทย์แม็กซ์พลังค์สืบค้นเมื่อ16 มิถุนายน 2026
  3. ^ Wilensky, Michael; Brown, Jordan; Hazelton, Bryna (มิถุนายน 2023). "เหตุใดและเมื่อใดจึงคาดหวังการกระจายข้อผิดพลาดแบบเกาส์เซียนในการวัดสเปกตรัมกำลัง 21 ซม. ในยุคการแตกตัวเป็นไอออนใหม่" . Monthly Notices of the Royal Astronomical Society . 521 (4): 5191– 5206. arXiv : 2211.13576 . doi : 10.1093/mnras/stad863 .
  4. ^ Weisstein, Eric W. "ฟังก์ชันเบสเซลทรงกลมชนิดที่สอง" . MathWorld .
  5. ^ a b Weisstein, Eric W. "ฟังก์ชันเบสเซลชนิดที่สอง" . MathWorld .
  6. อรรถ เป็นข อับราโมวิตซ์ และ สเตกุน, พี. 360, 9.1.10 .
  7. ^ Whittaker, Edmund Taylor ; Watson, George Neville (1927). A Course of Modern Analysis (ฉบับที่ 4). สำนักพิมพ์มหาวิทยาลัยเคมบริดจ์. หน้า 356.ตัวอย่างเช่น ฮันเซน (1843) และ ชโลมิลช์ (1857)
  8. ↑ อับราโมวิทซ์ และสเตกุน, p. 358, 9.1.5 .
  9. ^ a b Temme, Nico M. (1996). ฟังก์ชันพิเศษ: บทนำสู่ฟังก์ชันคลาสสิกของฟิสิกส์คณิตศาสตร์ (ฉบับพิมพ์ครั้งที่ 2). นิวยอร์ก: Wiley. หน้า  228–231 . ISBN 0471113131.
  10. ไวส์สไตน์, เอริก ดับเบิลยู. "สูตรแฮนเซน-เบสเซล" . แมทเวิลด์ .
  11. ^ Bessel, F. (1824). อินทิกรัลที่เกี่ยวข้องคือสมการที่ไม่มีหมายเลขระหว่างสมการที่ 28 และ 29 โปรดทราบว่าในปัจจุบัน Besselจะเขียนว่า.
  12. ^วัตสัน,หน้า 176
  13. ^ "คุณสมบัติของฟังก์ชัน Hankel และ Bessel" . เก็บถาวรจากต้นฉบับเมื่อ 2010-09-23 . เรียกดูเมื่อ2010-10-18 .
  14. ^ "การแสดงฟังก์ชันเบสเซลในรูปแบบอินทิกรัล" . www.nbi.dk . เก็บถาวรจากต้นฉบับเมื่อวันที่ 3 ตุลาคม 2022 . เรียกดูเมื่อวันที่ 25 มีนาคม 2018 .
  15. ^ Arfken & Weber, แบบฝึกหัด 11.1.17
  16. ↑ อับราโมวิทซ์ และสเตกุน, p. 362, 9.1.69 .
  17. เซเกโก, กาบอร์ (1975) พหุนามมุมฉาก (ฉบับที่ 4) พรอวิเดนซ์ โรดไอแลนด์: AMS
  18. ^ "ฟังก์ชันเบสเซลชนิดที่หนึ่งและชนิดที่สอง" (PDF) . mhtlab.uwaterloo.ca . หน้า 3. เก็บถาวร(PDF)จากต้นฉบับเมื่อ 2022-10-09 . เรียกดูเมื่อ24 พฤษภาคม 2022 .
  19. ^ห้องสมุดดิจิทัลของ NIST เกี่ยวกับฟังก์ชันทางคณิตศาสตร์ (10.8.1) เข้าถึงทางออนไลน์เมื่อวันที่ 25 ตุลาคม 2559
  20. ^ a b Watson, หน้า 178 .
  21. ↑ อับราโมวิทซ์ และสเตกุน, p. 358, 9.1.3, 9.1.4 .
  22. ↑ อับราโมวิทซ์ และสเตกุน, p. 358, 9.1.6 .
  23. ↑ อับราโมวิทซ์ และสเตกุน, p. 360, 9.1.25 .
  24. ↑ อับราโมวิทซ์ และสเตกุน, p. 375, 9.6.2, 9.6.10, 9.6.11 .
  25. ^ Dixon; Ferrar, WL (1930). "การพิสูจน์โดยตรงของอินทิกรัลของ Nicholson" The Quarterly Journal of Mathematics . Oxford: 236–238 . doi : 10.1093/qmath/os- 1.1.236
  26. ↑ อับราโมวิทซ์ และสเตกุน, p. 375, 9.6.3, 9.6.5 .
  27. ↑ อับราโมวิทซ์ และสเตกุน, p. 374, 9.6.1 .
  28. ^ Greiner, Walter; Reinhardt, Joachim (2009). Quantum Electrodynamics . Springer. หน้า 72. ISBN 978-3-540-87561-1.
  29. ^วัตสัน,หน้า 181
  30. ^ "ฟังก์ชันเบสเซลที่ดัดแปลง §10.32 การแสดงผลแบบอินทิกรัล"ห้องสมุดดิจิทัลของฟังก์ชันทางคณิตศาสตร์ NIST NIST สืบค้นเมื่อ2024-11-20
  31. ^ Khokonov, M. Kh. (2004). "กระบวนการเรียงลำดับของการสูญเสียพลังงานโดยการปล่อยโฟตอนแข็ง" วารสารฟิสิกส์เชิงทดลองและทฤษฎี99 (4): 690– 707. Bibcode : 2004JETP...99..690K . doi : 10.1134/1.1826160 . S2CID 122599440 . ดัดแปลงมาจากสูตรที่อ้างอิงจากIS Gradshteyn และ IM Ryzhikในหนังสือ Table of Integrals, Series, and Products (Fizmatgiz, Moscow, 1963; Academic Press, New York, 1980)
  32. ^อ้างอิงถึงใน: Teichroew, D. (1957). "The Mixture of Normal Distributions with Different Variances" (PDF) . The Annals of Mathematical Statistics . 28 (2): 510– 512. doi : 10.1214/aoms/1177706981 .
  33. ↑ อับราโมวิทซ์ และสเตกุน, p. 437, 10.1.1 .
  34. ↑ อับราโมวิทซ์ และสเตกุน, p. 439, 10.1.25, 10.1.26 .
  35. ↑ อับราโมวิทซ์ และสเตกุน, p. 438, 10.1.11 .
  36. ↑ อับราโมวิทซ์ และสเตกุน, p. 438, 10.1.12 .
  37. ↑ อับราโมวิทซ์ และสเตกุน, p. 439, 10.1.39 .
  38. ^ LV Babushkina, MK Kerimov, AI Nikitin, อัลกอริทึมสำหรับการคำนวณฟังก์ชันเบสเซลลำดับครึ่งจำนวนเต็มที่มีอาร์กิวเมนต์เชิงซ้อน,หน้า 110, หน้า 111
  39. ↑ อับราโมวิทซ์ และสเตกุน, p. 439, 10.1.23, 10.1.24 .
  40. ^ Griffiths. บทนำสู่กลศาสตร์ควอนตัม ฉบับพิมพ์ครั้งที่ 2 หน้า 154
  41. ^ Du, Hong (2004). "การคำนวณการกระเจิงของ Mie". Applied Optics . 43 (9): 1951– 1956. Bibcode : 2004ApOpt..43.1951D . doi : 10.1364/ao.43.001951 . PMID 15065726 . 
  42. ↑ อับราโมวิทซ์ และสเตกุน, p. 364, 9.2.1 .
  43. ^คลังข้อมูลดิจิทัลของฟังก์ชันทางคณิตศาสตร์ NISTส่วน ที่ 10.11
  44. ↑ อับราโมวิทซ์ และสเตกุน, p. 377, 9.7.1 .
  45. ↑ อับราโมวิทซ์ และสเตกุน, p. 378, 9.7.2 .
  46. ^ Fröhlich และ Spencer 1981 ภาคผนวก B
  47. ^ a b Edwards, C. Henry; Penney, David E. (1994). สมการเชิงอนุพันธ์เบื้องต้นพร้อมการประยุกต์ใช้ (ฉบับที่ 3). Englewood Cliffs, NJ: Prentice-Hall. หน้า  273–274 . ISBN 978-0-13-312075-2.
  48. ^ Sung, S.; Hovden, R. (2022). "เกี่ยวกับอนุกรมอนันต์ของฟังก์ชันเบสเซลชนิดแรก". arXiv : 2211.01148 [ math-ph ].
  49. ↑ อับราโมวิทซ์ และสเตกุน, p. 363, 9.1.82 ff.
  50. ^ Watson, GN (25 สิงหาคม 1995). ตำราว่าด้วยทฤษฎีฟังก์ชันเบสเซล . สำนักพิมพ์มหาวิทยาลัยเคมบริดจ์. ISBN 9780521483919สืบค้นข้อมูลเมื่อวันที่ 25 มีนาคม 2561ผ่านทาง Google Books
  51. กราดชเตน, อิซเรล โซโลโมโนวิช ; ริซิค, อิโอซิฟ มอยเซวิช ; เจโรนิมัส, ยูริ เวเนียมิโนวิช ; เซย์ตลิน, มิคาอิล ยูลิเยวิช ; เจฟฟรีย์, อลัน (2015) [ตุลาคม 2014] "8.411.10". ในซวิลลิงเจอร์, แดเนียล; โมลล์, วิคเตอร์ ฮูโก (บรรณาธิการ). ตารางอินทิกรัล อนุกรม และผลิตภัณฑ์ แปลโดย Scripta Technica, Inc. (ฉบับพิมพ์ 8) สำนักพิมพ์วิชาการ อิงค์ISBN 978-0-12-384933-5. ลคซีเอ็น 2014010276 .
  52. ^อาร์ฟเคนและเวเบอร์ ส่วนที่ 11.2
  53. ↑ อับราโมวิทซ์ และสเตกุน, p. 361, 9.1.27 .
  54. ↑ อับราโมวิทซ์ และสเตกุน, p. 361, 9.1.30 น .
  55. ซีเกล, คาร์ล แอล. (2014) "Über einige Anwendungen diophantischer Approximationen" . ในการใช้งานบางประการของการประมาณค่าไดโอแฟนไทน์: การแปลของ Über einige Anwendungen diophantischer Approximationen ของ Carl Ludwig Siegel โดย Clemens Fuchs พร้อมคำอธิบายและบทความ จุดอินทิกรัลบนเส้นโค้ง: ทฤษฎีบทของ Siegel หลังจากการพิสูจน์ของ Siegel โดย Clemens Fuchs และ Umberto Zannier (ในภาษาเยอรมัน) สกูโอลา นอร์มอลเล ซูพีเรียเร หน้า  81– 138. ดอย : 10.1007/978-88-7642-520-2_2 . ไอเอสบีเอ็น 978-88-7642-520-2.
  56. ^ a b James, RD (พฤศจิกายน 1950). "บทวิจารณ์: Carl Ludwig Siegel, จำนวนอดิศัย" . Bulletin of the American Mathematical Society . 56 (6): 523– 526. doi : 10.1090/S0002-9904-1950-09435-X .
  57. อรรถ เป็นข อับราโมวิตซ์ และ สเตกุน, พี. 363, 9.1.74 .
  58. ^ Truesdell, C. (1950). "เกี่ยวกับทฤษฎีบทการบวกและการคูณสำหรับฟังก์ชันพิเศษ" . Proceedings of the National Academy of Sciences . 1950 (12): 752– 757. Bibcode : 1950PNAS...36..752T . doi : 10.1073/pnas.36.12.752 . PMC 1063284 . PMID 16578355 .  
  59. ^เบสเซล, เอฟ. (1824), บทความที่ 14.
  60. ^วัตสัน, หน้า 484–485.
  61. ^ a b Lorch, Lee; Muldoon, Martin E. (1995). "Transcendentality of zeros of higher derivatives of functions involving Bessel functions" . International Journal of Mathematics and Mathematical Sciences . 18 (3): 551– 560. doi : 10.1155/S0161171295000706 .
  62. ^อับราโมวิทซ์และสเตกัน, หน้า 409
  • Lizorkin, PI (2001) [1994], "ฟังก์ชันเบสเซล" , สารานุกรมคณิตศาสตร์ , EMS Press.
  • Karmazina, LN; Prudnikov, AP (2001) [1994], "ฟังก์ชันทรงกระบอก" , สารานุกรมคณิตศาสตร์ , EMS Press.
  • Rozov, N. Kh. (2001) [1994], "สมการเบสเซล" , สารานุกรมคณิตศาสตร์ , EMS Press.
  • หน้าเว็บ Wolfram เกี่ยวกับฟังก์ชัน Bessel JและYรวมถึงฟังก์ชัน Bessel IและK ที่ปรับปรุงแล้ว หน้าเว็บเหล่านี้ประกอบด้วยสูตร เครื่องมือประเมินค่าฟังก์ชัน และเครื่องคำนวณการสร้างกราฟ
  • หน้าเว็บ Wolfram เกี่ยวกับฟังก์ชัน Riccati-Bessel
  • ไวส์สไตน์, เอริค ดับเบิลยู. "ฟังก์ชันเบสเซลชนิดแรก" . MathWorld .
  • ฟังก์ชัน Bessel J ν , Y ν , I νและK νในคู่มือฟังก์ชัน Librow
  • FWJ Olver, LC Maximon, ฟังก์ชันเบสเซล (บทที่ 10 ของห้องสมุดดิจิทัลของฟังก์ชันทางคณิตศาสตร์)
  • Moler, CB (2004). การคำนวณเชิงตัวเลขด้วย MATLAB (PDF) . สมาคมคณิตศาสตร์อุตสาหกรรมและประยุกต์. เก็บถาวรจากต้นฉบับ(PDF)เมื่อ 2017-08-08.
ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Bessel_function&oldid=1359607442 "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ ฟังก์ชันเบสเซล

ฟังก์ชันเบสเซลเป็นคลาสของฟังก์ชันพิเศษที่มักปรากฏในปัญหาที่เกี่ยวข้องกับการเคลื่อนที่ของคลื่นการนำความร้อนและปรากฏการณ์ทางกายภาพอื่นๆ ที่มี สมมาตร แบบวงกลมหรือทรงกระบอก

แอปพลิเคชัน

สมการของเบสเซลเกิดขึ้นเมื่อต้องการหาคำตอบที่แยกตัวแปรได้ของ สมการลาปลาส และ สมการเฮล์มโฮลทซ์ ในระบบพิกัดทรงกระบอกหรือ ทรงกลม ฟังก์ชัน เบสเซลจึงมีความสำคัญอย่างยิ่งสำหรับปัญหา การแพร่กระจายคลื่น และศักย์สถิตหลายๆ ปัญหา ในการแก้ปัญหาในระบบพิกัดทรงกระบอก...

คำจำกัดความ

เนื่องจากสมการนี้เป็น สมการเชิงอนุพันธ์เชิงเส้น จึงสามารถปรับขนาดคำตอบให้มีขนาดใดก็ได้ ขนาดที่เลือกใช้สำหรับฟังก์ชันนั้นมาจากงานวิจัยในยุคแรกๆ ที่ฟังก์ชันเหล่านี้ปรากฏเป็น คำตอบ ของอินทิกรัลจำกัด แทนที่จะเป็นคำตอบของสมการเชิงอนุพันธ์...

ฟังก์ชันเบสเซลชนิดแรก: J α

ฟังก์ชันเบสเซลชนิดแรก ซึ่งแสดงด้วย J α ( x ) คือคำตอบของสมการเชิงอนุพันธ์ของเบสเซล สำหรับ α ที่เป็นจำนวนเต็มหรือบวก ฟังก์ชันเบสเซลชนิดแรกจะมีค่าจำกัดที่จุดกำเนิด ( x = 0 ) ในขณะที่สำหรับ α ที่เป็นลบและไม่ใช่จำนวนเต็ม...