关于二次曲面拟合，matlab，这个程序哪位大神给我讲讲啊，GPS高程拟合 15
[P1,F1]=uigetfile('*.txt','打开GPS水准点平面坐标文件');fnal1=strcat(F1,P1);fpath=P1;fid1=fopen(f...
[P1,F1]=uigetfile('*.txt','打开GPS水准点平面坐标文件');
fnal1=strcat(F1,P1);fpath=P1;
fid1=fopen( fnal1,'r');
lineNum=0;
zbjz=[];
format long e
while (~feof(fid1))
curline = fgetl(fid1);
lineNum = lineNum + 1;
zbjz(lineNum,1)=1;
zbjz(lineNum,2)=str2num(curline(5:16));
zbjz(lineNum,3)=str2num(curline(18:28));
zbjz(lineNum,4)=str2num(curline(30:36));
end
zbjz;
Y=zbjz(:,2);
X=zbjz(:,3);
Z=zbjz(:,4);
o=zbjz(:,1);
X2=X.*X;
Y2=Y.*Y;
XY=X.*Y;
B=[o,X,Y,X2,Y2,XY];
[P4,F4]=uigetfile('*.txt','打开GPS水准点大地高文件');
fnal4=strcat(F4,P4);fpath4=P4;
fid4=fopen( fnal4,'r');
LineNum=0;
ddg1=[];
while (~feof(fid4))
LineNum = LineNum + 1;
urline = fgetl(fid4);
ddg1(LineNum,1)=str2num(urline(5:10));
end
ddg1;
gcyc1=ddg1-Z;
%---------------求取模型参数-----------%
%a0,a1,a2,a3 ,a4,a5为模型参数
A=pinv(B'*B)*B'*gcyc1
[miu,M,cc]=pdjd(A,ddg1); 展开
fnal1=strcat(F1,P1);fpath=P1;
fid1=fopen( fnal1,'r');
lineNum=0;
zbjz=[];
format long e
while (~feof(fid1))
curline = fgetl(fid1);
lineNum = lineNum + 1;
zbjz(lineNum,1)=1;
zbjz(lineNum,2)=str2num(curline(5:16));
zbjz(lineNum,3)=str2num(curline(18:28));
zbjz(lineNum,4)=str2num(curline(30:36));
end
zbjz;
Y=zbjz(:,2);
X=zbjz(:,3);
Z=zbjz(:,4);
o=zbjz(:,1);
X2=X.*X;
Y2=Y.*Y;
XY=X.*Y;
B=[o,X,Y,X2,Y2,XY];
[P4,F4]=uigetfile('*.txt','打开GPS水准点大地高文件');
fnal4=strcat(F4,P4);fpath4=P4;
fid4=fopen( fnal4,'r');
LineNum=0;
ddg1=[];
while (~feof(fid4))
LineNum = LineNum + 1;
urline = fgetl(fid4);
ddg1(LineNum,1)=str2num(urline(5:10));
end
ddg1;
gcyc1=ddg1-Z;
%---------------求取模型参数-----------%
%a0,a1,a2,a3 ,a4,a5为模型参数
A=pinv(B'*B)*B'*gcyc1
[miu,M,cc]=pdjd(A,ddg1); 展开
2个回答
展开全部
请问，这段程序是采用的何种方法？三次样条？还是？
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
Sievers分析仪
2025-07-02 广告
是的。传统上，对于符合要求的内毒素检测，最终用户必须从标准内毒素库存瓶中构建至少一式两份三点标准曲线；必须有重复的阴性控制；每个样品和PPC必须一式两份。有了Sievers Eclipse内毒素检测仪，这些步骤可以通过使用预嵌入的内毒素标准...
点击进入详情页
本回答由Sievers分析仪提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
