matlab 用命令orth把向量组正交规范化

1、在命令行窗口中输入“r = [ 1 2 3 4 ]tr = r'”，其中单引号（'）是用于转置 *** 作的。

2、按回车键之后，可以看到将行向量r转置成为了列向量tr。

3、也可以将列向量转置成为行向量，在命令行窗口中输入“v = [1234]tv = v'”。

4、按回车键之后，可以看到将列向量v转置成为了行向量tv。

5、矩阵转置的 *** 作也是通过单引号（'）来实现的。

6、按回车键之后，可以看到将a矩阵的每列转置成了b矩阵的每行，a矩阵的每行转置成了b矩阵的每列。

wp=35

ws=60

rp=0.5

rs=40

Fs=360%%设置各个参数

[N,Wn] = buttord(wp/(Fs/2),ws/(Fs/2),rp,rs)%%wp/(Fs/2)通带拐点频率；ws/(Fs/2)阻带拐点频率；rp通带内最大允许衰减度，以分贝表示；rs从通带到阻带的最小衰减度；N为最小阶数；Wn为截止频率；

[b,a]=butter(N,Wn)%%设计滤波器，阶数为N，截止频率为Wn，b，a为最终滤波器的系数

I=filter(b,a,I)%%得到滤波器的系数b和a，并对数据I进行滤波

for i=1:COUNT%%COUNT应该是个数吧，还有STEP

[maxx lag]=max(abs(I((i-1)*STEP+1:i*STEP)))%%在I中取STEP大小的窗口，并计算窗口内的最大值，maxx是最大值，lag是最大值所在的列

if maxx>=0.4

n=n+1%%如果最大值大于等于0.4

maxium(n)=maxx%%把最大值赋给maxium的第n个元素

t(n)=(i-1)*STEP+lag%%把最大值在I中所在的列赋给t的第n个元素

f(n)=I((i-1)*STEP+lag)%%在I中取窗口，并把窗口内的最大值赋给f的第n个元素

end

end

m=1

tt(1)=t(1)

ff(1)=f(1)%%将f的第一个元素赋给ff的第一个元素

for j=2:n%%n为上面的循环结束后最终的n

if t(j)-t(j-1)<=80%%t中存的是最大值对应的列，如果相邻两列的差值小于等于80

if maxium(j)>=maxium(j-1)

tt(m)=t(j)

ff(m)=f(j) %%maxium中是最大值，如果第j个比第j-1个大，那么就把第j个最大值赋给tt(1),ff(1)

end

else%%如果t中存的是最大值对应的列，相邻两列的差值不小于等于80

m=m+1%%m更新

tt(m)=t(j)

ff(m)=f(j)%%将第j个最大值赋给tt(2),ff(2)

end

end

plot(I)

hold on

stem(tt,ff,'r*')%%画tt和ff的离散图

orth函数可以把一个矩阵正交化。但不知道它用的什么算法

另外，可以用QR分解直线正交化，函数是qr，可参考：

http://www.ilovematlab.cn/thread-168827-1-1.html

---