摘要:
function N=wepnewton(x,y,xi)%Newton基本插值方法n=length(x);m=length(y);if m~=n error('x or y 输入有误,再来');endA=zeros(n);Z=1.0;A(:,1)=y;N=A(1,1);for k=2:n % k 为列标 for i=k:n % i 为行...
阅读全文
posted @ 2008-12-13 21:32
秋波渡
阅读(1413)
推荐(0)
编辑
摘要:
每种迭代方法的程序都返回迭代步数,可以很方便的进行对比雅克比方法function [x,n]=jacob(A,b,x0,eps,M)%雅可比迭代%x 解 n 达到所需精度实际用的步数%对输入单数进行默认设置 x0初始值 eps精度 M限制步数if nargin==3 eps=1.0e-6; M=200;elseif nargin==4 M=200;elseif nargin=e...
阅读全文
posted @ 2008-12-13 21:28
秋波渡
阅读(1689)
推荐(0)
编辑
摘要:
平方根法function [x]=pingfg(A,b)%乔累斯基分解[n,n]=size(A);L=zeros(n,n);%实际上不用为 L 申请空间,使用 A 即可L(1,1)=sqrt(A(1,1));for k=2:n L(k,1)=A(k,1)/L(1,1);endfor k=2:n-1 L(k,k)=sqrt(A(k,k)-sum(L(k,1:k-1).^2)); f...
阅读全文
posted @ 2008-12-13 21:24
秋波渡
阅读(6083)
推荐(0)
编辑
摘要:
下面的程序代码相当规范,不是我写的,记得好像是在MATLAB网站的交流区下载的。 每个*.img文件对应一个*.hdr的头文件,程序首先从头文件里得到影像的行列和波段数,然后按照二进制方式读取,并按照行列波段数重排为多维矩阵。 function [image,p,t]=freadenvi(fname); % freadenvi - read envi image (V. Gui...
阅读全文
posted @ 2008-12-13 15:43
秋波渡
阅读(6477)
推荐(0)
编辑
摘要:
MNF变化在遥感图像处理中应该很有名气吧,这里稍微偷了下懒,没进行主成分变换的标准化,如果有朋友看的时候仔细了。之所以没有进行标准化是因为看到某老师的PPT说其实无所谓,在数学上是等价的。如果标准化就可以用相关矩阵代替协方差矩阵。我概率学的真是不好,这里有点儿晕。没有标准化,直接用协方差矩阵计算了。 function T=mnf_plus(n1,n2)%————————————————第一步% n...
阅读全文
posted @ 2008-12-13 15:21
秋波渡
阅读(2918)
推荐(0)
编辑