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

อ่าน 11 นาที

บอร์ไวน์อินทิกรัล

ปริพันธ์

ในทางคณิตศาสตร์อินทิกรัลบอร์เวนเป็นอินทิกรัล ที่มีคุณสมบัติพิเศษซึ่งนักคณิตศาสตร์ David BorweinและJonathan Borweinเป็นผู้นำเสนอเป็นครั้งแรกในปี 2001...

บอร์ไวน์อินทิกรัล

ในทางคณิตศาสตร์อินทิกรัลบอร์เวนเป็นอินทิกรัล ที่มีคุณสมบัติพิเศษซึ่งนักคณิตศาสตร์ David BorweinและJonathan Borweinเป็นผู้นำเสนอเป็นครั้งแรกในปี 2001 [ 1 ]อินทิกรัลบอร์เวนเกี่ยวข้องกับผลคูณของซินซ์(เอx){\displaystyle \operatorname {sinc} (ax)}โดยที่ฟังก์ชัน sincกำหนดโดยซินซ์(x)=บาป(x)/x{\displaystyle \operatorname {sinc} (x)=\sin(x)/x}สำหรับx{\displaystyle x}ไม่เท่ากับ 0 และซินซ์(0)=1{\displaystyle \operatorname {sinc} (0)=1}[ 1 ] [ 2 ]

อินทิกรัลเหล่านี้มีความโดดเด่นตรงที่แสดงให้เห็นรูปแบบที่ชัดเจน แต่ในที่สุดก็จะแตกสลายไป ตัวอย่างต่อไปนี้

0บาป(x)xx=π20บาป(x)xบาป(x/3)x/3x=π20บาป(x)xบาป(x/3)x/3บาป(x/5)x/5x=π2{\displaystyle {\begin{aligned}&\int _{0}^{\infty }{\frac {\sin(x)}{x}}\,dx={\frac {\pi }{2}}\\[10pt]&\int _{0}^{\infty }{\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\,dx={\frac {\pi }{2}}\\[10pt]&\int _{0}^{\infty }{\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}{\frac {\sin(x/5)}{x/5}}\,dx={\frac {\pi }{2}}\end{aligned}}}

รูปแบบนี้ยังคงดำเนินต่อไปจนถึง

0บาป(x)xบาป(x/3)x/3บาป(x/13)x/13x=π2.{\displaystyle \int _{0}^{\infty }{\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\cdots {\frac {\sin(x/13)}{x/13}}\,dx={\frac {\pi }{2}}.}

ในขั้นตอนถัดไป รูปแบบดังกล่าวเกิดความผิดพลาด

0บาป(x)xบาป(x/3)x/3บาป(x/15)x/15x=467807924713440738696537864469935615849440640907310521750000 π0.499999999992646859 ππ22.31×1011.{\displaystyle {\begin{aligned}\int _{0}^{\infty }{\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\cdots {\frac {\sin(x/15)}{x/15}}\,dx&={\frac {467807924713440738696537864469}{935615849440640907310521750000}}~\pi \\[5pt]&\approx 0.499999999992646859~\pi \\[5pt]&\approx {\frac {\pi }{2}}-2.31\times 10^{-11}.\end{aligned}}}

โดยทั่วไป อินทิกรัลที่คล้ายกันจะมีค่าเท่ากับ ⁠π / 2⁠ เมื่อ ใดก็ตามที่ แทนที่ตัวเลข3, 5, 7… ด้วย จำนวนจริง บวกที่ผลรวมของส่วนกลับของจำนวนเหล่านั้นมี ค่าน้อยกว่า 1

ในตัวอย่างข้างต้น1 / 3 + 1 / 5 + … + 1 / 13 < 1แต่1 / 3 + 1 / 5 + … + 1 / 15 > 1

โดยการรวมปัจจัยเพิ่มเติมเข้าไปด้วย2คอส(x){\displaystyle 2\cos(x)}รูปแบบดังกล่าวยังคงใช้ได้ผลดีในระยะยาว[ 3 ]

02คอส(x)บาป(x)xบาป(x/3)x/3บาป(x/111)x/111x=π2,{\displaystyle \int _{0}^{\infty }2\cos(x){\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\cdots {\frac {\sin(x/111)}{x/111}}\,dx={\frac {\pi }{2}},}

แต่

02คอส(x)บาป(x)xบาป(x/3)x/3บาป(x/111)x/111บาป(x/113)x/113xπ22.3324×10138.{\displaystyle \int _{0}^{\infty }2\cos(x){\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\cdots {\frac {\sin(x/111)}{x/111}}{\frac {\sin(x/113)}{x/113}}\,dx\approx {\frac {\pi }{2}}-2.3324\times 10^{-138}.}

ในกรณีนี้1/3 + 1/5 + + 1/111 < 2แต่1/3 + 1/5 + + 1/113 > 2 คำ ตอบ ที่ ถูกต้องสามารถคำนวณได้ โดยใช้สูตรทั่วไป ที่ ให้ไว้ ในส่วนถัด ไป และแสดงไว้ด้านล่าง เมื่อ ขยาย ค่า นี้ อย่างสมบูรณ์จะได้เป็นเศษส่วนที่มีจำนวนเต็ม 2736 หลักสองจำนวน

π2(135113(1/3+1/5++1/1132)5625556!){\displaystyle {\frac {\pi }{2}}\left(1-{\frac {3\cdot 5\cdots 113\cdot (1/3+1/5+\dots +1/113-2)^{56}}{2^{55}\cdot 56!}}\right)}

เหตุผลที่ชุดดั้งเดิมและชุดขยายล้มเหลวได้รับการอธิบายด้วยคำอธิบายทางคณิตศาสตร์ที่เข้าใจง่าย[ 4 ​​] [ 5 ]โดยเฉพาะอย่างยิ่ง การปรับปรุง การเดินแบบสุ่มด้วยข้อโต้แย้งเชิงสาเหตุช่วยให้เข้าใจถึงการแตกของรูปแบบและเปิดทางให้กับการสรุปทั่วไปหลายประการ[ 6 ]

สูตรทั่วไป

กำหนดลำดับของจำนวนจริงที่ไม่เป็นศูนย์เอ0,เอ1,เอ2,{\displaystyle a_{0},a_{1},a_{2},\ldots }สูตรทั่วไปสำหรับอินทิกรัล

0เค=0nบาป(เอเคx)เอเคxx{\displaystyle \int _{0}^{\infty }\prod _{k=0}^{n}{\frac {\sin(a_{k}x)}{a_{k}x}}\,dx}

สามารถให้ได้[ 1 ]ในการระบุสูตร จะต้องพิจารณาผลรวมที่เกี่ยวข้องกับเอเค{\displaystyle a_{k}}โดยเฉพาะอย่างยิ่ง ถ้าγ=(γ1,γ2,,γn){±1}n{\displaystyle \gamma =(\gamma _{1},\gamma _{2},\ldots ,\gamma _{n})\in \{\pm 1\}^{n}}เป็นn{\displaystyle n}- ทูเปิลที่แต่ละรายการคือ±1{\displaystyle \pm 1}จากนั้นเราจึงเขียนγ=เอ0+γ1เอ1+γ2เอ2++γnเอn{\displaystyle b_{\gamma }=a_{0}+\gamma _{1}a_{1}+\gamma _{2}a_{2}+\cdots +\gamma _{n}a_{n}}ซึ่งเป็นผลรวมสลับกันของค่าแรกๆ ไม่กี่ค่าเอเค{\displaystyle a_{k}}และเราได้ตั้งεγ=γ1γ2γn{\displaystyle \varepsilon _{\gamma }=\gamma _{1}\gamma _{2}\cdots \gamma _{n}}ซึ่งก็คืออย่างใดอย่างหนึ่ง±1{\displaystyle \pm 1}ด้วยสัญลักษณ์นี้ ค่าของอินทิกรัลข้างต้นคือ

0เค=0nบาป(เอเคx)เอเคxx=π2เอ0ซีn{\displaystyle \int _{0}^{\infty }\prod _{k=0}^{n}{\frac {\sin(a_{k}x)}{a_{k}x}}\,dx={\frac {\pi }{2a_{0}}}C_{n}}

ที่ไหน

ซีn=12nn!เค=1nเอเคγ{±1}nεγγnsgn(γ){\displaystyle C_{n}={\frac {1}{2^{n}n!\prod _{k=1}^{n}a_{k}}}\sum _{\gamma \in \{\pm 1\}^{n}}\varepsilon _{\gamma }b_{\gamma }^{n}\operatorname {sgn}(b_{\gamma })}

ในกรณีที่เอ0>|เอ1|+|เอ2|++|เอn|{\displaystyle a_{0}>|a_{1}|+|a_{2}|+\cdots +|a_{n}|}เรามีซีn=1{\displaystyle C_{n}=1}.

นอกจากนี้ หากมีn{\displaystyle n}โดยที่สำหรับแต่ละเค=0,,n1{\displaystyle k=0,\ldots ,n-1}เรามี0<เอn<2เอเค{\displaystyle 0<a_{n}<2a_{k}}และเอ1+เอ2++เอn1<เอ0<เอ1+เอ2++เอn1+เอn{\displaystyle a_{1}+a_{2}+\cdots +a_{n-1}<a_{0}<a_{1}+a_{2}+\cdots +a_{n-1}+a_{n}}ซึ่งหมายความว่าn{\displaystyle n}คือค่าแรกเมื่อผลรวมย่อยของค่าแรกn{\displaystyle n}องค์ประกอบของลำดับเกินกว่าเอ0{\displaystyle a_{0}}, แล้วซีเค=1{\displaystyle C_{k}=1}สำหรับแต่ละคนเค=0,,n1{\displaystyle k=0,\ldots ,n-1}แต่

ซีn=1(เอ1+เอ2++เอnเอ0)n2n1n!เค=1nเอเค{\displaystyle C_{n}=1-{\frac {(a_{1}+a_{2}+\cdots +a_{n}-a_{0})^{n}}{2^{n-1}n!\prod _{k=1}^{n}a_{k}}}}

ตัวอย่างแรกคือกรณีที่เอเค=12เค+1{\displaystyle a_{k}={\frac {1}{2k+1}}}.

โปรดทราบว่าหากn=7{\displaystyle n=7}แล้วเอ7=115{\displaystyle a_{7}={\frac {1}{15}}}และ13+15+17+19+111+1130.955{\displaystyle {\frac {1}{3}}+{\frac {1}{5}}+{\frac {1}{7}}+{\frac {1}{9}}+{\frac {1}{11}}+{\frac {1}{13}}\approx 0.955}แต่13+15+17+19+111+113+1151.02{\displaystyle {\frac {1}{3}}+{\frac {1}{5}}+{\frac {1}{7}}+{\frac {1}{9}}+{\frac {1}{11}}+{\frac {1}{13}}+{\frac {1}{15}}\approx 1.02}ดังนั้นเพราะเอ0=1{\displaystyle a_{0}=1}เราเข้าใจแล้ว

0บาป(x)xบาป(x/3)x/3บาป(x/13)x/13x=π2{\displaystyle \int _{0}^{\infty }{\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\cdots {\frac {\sin(x/13)}{x/13}}\,dx={\frac {\pi }{2}}}

ซึ่งยังคงเป็นความจริงแม้ว่าเราจะนำผลิตภัณฑ์ใดๆ ออกไปก็ตาม แต่ว่า

0บาป(x)xบาป(x/3)x/3บาป(x/15)x/15x=π2(1(31+51+71+91+111+131+1511)7267!(1/31/51/71/91/111/131/15)),{\displaystyle {\begin{aligned}&\int _{0}^{\infty }{\frac {\sin(x)}{x}}{\frac {\sin(x/3)}{x/3}}\cdots {\frac {\sin(x/15)}{x/15}}\,dx\\[5pt]={}&{\frac {\pi }{2}}\left(1-{\frac {(3^{-1}+5^{-1}+7^{-1}+9^{-1}+11^{-1}+13^{-1}+15^{-1}-1)^{7}}{2^{6}\cdot 7!\cdot (1/3\cdot 1/5\cdot 1/7\cdot 1/9\cdot 1/11\cdot 1/13\cdot 1/15)}}\right),\end{aligned}}}

ซึ่งเท่ากับค่าที่ระบุไว้ก่อนหน้านี้

วิธีการแก้สมการอินทิกรัลของบอร์เวน

มีการกล่าวถึงวิธีการอินทิเกรตที่แม่นยำซึ่งมีประสิทธิภาพสำหรับการประเมินอินทิกรัลแบบ Borwein [ 7 ]วิธีการอินทิเกรตนี้ทำงานโดยการปรับปรุงการอินทิเกรตใหม่ในแง่ของชุดของการอนุพันธ์ และให้สัญชาตญาณเกี่ยวกับพฤติกรรมที่ผิดปกติของอินทิกรัล Borwein วิธีการอินทิเกรตโดยการอนุพันธ์นั้นใช้ได้กับอินทิกรัลทั่วไป รวมถึงการแปลงฟูริเยร์และลาปลาส มันถูกใช้ในเอนจินการอินทิเกรตของ Maple ตั้งแต่ปี 2019 วิธีการอินทิเกรตโดยการอนุพันธ์นั้นเป็นอิสระจากวิธีการของ Feynman ซึ่งใช้การอนุพันธ์ในการอินทิเกรตเช่นกัน

ผลิตภัณฑ์ไม่จำกัด

ในขณะที่อินทิกรัล

0เค=0nบาป(x/(2เค+1))x/(2เค+1)x{\displaystyle {\begin{aligned}\int _{0}^{\infty }\prod _{k=0}^{n}{\frac {\sin(x/(2k+1))}{x/(2k+1)}}\,dx\end{aligned}}}

กลายเป็นน้อยกว่าπ2{\displaystyle {\frac {\pi }{2}}}เมื่อไรn{\displaystyle n}เมื่อเกิน 6 มันจะไม่มีวันน้อยลงมากนัก และในความเป็นจริง Borwein และ Bailey [ 8 ]ได้แสดงให้เห็นแล้ว

0เค=0บาป(x/(2เค+1))x/(2เค+1)x=0ลิมnเค=0nบาป(x/(2เค+1))x/(2เค+1)x=ลิมn0เค=0nบาป(x/(2เค+1))x/(2เค+1)xπ20.0000352{\displaystyle {\begin{aligned}\int _{0}^{\infty }\prod _{k=0}^{\infty }{\frac {\sin(x/(2k+1))}{x/(2k+1)}}\,dx&=\int _{0}^{\infty }\lim _{n\to \infty }\prod _{k=0}^{n}{\frac {\sin(x/(2k+1))}{x/(2k+1)}}\,dx\\[5pt]&=\lim _{n\to \infty }\int _{0}^{\infty }\prod _{k=0}^{n}{\frac {\sin(x/(2k+1))}{x/(2k+1)}}\,dx\\[5pt]&\approx {\frac {\pi }{2}}-0.0000352\end{aligned}}}

โดยที่เราสามารถดึงลิมิตออกมาจากอินทิกรัลได้ด้วยทฤษฎีบทการลู่เข้าแบบครอบงำ ในทำนองเดียวกัน ในขณะที่

02คอสxเค=0nบาป(x/(2เค+1))x/(2เค+1)x{\displaystyle \int _{0}^{\infty }2\cos x\prod _{k=0}^{n}{\frac {\sin(x/(2k+1))}{x/(2k+1)}}\,dx}

กลายเป็นน้อยกว่าπ2{\displaystyle {\frac {\pi }{2}}}เมื่อไรn{\displaystyle n}เกิน 55 เรามี

02คอสxเค=0nบาป(x/(2เค+1))x/(2เค+1)xπ22.96291042{\displaystyle \int _{0}^{\infty }2\cos x\prod _{k=0}^{n}{\frac {\sin(x/(2k+1))}{x/(2k+1)}}\,dx\approx {\frac {\pi }{2}}-2.9629\cdot 10^{-42}}

นอกจากนี้ การใช้การแยกตัวประกอบของไวเออร์สตรัส

บาปxx=n=1(1x2π2n2)คอสx=n=0(14x2π2(2n+1)2){\displaystyle {\frac {\sin x}{x}}=\prod _{n=1}^{\infty }\left(1-{\frac {x^{2}}{\pi ^{2}n^{2}}}\right)\qquad \cos x=\prod _{n=0}^{\infty }\left(1-{\frac {4x^{2}}{\pi ^{2}(2n+1)^{2}}}\right)}

สามารถแสดงให้เห็นได้

n=0บาป(2x/(2n+1))2x/(2n+1)=n=1คอส(xn){\displaystyle \prod _{n=0}^{\infty }{\frac {\sin(2x/(2n+1))}{2x/(2n+1)}}=\prod _{n=1}^{\infty }\cos \left({\frac {x}{n}}\right)}

และด้วยการเปลี่ยนตัวแปรจะได้[ 9 ]

0n=1คอส(xn)x=120n=0บาป(x/(2n+1))x/(2n+1)xπ40.0000176{\displaystyle \int _{0}^{\infty }\prod _{n=1}^{\infty }\cos \left({\frac {x}{n}}\right)\,dx={\frac {1}{2}}\int _{0}^{\infty }\prod _{n=0}^{\infty }{\frac {\sin(x/(2n+1))}{x/(2n+1)}}\,dx\approx {\frac {\pi }{4}}-0.0000176}

และ[ 8 ] [ 10 ]

0คอส(2x)n=1คอส(xn)x=120คอส(x)n=0บาป(x/(2n+1))x/(2n+1)xπ87.40731043{\displaystyle \int _{0}^{\infty }\cos(2x)\prod _{n=1}^{\infty }\cos \left({\frac {x}{n}}\right)\,dx={\frac {1}{2}}\int _{0}^{\infty }\cos(x)\prod _{n=0}^{\infty }{\frac {\sin(x/(2n+1))}{x/(2n+1)}}\,dx\approx {\frac {\pi }{8}}-7.4073\cdot 10^{-43}}

การกำหนดสูตรความน่าจะเป็น

Schmuland [ 11 ] ได้นำเสนอสูตรความน่าจะเป็นที่น่าสนใจของอินทิกรัล Borwein ผลคูณอนันต์ ตัวอย่างเช่น พิจารณาอนุกรมฮาร์มอนิกแบบสุ่ม

±1±12±13±14±15±{\displaystyle \pm 1\pm {\frac {1}{2}}\pm {\frac {1}{3}}\pm {\frac {1}{4}}\pm {\frac {1}{5}}\pm \cdots }

โดยการโยนเหรียญที่ยุติธรรมอย่างอิสระเพื่อเลือกเครื่องหมาย อนุกรมนี้ลู่เข้าเกือบแน่นอนนั่นคือด้วยความน่าจะเป็น 1 ฟังก์ชันความหนาแน่นความน่าจะเป็นของผลลัพธ์เป็นฟังก์ชันที่กำหนดไว้อย่างดี และค่าของฟังก์ชันนี้ที่ 2 ใกล้เคียงกับ 1/8 อย่างไรก็ตาม มันใกล้เคียงกับ

0.124999999999999999999999999999999999999999764{\displaystyle 0.124999999999999999999999999999999999999999764\ldots }

คำอธิบายของ Schmuland คือปริมาณนี้คือ1/π{\displaystyle 1/\pi }ครั้ง

0คอส(2x)n=1คอส(xn)xπ87.40731043{\displaystyle \int _{0}^{\infty }\cos(2x)\prod _{n=1}^{\infty }\cos \left({\frac {x}{n}}\right)\,dx\approx {\frac {\pi }{8}}-7.4073\cdot 10^{-43}}
  • ไวส์สไตน์, เอริค ดับเบิลยู. "อินทิกรัลผลคูณโคไซน์อนันต์" . MathWorld . สืบค้นเมื่อ10 มกราคม 2023 .
ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Borwein_integral&oldid=1358355713 "

สรุปเนื้อหา

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

ข้อมูลสำคัญเกี่ยวกับ บอร์ไวน์อินทิกรัล

ในทางคณิตศาสตร์อินทิกรัลบอร์เวนเป็นอินทิกรัล ที่มีคุณสมบัติพิเศษซึ่งนักคณิตศาสตร์ David BorweinและJonathan Borweinเป็นผู้นำเสนอเป็นครั้งแรกในปี 2001...

สูตรทั่วไป

กำหนดลำดับของจำนวนจริงที่ไม่เป็นศูนย์ เอ 0 , เอ 1 , เอ 2 , … {\displaystyle a_{0},a_{1},a_{2},\ldots } สูตรทั่วไปสำหรับอินทิกรัล

วิธีการแก้สมการอินทิกรัลของบอร์เวน

มีการกล่าวถึงวิธีการอินทิเกรตที่แม่นยำซึ่งมีประสิทธิภาพสำหรับการประเมินอินทิกรัลแบบ Borwein [ 7 ] วิธีการอินทิเกรตนี้ทำงานโดยการปรับปรุงการอินทิเกรตใหม่ในแง่ของชุดของการอนุพันธ์ และให้สัญชาตญาณเกี่ยวกับพฤติกรรมที่ผิดปกติของอินทิกรัล Borwein...

การกำหนดสูตรความน่าจะเป็น

Schmuland [ 11 ] ได้นำเสนอสูตรความน่าจะเป็นที่น่าสนใจของ อินทิกรัล Borwein ผล คูณอนันต์ ตัวอย่างเช่น พิจารณาอนุกรมฮาร์มอนิกแบบสุ่ม