linefunction.m

这是用于工程测量进行大坝变形分析的软件

function result=linefunction(x0,y0,temp) 
global edit2 edit3 edit4 edit5 fit2;
[a,b,cigema_squre_line]=line_nihe(x0,y0);
plotx=[x0(1):0.1:x0(length(x0))];
plot(plotx,a.*plotx+b,'r');
hold on
plot(x0,y0,'b*');
if nargin==2
	title('这是整个滑坡在所求方向的直线拟合图')
elseif nargin==3
	title(['这是第' num2str(temp) '个点在所求方向的直线拟合图'])
end
xlabel(['从' get(edit2,'string') '年' get(edit4,'string')  '月到' get(edit3,'string') '年' get(edit5,'string')  '月']);
ylabel('在滑坡方向的位移量/米')
hold off
%写函数参数------------------------------------------------------------------------------------------------------        
if nargin==2 
    fprintf(fit2,'以下是整个滑坡直线拟合结果\n') ;
elseif nargin==3
    fprintf(fit2,'以下是第%d点直线拟合结果\n',temp) ;
end
fprintf(fit2,'其方差σ=%f\n',cigema_squre_line);
fprintf(fit2,[num2str(a) 'x+'  num2str(b) '\n']);
fprintf(fit2,'各参与拟合点的信号及噪声为\n') ;
for i=1:length(x0)      
    y=a*x0(i)+b;
    fprintf(fit2,['y' num2str(i) '的信号为' num2str(y*1000) 'MM\t' '噪声为' num2str((y0(i)-y)*1000) 'MM\n']);
end
return