matlab闭合导线代码,【求助】闭合导线简易平差Matlab代码需要大神帮忙
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
%观测数据的输入
n=input('请输入闭合导线的点数:');
for i=1:n %按边数进行循环
a=num2str(i); %数字转换为字符串
string1=(char(('请输入第')',a,('条导线的边长:')'))';
bian(i)=input(string1);
b=num2str(i);
string2=(char(('请输入第')',b,('个转折角:')'))';
jiao1=input(string2,'s');
c=jiao1;
du=str2num(c(1:3));
fen=str2num(c(4:5));
miao=str2num(c(6:7));
jiao2(i)=du+fen/60+miao/3600;
end
%起算数据的输入
x=input('请输入起始已知点x的坐标:');
y=input('请输入起始已知点y的坐标:');
string3=(char('请输入起始已知方向的起算方位角:'));
fwj1=input(string3,'s');
c=fwj1;
du=str2num(c(1:3));
fen=str2num(c(4:5));
miao=str2num(c(6:7));
fwj2=du+fen/60+miao/3600;
%角度闭合差的计算
fb=sum(jiao2)-(n-2)*180;
for i=1:n
v(i)=(-1)*fb/n;
end
%真方位角计算
for i=1:n
if i==1
fwjz(1)=fwj2+jiao2(1)+v(1)-180;
end
if i>1
fwjz(i)=fwjz(i-1)+jiao2(i)+v(i)-180;
end
if fwjz(i)>0
fwjz1(i)=fwjz(i);
end
if fwjz(i)<0
fwjz1(i)=fwjz(i)+360;
end
end
%坐标闭合差计算
xzq=cosd(fwj2)*bian(1);
yzq=sind(fwj2)*bian(1);
for i=1:(n-1)
xz(i)=cosd(fwjz1(i))*bian(i+1);
yz(i)=sind(fwjz1(i))*bian(i+1);
end
fx=xzq+sum(xz);
fy=yzq+sum(yz);
fs=sqrt(fx^2+fy^2);
k=fs/sum(bian);
%坐标计算
for i=1:n
vx=(-1)*fx/sum(bian)*bian;
vy=(-1)*fy/sum(bian)*bian;
end
for i=1:n
if i==1
X(1)=x+xzq+vx(1);
Y(1)=y+yzq+vy(1);
end
if i>1
X(i)=X(i-1)+xz(i-1)+vx(i);
Y(i)=Y(i-1)+yz(i-1)+vy(i);
end
end
matlab闭合导线代码,【求助】闭合导线简易平差Matlab代码需要大神帮忙相关推荐
- 求大神帮忙,想要一元20次方程的解matlab代码,在线等,(解一元20次方程的解matlab代码)
求大神帮忙,想要一元20次方程的解matlab代码,在线等,(解一元20次方程的解matlab代码) 已知y是一个1*101的数组如下,函数如下,想求得x 135 130 134 128 129 13 ...
- ubuntu14.04扩展屏幕后,打开matlab就死机。求大神帮忙解决!!!万分感谢
ubuntu14.04扩展屏幕后,打开matlab就死机.求大神帮忙解决!!!万分感谢
- java主界面设置背景图片_java 窗体设置背景图片问题?(附上登陆界面代码,我想加个背景图片,求大神帮忙改改)...
java 窗体设置背景图片问题?(附上登陆界面代码,我想加个背景图片,求大神帮忙改改) 关注:223 答案:4 mip版 解决时间 2021-01-26 22:09 提问者非莪莫属 2021-01 ...
- 如何修改matlab中的语句,求大神帮忙看一下这个语句怎么改!!!!!
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 clc clear W=imread('e:\matlab\aaa\1.jpg'); imshow (W); I = rgb2gray(W); % 灰度处 ...
- 清华姚班毕业生开发新特效编程语言,99 行代码实现《冰雪奇缘》,网友:大神碉堡!创世的快乐...
公众号关注 "GitHubDaily" 设为 "星标",每天带你逛 GitHub! 转自量子位,作者边策.鱼羊 只用 99 行代码,你也可以像<冰雪奇缘& ...
- matlab程序闭合导线网程序,【求助】闭合导线简易平差Matlab代码需要大神帮忙
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 %观测数据的输入 n=input('请输入闭合导线的点数:'); for i=1:n %按边数进行循环 a=num2str(i); %数字转换为字符串 s ...
- matlab根号图像,求助高手,用matlab求两幅图像平方和再开根号公式怎样表达?
不用那么麻烦,matlab有现成的函数:均值:mean2(img)标准差:std2(img)方差是标准差的平方,把上面的结果平方一下www.mh456.com防采集. ^p1=imread('data ...
- matlab 降低维度,求助。。。matlab索引超出维度要怎么修改。。。谢谢
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 程序如下 %%六种板形缺陷 y=-1:0.2:1;%残余应力数目(横向) p1=y;%左边浪 p2=-p1;%右边浪 p3=3*(y.^2)/2-0.5; ...
- matlab超出维度,求助。。。matlab索引超出维度要怎么修改。。。谢谢
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 程序如下 %%六种板形缺陷 y=-1:0.2:1;%残余应力数目(横向) p1=y;%左边浪 p2=-p1;%右边浪 p3=3*(y.^2)/2-0.5; ...
- 断点回归matlab,断点回归求助 - Stata专版 - 经管之家(原人大经济论坛)
你的看rd的帮助啊 你需要设置cutoff点在哪里 比如随便设置一个cutoff在2007 . rd si d, gr mbw(100) z0(2007) Two variables specifie ...
最新文章
- BERT的通俗理解 预训练模型 微调
- CSS捡屎记 // Web开发之精通CSS
- 按钮 小程序 弹出菜单_公众号怎么关联小程序
- R语言观察日志(part20)--包的组件之R代码
- asp.net core结合NLog搭建ELK实时日志分析平台
- Windows Phone开发(8):关于导航的小技巧
- 为何我中断执行的线程不起作用,Why
- KG—Linux添加新的环境变量以及对PATH环境变量的修改
- SCO UNIX环境下自动增加网关的两种方法
- Android接入融云即时通讯SDK,实现单聊和发送自定义消息。
- 移动硬盘格式化了?这样恢复数据
- 利用OpenGL模拟太阳系
- python的学习笔记案例1--汇率兑换3.0
- VMware跨电脑移动Linux虚拟机
- 16岁男孩发布病毒“小兔子” 与熊猫烧香类似
- stm32f103r6最小系统原理图_PCB stm32的最小系统原理图 基于STM32F103VBT6的最小系统板 - 下载 - 搜珍网...
- 威世 | 新款 LLC 谐振变压器/电感器,优化 PCB 空间布局
- HTTPS能否避免流量劫持?
- n1 盒子安装armbian 5.6
- Unity的Handles类