最新matlab导线网程序
m a t l a b导线网程序
课程编号：课程性质：
MATLAB及其应用
课程论文
学院：测绘学院
专业：测绘工程
姓名：
学号：
目录
一、题目内容 (3)
二、程序编写思路 (4)
1、数学模型 (4)
2、程序解析 (4)
三、程序运行结果 (5)
四、程序源代码 (6)
注：这是《误差理论与测量平差基础》这一课程在学习其基本原理之后为解决实际问题而编写的一个程序。
一、题目内容
要求：对于给出的导线控制网图2 ，根据已知条件和观测数据，设计平差方案，编制平差程序，给出平差结果，评定精度。
题目：A 、B 为已知点，已知点坐标已给出如下。 为待定点。同精度观测了24个角度 ，测角中误差 。测量了17条边长 ，其观测结果及中误差列于表4中，试按间接平差法求： （1）平差后单位权中误差；
（2） 14个待定点的的坐标平差值及中误差；
（3）平差后 边的相对中误差。
已知点坐标：X A =730.024 ，Y A =126.040 ，111.855=X B ，
232.172=Y
B
表一：导线网角度观测值
表二：导线网边长观测值及精度 14
1~P P 241~L L "
10β
σ=171~S S 17
S
2 195.782 14 8 136.578 11.7 14 119.1 10.9
3 184.98
4 13.6 9 112.502 10.6 1
5 123.94 11.1 4 110.03
6 10.5 10 131.074 11.5 16 122.84 11.1 5 163.14
7 12.
8 11 140.88 11.
9 17 139.46 11.8 6 219.482 14.9 12 132.608 11.5
图一：
二、程序编写思路
1、数学模型
◆函数模型 观测方程 L=BX+d 误差方程 V B l =x
◆随机模型
2
X X X X
D =Q 0σ
◆方程的解与精度评定
令 PB B N T BB = pl W B T
= 法方程即为 0?=-W x
N BB 可得解为：W BB x
N 1
?-= 故可得平差结果：x
X V L L
X
??,?0
+=
+= 从而解得单位权中误差为
T
2
0V PV /R
σ 2
X X
X X D
=Q 0σ
2、程序解析
◆计算框图
三、程序运行结果
◆平差后单位权中误差： = 5.30 (由于程序取的先验单位权中误差为10，两
者相比较，故可以认为这个结果是可以接受的)
◆14个待定点的坐标平差值及中误差见下表： 表三：
表七：
◆平差后S17边的相对中误差为：253451
四、程序源代码
这是要输入的文件
软件无线电实验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;
基于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 % 保存旧的聚类中心
基于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
软件无线电实验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;
实验三-利用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
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;
测绘 matlab 编程 程序
导线网严密平差 完整程序代码及注释 已知数据 已知点坐标： X A =730.024 ， Y A =126.040 ， 111.855=X B ，232.172=Y B
应得结果 = 5.30 表七：
实际结果 点号X(m) Y(m) 中误差X(mm) 中误差Y(mm) 点位中误差 1 678.1641 287.3411 3.3868 5.8737 7.6078 2 564.6919 446.8827 6.0102 7.5137 9.1937 3 533.8468 629.2749 8.7985 7.3927 10.0596 4 475.6776 535.8728 7.3949 8.4987 9.9667 5 499.1408 374.4231 5.807 6 7.6840 9.1827 6 594.6190 176.7976 5.2308 3.4068 7.3475 7 826.3848 305.7569 2.7439 5.4473 7.1551 8 767.8778 401.8510 4.8371 6.3810 8.3734 9 745.7230 531.0410 6.8992 6.6276 9.1947 10 671.6256 650.8637 9.0695 6.0533 9.7220 11 794.2833 600.4688 8.4994 6.1800 9.5784 12 859.0528 511.5349 6.7399 6.6874 9.1608 13 898.2695 399.0763 4.7580 6.4207 8.3586 14 919.3459 276.9409 3.1903 4.8006 7.0670 单位权中误差： 5.3393
常用的MATLAB程序和函数
==等于 <小于 >大于 <=小于或等于>=大于或等于~=不等于 4 常用内部数学函数 指数函数exp(x)以e为底数 对数函数log(x)自然对数，即以e为底数的对数log10(x)常用对数，即以10为底数的对数log2(x)以2为底数的x的对数 开方函数sqrt(x)表示x的算术平方根
绝对值函数abs(x)表示实数的绝对值以及复数的模 三角函数（自变量的单位为弧度）sin(x)正弦函数cos(x)余弦函数tan(x)正切函数cot(x)余切函数sec(x)正割函数csc(x)余割函数 反三角函数asin(x)反正弦函数acos(x)反余弦函数atan(x)反正切函数acot(x)反余切函数asec(x)反正割函数acsc(x)反余割函数 双曲函数sinh(x)双曲正弦函数cosh(x)双曲余弦函数tanh(x)双曲正切函数coth(x)双曲余切函数sech(x)双曲正割函数csch(x)双曲余割函数 反双曲函数asinh(x)反双曲正弦函数acosh(x)反双曲余弦函数atanh(x)反双曲正切函数acoth(x)反双曲余切函数asech(x)反双曲正割函数acsch(x)反双曲余割函数 求角度函数atan2(y,x)以坐标原点为顶点，x轴正半轴为始边，从原点到点（x，y）的射线为终边的角，其单位为弧度，范围为（ ， ] 数论函数gcd(a,b)两个整数的最大公约数lcm(a,b)两个整数的最小公倍数
排列组合函数factorial(n)阶乘函数，表示n的阶乘 复数函数real(z)实部函数 imag(z)虚部函数 abs(z)求复数z的模 angle(z) 求复数z的辐角，其范围是（ ， ] conj(z)求复数z的共轭复数 求整函数与截尾函数ceil(x)表示大于或等于实数x的最小整数floor(x)表示小于或等于实数x的最大整数round(x)最接近x的整数 最大、最小函数max([a，b，c，．．．]) 求最大数min([a，b， c，．．]) 求最小数 符号函数 sign(x)
基于matlab的测量导线描绘
测绘学院学号:2009301610208 夏少波手机: 基于MATLAB的测量导线描绘 摘要: 本文从导线的描绘出发，提出了该问题的MATLAB解决方法。给出了数据的处理结果和图表。并且，设计了二维和三维的代码描述。其中对待处理的数据有一些优化，但计算方法和计算的逻辑是没有问题的。最后，MATLAB对这个问题的解决也较为令人满意。 关键字: 控制测量导线代码运行结果MATLAB 任何测量工作都会产生误差，而误差在绝对意义上是不可避免的，所以在测量工作中必须采用一定的方法来控制误差的范围，以期达到工作需求。测量工作中遵循“先整体后局部，先控制后碎部”的原则。简而言之就是在工作区域内建立一定精度的由一定数目的已知坐标点构成的几何图形，然后在几何图形内进行测量。在测量过程中，我们可以通过各个已知点的坐标来验算测得的未知点的坐标。这样通过不断地检核、测量、检核，成果的误差就可以控制在一定范围内。其中，已知点被称作“控制点”，几何图形被称为“控制网”，这种测量方法被称为“控制测量”。 一般意义上的控制测量可以分为“平面控制测量”和“高程控制测量”。平面控制测量是基平面坐标系X-Y，而高程控制测量是基于高程H的。目前，随着测绘技术的发展，也出现了把两种控制测量结合起来的“三维控制测量”。 控制网根据其精度从高到低大致可以分为一、二、三、四等，精度越高的点数目越少。点数不够自然无法形成有效地控制网。因此，测量工作中需要通过已知的几个点，“一生二，二生三”，绘出一个可用于实际工程的控制网。“一生二，二生三”的过程需要已知点前后通视（即互相看得见），然而，在城市中，由于建筑物较多，视线不好，往往会出现已知点之间无法通视的现象。所以就引入了导线测量（traverse survey）。 导线测量又有很多分类，在此就不在赘述了。下面来介绍导线测量中的附和导线。所谓附和导线就是：“导线起始于一个已知控制点而终至于另一个已知控制点”。通过测量导线上未知点与已知控制点坐标的对比、计算，进行平差（控制误差的一种手段），得出一条导线上的各个待测点的坐标。这个过程就是导线测量的主体步骤了。导线示意图如下
简单的matlab程序
实验报告 （2010学年第一学期） 课程名称Matlab 教学院（部）工学院 任课教师及职称许芹 专业、班级电子信息工程 学号1104080129 姓名赵超
实验一Matlab基础入门四：实例练习 1，设两个复数a=2+3i,b=4+5i,计算a+b,a-b,a*b,a/b function one_1(a,b) a1=a+b a2=a-b a3=a*b a4=a/b >> a=complex(2,3); >> b=complex(4,5); >> one_1(a,b) a1 = 6.0000 + 8.0000i a2 = -2.0000 - 2.0000i a3 = -7.0000 +22.0000i a4 = 0.5610 + 0.0488i 2,计算下列函数的结果，其中x=-3.5?,y=8?, F= function one_2(x,y) x=x*pi/180 y=y*pi/180 f=(sin(abs(x)+abs(y)))/(sqrt(cos(abs(x+y)))) >> one_2(-3.5,8) x = -0.0611 y = 0.1396 f = 0.1997
3.我国人口按2000年第五次人口普查的结果为12.9533亿，如果人口增长率为1.07%，求公元2010年的人口数。 function one_3(n) a=12.9533 y=a*(1+0.00107)^(n-2000) >>one_3(2010) a = 12.9533 y = 13.0926 4，求解ax2+bx+c=0方程的根，其中a=1,b=2,c=3 function one_4(a,b,c) p1=[a,b,c] roots(p1) >> one_4(1,2,3) p1 = 1 2 3 ans = -1.0000 + 1.4142i -1.0000 - 1.4142i 5,已知三角形三边a=8.5,b=14.6.c=18.4,求三角形面积 function one_5(a,b,c) s=(a+b+c)/2; area=sqrt(s*(s-a)*(s-b)*(s-c)) >> one_5(8.5,14.6,18.4) area = 60.6106 实验二Matlab 符号运算 2，微积分 1）求极限：,,
电机设计matlab程序
%电机设计程序 clear all format short e m1=3;p=2;f=50 %1.额定功率 PN=5.5*10^3 ; %2.额定电压（单位V，三角形接法） UN=380;UN0=380; %3.功电流(单位A) IKW=PN/(m1*UN0) %4.效率eta按照技术条件的规定eta=0.875 eta=0.855 ; %5.功率因数cos(phi) =0.84,按照技术条件的规定cos(phi)=0.84 phi=acos(0.84); cos(phi); %6.极对数p=2 p=2; %7.定转子槽数：每极每相槽数取整数。参考类似规格电机取q1=3,则Z1=2m1pq1，再查表10-8选Z2=32，并采用转子斜槽。 q1=3; Z1=2*m1*p*q1 Z2=32 ; %8.定转子每极槽数 Zp1=Z1/(2*p) Zp2=Z2/(2*p) %9.确定电机的主要尺寸；一般可参考类似电机的主要尺寸来确定Di1和lef.现按10-2中的 KE1=0.0108*log(PN/1000)-0.013*p+0.931 P1=KE1*PN/(eta*cos(phi)) alphap1=0.68;KNm1=1.10;Kdp1=0.96;A1=25000; Bdelta1=0.69;n1=1450; V=(6.1/(alphap1*KNm1*Kdp1))*(1/(A1*Bdelta1 ))*(P1/n1) D1=0.21; %铁心的有效长度 Di1=0.136; lef =V/((Di1)^2) %气隙的确定 %参考类似产品或由经验公式（10-10a），得 lt=0.115; delta =0.0004 lef=lt + 2*delta D2=Di1-2*delta %转子内径先按转轴直径决定（以后再校验转子轭部磁密） Di2=0.048 ;
最新列车运行图MATLAB程序
%第一个 %设6:42时间为0时刻，转换如下。 a1=[6 7 7 7 7 7 7 8 8 9 9 9 9 9 10 10 10 10 10 10 11 11 12 12 12 12 13 13 14 14 15 15 15 16 16 16 16 17 17 17 17 18 18 18 18 19 20 21]; b1=[42 0 10 22 33 40 52 15 30 12 17 22 27 40 5 10 15 28 45 50 33 43 7 13 34 43 5 18 25 43 6 40 45 0 18 35 57 3 19 32 37 2 12 36 42 22 47 6]; c1=7*ones(1,length(a1)); d1=60*(a1-c1)+b1+18; e1=[7 7 7 8 8 8 8 8 9 9 9 10 10 10 10 10 11 11 11 11 12 12 12 12 13 13 13 14 15 15 15 16 16 16 16 17 17 18 18 18 18 18 19 19 19 20 21 21]; f1=[23 41 58 3 15 28 37 56 11 53 58 10 15 21 46 51 3 9 26 31 14 24 48 54 15 24 46 15 6 24 47 21 26 48 59 23 38 0 7 13 18 43 0 17 30 17 35 47]; g=7*ones(1,length(a1)); h1=60*(e1-g)+f1+18; %输入经历时间，求众数，得到一个差值矩阵f。 x1=[41 41 48 41 42 48 45 41 41 41 41 48 48 41 41 41 48 41 41 41 41 41 41 41 41 41 41 57 41
波束形成Matlab程序
1.均匀线阵方向图 %8阵元均匀线阵方向图，来波方向为0度 clc; clear all; close all; imag=sqrt(-1); element_num=8;%阵元数为8 d_lamda=1/2;%阵元间距d与波长lamda的关系 theta=linspace(-pi/2,pi/2,200); theta0=45/180*pi;%来波方向(我觉得应该是天线阵的指向) %theta0=0;%来波方向 w=exp(imag*2*pi*d_lamda*sin(theta0)*[0:element_num-1]'); for j=1:length(theta) %(我认为是入射角度，即来波方向,计算阵列流形矩阵A) a=exp(imag*2*pi*d_lamda*sin(theta(j))*[0:element_num-1]'); p(j)=w'*a; %(matlab中的'默认为共轭转置，如果要计算转置为w.'*a) end figure; plot(theta,abs(p)),grid on xlabel('theta/radian') ylabel('amplitude') title('8阵元均匀线阵方向图') 见张小飞的书《阵列信号处理的理论和应用2.3.4节阵列的方向图》
当来波方向为45度时，仿真图如下： 8阵元均匀线阵方向图如下，来波方向为0度，20log（dB）
随着阵元数的增加，波束宽度变窄，分辨力提高：仿真图如下：
2.波束宽度与波达方向及阵元数的关系 clc clear all close all ima=sqrt(-1); element_num1=16; %阵元数 element_num2=128; element_num3=1024; lamda=0.03; %波长为0.03米 d=1/2*lamda; %阵元间距与波长的关系 theta=0:0.5:90; for j=1:length(theta); fai(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_num1*d)); psi(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_num2*d)); beta(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_num3*d)); end figure; plot(theta,fai,'r',theta,psi,'b',theta,beta,'g'),grid on xlabel('theta'); ylabel('Width in radians') title('波束宽度与波达方向及阵元数的关系') 仿真图如下：
matlab导线平差代码
function [zjbl1 kp ukp sp sa ss pn pn1 rlf12 x y x1 y1 k1 k2 k3 ma ms1 ms2 angles zjbl2 k11 k22 s]=readdata; global zjbl1 kp ukp sp sa ss pn pn1 rlf12 x y x1 y1 k1 k2 k3 ma ms1 ms2 angles zjbl2 k11 k22 s; [filename,pathname]=uigetfile('*.dat;*.txt','?áè?êy?Y'); if filename~=0 [fid1 msg]=fopen(strcat(pathname,filename),'rt'); if fid1>0 kp=fscanf(fid1,'%f',1); ukp=fscanf(fid1,'%f',1); sp=kp+ukp; sa=fscanf(fid1,'%f',1); ss=fscanf(fid1,'%f',1); pn=fscanf(fid1,'%s',[1 sp]); pn1=pn; rlf12=fscanf(fid1,'%f',2); rlf12=de_to_rad(rlf12); x=fscanf(fid1,'%f',1); y=fscanf(fid1,'%f',1); x1=zeros(1,2); x1(2)=x(1); y1=zeros(1,2); y1(2)=y(1); ma=fscanf(fid1,'%f',1); ms1=fscanf(fid1,'%f',1); ms2=fscanf(fid1,'%f',1); zjbl1=fscanf(fid1,'%s%s%s%f',[4 sa]); %?D??±?á? zjbl1=zjbl1'; k1=zjbl1(:,1); k2=zjbl1(:,2); k3=zjbl1(:,3); angles=zjbl1(:,4); angles=de_to_rad(angles); zjbl2=fscanf(fid1,'%s%s%f',[3 ss]); %?D??±?á? zjbl2=zjbl2'; k11=zjbl2(:,1); k22=zjbl2(:,2); s=zjbl2(:,3); fclose(fid1);
matlab 常见经典平差 程序
希望本人编写的经典平差可以对matlab初学者以及测绘专业的学生有用 By cowboy function void() sprintf('请选择平差类型') sprintf('1: 参数平差\n2: 条件平差\n3：具有条件的参数平差\n4: 具有参数的条件平差\n') a=input('请输入对应号码\n'); switch (a) case 1 % 参数平差V=AX-L sprintf('1:参数平差V=AX-L') n=input('请输入n='); t=input('请输入t='); A=input('请输入系数矩阵A='); L=input('请输入常数项矩阵L='); P=diag(input('请输入权阵P=')) %P=input('请输入权阵P='); N=A'*P*A; U=A'*P*L; sprintf('开始计算') X=inv(N)*(U); V=A*X-L; Qx=inv(N); u=sqrt((V'*P*V)/(n-t)); fid=fopen('data.txt','wt'); %写入文件路径，文件输出 fprintf(fid,'=====================参数平差: V=AX-L =====================\n'); fprintf(fid,'输出n,t\n') ; fprintf(fid,'%12.5f %12.5f\n',n ,t) ;%n,t fprintf(fid,'输出矩阵A\n') [m,n]=size(A); %A for i=1:1:m for j=1:1:n if j==n fprintf(fid,'%12.5f\n',A(i,j)); else fprintf(fid,'%12.5f\t',A(i,j)); end end end fprintf(fid,'\n');
常用离散信号MATLAB程序及波形
常用离散信号MATLAB 表示 单位抽样序列 n=-10:10; x=[n==0]; stem(n,x); title('单位抽样信号'); xlabel('n');ylabel('x');grid on ; -10 -8 -6 -4 -2 02 4 6 8 10 00.10.20.30.40.50.60.70.80.91单位抽样信号 n x 延迟k=4个单位后的图像为： -10 -8 -6 -4 -2 02 4 6 8 10 00.10.20.30.40.50.60.70.80.91单位抽样信号 n x
单位阶跃信号 n=-5:10; x=[zeros(1,5),ones(1,11)]; stem(n,x,'m','p'); axis([-5,10,-0.5,1.5]); title('单位阶跃信号'); xlabel('n');ylabel('幅度') grid on ; -5 510 -0.50 0.5 1 1.5 单位阶跃信号 n 幅度 矩形序列 n=-5:15; x=[zeros(1,5),ones(1,11),zeros(1,5)]; stem(n,x,'r','h') axis([-5,15,-0.5,1.5]); title('矩形序列'); xlabel('时间');ylabel('幅度') grid on ;
-5 51015 -0.50 0.5 1 1.5 矩形序列 时间 幅度 正弦序列 n=-5:0.1:5; xn=2*sin(0.5*pi*n+pi/3); stem(n,xn) axis([-5,5,-3,3]); title('正弦序列'); xlabel('时间');ylabel('幅度'); grid on ; -5 -4-3-2-1 012345 -3-2 -1 1 2 3 正弦序列 时间 幅度
推荐-Matlab的导线网坐标计算 精品
东华理工大学长江学院 题目：基于Matlab的导线网坐标计算 English Title：Traverse Network Coordinate Calculation Based On Matlab 学生姓名：闫晓天申请学位门类：工学学士 学号：09351122 专业：测绘工程 系别：测绘工程系 指导老师：刘波职称：讲师 二零一三年六月
摘要 导线计算是在所有测量工作中经常遇见的问题之一，同时导线计算的方法也有很多种，本文主要是利用简单易懂的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 puter 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; puter Programming