摘要： 非常类似前面的三弯矩法，这里的sanzhj函数和intersanzhj作用相当于前面的sanwanj和intersanwj，追赶法程序通用，代码如下。%%%%%%%%%%%%%%%%%%%function [newu,w,newv,d]=sanzhj(x,y,x0,y0,y1a,y1b)% 三转角样条插值% 将插值点分两次输入，x0 y0 单独输入% 边值条件a的一阶导数 y1a 和b的一阶导数... 阅读全文

摘要： 首先说这个程序并不完善，为了实现通用（1，2，…，n）格式解题，以及为调用追赶法程序，没有针对节点数在三个以下的情况进行分类讨论。希望能有朋友给出更好的方法。首先，通过函数 sanwanj得到方程的系数矩阵，即追赶法方程的四个向量参数，接下来调用追赶法（在intersanwj函数中），得到三次样条分段函数系数因子，然后进行多项式合并得到分段函数的解析式，程序最后部分通过判断输入值的区间自动选择对应... 阅读全文

摘要： 从RTK到网络RTK 一、 引言随着全球定位系统(GPS)技术的快速发展，RTK技术也日益成熟。实时动态（Real Time Kinematic） 测量系统，是GPS测量技术与数据传输技术的结合，是GPS测量技术中的一个新突破。RTK的出现使得GPS工程应用成为可能，其速度快，成本相对低廉。然而，它对基站无线电信号的依赖，导致有效距离较小，且容易受到地形和建筑物的限制，... 阅读全文

摘要： 通过实例验证加速对简单迭代的改进效果，注意到加速法使在简单迭代情况下不收敛的迭代方程也能收敛，这也是迭代法的特色哦。 实例： f(x)=x^3+10*x-20 精确到6位小数迭代方程： φ（1） 20/(x.^2+10)；φ（2） x.^3+11*x-20；φ（3） 2-x.^3/10 简单迭代：%简单迭代法 use structure i=2;N=100... 阅读全文

摘要： 文件用UltraEdith32打开如下： 2.11 NAVIGATION DATA G (GPS) RINEX VERSION / TYPEDAT2RINW 3.10 001 1 16APR08 9:51:13 PGM / RUN BY / DATE ... 阅读全文

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