1、精選優(yōu)質(zhì)文檔-傾情為你奉上 GPS衛(wèi)星位置計(jì)算程序設(shè)計(jì)報(bào)告學(xué) 院： 專業(yè)班級(jí)： 學(xué)生姓名： 學(xué)生學(xué)號(hào)： 指導(dǎo)老師： 設(shè)計(jì)目的 運(yùn)用計(jì)算機(jī)語(yǔ)言和課本內(nèi)容結(jié)合進(jìn)行編程，更多熟悉專業(yè)知識(shí)。設(shè)計(jì)工具M(jìn)atlab R2010b 設(shè)計(jì)代碼 %讀取數(shù)據(jù). fid,msg = fopen('單點(diǎn)定位實(shí)驗(yàn)數(shù)據(jù).txt','r'); arr_skip,count = fscanf(fid,'%s',5); arr_data = textscan(fid,'%d %f %f %f %f ',-1);satellite_num = arr_data1;d

2、is = arr_data2;x_coor = arr_data3;y_coor = arr_data4;z_coor = arr_data5;dis_app = 1:6;L = 0;0;0;0;0;0;B = L,L,L,L; X=0;Y=0;Z=0;for t = 1:10 for ii = 1:6 t1 = (x_coor(ii)-X)2; t2 = (y_coor(ii)-Y)2; t3 = (z_coor(ii)-Z)2; dis_app(ii) = sqrt(t1+t2+t3);end for ii = 1:6 B(ii,1) = (x_coor(ii)-X) / dis_app(

3、ii); B(ii,2) = (y_coor(ii)-Y) / dis_app(ii); B(ii,3) = (z_coor(ii)-Z) / dis_app(ii); B(ii,4) = -1; L(ii) = dis(ii)-dis_app(ii); enddx = -inv(B'*B)*B'*L; X= X+dx(1);Y = Y+dx(2);Z = Z+dx(3);end% 結(jié)果顯示str1=X;str2=Y;str3=Z; % 計(jì)算 L% L_l = atan(Y/X)L_l=asin(Y/sqrt(X2+Y2); L_b = atan(Z/sqrt(X2+Y2);%

4、 計(jì)算橢球的相關(guān)參數(shù)sqr_e = 0.013;for ii = 1:5 L_b1 = L_b; L_b = atan(Z+(X/cos(L_l)/cos(L_b1)*sqr_e*sin(L_b1)/sqrt(X2+Y2);end% 精度評(píng)定v = B*dx+L;sigma0 = sqrt(v'*v/2)t1 = sin(L_b);t2 = cos(L_b);t3 = sin(L_l);t4 = cos(L_l);r = -t1*t4,-t1*t3,t2;-t3,t4,0;t2*t4,t2*t3,t1;Qx = inv(B'*B);qx = Qx(:,1:3);qx = qx(

5、1:3,:);qb = r*qx*r;% 平面位置精度因子HDOP = sqrt(qb(1,1)/L_l+qb(2,2);Mh = HDOP*sigma0;% 高程精度因子VDOP = sqrt(qb(3,3);Mv = VDOP*sigma0;% 空間位置精度因子PDOP = sqrt(qx(1,1)+qx(2,2)+qx(3,3);Mp = PDOP*sigma0;% 接收機(jī)鐘差精度因子TDOP = sqrt(Qx(4,4);Mt = TDOP*sigma0;% 幾何精度因子GDOP = sqrt(Qx(1,1)+Qx(2,2)+Qx(3,3)+Qx(4,4);Mg = GDOP*sigma0;% 輸出結(jié)果str4=HDOP;str5=VDOP;str6=PDOP;str7=TDOP;str8=GDOP;str0=sigma0;str=str1,str2,str3,str4,str5,str6,str7,str8,str0;set(handles.edit1,'string',str);運(yùn)行界面課程設(shè)計(jì)感想 這次專業(yè)的最后一次實(shí)驗(yàn)所以拖得有點(diǎn)久，因?yàn)檫@次有矩陣計(jì)算的代碼matlab對(duì)矩陣的計(jì)算很方便，所以考慮使用了matlab來編程，對(duì)于matlab不