空间后交-前交程序设计实验报告
上机实验报告(精选11篇)
上机实验报告篇1用户名se××××学号姓名学院①实验名称：②实验目的：③算法描述（可用文字描述，也可用流程图）：④源代码：（.c的文件）⑤用户屏幕（即程序运行时出现在机器上的画面）：2．对c文件的要求：程序应具有以下特点：a可读性：有注释。
b交互性：有输入提示。
c结构化程序设计风格：分层缩进、隔行书写。
3.上交时间：12月26日下午1点-6点，工程设计中心三楼教学组。
请注意：过时不候哟！四、实验报告内容0．顺序表的插入。
1.顺序表的删除。
2．带头结点的单链表的\'插入。
3.带头结点的单链表的删除。
注意：1.每个人只需在实验报告中完成上述4个项目中的一个，具体安排为：将自己的序号对4求余，得到的数即为应完成的项目的序号。
例如：序号为85的同学，85%4=1，即在实验报告中应完成顺序表的删除。
2.实验报告中的源代码应是通过编译链接即可运行的。
3.提交到个人空间中的内容应是上机实验中的全部内容。
上机实验报告篇2一、《软件技术基础》上机实验内容1．顺序表的建立、插入、删除。
2．带头结点的单链表的建立（用尾插法）、插入、删除。
二、提交到个人10m硬盘空间的内容及截止时间1．分别建立二个文件夹，取名为顺序表和单链表。
2．在这二个文件夹中，分别存放上述二个实验的相关文件。
每个文件夹中应有三个文件(.c文件、.obj文件和.exe文件)。
3. 截止时间：12月28日（18周周日）晚上关机时为止，届时服务器将关闭。
三、实验报告要求及上交时间（用a4纸打印）1．格式：《计算机软件技术基础》上机实验报告用户名se××××学号姓名学院①实验名称：②实验目的：③算法描述（可用文字描述，也可用流程图）：④源代码：（.c的文件）⑤用户屏幕（即程序运行时出现在机器上的画面）：2．对c文件的要求：程序应具有以下特点：a 可读性：有注释。
b 交互性：有输入提示。
顺序结构程序设计-C语言程序实验报告
安徽工程大学高级语言程序设计实验报告班级姓名同组者/ 成绩日期2019.09.30 指导教师实验名称顺序结构程序设计一、实验目的1.掌握数据的输入/输出方法，能正确使用有关格式转换符。
2.掌握顺序结构程序中语句的执行过程。
3.掌握顺序结构程序的设计方法。
二、实验内容1.P47页第一个程序的作用是依次输入2个整数，计算并输出这2个整数之差。
（1）分析程序，若运行时输出：200,160<回车>，预期结果是多少？（2）上机运行该程序，查看程序运行结果是否符合题目要求。
如果不符合，请分析原因并修改程序，直至符合要求为止。
2.P47页第二个程序用于实现按下列公式计算并输出s1和s2的值：s1=3/(a+b)2,s2=ab/(a+b)3,其中a,b为整型数据。
（1）根据题意修改上述程序，并进行调试，直到正确为止。
（2）在（1）的基础上，将“scanf("%d,%d",&a,&b);”改为“scanf("%d%d",&a,&b);”后再编译、连接、运行。
3.分析P47页第三个程序，写出预期结果，然后输入调试，查看运行结果与预期结果是否一致，并分析其原因。
4.编程实现下列功能并上机调试运行。
（1）设圆半径为r，求圆周长和面积。
要求用scanf函数输入数据，输出时取小数点后两位。
（2）输入一个3位十进制整数，分别输出百位、十位以及个位上的数。
（3）从键盘输入一个带两位小数的实数，将其整数部分和小数部分分离后输出。
（4）用getchar函数读入两个字符，然后分别用putchar和printf函数输出这两个字符。
*思考题：5. 若实验内容1的程序改为P48页第一个程序段，运行该程序，输入5,3，查看程序运行结果是否与自己预测的结果一致，并分析原因。
6.若实验内容1的程序改为P48页第二个程序段，运行该程序，输入5,3，查看程序运行结果是否与自己预测的结果一致，并分析原因。
空间后方交会程序设计
单片空间后方交会程序设计实验报告专业：测绘工程班级：081姓名：张发伟学号：09[一]、实习任务：用C或VC++语言实现单片后方交汇的计算。
[二]、实习目的：1.深入理解单片空间后方交会的原理；2.在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程；3．通过上机调试程序加强动手能力的培养。
[三]、实习原理:以单幅影象为基础，从该影象所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，求解该影象在航空摄影时刻的像片外方位元素Xs，Ys，Zs,ф,ω,κ.共线条件方程如下：x-x0=-f*[a1(X-Xs)+b1(Y-Ys)+c1(Z-Zs)]/[a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)] y-y0=-f*[a2(X-Xs)+b2(Y-Ys)+c2(Z-Zs)]/[a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)]x,y为像点的像平面坐标；x0,y0,f为影像的外方位元素；Xs,Ys,Zs为摄站点的物方空间坐标；X,Y,Z为物方点的物方空间坐标；[四]、程序框图：输入原始数据归算像点坐标x,,y计算和确定初值X s 0, Y s 0, Z s 0, φ0,ω0,κ0组成旋转矩阵R计算（x ），（y ）和l x ，l y逐点组成误差方程并法化所有点完否？ 解法方程，求未知数改正数 计算改正后外方位元素未知数改正数<限差否？整理并输出计算结果正常输出迭代次数小于限差否？ 输出中间结果和出错信息非正常结束[五]、实验数据：4个地面控制点的地面坐标及其对应像点的像片坐标：[六]、结果输出：已知条件:像点坐标x,y:-53.4 82.21-14.78 -76.6310.46 64.43153.24 -86.15地点坐标Xa,Ya,Za:37631.1 31324.5 728.6939101 24935 2386.540426.5 30319.8 757.31-68.99 36589.4 25273.3摄影比例尺：1:50000内方位元素:x0=y0=0 f=153.24mm计算结果:旋转矩阵：0.997709 0.0675334 0.00398399-0.0675254 0.997715 -0.00211178-0.0041175 0.00183792 0.99999像点坐标位:(单位：mm)-86.15 -68.99-53.41 82.21-14.78 -76.6310.47 64.43外方位元素：Xs=39795.435 精度为:1.1254813Ys=27476.479 精度为:1.2437625Zs=7572.6929 精度为:0.48380521q=-0.0039840098 精度为:0.00018182003w=0.0021117837 精度为:0.00015959235k=-0.067576934 精度为:7.2440432e-005迭代次数：4Press any key to continue[七]、实习总结：掌握运用程序设计实现单片后方交汇的计算，了解在多余观测情况下，用最小二乘平差方法计算相片外方位元素，及其相关精度。
实验报告 实验报告总结万能(优秀7篇)
实验报告实验报告总结万能（优秀7篇）（经典版）编制人：__________________审核人：__________________审批人：__________________编制单位：__________________编制时间：____年____月____日序言下载提示：该文档是本店铺精心编制而成的，希望大家下载后，能够帮助大家解决实际问题。
文档下载后可定制修改，请根据实际需要进行调整和使用，谢谢!并且，本店铺为大家提供各种类型的经典范文，如工作总结、计划大全、策划方案、报告大全、心得体会、演讲致辞、条据文书、作文大全、教案资料、其他范文等等，想了解不同范文格式和写法，敬请关注!Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you!Moreover, this store provides various types of classic sample essays for everyone, such as work summaries, plan summaries, planning plans, report summaries, insights, speeches, written documents, essay summaries, lesson plan materials, and other sample essays. If you want to learn about different formats and writing methods of sample essays, please stay tuned!实验报告实验报告总结万能（优秀7篇）总结就是对一个时期的学习、工作或其完成情况进行一次全面系统的回顾和分析的书面材料，它能使我们及时找出错误并改正，因此好好准备一份总结吧。
程序设计基础实验报告
实验一 用C 语言编写简单程序一、实验目的1. 熟悉VC6.0的编程环境，掌握运行C 程序的基本步骤。
2. 了解C 程序的基本框架，模仿例题编写简单的C 语言程序。
3. 正确书写算术表达式、赋值表达式和关系表达式。
4. 掌握基本输入输出函数的使用，正确调用C 语言提供的数学库函数。
5. 掌握简单的单步调试方法。
二、实验环境Windows XP ；Visual C++ 6.0。
三、实验内容1．在屏幕上显示一个句子“What is a computer?”。
2．求华氏温度150F 对应的摄氏温度。
计算公式如下：其中，c 表示摄氏温度；f 表示华氏温度。
3．输入x ，计算下列分段函数的值（保留2位小数），请调用sqrt 函数求平方根，调用pow 函数求幂。
4. 输入一个正整数m (0100)m ≤≤，求100i m i =∑。
四、实验要求1. 将上机验证正确的源代码写到实验报告上。
2．根据自己的真实感受，认真填写实验分析和实验心得以及问题和建议。
3．按时提交实验报告。
553299c f =⨯-⨯21(1)2 x<02()0x x f x ⎧+++⎪=≥一、实验目的1．熟练掌握关系表达式的使用。
2．熟练掌握使用else-if 语句实现多分支结构程序设计。
二、实验环境Windows XP ；Visual C++ 6.0。
三、实验内容1．输入x ，计算并输出下列分段函数sign(x)的值。
2．输入月薪salary ，输出应交的个人所得税tax （保留2位小数）。
按照2011年开始实行的新的个人所得税法，计算公式为：tax = rate*(salary-3500)-deduction 。
当salary≤3500时，rate=0、deduction=0当3500<salary≤5000时，rate=3%、deduction=0当5000<salary≤8000时，rate=10%、deduction=105当8000<salary≤12500时，rate=20%、deduction=555当12500<salary≤38500时，rate=25%、deduction=1005当38500<salary≤58500时，rate=30%、deduction=2775当83500<salary 时，rate=45%、deduction=135053．输入一个正整数n ，再输入n 个学生的百分制成绩，统计各等级成绩的个数。
原版Javascript程序设计实验报告
20—20学年第学期
Javascript程序设计实验
报告
系别：
专业：
班级：
姓名：
学号：
指导教师：
教务处制
实验项目：javascript基础实验
实验要求
1.请认真阅读下面的项目描述。
2.按照课程要求，每个班级分成若干个项目小组，每组人数大约2~3名同学，每个小组选出一名负责的同学。
请负责同学做好小组内分工。
项目描述
1.实验项目是围绕javascript综合实验平台展开，涉及到平台的外围设备。
2.要求熟练掌握javascript编程的基本流程。
3.要求熟练使用javascript编程软件平台。
4.项目实验内容由简单到复杂，采用循环渐进的引导方式，使学生在轻松的氛围中掌握javascript开发及软件的使用。
实验一Javascript基本页面操作
实验二Javascript基本功能控件练习
实验三文本框与函数参数传递练习
实验四HTML页面交互设计练习
实验五Javascript中text控件的设计练习
实验六Javascript中函数的简单应用
实验七Javascript控件交互设计
实验八新函数的学习与应用训练
实验九div与text的显示应用练习
实验十Javascript程序修改练习
实验十一Javascript中的函数应用练习
实验十二数字计算程序设计练习
实验十三数组程序设计练习
实验十四Javascript综合应用练习
实验十五选择结构程序设计练习
实验十六Javascript控件综合应用练习
实验十七循环结构程序设计练习
实验十八for结构程序设计练习。
java_实验四
武汉工程大学计算机科学与工程学院实验报告[4]一、实验目的1）了解Java.awt包及Swing的新增特性。
2）了解Java系统图形用户界面的基本工作原理和界面设计的步骤。
3）了解图形用户界面的事件响应机制。
4）了解图形用户界面各种常用的基本组件的使用方法及其事件响应。
5）掌握Java常用组件的应用6）掌握图形用户界面各种布局策略的设计与使用二、实验要求1）创建一个简单的图形用户界面。
2）使用Java类库中提供的组件，并添加到图形用户界面中。
3）通过Java事件处理机制，实现与用户的交互信息。
4）练习使用布局管理器5）使用Java类库中提供的组件，并添加到图形用户界面中6）通过Java事件处理机制，实现与用户的交互信息三、实验内容对一个优秀的应用程序来说，良好的图形用户界面是必不可少的。
因为只有通过图形界面，用户才能够与程序友好的交互。
缺少良好的图形用户界面，将会给用户理解和应用程序带来很多不便。
创建一个完整的图形用户界面的应用程序大致流程如下：1）引用需要的包和类。
2）设置一个顶层的容器。
3）根据需要为容器设置布局管理器或使用默认布局管理器。
4）将组件添加到容器内，位置自行设计。
5）为响应事件的组件编写事件处理代码。
下面我们先建立一个简单的图形用户界面，然后通过合适的布局管理器，将各种常用组件添加到图形用户界面中，实现与用户的交互信息。
1. 创建简单的图形用户界面（1）创建简单的 Frame 框架在 AWT 库中，对应于框架的类是 Frame，在 Swing 库中，相应的类是 JFrame。
JFrame 类扩展了 Frame 类。
大部分 ATWT 组件在 Swing 中都有等价的组件，它们在表示形式上差一个“J”。
框架是一个容器，可以在这个容器中放入其他一些图形用户界面组件，如按钮、菜单等，从而组成一个功能完善的程序。
下面的例子是建立一个 Frame 框架，通过这个例子学习如何创建一个简单的图形用户界面。
《Java程序设计》课程实验指导书2
《Java程序设计》课程实验指导书2《Java程序设计》课程实验指导书合肥学院计算机系实验⼀ Java开发环境的安装与配置,熟悉Java程序结构⼀、实验⽬的：1. 掌握JDK的安装步骤。
2. 理解环境变量PATH, CLASSPATH的作⽤，以及它们的设置⽅法。
3. 熟悉Editplus（或TextPad,JCreator）编辑环境，编写简单的Application程序和Applet程序，并编译和执⾏。
⼆、实验内容：熟悉JDK的安装和配置，学习如何编写并运⾏简单的Application程序和Applet程序(能输出⼀条简单的问候信息)；三、实验要求：1. 能正确地安装JDK2. 熟悉环境变量PATH, CLASSPATH设置⽅法，熟悉编辑环境3. 调试程序、编译，运⾏后得到正确的结果4. 写出实验报告，要求记录编译和执⾏Java程序当中的系统错误信息提⽰，并给出解决办法。
四、实验步骤：1．从/doc/9ac8d1619b6648d7c1c74618.html 上下载最新版本的JDK，并安装。
2．设置环境变量PATH, CLASSPATH, 使得Java程序能正确编译和执⾏。
3．在Editplus（或JCreator ，Textpad）环境下编写⼀个HelloWorld.jav a程序，(1)在主⽅法static public void main(String[ ] args)中调⽤System.out.println()⽅法，使程序输出⼀条问候信息；(2) 编译运⾏程序，观察运⾏情况和输出结果。
（使⽤JDK环境，调⽤javac.exe和java.exe编译和执⾏程序）4.在Editplus（或JCreator ，Textpad）环境下编写⼀个HelloWorldApplet.java程序。
(1)在public void paint(Graphics g)⽅法中调⽤g.drawString()⽅法，使程序输出⼀条问候信息；(2) 编译运⾏程序，观察运⾏情况和输出结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
空间后交-前交程序设计(实验报告)姓名：班级：学号：时间：空间后交-前交程序设计一、实验目的用 C 、VB或MATLAB语言编写空间后方交会-空间前方交会程序⑴提交实习报告：程序框图、程序源代码、计算结果、体会⑵计算结果：像点坐标、地面坐标、单位权中误差、外方位元素及其精度二、实验数据f=150.000mm，x0=0，y0=0三、实验思路1.利用空间后方交会求左右像片的外方位元素(1).获取m(于像片中选取两点，于地面摄影测量坐标系中选取同点，分别计算距离，距离比值即为m)，x,y,f,X,Y,Z(2).确定未知数初始值Xs,Ys,Zs,q,w,k(3).计算旋转矩阵R(4).逐点计算像点坐标的近似值（x），（y）(5).组成误差方程式(6).组成法方程式(7).解求外方位元素(8).检查是否收敛，即将求得的外方位元素的改正数与规定限差比较，小于限差即终止；否则用新的近似值重复步骤（3）-（7）2.利用求出的外方位元素进行空间前交，求出待定点地面坐标(1).用各自像片的角元素计算出左、右像片的方向余弦值，组成旋转矩阵R1,R2(2).根据左、右像片的外方位元素，计算摄影基线分量Bx，By，Bz(3).计算像点的像空间辅助坐标（X1,Y1,Z1）和（X2,Y2,Z2）(4).计算点投影系数N1和N2(5).计算未知点的地面摄影测量坐标四、实验过程⑴程序框图函数AandL%求间接平差时需要的系数%%%已知%a=像点坐标x,b=像点坐标y，f内方位元素主距%φ=q,ψ=w,κ=k%像空间坐标系X,Y,Z%地面摄影测量坐标系Xs,Ys,Zsfunction [A1,L1,A2,L2]=AandL(a,b,f,q,w,k,X,Y,Z,Xs,Ys,Zs) %%%%%%%%%%%选择矩阵元素a1=cos(q)*cos(k)-sin(q)*sin(w)*sin(k);a2=-cos(q)*sin(k)-sin(q)*sin(w)*cos(k);a3=-sin(q)*cos(w);b1=cos(w)*sin(k);b2=cos(w)*cos(k);b3=-sin(w);c1=sin(q)*cos(k)+cos(q)*sin(w)*sin(k);c2=-sin(q)*sin(k)+cos(q)*sin(w)*cos(k);c3=cos(q)*cos(w);%%%%%%%共线方程的分子分母X_=a1*(X-Xs)+b1*(Y-Ys)+c1*(Z-Zs);Y_=a2*(X-Xs)+b2*(Y-Ys)+c2*(Z-Zs);Z_=a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs);%%%%%%%近似值x=-f*X_/Z_;y=-f*Y_/Z_;%%%%%%%A组成L组成a11=1/Z_*(a1*f+a3*x);a12=1/Z_*(b1*f+b3*x);a13=1/Z_*(c1*f+c3*x);a21=1/Z_*(a2*f+a3*y);a22=1/Z_*(b2*f+b3*y);a23=1/Z_*(c2*f+c3*y);a14=y*sin(w)-(x/f*(x*cos(k)-y*sin(k))+f*cos(k))*cos(w);a15=-f*sin(k)-x/f*(x*sin(k)+y*cos(k)); a16=y;a24=-x*sin(w)-(y/f*(x*cos(k)-y*sin(k))-f*sin(k))*cos(w);a25=-f*cos(k)-y/f*(x*sin(k)+y*cos(k));a26=-x;lx=a-x;ly=b-y;%%%%%%%%%组成一个矩阵，并返回A1=[a11,a12,a13,a14,a15,a16];A2=[a21,a22,a23,a24,a25,a26];L1=lx;L2=ly;函数deg2dms%%%%%%%%角度转度分秒function y=deg2dms(x)a=floor(x);b=floor((x-a)*60);c=(x-a-b/60)*3600;y=a+(b/100)+(c/10000);函数dms2deg%%%%%度分秒转度function y=dms2deg(x)a=floor(x);b=floor((x-a)*100);c=(x-a-b/100)*10000;y=a+b/60+c/3600;函数ok%%%%%%%%%%%%%%目的是为了保证各取的值的有效值%%xy为n*1,a为1*nfunction result=ok(xy,a)format short gi=size(xy,1);for n=1:io=xy(n)-floor(xy(n,1));o=round(o*(10^a(n)))/(10^a(n));xy(n,1)=floor(xy(n,1))+o;endformat long gresult=xy;函数rad2dmsxy%%%%求度分秒表现形式的三个外方位元素，三个角度function xydms=rad2dmsxy(xy)[a,b,c,d,e,f]=testvar(xy);d=deg2dms(rad2deg(d));e=deg2dms(rad2deg(e));f=deg2dms(rad2deg(f));xydms=[a,b,c,d,e,f]';函数spacehoujiao%%%%%%%空间后交%%% f%%输入p（2*n,1）%%像点坐标x,y,X,Y,Z,均为（n,1）function [xy,m,R]=spacehoujiao(p,x,y,f,X,Y,Z)format long;%%%%%权的矢量化,这是等精度时的，如果非，将函数参数改为P P=diag(p);%%求nj=size(X,2);%%初始化Xs=0;Ys=0;Zs=0;for n=1:jXs=Xs+X(n);Ys=Ys+Y(n);Zs=Zs+Z(n);endSx=sqrt((x(2)-x(1))^2+(y(2)-y(1))^2);%%%%两像点之间距离Sd=sqrt((X(2)-X(1))^2+(Y(2)-Y(1))^2);%%%%两地面控制点之间距离m=Sd/Sx; %%%%图像比例系数Xs=Xs/j;Ys=Ys/j;Zs=m*f+Zs/j;m0=0;q=0;w=0;k=0;i=0;a=rand(2*j,6);l=rand(2*j,1);%%%%for n=1:j[a(2*n-1,:),l(2*n-1,1),a(2*n,:),l(2*n,1)]=AandL(x(n),y(n),f,q,w,k,X(n),Y(n),Z(n ),Xs,Ys,Zs);enddet=inv(a'*P*a)*transpose(a)*P*l;%%%%%%%%%循环体while 1%%%%%%%%%%%%%%%%[dXs,dYs,dZs,dq,dw,dk]=testvar(det);detXs=abs(dXs);detYs=abs(dYs);detZs=abs(dZs);detq=abs(dq);detw=abs(dw);detk=abs(dk);%%%%%%%%%if((detXs<0.01)&&(detYs<0.01)&&(detZs<0.01)&&(detq<pi/648000)&&(detw<pi/648000)&& (detk<pi/648000))break;elseV=(a*det-l);Q=inv(a'*P*a);m0=m0+sqrt((V'*P*V)/(2*j-6));%%m0需要每次的改正数算出来相加%%%Xs=Xs+dXs;Ys=Ys+dYs;Zs=Zs+dZs;q=q+dq;w=w+dw;k=k+dk;%%%for n=1:j[a(2*n-1,:),l(2*n-1,1),a(2*n,:),l(2*n,1)]=AandL(x(n),y(n),f,q,w,k,X(n),Y(n),Z(n ),Xs,Ys,Zs);enddet=inv(a'*P*a)*transpose(a)*P*l;i=i+1;%%%%end%%%end[dXs,dYs,dZs,dq,dw,dk]=testvar(det);detXs=abs(dXs);detYs=abs(dYs);detZs=abs(dZs);detq=abs(dq);detw=abs(dw);detk=abs(dk);V=(a*det-l);Q=inv(a'*P*a);m0=m0+sqrt((V'*P*V)/(2*n-6));%%%Xs=Xs+dXs;Ys=Ys+dYs;Zs=Zs+dZs;q=q+dq;w=w+dw;k=k+dk;%%%%%%%%%%%%%可以输出迭代次数的i%%%%%%%%%%%%Xs,Ys,Zs,q,w,k,i,dXs,dYs,dZs,dq,dw,dk,detXs,detYs,detZs %%%%%%%%%%%精度mo=m0*sqrt(Q);m=[mo(1,1),mo(2,2),mo(3,3),mo(4,4),mo(5,5),mo(6,6)]';[mXs,mYs,mZs,mq,mw,mk]=testvar(m);%%%%%%%%%输出xy=[Xs,Ys,Zs,q,w,k]';%%输出（6,1）的外方位元素m=[m0,mXs,mYs,mZs,mq,mw,mk]';%%单位误差，各元素中误差R=xyR(xy);%%旋转矩阵函数spaceqianjiao%%空间前交%输入f%输入x1,y1,x2,y2,R1,R2,xy1,xy2 （n,1）%输出X,Y,Z （n,1）function [X,Y,Z]=spaceqianjiao(x1,y1,x2,y2,f,R1,R2,xy1,xy2) i=size(x1,2);[Xs1,Ys1,Zs1,q1,w1,k1]=testvar(xy1);[Xs2,Ys2,Zs2,q2,w2,k2]=testvar(xy2);for n=1:i[X1(n),Y1(n),Z1(n)]=testvar(R1*[x1(n),y1(n),-f]');[X2(n),Y2(n),Z2(n)]=testvar(R2*[x2(n),y2(n),-f]');Bx=Xs2-Xs1;By=Ys2-Ys1;Bz=Zs2-Zs1;N1=(Bx*Z2(n)-Bz*X2(n))/(X1(n)*Z2(n)-X2(n)*Z1(n));N2=(Bx*Z1(n)-Bz*X1(n))/(X1(n)*Z2(n)-X2(n)*Z1(n));X(n)=Xs1+N1*X1(n);Z(n)=Zs1+N1*Z1(n);Y(n)=0.5*((Ys1+N1*Y1(n))+(Ys2+N2*Y2(n)));end函数testvar%分割矩阵。