摘要： 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 为行... 阅读全文

摘要： 每种迭代方法的程序都返回迭代步数，可以很方便的进行对比雅克比方法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... 阅读全文

摘要： 平方根法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... 阅读全文

摘要： 下面的程序代码相当规范，不是我写的，记得好像是在MATLAB网站的交流区下载的。 每个*.img文件对应一个*.hdr的头文件，程序首先从头文件里得到影像的行列和波段数，然后按照二进制方式读取，并按照行列波段数重排为多维矩阵。 function [image,p,t]=freadenvi(fname); % freadenvi - read envi image (V. Gui... 阅读全文

摘要： MNF变化在遥感图像处理中应该很有名气吧，这里稍微偷了下懒，没进行主成分变换的标准化，如果有朋友看的时候仔细了。之所以没有进行标准化是因为看到某老师的PPT说其实无所谓，在数学上是等价的。如果标准化就可以用相关矩阵代替协方差矩阵。我概率学的真是不好，这里有点儿晕。没有标准化，直接用协方差矩阵计算了。 function T=mnf_plus(n1,n2)%————————————————第一步% n... 阅读全文