高教杯比賽論文-古塔變形預(yù)測(cè)模型-資料下載頁(yè)
2025-06-07 12:51本頁(yè)面
【正文】 2)。 z=sj(:,3)。 A=[x,y,ones(size(x))]。 a=A\z 附錄二 : 求解中心點(diǎn)坐標(biāo) %求中心點(diǎn)坐標(biāo)（以 1986 第一層中心為例） clear all clc % 球面參數(shù)為 1，平面參數(shù)為 2 %數(shù)據(jù)輸入： sj=[ ]。 x=sj(:,1)。 y=sj(:,2)。 z=sj(:,3)。 %擬合平面方程： A=[x,y,ones(size(x))]。 pingmian=A\z A2=pingmian(1,1)/pingmian(3,1) B2=pingmian(2,1)/pingmian(3,1) C2=1/pingmian(3,1) 2 %擬合球面方程 ： a=[sum(x.^2) sum(x.*y) sum(x.*z) sum(x) sum(x.*y) sum(y.^2) sum(y.*z) sum(y) sum(x.*z) sum(y.*z) sum(z.^2) sum(z) sum(x) sum(y) sum(z) 8] b=[sum(x.*(x.^2+y.^2+z.^2)) sum(y.*(x.^2+y.^2+z.^2)) sum(z.*(x.^2+y.^2+z.^2)) sum(x.^2+y.^2+z.^2)] p=inv(a)*b t=vpa(p) A1=t(1,1) B1=t(2,1) C1=t(3,1) D1=t(4,1) a=A1/2 b=B1/2 c=C1/2 R=(a^2+b^2+c^2D1)^ x1=a。 y1=b。 z1=c。 %求平面圓圓心 （ 即中心點(diǎn)坐標(biāo) ?。? M0=(2+A1*A2+B1*B2+C1*C2)/(2*(A2^2+B2^2+C2^2)) x2=A2*M0+A1/2 y2=B2*M0+B1/2 z2=C2*M0+C1/2 附錄三 ： 空間直線擬合 1996年 matlab 求解 %擬合空間直線（ 1996 年 為例） x=[ 3 ]39。 y=[ ]39。 z=[ ]39。 F=[z。1 1 1 1 1 1 1 1 1 1 1 1 1]。 M=F*F39。 N=F*x39。 O=F*y39。 A=(M\N)39。 B=(M\O)39。 x1=*z+。 y1=*z+。 z1=z。 plot3(x1,y1,z1,39。o39。,x,y,z,39。39。) 4 grid 方程為??? ??? ?? 74 26 845800 82 4375 66 87 826100 87 3194 zy zx 2021 年方程??? ??? ?? zy zx 2021 年的方程??? ??? ?? zy zx 附錄四 傾斜度計(jì)算 %1986 年直線的傾斜度 clear all clc n1=[1,0,]。 n2=[0 1 ]。 5 r=[0 0 1] s=cross(n1,n2) cita=acos(dot(s,r)/(norm(s)*norm(r)))*180/pi %1996 年直線的傾斜度 clear all clc n1=[1,0, ]。 n2=[0 1 ]。 r=[0 0 1] s=cross(n1,n2) cita=acos(dot(s,r)/(norm(s)*norm(r)))*180/pi %2021 年直線的傾斜度 clear all clc n1=[1,0, ]。 n2=[0 1 ]。 r=[0 0 1] s=cross(n1,n2) cita=acos(dot(s,r)/(norm(s)*norm(r)))*180/pi %2021 年直線的傾斜度 clear all clc n1=[1,0, ]。 n2=[0 1 ]。 r=[0 0 1] s=cross(n1,n2) cita=acos(dot(s,r)/(norm(s)*norm(r)))*180/pi 附錄五 ： 計(jì)算曲線的曲率 %1986 年 曲線曲率 clear y=[ 6 ]。 z=abs(66960.*y.^.*y+)./(1+(22320.*y.^.*y.^2+9.*).^2).^ 附錄六 : 1996 年第一層 8 個(gè)投影點(diǎn)移動(dòng)后的坐標(biāo) %1996 年第一層 8 個(gè)投影點(diǎn)移動(dòng)后的坐標(biāo) XN1=[ ]。 YN1=[ ]。 xm= ym= XOM=ones(8,1)*xm YOM=ones(8,1)*ym XN2=XN1+XOM YN2=YN1+YOM 7 ZUOBIAO=[XN2 YN2] %1996 年與 1986 年第一層第四個(gè)點(diǎn)與圓點(diǎn)連線的夾角 arfa OT1=[ ]。 OT2=[ ]。 s=OT1 r=OT2 arfa=acos(dot(s,r)/(norm(s)*norm(r)))*180/pi 附錄 7: 傾斜角角與年份的擬合 t=[1986 1996 2021 2021 ]39。 z=[ ]39。 a= b= c= 2062 x=1986:2021 y=b/(2*a)+(b^24*a.*(cx)).^(2*a) plot(t,z,39。o39。,x,y)
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1