基于Matlab的导线网平差程序设计

软件无线电实验matlab程序

附录 (1) 、SDR 低通采样理论 %parameters fs = 5e4;%采样频率注意 或 2 倍以上 f = 5e3;%信号的频率 N = 1024; n = 1:1024; t= n/fs; deta_t = 1/fs;% 采样间隔 %signals sn = sin(2*pi*f*t);% 时域采样后的信号w = hanning(N);% 加汉宁窗sn1 = sn.*w';% 加窗减少频率泄露 Sn_fft = fft(sn1);% 频域信号 Sn = abs(Sn_fft); %figures figure(1); plot(t,sn); title('s(n)的时域波形'); grid on; xlabel('t'); xlim([0 0.021]); ylabel('幅度'); figure(2); xax = [-N/2:N/2-1]/N*fs/1000;% 将f 轴单位变为kHz plot(xax,20*log10(fftshift(Sn))); title('s(n)的频域波形'); grid on; xlabel('f(kHz)'); ylabel('幅度(dB)'); ylim([-200 100]); (2) 、SDR 带通采样理论 %parameters f = 1.3e3;% 基频 fo = 100e6;% 载波频率 N = 1024; fs = 4e3;%采样频率 t= 0:0.001:2; st = cos(2*pi*(f+fo)*n/fs);% 带通抽样信号 St1 = fft(st); St = abs(St1); st0 = cos(2*pi*(f+fo)*t); St0 = fft(st0); St2 = abs(St0); figure(1); plot(n,st); title('带通抽样信号时域图像’)； xlabel('n'); xlim([0 1026]); ylabel('幅度'); grid on; figure(2); plot(t,st0); title('带通原始信号时域图像'); xlabel('t'); ylabel('幅度'); grid on; figure(3); xax =[-N/2:N/2-1]/N*fs/1000; plot(xax,20*log10(fftshift(St))); title('带通信号抽样后频域图像’)； xlabel('f(kHz)'); ylabel('幅度(dB)'); grid on; (3) 、频率调制信号 %频率调制 %parameters fs = 1e3;%抽样频率 ts = 1/fs;%采样率 N = 1024; n = 0:1:N-1; t = n/fs; kf = 100.1; A = 100; Aw = 10; fc = 10000;% 载波频率fm = 10;% 调制信号频率mf = kf*Aw/(2*pi*fm); %signals s1 A*cos(2*pi*fc*t).*cos(mf*sin(2*pi*fm*t)); s2 fs 至少是f 的2 倍n = 1:1024;

闭合导线平差计算步骤

闭合导线平差计算步骤： 1、绘制计算草图。在图上填写已知数据和观测数据。 2、角度闭合差的计算与调整 （1）计算闭合差： （2）计算限差：（图根级） （3）若在限差内，则按平均分配原则，计算改正数： （4）计算改正后新的角值： 3、按新的角值，推算各边坐标方位角。 4、按坐标正算公式，计算各边坐标增量。 5、坐标增量闭合差的计算与调整 （1）计算坐标增量闭合差。有： 导线全长闭合差： 导线全长相对闭合差: （2）分配坐标增量闭合差 若 K<1/2000 （图根级），则将、以相反符号，按边长成正比分配到各坐标增量上去。并计算改正后的坐标增量。

6、坐标计算 根据起始点的已知坐标和经改正的新的坐标增量，来依次计算各导线点的坐标。 [ 例题 ] 如图所示闭合导线，试计算各导线点的坐标。 计算表格见下图：

闭合水准路线内业计算的步骤： (1) 填写观测数据 (2) 计算高差闭合差 h f =∑h ，若h f ≤容h f 时，说明符合精度要求，可以进行高差闭合差的调整；否则，将重新进行观测。 (3) 调整高差闭合差 各段高差改正数： i h i i h i L L f V n n f V ·· ∑-= ∑-= 或 各段改正高差： i i i V h h +=改 (4) 计算待定点的高程 闭合差(ｆh ) 水准路线中各点间高差的代数和应等于两已知水准点间的高差。若不等两者之差称为闭合差 高差闭合差的计算 .支水准路线闭合差的计算方法 .附合水准路线闭合差的计算方法 .闭合水准路线闭合差的计算方法 高差闭合差容许值 （n 为测站数，适合山地） (L 为测段长度，以公里为单位，适合平地) 水准测量中，消除闭合差的原则一般按距离或测站数成正比地改正各段的观测高差

附合闭合导线平差DOC

测量班冬季培训复习题纲 一．填空（每空0分，共计0分） 1．水准仪的结构是由（望远镜、水准器、基座）构成的。 2. 测量工作中常用（方位角）来表示直线的方向。 3. 地面点到大地水准面的铅垂距离叫(绝对高程). 4. 经纬仪导线外业测量工作的主要内容是（选点并设立标志，距离测量，角度测量）。 5. 水准测量中,水准尺的竖立应( 垂直) 6. 设A为后视点，B为前视点，A点高程为35.712m,后视读数为0.983 m ,前视读数为1.149 m ,则 A.B两点的高差是(—0.166 m ). 7. 在测区内,选定一些具有控制意义和作用的点子称为(控制点). 8. 转点在水准测量中起传递(高程)的作用. 9. 距离测量的基本单位是(米). 10. 方位角的角值范围是(0°--360°). 11．将经纬仪安置在O点，盘左照准左测目标A点，水平盘读数为0°01′30″，顺时针方向瞄准B点，水平盘读数为68°07′12″，则水平夹角为（68°05′42″） 12．地面点的空间位置是由坐标和高程决定的． 13．水准测量的检核方法有测站检核和路线长度检核等． 14. 测量的基本工作包括（.测角、测边、测高差）。 15. 水准仪的粗略整平是通过调节(脚螺旋)来实现的. 16. 竖直角的角值范围是(0°--90°). 17. 在调节水准仪粗平时,要求气泡移动的方向与左手大拇指转动脚螺旋的方向(相同). 18．控制导线的布设形式一般为闭合导线和附合导线。 19．水准测量测站校核的方法有改变仪器高法、双面尺法。 20．调节目镜螺旋，可以调节十字丝清晰。 21．从基本方向的北端起，顺时针方向到某一直线的水平夹角，称为该直线的（方位角） 22．经纬仪导线最终计算的目的是得到控制点的（坐标） 23．水准仪的圆水准器气泡居中，表明仪器达到粗平，长水准管气泡居中，表明仪器达到精平。 24．地面上一条直线的正反方位角之间相差（ 180 ）度。 25．测量地面两点高差时，尺垫安置在（转点）上。 26．从一个已知的水准点出发，沿途经过各点，最后到达另一个已知的水准点上，这样的水准路线是（附合水准路线） 27．当已知两点的坐标，求两点的平距和方位角的方法称为（坐标反算问题） 28．用经纬仪测定水平角时，应尽量瞄准目标的（底部） 29．水准面的特性是(处处与铅垂线相垂直) 30．我国常用的两个大地坐标系是(1954年北京坐标系)和(1980年国家大地坐标系) 31．测量工作在实地要测量的三个基本要素都是：(距离,)( 角度,) ( 高程) 32．测量工作必须遵循的第一条基本原则是(“从整体到局部”，“先控制后碎部”)的原则 33．进行高程测量的主要方法有(水准测量)和(三角高程测量) 34．常用的导线的形式(闭合导线)(附合导线)(.支导线) 二．简答题(共计0分) 1．水准测量的原理是什么?（0分） 答:利用水准仪提供的一条水平视线,借助于水准尺来测定地面点的高差,然后再根据已知点的高程去推求待定点的高程. 2.在测定地面上两点高差时，为何尽量将水准仪安置在地面两点的中点处？ 在测定地面上两点高差时，将水准仪安置在地面两点的中点处，可以消除水准管轴与视准轴不平行的

基于MATLAB的径向基网络源程序

%一维输入，一维输出，逼近效果很好！ 1.基于聚类的RBF 网设计算法 SamNum = 100; % 总样本数 TestSamNum = 101; % 测试样本数 InDim = 1; % 样本输入维数 ClusterNum = 10; % 隐节点数，即聚类样本数 Overlap = 1.0; % 隐节点重叠系数 % 根据目标函数获得样本输入输出 rand('state',sum(100*clock)) NoiseVar = 0.1; Noise = NoiseVar*randn(1,SamNum); SamIn = 8*rand(1,SamNum)-4; SamOutNoNoise = 1.1*(1-SamIn+2*SamIn.^2).*exp(-SamIn.^2/2); SamOut = SamOutNoNoise + Noise; TestSamIn = -4:0.08:4; TestSamOut = 1.1*(1-TestSamIn+2*TestSamIn.^2).*exp(-TestSamIn.^2/2); figure hold on grid plot(SamIn,SamOut,'k+') plot(TestSamIn,TestSamOut,'k--') xlabel('Input x'); ylabel('Output y'); Centers = SamIn(:,1:ClusterNum); NumberInClusters = zeros(ClusterNum,1); % 各类中的样本数，初始化为零IndexInClusters = zeros(ClusterNum,SamNum); % 各类所含样本的索引号while 1, NumberInClusters = zeros(ClusterNum,1); % 各类中的样本数，初始化为零IndexInClusters = zeros(ClusterNum,SamNum); % 各类所含样本的索引号 % 按最小距离原则对所有样本进行分类 for i = 1:SamNum AllDistance = dist(Centers',SamIn(:,i)); [MinDist,Pos] = min(AllDistance); NumberInClusters(Pos) = NumberInClusters(Pos) + 1; IndexInClusters(Pos,NumberInClusters(Pos)) = i; end % 保存旧的聚类中心

5800导线平差程序

5800计算器导线平差程序 一、程序用途及使用范围 本程序适用于一般导线复测平差计算，利用左角复测复合导线、闭合导线的平差时可直接使用。复合导线平差时输入起始和终止边两个方位角，如果将终止边方位角输入为0，程序会自动转入闭合导线平差界面。如果想用观测右角平差时，只需将程序中带下划线的地方按使用说明稍加改动即可。 二、源程序清单：1、JDPCA（文件名称） 0→M:0→I:0→K:“AFWJ”?A:“BFWJ”?B:“CZS”?→N:“JDRXBHCA”:24√ˉ(N)→P◣“∑(ZJ)”?→C:If B≠0 Then A-B-180N+C→D: EIseC-180(N-2)→D:IfEnd: “JDBHCA.f=”:D ?DMS◣PGoto 4:“JDGZ=”:-D÷N→F:F?DMS◣LbI 1:“CJ”?→P:B=0 And P=0＝>Goto3:“JGH=”:P+F→J:J?DMS ◣A+J-180→E:If E>0 And E<360:ThenE:IfEnd:If E>360:ThenE-360→E:IfEnd: :If E<360:ThenE+360→E:IfEnd: “FWJ=”:E?DMS◣“L”?→L:L=0 And B≠0＝>Goto3:“XO=”:Lcos(E)→X◣“YO=”:Lsin(E)→Y◣M+L→M:I+X→I:K+Y→K:E→A:Goto1:LbI3:“∑(L)=”:M◣“∑(XO)=”:I◣“∑(YO)=”:K◣Prog“DXPCA”：LbI4 DXPCA（文件名称，可单独运行） “AX”?→A:“A Y”?→Z:If B≠0:Then “BX”?→C:“BY”?→D:“FX=”:I-C+A→F◣“FY=”:K-D+Z→W◣EIse “FX=”:I

基于Matlab的导线网坐标计算设计

基于Matlab的导线网坐标计算设计

东华理工大学长江学院毕业设计论文题目：基于Matlab的导线网坐标计算 English Title：Traverse Network Coordinate Calculation Based On Matlab 学生姓名： 申请学位门类：工学学士 专业：测绘工程 系别：测绘工程系

摘要 导线计算是在所有测量工作中经常遇见的问题之一，同时导线计算的方法也有很多种，本文主要是利用简单易懂的Matlab对附合导线、闭合导线和支导线进行相应的平差计算。文章首先介绍了附合导线、闭合导线、支导线基本概念和计算方法，其次利用Matlab计算机编程语言对三种导线的计算进行编程实现；最后通过实例验证，本文利用Matlab编写的程序正确，通过输入边长和角度，能够快速的得到各控制点的准确坐标。 关键词：Matlab；导线计算；精度评价；计算机编程

ABSTRACT Traverse calculated in all measurements often met one of the problems,Also there are many kinds of traverse the way.This article mainly is to use simple Matlab to connecting traverse、closed traverse、spur traverse to the corresponding adjustment calculation The article first introduces the connecting traverse、closed traverse、spur traverse is basic concept and calculation method.Secondly using Matlab computer programming language was realized by programming calculation of three traverse.At last,through example validation, In this paper, using the Matlab program, right through the input variable length and Angle, able to quickly get the accurate coordinates of each control point. Key words：Matlab; Traverse Calculated ; Precision Evaluation; Computer Programming

导线平差计算

导线平差计算 1 简介 闭合导线和附合导线是长输管道站场和穿跨越测量常用的控制手段，其优点是可以同时完成平面和高程控制测量。导线平差原理请查阅相关文献。不同平差软件的平差方法步骤基本相同，本文件基于南方平差易软件平台介绍导线（闭合导线、附合导线是最简单的导线控制网）平差的操作方法。 2 规范性引用文件 下列文件中的条款通过本标准的引用而成为本标准的条款。 《长距离输油输气管道测量规范》（SY/T 0055-2003） 《工程测量规范》（GB 50026-2007） 3 操作步骤 （1）录入数据 录入数据是将导线测量数据录入平差软件。可以采用手工或文件方式录入（建议采用后者，选菜单“文件/打开”）。其数据格式如下： [NET] 控制网信息 [PARA] 控制网参数 [STATION]坐标和高程信息（11表示高程已知，如果无坐标则无法在平差易中看到和输出地图）[OBSER] 观测的转角、平距、高差等信息 下图为导入数据窗口： 图3-1 导入数据窗口 （2）坐标推算（F3）

选菜单“平差/推算坐标”，根据已知条件（测站点信息和观测信息）推算出待测点的近似坐标。为构建动态网图和导线平差作基础。 （3）概算 选菜单“平差/选择概算”→配置概算参数→输出概算结果。下图为“选择概算”的配置参数窗口： 图3-2 配置概算参数 （4）调整观测数据 将概算结果调整到输入的观测数据中，重新导入。 （5）计算方案的选择 对于同时包含了平面数据和高程数据的导线, 一般处理过程应为：先进行平面处理, 然后在高程处理时软件会使用已经较为准确的平面数据（如距离等）来处理高程数据。对精度要求很高的平面高程混合平差，您也可以在平面和高程处理间多次切换，迭代出精确的结果（但建议平面和高程分开了平差）。 针对导线平差，需要设置中误差及仪器参数、高程平差参数、限差及等级内容。 选菜单“平差/平差方案”即可进行参数的设置，如下图：

软件无线电实验matlab程序

附录 (1)、SDR低通采样理论 %parameters fs = 5e4;%采样频率注意fs至少是f的2倍或2倍以上 f = 5e3;%信号的频率 N = 1024; n = 1:1024; t= n/fs; deta_t = 1/fs;%采样间隔 %signals sn = sin(2*pi*f*t);%时域采样后的信号 w = hanning(N);%加汉宁窗 sn1 = sn.*w';%加窗减少频率泄露 Sn_fft = fft(sn1);%频域信号 Sn = abs(Sn_fft); %figures figure(1); plot(t,sn); title('s(n)的时域波形'); grid on; xlabel('t'); xlim([0 0.021]); ylabel('幅度'); figure(2); xax = [-N/2:N/2-1]/N*fs/1000;%将f轴单位变为kHz plot(xax,20*log10(fftshift(Sn))); title('s(n)的频域波形'); grid on; xlabel('f(kHz)'); ylabel('幅度(dB)'); ylim([-200 100]); (2)、SDR带通采样理论 %parameters f = 1.3e3;%基频 fo = 100e6;%载波频率 N = 1024; n = 1:1024; fs = 4e3;%采样频率 t= 0:0.001:2; st = cos(2*pi*(f+fo)*n/fs);%带通抽样信号St1 = fft(st); St = abs(St1); st0 = cos(2*pi*(f+fo)*t); St0 = fft(st0); St2 = abs(St0); figure(1); plot(n,st); title('带通抽样信号时域图像'); xlabel('n'); xlim([0 1026]); ylabel('幅度'); grid on; figure(2); plot(t,st0); title('带通原始信号时域图像'); xlabel('t'); ylabel('幅度'); grid on; figure(3); xax =[-N/2:N/2-1]/N*fs/1000; plot(xax,20*log10(fftshift(St))); title('带通信号抽样后频域图像'); xlabel('f(kHz)'); ylabel('幅度(dB)'); grid on; (3)、频率调制信号 %频率调制 %parameters fs = 1e3;%抽样频率 ts = 1/fs;%采样率 N = 1024; n = 0:1:N-1; t = n/fs; kf = 100.1; A = 100; Aw = 10;

闭合导线平差计算步骤

闭合导线平差计算步骤: 1、绘制计算草图。在图上填写已知数据和观测数据。 2、角度闭合差的计算与调整 Z 鬲-ZAe =（角 + 腐+、，+ 岗）- （〃-2）180°（1）计算闭合差：二 （2）计算限差： @二±4°扁（图根级） （3）若在限差，则按平均分配原则，计算改正数：尸 （4）计算改正后新的角值：禹=小弓 3、按新的角值，推算各边坐标方位角。 4、按坐标正算公式，计算各边坐标增量。 5、坐标增量闭合差的计算与调整 （1）计算坐标增量闭合差。有： 导线全长闭合差：f = 导线全长相对闭合差：# （2）分配坐标增量闭合差 若K< 1/2000 （图根级），则将兀、5以相反符号，按边长成正比分配到各坐标增量上去。并计算改正后的坐标增量。

6、坐标计算 根据起始点的已知坐标和经改正的新的坐标增量，来依次计算各导线点的坐标。［例题］如图所示闭合导线，试计算各导线点的坐标。 计算表格见下图:

(1)填写观测数据 (2) 计算高差闭合差 fh = ￡h,若hwh 容 时，说明符合精度要求，可以进行高差 闭合差的调整；否则，将重新进行观测。 (3) 调整高差闭合差 各段高差改正数: V,- 〃 或 y 乙 ， XL 1 各段改正高差： 九改=h i + v i (4) 计算待定点的高程 闭合差(fh ) 水准路线中各点间高差的代数和应等于两已知水准点间的高差。若不等两者之差称为闭 合差 高差闭合差的计算 fh = 22/ +2Z 知 .支水准路线闭合差的计算方法(或名= .附合水准路线闭合差的计算方法 九- ￡编 .闭合水准路线闭合差的计算方法其=￡% 高差闭合差容许值 义￡±12诚(所0 (n 为测站数，适合山地) 丁磅M±40JZ(*^)(L 为测段长度，以公里为单位，适合平地) 水准测量中，消除闭合差的原则一般按距离或测站数成正比地改正各段的观测高差 :

闭合导线平差公式

闭合导线坐标计算 闭合导线计算式根据外业观测的边长、夹角和方位角以及其中一个导线点的坐标，结合平差计算，来推算其余各导线点的坐标。 设对闭合导线n 个内角分别进行了观测，各个符号精度要求的观测值为β i 测，并对 闭合多边形的n 个边长分别进行了测量，各个符号精度要求的观测值为L i ；其中一个导 线点的坐标为 x i y i ；确定其余各个导线点的坐标x x i 1 +,y i 1 + 1 角度闭合差的计算也调整 （1）实测角度闭合差的计算 闭合导线n 个实测内角的和∑ 测 β不等于其理论值（n-2）*180,其差称为角度闭合差以 f β 表示： ?--=∑180*2）（测 n f β β (2)实测角度闭合差检核 角度闭合差校核是将实测角度闭合差也同级导线角度闭合差的容许值 f 容 β ，按各级导线测 量主要技术要求比较，以确定角度综合限差是否满足要求。这里角度综合限差采用图根导线数据，即 f 容β =40 ' 'n 。 （3）角度闭合差的调整 若 f β≤ f 容β ，则可以进行角度闭合的调整，否则，应分析情况重测。角度闭 合差的调整原则是，将 f β 以相反的符号平均分配到各个观测角中，即各点改正数为式 v β= f β／n 计算时，根据角度的取位的要求，改正数可凑整到1″、6″、10″.若不能均分，一般情况 下，因短边角引起的误差较大，因此给短边角的夹角多分配一点，使各角改正数的总和也反号的闭合差相等，即 f v β β-=∑ 2、推算各边的坐标方位角 推算各边的坐标方位角目的是为了计算坐标增量。推算方法根据起始方位角及改正后的

转折角，按式依次推算出各边的坐标方位角。 或 βα α右 -+=+1801 i i 1801 -+=+βα α左 i i 式中：αi ----------第i 条边的正方位角 α1 -i ---------第i+1条的正方位角 ββ右 左 --------分别为第i-1条边与第i 条边间所夹的左右角。 在推算过程中，如果算出 αi ＞360°，则应减去360°如果算出的α i ＜0°，则应加 上360° 为了发现推算过程中的差错，最后必须推算至起始边的坐标方位角，看其是否与已知值相等，以此作为计算校核。 3 坐标增量闭合差的计算和调整 （1）计算实测各边的坐标增量 设第i 条实测边的终、横坐标增量分别为 α i i i L X cos .=?测 αi i i L Y sin .=?测 (2)确定理论纵、横坐标增量∑△Xi 理、∑△Yi 理 闭合导线的纵横坐标增量总和的理论值应为零，则有 ∑△Xi 理=0 ∑△Yi 理=0 (3)计算坐标增量闭合差fx.fy 由于测量误差，改正后的角度仍有残余误差，坐标增量总和的测量计算值∑△X 测与∑△Y 测一般都不为零，其值称为坐标增量闭合差，fx.fy 表示，则 fx=∑△Xi 测-∑△Xi 理=∑△Xi 测 fy=∑△Yi 测-∑△Yi 理=∑△Yi 测 （4）计算导线全长闭合差f 并检核全长相对闭合差K 因计算的闭合导线并不闭合，而存在一个缺口，这个缺口的长度称为导线闭合差f fy fx f 2 2 + = 导线越长，全长闭合差也越大。通常用相对闭合差来衡量导线测量的精度，导线的全长相对

实验三-利用matlab程序设计语言完成某工程导线网平差计算

实验三利用matlab程序设计语言完成某工 程导线网平差计算 实验数据； 某工程项目按城市测量规范（CJJ8-99）不设一个二级导线网作为首级平面控制网，主要技术要求为：平均边长200cm，测角中误差±8,导线全长相对闭合差<1/10000,最弱点的点位中误差不得大于5cm，经过测量得到观测数据，设角度为等精度观测值、测角中误差为m=±8秒，鞭长光电测距、测距中误差为m=±0.8√smm，根据所学的‘误差理论与测量平差基础’提出一个最佳的平差方案，利用matlab完成该网的严密平差级精度评定计算; 平差程序设计思路: 1采用间接平差方法,12个点的坐标的平差值作为参数.利用matlab进行坐标反算，求出已知坐标方位角；根据已知图形各观测方向方位角； 2计算各待定点的近似坐标,然后反算出近似方位角，近似边.计算各边坐标方位角改正数系数； 3确定角和边的权,角度权Pj=1；边长权Ps=100/S； 4计算角度和边长的误差方程系数和常数项，列出误差方程系数矩阵B，算出Nbb=B’PB,W=B’Pl,参数改正数x=inv(Nbb)*W;角

度和边长改正数V=Bx-l； 6 建立法方程和解算x,计算坐标平差值, 精度计算；程序代码以及说明: s10=238.619;s20=170.759; s30=217.869;s40=318.173; s50=245.635;s60=215.514; s70=273.829;s80=241.560; s90=224.996;s100=261.826; s110=279.840;s120=346.443; s130=312.109;s140=197.637; %已知点间距离 Xa=5256.953;Ya=4520.068; Xb=5163.752;Yb=4281.277; Xc=3659.371;Yc=3621.210; Xd=4119.879;Yd=3891.607; Xe=4581.150;Ye=5345.292; Xf=4851.554;Yf=5316.953; %已知点坐标值 a0=atand((Yb-Ya)/(Xb-Xa))+180; d0=atand((Yd-Yc)/(Xd-Xc)); f0=atand((Yf-Ye)/(Xf-Xe))+360; %坐标反算方位角a1=a0+(163+45/60+4/3600)-180 a2=a1+(64+58/60+37/3600)-180; a3=a2+(250+18/60+11/3600)-180;

Matlab绘制频散曲线程序代码(20210119130722)

Matlab绘制频散曲线 程序代码 -CAL-FENGHAI-(2020YEAR-YI function disper %绘制平板频散曲线 %tic

clc;clear; cl=5790;%材料纵波波速(钢板) cs=3200;%材料横波波速(钢板) dfd=*le3; fdO=:dfd/le3:2O)*le3;%频厚积(MHz*mm) d_Q235二6; cps_mi n二2700; cpa_min=100; cp_max=10000; mode=3;%绘制的模式数 precision=le-8; cpa=zeros(length(fdO),mode); cps=zeros(le ng th(fdO),mode); for i=l:length(fdO) fd=fdO(i); [cpl2 n]=ss(cps_min/cp_max/fd/cl,cs,mode); for j=l:n cpl=cpl2(j,l); cp2=cpl2(j,2); cps(i,j)=serfe n(cpl,cp2,fctcl￡S'precisi on); end [cpl2 n]=aa(cpa_min,cp_max/fd/cl/cs,mode); for j=l:n cpl=cpl2(j,l); cp2=cpl2(j,2); cpa(ij)=aerfe n(cpbcp2,fd￡l‘cs,precisi on); end end h=zeros(mode,2); %相速度 figure(l) for j=l:2 ifj==l cp=cps; color=,b,; else cp=cpa; color二T; end for i=l:mode cpp=cp(:,i); in d=fi nd(cpp==0); if ^isempty(ind) h(i/j)=plot((fdO(ind(end)+l:end))/d_Q235/cpp(ind(end)+l:end),color); else h(i/j)=plot(fdO/d_Q235,cpp/color); end hold on end ifj==2 xlabel('f/(KHz)') ylabel('C_{p}/(km-sA{-l})')

matlab在电路中的全响应程序

A=[1 0 0 1 1 0;0 0 1 0 -1 1;0 1 0 -1 0 -1] %输入相关矩阵A Is=[2;0;-5;0;0;0] %输入电流源列向量 Us=[0;3;0;-2;0;0] %输入电压源列向量 Y=[1 0 0 0 0 0 ;-2 1/3 0 0 0 0 ;0 0 1/4 0 0 0;0 0 0 1/5 0 1/3 ;0 0 0 0 1/2 0;0 0 0 0 0 1/8] %输入导纳矩阵Y Z=inv(Y)%求解阻抗矩阵Z Yn=A*Y*A' %求解节点导纳矩阵Yn Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un U=A'*Un %求解之路电压U I=Y*(U+Us)-Is %求解支路电流I 结果如下： Z = 1.0000 0 0 0 0 0 6.0000 3.0000 0 0 0 0 0 0 4.0000 0 0 0 0 0 0 5.0000 0 -13.3333 0 0 0 0 2.0000 0 0 0 0 0 0 8.0000 Yn = 1.7000 -0.1667 -0.5333 -0.5000 0.8750 -0.1250 -2.2000 -0.4583 0.9917 Un = -7.8339 -13.7855 -25.1626 U = -7.8339 -25.1626 -13.7855 17.3287 5.9516 11.3772 I = -9.8339 8.2803 1.5536 6.8581 2.9758 1.4221

A=input('相关矩阵A=:') Is=input('电流源Is= :') Us=input('电压源Us=:' ) display ('是否含有受控源？是（a=1)否（a=0)') a=input('输入a=:' ) if(a==0) %判断语句 Z=input('阻抗矩阵Z=:') Y=inv(Z) %求解导纳矩阵Y else (a==1) Y=input('导纳矩阵Y=:') Z=inv(Y) %求解阻抗矩阵Z end Yn=A*Y*A' %求解节点导纳矩阵Yn Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un U=A'*Un %求解之路电压U I=Y*(U+Us)-Is %求解支路电流I 含受控源的运行结果如下 相关矩阵A=:[1 0 0 1 1 0;0 0 1 0 -1 1;0 1 0 -1 0 -1] 电流源Is= :[2;0;-5;0;0;0] 电压源Us=:[0;3;0;-2;0;0] 是否含有受控源？是（a=1)否（a=0) 输入a=:1 导纳矩阵Y=:[1 0 0 0 0 0 ;-2 1/3 0 0 0 0 ;0 0 1/4 0 0 0;0 0 0 1/5 0 1/3;0 0 0 0 1/2 0;0 0 0 0 0 1/8] Z = 1.0000 0 0 0 0 0 6.0000 3.0000 0 0 0 0 0 0 4.0000 0 0 0 0 0 0 5.0000 0 -13.3333 0 0 0 0 2.0000 0 0 0 0 0 0 8.0000 Yn = 1.7000 -0.1667 -0.5333 -0.5000 0.8750 -0.1250 -2.2000 -0.4583 0.9917 Un = -7.8339 -13.7855 -25.1626 U = -7.8339 -25.1626 -13.7855

闭合导线的内业计算

导线的内业计算 导线闭合平差计算步骤： （绘制草图，在图上填写已知数据和观测数据。） 一、角度闭合差的计算和角度的调整 （1)计算角度闭合差： f β =Σβ测－Σβ理=Σβ测－(n －2)×180°（n 为内角的个数） (2)计算限差： f β允许=±40"√n (n 为角的个数） （3）若在限差内，则平均分配原则，计算改正数: V=－f β/n 计算改正后新的角值： βi'=βi +v β 注：当角度合差不能整除时，可将余数再分配到含有短边的角上。（原因：由于仪器对中和目标偏心的原 因。含有短边的角可能产生较大的误差） 二、方位角的计算:（按新的角值，推算各边坐标方位角）导线各边的坐标方位角，是按各边的已知坐标方位角和导线的 转折角依次推算出来的。 注意：(①当计算结果出现负值时，则加上360? ②当计算结果出现大于360?时，则减360?) αbc =αab +180?-β右 αbc =αab -180?+β左 三、坐标增量的计算和坐标增量的调整。 按坐标正算公式，计算各边坐标增量 1.坐标增量的计算：△x ab = x b -x a =dcos αab △y ab =y b -y a =dsin αab 2.坐标增量闭合差的计算 f x=Σ△x 测 f y=Σ△y 测 3.由于坐标增量闭合差的存在，使闭合导线不能闭合 故：导线长相对闭合差k: =K T D f 1 = ∑ 5.分配坐标增量闭合差。 若k <1/2000（图根级），则将f x,f y 以相反符号，按边长成正比分配到各坐标增量上去，并计算改正后 的坐标增量。 v △x i =- ∑D f x D i △ x i ' =△x +v △x i v △y i =- ∑D f x D i △ y i '=△x +v △y i 6.坐标计算。根据起始点的已知坐标和经改正的新的坐标增量，来依次计算各导线点的坐标 图根导线测量的主要技术要求 测角中误差（"） 方位角闭合差 导线长度（m ） 相对闭合差 边长 一般 首级控制 一般 首级控制 ≤1.0m ≤1/2000 ≤1.5倍测图最大视距 30 20 60√n 40√n

全站仪测量闭合导线如何平差计算出各点坐标

如果你想学习导线(闭合、符合、支导线），我可以传份学习资料给你， 如果要严密平差建议用清华山维测量平差软件 如果简单平差可以先推算方位角闭合差，然后将闭合差平均分配到每站测的角度上进行角度平差。然后用平差后的角度推算坐标闭合差，得到的x和y的闭合差平均分配到每一站的坐标上即可 ，求得导线绝对闭合差，在除以导线全长得到导线全长相对闭合差 导线平差主要是看方向中误差和导线全长闭合差及导线全长相对闭合差 工程测量闭合导线差怎么计算的？？？ 闭合导线平差手算简单来讲分两步： 1.先计算出导线闭合环内角和，它与理论值（n-2）x180相减产生的闭合差平均分配到各个转角。使修正后的内角和等于理论值。 2.根据已知坐标方位角（已知两点坐标可求得坐标方位角）与修正后的各个转角值求出导线边坐标方位角。再通过方位角、导线边的长度计算出各个导线边产生的坐标增量。算到起算点后，X、Y 增量和的理论值均应为0。但因观测误差，坐标增量和往往不等于0。将偏差值平均分配到各个点位上，以消除偏差。 最后，用起算坐标依次加上修正后的坐标增量，就可以得到平差后各点的坐标值了。 从一个已知点出发，在连续测量多个点位后再到原出发点，这就叫导线的闭合，因为误差的存在，在闭合时有可能产生误差，这时你需要将误差以每根导线长度为权重进行平差。 看到百笑狂生的回答了，忍不住想说几句，这个什么“原位往复闭合”、“开路测量闭合”等等，你是从哪儿学来的概念？反正在测量这个专业里是没有这个概念的，楼主所说的导线测量，是建立平面控制网的一种最常见的测量方法，下面我简单介绍一下： 一、导线进行测量共有三种方法： 1、从一个已知点出发，依次对各个目标点进行测量，这种测量方式，因为累积误差的原因，在精度要求较高的场合一般不采用，也谈不上什么闭合不闭合的问题； 2、从一个已知点出发，依次对各个目标点进行测量后，再回到这个已知点，也就是说，将已知点做为测量的最后一个点也进行测量。如果没有误差出现，那么最后一个点的测量结果应与已知点相同，这就叫导线的闭合。如果有误差你再根据规范要求进行平差； 3、还有一种情况便是从一个已知点出发，依次对各个目标点进行测量后，最后回到另一个已知点，这种情况的闭合叫“导线闭合于某已知点”，由此可见，这个测量结果还包含了原已知点的误差，因此其精度不及前者，但有时限于测量对象分布条件的限制，可能也得采用。 上述测量方法，同时适用于水准测量，水准测量同样也采用导线法，但所采用的仪器不同（以前测量时是分开的，现在因为全站仪的出现已经同步进行了）。二、百笑狂生所谓的第一种测量方法，类似于导线中某个点的测量程序，对某个

附合导线平差程序设计报告

《测量平差程序》课程设计 （报告） 学生姓名：罗正材 学号：1108030128 专业：2011级测绘工程 指导教师：肖东升

目录 一、前言 (3) 二、平差程序的基本要求 (3) 三、平差程序模块化 (3)

图1 四、平差中的重要函数 （一）、角度制与弧度制的相互转化 C/C++程序设计中，关于角度的计算以弧度制为单位，而在测量以及具体工作中我们通常习惯以角度制为单位。这样，在数据处理中，经常需要在角度制与弧度制之间进行相互转化。这里，我们利用C/C++数学函数库math.h中的相关函数完成这两种功能。 这里，我们使用double类型数据表示角度制数和弧度制数。例如：123度44分58.445秒，用double类型表示为123.4458445，其中分、秒根据小数位确定。 在角度制与弧度制的转化中，涉及如下图2所示的两个环节。 度.分秒度弧度 图2 1.角度化弧度函数 double d_h(double angle) //角度化弧度 { double a,b; angle=modf(angle,&a);//a为提取的度值（int类型），angle为分秒值（小数） angle=modf(angle*100.0,&b); // b为提取的分值（int类型），angle为秒值（小数） return (a+b/60.0+angle/36.0)*(PI+3.0E-16)/180.0; } 2.弧度化角度函数 double h_d(double angle) //弧度化角度

{ double a,b,c; angle=modf(angle*180.0/(PI-3.0E-16),&a); angle=modf(angle*60.0,&b); angle=modf(angle*60.0,&c); return a+b*0.01+c*0.0001+angle*0.0001; } 其中，函数modf(angle,&a)为C语言数学库函数，返回值有两个，以引用类型定义的a 返回angle的整数部分，函数直接返回值为angle的小数部分。 （二）近似坐标计算 在平面网间接平差计算中，近似坐标计算是非常重要的一项基础工作。近似坐标是否计算成功是间接平差是否可以进行的必要条件。 1.两方向交会 已知条件：两个点的近似坐标，这两个点到未知点的方位角，如图3所示 图3两方向交会 根据图4.2，设 1 1 α tg k=， 2 2 α tg k=，则很容易写出 ? ? ? ? ? ? - = - - = B P B P A P A P y y k x x y y k 2 1 整理该式，得两方向交会的的计算公式 ?? ? ? ? ? - - = ?? ? ? ? ? ?? ? ? ? ? - - B B A A P P y x k y x k y x k k 2 1 2 1 1 1 （4.1）对（4.1）式计算，即可得到未知点的近似坐标。应用中需要注意的是，若两方向值相同或相反，则该式无解。 程序中，定义该问题的函数为：int xy0ang(obser &a1,obser &a2) 2.三边交会 如图4所示，为排除两边长交会的二义性，给出如下三边交会的模型，已知条件：三个

MATLAB水准网间接平差课程作业程序

条件平差 A=[1 -1 0 0 1 0 0;0 0 1 -1 1 0 0;0 0 1 0 0 1 1;0 1 0 -1 0 0 0]; Q=diag(s); W=zeros(4,1); W(1)=h(1)-h(2)+h(5); W(2)=h(3)-h(4)+h(5); W(3)=h(3)+h(6)+h(7); W(4)=h(2)-h(4)+H(1)-H(2); W=W.*1000; Naa=A*Q*(A'); K=-1.*inv(Naa)*W; V=Q*(A')*K; L=h'+V./1000; X=zeros(1,3); X(1)=H(1)+L(1); X(2)=H(1)+L(2); X(3)=H(2)-L(7); ZWC=sqrt(V'*inv(Q)*V/4); QLL=Q-Q*A'*inv(Naa)*A*Q; ZWC_h5=ZWC*sqrt(QLL(5,5)); 间接平差 h1=1.359; h2=2.009; h3=0.363; h4=1.012; h5=0.657; h6=0.238; h7=-0.595; H1=5.016 H2=6.016 h=[h1 h2 h3 h4 h5 h6 h7]' s=[1.1 1.7 2.3 2.7 2.4 1.4 2.6]' B=[1 0 0 ;0 1 0; 1 0 0;0 1 0 ; -1 1 0 ; -1 0 1 ;0 0 -1 ] p=diag(1./s) l=[0;0;4;3;7;2;0] W=B'*p*l Nbb=B'*p*B x=inv(Nbb)*W V=(B*x-l) H=h+V/1000 Q=inv(Nbb) n=7;