4.6.19 · D3 · Maths › Ordinary Differential Equations › Bessel's equation and Bessel functions (intro, physical rele
Intuition Yeh page kis liye hai
Parent note ne bataya tha ki Bessel's equation kya hai aur kahan se aati hai. Yeh page ek drill sheet hai. Hum har us situation ki list banate hain jo ek problem aapke saamne rakh sakti hai — alag-alag orders ν , tricky centre point x = 0 , "disguised" equations, drum frequencies, half-integer shortcuts, aur ek exam-style twist — aur phir hum har cell ke liye ek example work karte hain taaki aap kabhi aisa scenario na dekho jo aapne pehle dekha hi na ho.
Agar neeche koi symbol unfamiliar lage, toh sab kuch the parent note mein rebuilt hai, aur machinery Frobenius method aur Regular singular points se aati hai.
Definition Is page par "Verify" aur "VERIFY" ka matlab
Har example ek Verify line ke saath khatam hota hai: ek quick sanity check jo aap haath se kar sakte ho — answer wapas plug karo, units check karo, ya koi known limit lo. Alag se, har numeric answer bhi ek short computer-algebra snippet se machine-check hota hai jo is page ke saath stored hai (jis section ko site VERIFY kehti hai). Aap unhe khud reproduce kar sakte ho: same expression kisi bhi Python SymPy session mein type karo (ya decimals ke liye calculator mein bhi) aur confirm karo ki number match karta hai. Neeche "Checked in the automated checks" hamesha unhi reproducible snippets mein se kisi ek ki taraf point karta hai.
Pehle yeh table padho. Har cell ek alag cheez hai jo exam test kar sakta hai. Neeche ke examples us cell ke saath tagged hain jise woh cover karte hain.
Cell
Kya badalta hai
Kaunsa question puchha jaata hai
Example
A Integer order
ν = n (whole number)
Standard form pehchano, solution ka naam batao
Ex 1
B Half-integer order
ν = 2 1 , 2 3
Kya yeh sin/cos mein collapse ho jaata hai?
Ex 2
C Centre included (r = 0 )
region ek solid disc hai
Boundedness ke liye kaunsa solution bachta hai?
Ex 3
D Centre excluded (annulus)
beech mein hole hai
Kya hum dono solutions rakhte hain?
Ex 4
E Disguised equation
abhi x 2 y ′′ + x y ′ + … form mein nahi hai
ν reveal karne ke liye Rescale/multiply karo
Ex 5
F Rescaled argument
x = λ r
λ 's kahan jaate hain?
Ex 6
G Degenerate: ν = 0
order vanish ho jaata hai
J 0 ( 0 ) = 1 , frequencies
Ex 7
H Real-world word problem
physical drum, units
Physics ko ek number mein badlo
Ex 7
I Limiting behaviour
x → 0 aur x → ∞
Small/large-x shapes
Ex 8
J Exam twist
ek "− ν 2 " sign trap / wrong-order trap
Misread pakdo
Ex 9
Definition Is page par use hone wale Symbols (har ek use se pehle earned)
x — independent variable (ek number jo aap feed karte ho). Har graph ke horizontal axis par chalta hai.
r — radial coordinate : disc ya cylinder ke centre se us point tak seedhi-line ki doori jis par aap dhyan de rahe ho. Har physical example mein yeh "middle se kitna door" wala number hai, aur yeh x se ek rescaling jaise x = k r ya x = λ r ke through link hota hai.
ν (Greek "nu") — order , ek fixed constant jo har problem mein ek baar choose kiya jaata hai. Kabhi vary nahi karta.
J ν ( x ) — Bessel function of the first kind : wiggly, decaying "sine for circles" jo x = 0 par finite rehta hai.
Y ν ( x ) — Bessel function of the second kind : partner solution jo x = 0 par − ∞ ki taraf jaata hai .
c 1 , c 2 — do arbitrary constants (integration constants) jo koi bhi second-order linear ODE carry karta hai. Inhe baad mein boundary/finiteness conditions se fix kiya jaata hai; solid disc par unme se ek zero hone par majboor hota hai.
α ν , k — J ν ka k -wa positive zero (x ki ek value jahan curve axis cross karta hai).
Γ — Gamma function , factorial ka smooth version: Γ ( n + 1 ) = n ! .
x 2 y ′′ + x y ′ + ( x 2 − 25 ) y = 0 ko solve karo.
Forecast: Aage padhne se pehle ν ki value guess karo. Kaunsa number, square hone par, 25 deta hai?
Step 1 — Standard form se line up karo.
x 2 y ′′ + x y ′ + ( x 2 − ν 2 ) y = 0.
Yeh step kyun? Bessel's equation ka ek rigid skeleton hai — coefficients x 2 , x , ( x 2 − ν 2 ) . Term by term matching sabse fast legal move hai; pattern pehchan lene ke baad koi integration zaruri nahi.
Step 2 — Order read off karo. x 2 se subtract kiya gaya constant ν 2 hai, isliye ν 2 = 25 ⇒ ν = 5 (hum ν ≥ 0 lete hain).
Yeh step kyun? ν non-negative define hota hai; sign convention se fix hoti hai, aur J − 5 = J 5 waise bhi sign tak hota hai.
Step 3 — General solution likho.
y ( x ) = c 1 J 5 ( x ) + c 2 Y 5 ( x ) .
Yeh step kyun? Ek second-order linear ODE mein hamesha do independent solutions hote hain; J 5 aur Y 5 standard pair hain.
Verify: ν 2 = 5 2 = 25 ✓ "− 25 " se match karta hai. Equation mein y ′′ par koi extra multiplier nahi hai, isliye koi rescaling chupi hui nahi thi.
Dikhao ki series use karke J 1/2 ( x ) = π x 2 sin x hai, aur ise x = 2 π par evaluate karo.
Forecast: Half-integer orders secretly elementary hote hain. Kya aap x = π /2 par ek clean number expect karte ho?
Step 1 — ν = 2 1 ke saath series likho.
J 1/2 ( x ) = ∑ m = 0 ∞ m ! Γ ( m + 2 3 ) ( − 1 ) m ( 2 x ) 2 m + 1/2 .
Yeh step kyun? Frobenius method J ν ko is power series ke roop mein deta hai; ν = 2 1 plug karna ek direct substitution hai.
Step 2 — Gamma function half-integer values use karo. Γ ( m + 2 3 ) = 2 2 m + 1 m ! ( 2 m + 1 )! π ke saath, m ! cancel ho jaata hai aur sum sin x ki Taylor series mein rearrange ho jaata hai jise x se divide kiya gaya ho.
Yeh step kyun? Half-integer Gamma hi ek single bridge hai jo abstract series ko uss function se wapas jodte hai jise aap pehle se jaante ho. Yahi collapse ko hone deta hai.
Step 3 — Closed form padho.
J 1/2 ( x ) = π x 2 sin x .
Yeh step kyun? Jab rearranged sum literally 2/ ( π x ) times sin x Taylor series ho jaata hai, toh hum us series ko uske named function se replace karne ke haqdaar hain — yahi recognition sara payoff hai, ek infinite sum ko ek elementary expression mein turn karta hai jise hum graph aur differentiate kar sakte hain.
Step 4 — x = π /2 par evaluate karo. sin ( π /2 ) = 1 , isliye
J 1/2 ( 2 π ) = π ⋅ π /2 2 ⋅ 1 = π 2 4 = π 2 .
Yeh step kyun? Ek concrete number sabse best proof hai ki closed form sahi hai.
Verify: Numerically 2/ π ≈ 0.6366 . Ise kisi bhi calculator par reproduce kiya ja sakta hai, aur automated SymPy snippet symbolic identity 2/ ( π ⋅ π /2 ) sin ( π /2 ) = 2/ π confirm karta hai.
Ek solid cylinder (0 ≤ r ≤ a ) mein heat. Radial part obeys karta hai x 2 R ′′ + x R ′ + ( x 2 − 4 ) R = 0 jahan x = k r . Kaunsa solution bachta hai, aur kyun?
Forecast: General solution mein do constants c 1 , c 2 hain. Guess karo ki jab centre r = 0 region ka part ho toh kitne bachte hain.
Step 1 — Order identify karo aur general solution likho. ν 2 = 4 ⇒ ν = 2 , isliye
R ( r ) = c 1 J 2 ( k r ) + c 2 Y 2 ( k r ) .
Yeh step kyun? Standard form recognition — Ex 1 jaisa hi skeleton.
Step 2 — Centre par finiteness apply karo. Jab r → 0 , Y 2 ( k r ) → − ∞ (yeh x − 2 jaisa behave karta hai). Axis par physical temperature finite hai, isliye hume zaruri c 2 = 0 set karna hoga.
Yeh step kyun? Mathematics do solutions offer karta hai; physics (boundedness) extra condition hai jo ek select karta hai.
Neeche ki picture dono solutions ko same axes par plot karti hai. Magenta curve J 2 follow karo: yeh 0 par flat shuru hota hai, ek gentle hump tak rise karta hai, aur bounded rehta hai — ek legal temperature profile. Ab violet curve Y 2 follow karo: jab x right se 0 ki taraf approach karta hai yeh − ∞ ki taraf plunge karta hai (arrow us dive ko mark karta hai). Koi real cylinder apni axis par infinite temperature nahi rakh sakta, isliye violet solution physically forbidden hai aur uska coefficient c 2 vanish hona chahiye.
Step 3 — Final answer. R ( r ) = c 1 J 2 ( k r ) , disc par har jagah finite.
Verify: r = 0 par, J 2 ( 0 ) = 0 (finite ✓) jabki Y 2 ( k r ) → − ∞ (excluded ✓). Automated snippet confirm karta hai ki J 2 ke leading series term se J 2 ( 0 ) = 0 milta hai; aap dono facts seedhe figure se bhi padh sakte hain.
Ex 3 jaisa hi equation (ν = 2 ), lekin ab annulus b ≤ r ≤ a par (hole wala pipe). Mode condition setup karo. Kitne constants bachte hain?
Forecast: Centre r = 0 ab region mein nahi hai. Kya isse farq padta hai ki hum Y 2 drop karte hain ya nahi?
Step 1 — Poora general solution rakho.
R ( r ) = c 1 J 2 ( k r ) + c 2 Y 2 ( k r ) .
Yeh step kyun? Y 2 sirf r = 0 par misbehave karta hai, jo [ b , a ] ke bahar hai. Kuch bhi c 2 = 0 force nahi karta ab — classic "annulus keeps two constants" rule.
Step 2 — Do boundary conditions lagao. Dono walls par fixed: R ( b ) = 0 aur R ( a ) = 0 :
c 1 J 2 ( k b ) + c 2 Y 2 ( k b ) = 0 , c 1 J 2 ( k a ) + c 2 Y 2 ( k a ) = 0.
Yeh step kyun? Do unknowns ( c 1 , c 2 ) ko do conditions chahiye. Yeh ek 2 × 2 homogeneous system hai.
Step 3 — Non-trivial modes ke liye vanishing determinant zaruri hai.
J 2 ( k b ) Y 2 ( k a ) − J 2 ( k a ) Y 2 ( k b ) = 0.
Yeh step kyun? Ek homogeneous system mein non-zero ( c 1 , c 2 ) tabhi hote hain jab determinant zero ho — yeh transcendental equation k ko quantise karta hai.
Verify: b → 0 limit mein term Y 2 ( k b ) → − ∞ dominate karta hai, J 2 ( k a ) = 0 force karta hai — exactly Ex 3 ka solid-disc rule. Annulus condition sahi se degenerate hoti hai. (Consistency check, koi single number nahi.)
y ′′ + x 1 y ′ + ( 1 − 4 x 2 9 ) y = 0 ko Bessel form mein reduce karo aur order ka naam batao.
Forecast: 4 x 2 9 mein ν 2 chupi hai. Kaam karne se pehle ν guess karo.
Step 1 — Denominators clear karo. Pure equation ko x 2 se multiply karo:
x 2 y ′′ + x y ′ + ( x 2 − 4 9 ) y = 0.
Yeh step kyun? Bessel ke skeleton ko y ′′ par x 2 aur y ′ par x chahiye. x 2 se multiply karna woh ek move hai jo exactly woh skeleton produce karta hai.
Step 2 — ν read off karo. ν 2 = 4 9 ⇒ ν = 2 3 .
Yeh step kyun? Ex 1 jaisa hi rule — x 2 se subtract kiya gaya constant ν 2 hai.
Step 3 — Solution ka naam batao. y = c 1 J 3/2 ( x ) + c 2 Y 3/2 ( x ) , aur kyunki 2 3 half-integer hai yeh bhi elementary hai:
J 3/2 ( x ) = π x 2 ( x s i n x − cos x ) .
Yeh step kyun? Ex 2 mein ν = 2 1 ki tarah, har half-integer order sines aur cosines mein collapse ho jaata hai (usi Gamma function identity ke through), isliye answer ko abstract J 3/2 ke roop mein chhodne ki jagah hum uska elementary closed form quote karte hain — isse solution kuch aisa ban jaata hai jo hum evaluate, plot, aur sanity check ke roop mein use kar sakein.
Verify: ( 2 3 ) 2 = 4 9 ✓. J 3/2 ko x = π par evaluate karo: sin π = 0 , cos π = − 1 , jo π 2 2 ⋅ ( 0 − ( − 1 )) = 2 / π ≈ 0.4502 deta hai — automated snippet se confirm hota hai, aur decimal calculator par reproduce kar sakte ho.
Radial equation r 2 R ′′ + r R ′ + ( λ 2 r 2 − ν 2 ) R = 0 se shuru karke, x = λ r substitute karo aur dikhao ki yeh plain Bessel's equation ban jaata hai.
λ kya hai? Yahan λ (Greek "lambda") ek fixed positive constant hai — physically yeh wavenumber k (spatial frequency) hai jo Separation of variables se separation constant ke roop mein nikla; yeh batata hai ki per unit radius kitne wavelengths fit hote hain. Yeh ek variable nahi hai: ν ki tarah, yeh har problem mein ek baar choose hota hai.
Forecast: Kya λ 's final equation mein bachenge, ya cancel ho jayenge?
Step 1 — Substitution ke liye chain rule. x = λ r ke saath, d r d = λ d x d , isliye
R ′ = λ d x d R , R ′′ = λ 2 d x 2 d 2 R .
Yeh step kyun? Hum "ruler" ko r se x mein change kar rahe hain; chain rule hi derivatives ko rulers ke beech convert karne ka ek sahi tarika hai.
Step 2 — Substitute karo aur simplify karo. r = x / λ use karke:
λ 2 x 2 λ 2 R xx + λ x λ R x + ( λ 2 λ 2 x 2 − ν 2 ) R = 0.
Yeh step kyun? Har r x / λ ban jaata hai aur har λ 2 r 2 x 2 ban jaata hai; hum har factor track karte hain.
Step 3 — Dekho λ 's cancel ho jaate hain.
x 2 R xx + x R x + ( x 2 − ν 2 ) R = 0.
Yeh step kyun? Cancellation confirm karna hi point hai: frequency λ sirf argument x = λ r mein rehta hai, J ν ki shape mein kabhi nahi. Isliye R ( r ) = J ν ( λ r ) hota hai.
Verify: Coefficient audit — λ 2 x 2 λ 2 = x 2 ✓, λ x λ = x ✓, λ 2 λ 2 x 2 = x 2 ✓. Saare λ chale gaye. Symbolic snippet confirm karta hai ki har coefficient exactly reduce hota hai.
Radius a = 0.30 m wale circular drumhead ki wave speed c = 100 m/s hai. Symmetric mode ki lowest (fundamental) frequency nikalo. J 0 ka pehla zero use karo: α 0 , 1 = 2.4048 .
Forecast: Guitar string ka fundamental f = c / ( 2 L ) hota hai. Guess karo ki drum ka J 0 -zero 2.4048 iske fundamental ko us naive value se zyada banata hai ya kam.
Step 1 — Bounded solution setup karo. Symmetric mode ⇒ ν = 0 , centre included ⇒ sirf J 0 rakho:
R ( r ) = c 1 J 0 ( k r ) , k = c 2 π f .
Yeh step kyun? ν = 0 isliye kyunki koi angular variation nahi hai; Y 0 boundedness se drop ho jaata hai (Cell C logic). Yahan J 0 ( 0 ) = 1 , isliye centre sabse zyada move karta hai — yeh degenerate lekin bahut physical case hai.
Step 2 — Fixed rim quantisation deta hai. R ( a ) = 0 ⇒ J 0 ( k a ) = 0 ⇒ k a = α 0 , 1 = 2.4048.
Yeh step kyun? Rim clamped hai, isliye wahan displacement zero hai; sirf special k jo J 0 ke zero par hit karte hain allowed hain.
Step 3 — Frequency ke liye solve karo.
f = 2 π a c α 0 , 1 = 2 π × 0.30 100 × 2.4048 .
Yeh step kyun? k = 2 π f / c ko k a = α 0 , 1 mein substitute karo aur f isolate karo; yahi ek rearrangement hai jo quantisation condition ko ek audible number mein turn karta hai.
Step 4 — Number.
f = 1.88496 240.48 ≈ 127.6 Hz .
Yeh step kyun? Arithmetic ko ek single decimal value tak le jaana hi result ko checkable aur physically meaningful banata hai — middle C ke thoda neeche wala B note.
Verify: Units: m ( m/s ) = s − 1 = Hz ✓. Numeric value ≈ 127.6 Hz — aap ise calculator par reproduce kar sakte ho (100 × 2.4048/ ( 2 π × 0.30 ) ), aur automated snippet se confirm hota hai.
Neeche ki picture symmetric mode shape dikhati hai. Orange curve J 0 ( x ) hai; note karo ki yeh centre par height 1 se shuru hota hai (J 0 ( 0 ) = 1 , violet arrow se mark kiya gaya — drum ka middle sabse zyada move karta hai) aur apne first zero x = α 0 , 1 = 2.4048 (magenta dot) tak aata hai. Wahi zero exactly wahan hai jahan clamped rim hota hai, isliye k a = 2.4048 ; Step 2 ki saari cheez bas is crossing ko graph se padhna hai.
J 0 ke liye: (a) leading small-x value nikalo, (b) large-x asymptotic se J 0 ( 20 ) estimate karo aur true value ≈ 0.1670 se compare karo.
Forecast: Small x → 1 ke paas? Large x → tiny aur oscillating? J 0 ( 20 ) ka sign guess karo.
Step 1 — Small-x limit. Series J 0 ( x ) = 1 − 4 x 2 + ⋯ se, isliye J 0 ( 0 ) = 1 .
Yeh step kyun? Frobenius method series ka lowest term x → 0 par dominate karta hai; har higher term ek extra x 2 carry karta hai.
Step 2 — Large-x asymptotic. Bade x ke liye,
J 0 ( x ) ≈ π x 2 cos ( x − 4 π ) .
Yeh step kyun? 1/ x envelope encode karta hai energy ka spreading growing circle par; yeh formula jawab deta hai "door tak wiggle kaisa dikhta hai?"
Step 3 — x = 20 plug karo. x − π /4 = 20 − 0.7854 = 19.2146 rad. 2 π modulo reduce karte hue, 19.2146 − 6 π = 19.2146 − 18.8496 = 0.3650 rad, isliye cos ( 19.2146 ) = cos ( 0.3650 ) ≈ 0.9341 . Amplitude 2/ ( 20 π ) = 0.031831 ≈ 0.17841 . Product ≈ 0.17841 × 0.9341 ≈ 0.1667 .
Yeh step kyun? Direct evaluation dikhata hai ki asymptotic remarkably close hai (true ≈ 0.1670 ) moderate x par bhi — yahan 1% se better agreement.
Verify: Asymptotic ≈ 0.1667 vs true ≈ 0.1670 ; same sign, essentially same magnitude. Dono numbers automated snippet se confirm hote hain.
Neeche ki picture J 0 ( x ) (magenta) ko do dashed violet envelope curves ± 2/ ( π x ) ke beech trapped plot karti hai. Dekho kaise har hump envelope ko touch karta hai aur envelope kaise x badh ne ke saath zero ki taraf squeeze hoti hai — woh squeeze hi 1/ x decay hai. Orange dot x = 20 par us value ko mark karta hai jo humne abhi estimate ki, J 0 ( 20 ) ≈ 0.167 , envelope ke andar aaram se baitha hua.
Ek student likhta hai: "x 2 y ′′ + x y ′ + ( x 2 + 16 ) y = 0 mein ν = 4 hai, isliye y = c 1 J 4 + c 2 Y 4 ." Unki galti dhundho aur sahi solution do.
Forecast: Equation mein + 16 hai, − 16 nahi. Kya yeh phir bhi order 4 ka Bessel hai?
Step 1 — Signs dhyan se compare karo. Standard Bessel mein ( x 2 − ν 2 ) hota hai. Yahan hamein ( x 2 + 16 ) = ( x 2 − ( − 16 )) milta hai, isliye ν 2 = − 16 .
Yeh step kyun? Trap yeh hai ki + 16 ko − 16 jaisa padh lo. Constant ke aage ka sign sab kuch hai.
Step 2 — ν 2 = − 16 interpret karo. Tab ν = ± 4 i (imaginary order), kyunki − 16 = 4 i . Yeh order 4 nahi hai — x 2 mein added constant positive tha, jo koi real order produce nahi kar sakta.
Yeh step kyun? ν ko seedha ν 2 = − ( constant added to x 2 ) se padhna chahiye; blindly 16 lena sign ko ignore karta hai aur galat function galat zeros ke saath deta hai.
Step 3 — Sahi label batao. Solution order 4 i Bessel functions use karta hai:
y ( x ) = c 1 J 4 i ( x ) + c 2 Y 4 i ( x ) ,
equivalently inka ek real oscillatory combination — J 4 , Y 4 nahi . Isliye student ka answer do jagah galat hai: order imaginary hai (4 i , na 4 ), aur isliye named functions aur unka zero-spacing bilkul alag hain.
Yeh step kyun? J 4 naam lene se galat zeros ka set predict hoga aur isliye galat drum spectrum; sahi label us downstream error ko rokta hai.
Step 4 — Original ka sahi reading. Standard form mein likha, x 2 y ′′ + x y ′ + ( x 2 − ( − 16 )) y = 0 ka order ν = 4 i hai. Student ko jo fix apply karna chahiye : hamesha constant ko "x 2 minus something" ke roop mein rewrite karo aur ν 2 = woh something apne sign ke saath lo.
Verify: Sanity: ν 2 = 16 ke liye equation ko ( x 2 − 16 ) padhna hoga. Kyunki yeh actually ( x 2 + 16 ) padhti hai, isliye ν 2 = − 16 = 16 — student ka ν = 4 refute ho jaata hai. Sign check automated snippet se confirm hota hai.
Common mistake Sabhi cells mein #1 recurring error
Constant ka sign galat padhna. Hamesha "x 2 minus something" ke roop mein rewrite karo: ( x 2 − ν 2 ) . Agar aap + 16 dekhte ho, toh "something" − 16 hai, aur ν imaginary hai — bilkul alag problem (Ex 9).
Recall Quick self-test (kaunsa cell?)
Diya gaya x 2 y ′′ + x y ′ + ( x 2 − 4 49 ) y = 0 solid disc par — kaunse examples ki logic apply hoti hai? ::: Cell B/E (half-integer ν = 7/2 ) order ke liye, aur Cell C (drop Y ) disc ke liye.
Annulus par, kitne arbitrary constants bachte hain? ::: Do — aap dono J ν aur Y ν rakhte ho (Ex 4).
Solution mein frequency λ kahan chupta hai? ::: Sirf argument ke andar, J ν ( λ r ) ke roop mein; equation ki shape mein koi λ nahi hota (Ex 6).
Bessel's equation and Bessel functions (intro, physical relevance) — woh parent jise yeh drill sheet expand karti hai.
Frobenius method — J ν series produce ki jo Ex 2, 5, 8 mein use hui.
Regular singular points — kyun Y ν x = 0 par diverge karta hai (Ex 3, 4).
Gamma function — half-integer values jo Ex 2 aur Ex 5 ko collapse karte hain.
Separation of variables & Laplacian in polar and cylindrical coordinates — jahan Ex 7 ka drum equation aata hai, aur separation constant λ ka source.
Sturm-Liouville theory — woh framework jo guarantee karta hai ki zeros α ν , k ek complete orthogonal set dete hain.