數值積分(近似值)

2009-01-31 3:26 am
設∫[-1, 1] f(x) dx = w(1)f(x1)+w(2)f(x2)+...+w(n)f(xn)
對所有多項式f(x), deg(f(x))< 2n 均成立,
其中 n 為正整數, w(1),..,w(n), x1,..,xn為與 f(x) 無關之常數
1. 證明: x1, ... , xn 滿足 (x^2-1)^n之 n 階導數 = 0
2. 若 n= 4, 試求 x1,.., x4=?
更新1:

註:與棣之前問題有關!

回答 (5)

2009-01-31 11:13 am
✔ 最佳答案
已知 (x^2-1)^n之 n 階導數為Legendre polynomial Pn(x)。
需應用Legendre polynomials關系式:
(1 - x2) P"n(x) -2x P'n(x) + n(n + 1) Pn(x) = 0

n(n + 1) Pn(x) = 2x P'n(x) + (x^2 -1) P"n(x)
考慮式右,由-1至1積分之
∫ 2x P'n(x) dx + ∫ (x^2 -1) P"n(x) dx
= ∫ P'n(x) d(x^2) + ∫ (x^2 -1) P"n(x) dx
= { x^2 P'n(x)}_{-1 to 1} - ∫x^2 P"n(x) dx + ∫ (x^2 -1) P"n(x) dx
= { x^2 P'n(x)}_{-1 to 1} + ∫ ( -1) P"n(x) dx
= { P'n(1) - P'n(-1)} + ∫ ( -1) P"n(x) dx
= { P'n(1) - P'n(-1)} + { ( -1) P'n(x)}_{-1 to 1}
= 0
這表示式左由-1至1積分之也是零,因此
∫Pn(x) dx = 0 .... (*)
顯見 w(1), w(2), ..., w(n) 不會全等於零。例如:
g(x) = (x -x2)^2 (x -x3)^2 ... (x -xn)^2
我們有 g(x) >=0 for all x 及 deg( g) < 2n,故
0 < ∫g(x) dx = w(1) g(x1)
由於 g(x1) >0 ,所以 w(1) >0
同理,我們有 w(k) >0 , k=1,2, 3, ..., n .... (**)
由 (*)及(**),我們必有 Pn(xk) = 0 , k= 1, 2, ..., n


用軟件 MAPLE,我們有
P4(x) = 1680*x^4 - 1440*x^2 + 144
再用 MAPLE 得
x1= -(1/35)*sqrt(525+70*sqrt(30))
x2 = (1/35)*sqrt(525+70*sqrt(30))
x3 = -(1/35)*sqrt(525-70*sqrt(30))
x4 = (1/35)*sqrt(525-70*sqrt(30))

2009-01-31 18:03:09 補充:
犯了致命錯誤,不過今次夠時間可以敗者復活。
....
同理,我們有 w(k) >0 , k=1,2, 3, ..., n .... (**)
另一方面,Legendre polynomials 有正交性質,即
∫{-1,1} Pm(x)*Pn(x) dx = 0 if m≠ n

2009-01-31 18:04:48 補充:
又由於 deg(Pn(x) )=n,所以 { Pk(x)} , k=1,2,.... 是多項式函數空間的線性基底。特別地,對於任何多項式 g , deg(g)

2009-01-31 18:05:39 補充:
特別地,對於任何多項式 g , deg(g)




2009-01-31 18:06:24 補充:
特別地,對於任何多項式 g
deg(g)


2009-01-31 18:07:27 補充:
特別地,對於任何多項式 g , deg(g)少於n,g 是deg




2009-01-31 18:08:41 補充:
特別地,對於任何多項式 g , deg(g)少於n,g 是少於n次的Legendre polynomials的線性組合,所以當 deg(g)




2009-01-31 18:09:29 補充:
所以當 deg(g)少於n,有
∫{-1,1} Pn(x)*g(x) dx
= ∫{-1,1} Pn(x)*{deg


2009-01-31 18:12:09 補充:
特別地,對於任何多項式 g , deg(g)少於n,g 是deg少於n的Legendre polynomials的線性組合,所以當 deg(g)少於n,有
∫{-1,1} Pn(x)*g(x) dx
= ∫{-1,1} Pn(x)*{deg少於n的Legendre polynomials的線性組合} dx
= 0

2009-01-31 18:13:15 補充:
現在假設 g1(x) = (x -x2)(x -x3) ... (x -xn),有 deg(g1)少於n,所以
0 = ∫{-1,1} Pn(x)*g1(x) dx = w(1)(x1 -x2)(x1 -x3) ... (x1 -xn)Pn(x1)
由此可見 Pn(x1) = 0
同理,Pn(xk) = 0 , k= 1, 2, ..., n

2009-01-31 18:16:30 補充:
亂碼食字好討厭。
2009-01-31 6:33 am
這題鐵定是二月份的三條推薦問題之一!!!

2009-01-30 22:36:56 補充:
EMK, 加油, 期待你升上知識長的一天.
菩提兄與 maximal_ideal_space 亦如是.
And 我希望官方可以招請多兩個數學分類管理員, 負責不同範疇的數學知識.
2009-01-31 4:28 am
Gaussian Quadrature...
2009-01-31 3:45 am
設L(x)=(x^2 - 1)^n , 在 x1, x2, ..., xn處, L(x)之 n 階導數值=0

2009-01-30 20:31:25 補充:
To: EMK
Yes!
Plz prove it.

2009-01-31 13:52:50 補充:
∫[-1,1] P(x) dx = 0, and w(1),..,w(n)>0
並沒有 implies P(x1),..,P(xn)= 0
P(x)值可為負數喔!

2009-01-31 22:18:31 補充:
Very good!
亂碼食字已經習慣了!
這題也請幫忙一下:
http://hk.knowledge.yahoo.com/question/question?qid=7009011701265
2009-01-31 3:39 am
證明: x1, ... , xn 滿足 (x^2-1)^n之 n 階導數 = 0 MEANING?

2009-02-01 01:32:28 補充:
及時而出:

你段自我介紹和你的數學水平不成比例既


收錄日期: 2021-04-26 13:04:27
原文連結 [永久失效]:
https://hk.answers.yahoo.com/question/index?qid=20090130000051KK01778

檢視 Wayback Machine 備份