lzth.net
当前位置:首页 >> mAtlAB函数orth()翻译成C或C++ >>

mAtlAB函数orth()翻译成C或C++

这种应该是迭代求解的,你可以去翻翻矩阵论。好久没弄矩阵了,忘光了。。。

B = orth(A)返回A的正交基,B的每列正交, B'*B = eye(rank(A)). B的列数等于A的秩.

例子如下: >> s=[1,1,0;0,1,1;1,0,1] s = 1 1 0 0 1 1 1 0 1 >> [Q,R]=qr(s) Q = -0.7071 -0.4082 -0.5774 0 -0.8165 0.5774 -0.7071 0.4082 0.5774 R = -1.4142 -0.7071 -0.7071 0 -1.2247 -0.4082 0 0 1.1547 Q即为所求.

假设要生成n阶的单位矩阵,先用A=rand(n,n)生成n*n的随机矩阵(一般情况下都是满秩的),然后用B=orth(A)即可得到n*n的单位矩阵。命令就三行 n=5 %假定生成5阶单位矩阵 A=rand(n,n) B=orth(A)

正义命令orth处理后,如a=orth(b),a的列数是等于b的秩数,所以,当b缺秩时,a的列就与b的列数不相等了。从此可以看出,你程序中就出现了这种情况。

用函数eig可以直接求出一个矩阵的特征值和特征向量;按照问题的意思A*x=c*x,有:(cE-A)*x=0,也就是求该方程的一个基础解系,matlab中有一个函数orth可以求出一组标准正交基。至于有没有其他函数,我也不是很清楚,毕竟matlab的函数实在太多了。

我来回答你吧,我也是今天才明白的,搞了一个多星期。dbwavef('db3’)求的是db3的尺度滤波器。wfilter和orthfilt函数都可以求正交小波滤波器组。只是前者功能更强,还可以求双正交的。你把代码换成orthfilt(dbwavef('db3‘))。是不是发现和wfilt...

A a abs 绝对值, 模 acos 反余弦 acosh 反双曲余弦 acot 反余切 acoth 反双曲余切 acsc 反余割 acsch 反双曲余割 all 所有元素均非零则为真 alpha 透明控制 angle 相角 ans 最新表达式的运算结果 any 有非零元则为真 area 面域图 asec 反正割 as...

A a abs 绝对值、模、字符的ASCII码值 acos 反余弦 acosh 反双曲余弦 acot 反余切 acoth 反双曲余切 acsc 反余割 acsch 反双曲余割 align 启动图形对象几何位置排列工具 all 所有元素非零为真 angle 相角 ans 表达式计算结果的缺省变量名 any 所...

你好!目前是没有的,估计将来也不会有。因为如果一个计算可以被其他运算代替的简单运算而且又不是最基本的函数运算,那么是不可能为它设置特定的函数的。 你可以尝试e=a/norm(a) 或者e=a/(a*a')

网站首页 | 网站地图
All rights reserved Powered by www.lzth.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com