PDA

查看完整版本 : 求教各位Matlab达人


dong81206
2008-09-14, 10:07
:lovely:

我的问题是这样的:

给出
x1=0, x2=1/N ,..., xN=(N-1)/N

Kx1=exp{-x^2/a^2},
Kx2=exp{-(x-1/N)^2/a^2},
Kx3=exp{-(x-2/N)^2/a^2},
Kxi=exp{-(x-(i-1))^2/a^2}
...


内积
<Kxi,Kxj>=exp{-(i/N-j/N)^2/a^2}


Gram-Smicth正交化过程

g1=Kx1

g2=Kx2-<Kx2,g1> g1/<g1,g1>
.
.
.

gn=Kxn-<Kxn,g1> g1/<g1,g1>-......-<Kxn,gn-1> gn-1/<gn-1,gn-1>


要求计算 <gi,gi>, i=1,2,3,....

请教各位该怎么计算呢,正交化程序还好弄,但是那个内积怎么弄呢?

急盼回复!!!感谢!

mathjiang
2008-09-14, 14:56
最好用matlab语句写出来,比如,x=[0,1/N,...,(N-1)/N]
这样更好些

剑如虹
2008-09-14, 15:33
我以为这样好些:

x=0:1/N:(N-1)/N;

hitzhjtopku
2008-09-15, 12:02
循环累加不就成了

dong81206
2008-09-15, 15:48
但是内积
<Kxi,Kxj>=exp{-(i/N-j/N)^2/a^2}
必须编写一个函数的吧
怎么写出这个函数?

求教各位大侠!

正交化的过程的语句倒是好写的。

mathjiang
2008-09-16, 17:58
两个同维向量的内积:x*y'
x,y同为列向量或同为行向量