24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cosmology

金虫 (正式写手)

应助: 15 (小学生)
金币: 1411.3
红花: 1
帖子: 418
在线: 48.6小时
虫号: 230145
注册: 2006-03-26
性别: GG
专业: 基础物理学

[交流] 60金币请教高人一个正交性问题

见附件,愿献60金币求证明过程

[ Last edited by cosmology on 2010-1-21 at 08:28 ]

回复此楼

» 猜你喜欢

1楼 2010-01-20 15:46:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9683.7
散金: 89
红花: 19
帖子: 3577
在线: 372.3小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
Doctorcbw(金币+40):按照楼主的要求还有你的热情应答，奖励40bb给你！ 1-24 20:31
Doctorcbw(金币+20):按照楼主的要求还有你的热情应答，奖励20bb给你！ 1-24 20:31

好吧，我试试。
结果
>> 2*int('sin(2*n*pi*x)*cos((2*m-1)*pi*x)',0,1)
ans =
-8*(n-sin(n*pi)*cos(n*pi)*sin(pi*m)*cos(pi*m)+2*n*cos(n*pi)^2*cos(pi*m)^2-n*cos(n*pi)^2-n*cos(pi*m)^2+2*m*sin(n*pi)*cos(n*pi)*sin(pi*m)*cos(pi*m))/pi/(2*n+2*m-1)/(-2*n+2*m-1)
整理一下（因n,m是整数,sin(k*pi)=0,cos(k*pi)=(-1)^k）得到

A(n,m)=-8*n/pi/(2*n+2*m-1)/(-2*n+2*m-1)
                                 8*n
      =  ——————————————————
                  [(2n)^2-(2*m-1)^2]*pi

从而sum(A(n,m)*A(l,m)

      64*l*n                         1
  =  —————— *sum( —————————————————————————  )
   (pi)^2                [(2n)^2-(2*m-1)^2] * [(2n)^2-(2*m-1)^2]


当n等于l 时，结果显然非零。（是否等于1无关紧要，只要将其归一化就可以了，验证了一下几个值，貌似还真是1 ！主要讨论n和l 不相等的情况）
当n不等于l时，可以将其化简
sum(A(n,m)*A(l,m)

      16*l*n                               1                            1
= ———————— *sum( —————————— -  ————————————  )
(pi)^2*(n^2-l^2)    [(2*m-1)^2-(2n)^2]       [(2*m-1)^2-(2l)^2]

而
                  1
sum( ——————————)=0          (对任意k)
   [(2*m-1)^2-(2k)^2]

从而得到正交性

[ Last edited by bluesine on 2010-1-23 at 14:15 ]