求教各位Matlab达人

[dong81206]

: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]

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

[剑如虹]

我以为这样好些:

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

[hitzhjtopku]

循环累加不就成了

[dong81206]

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

求教各位大侠！

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

[mathjiang]

两个同维向量的内积：x*y'
x,y同为列向量或同为行向量