matlab上机报告

课程名称：测绘程序设计 课程性质：必修

测绘程序设计 实习报告

学 院： 测绘工程学院 班 级： 实习地点：

1220501

核工楼410

班级组号： 学生姓名：

赵建新

学生学号： 201220050163 指导教师： 王建强

日期： 年 月 日

目 录

实验1 matlab软件练习 ................................................................... 2 实验2 matlab解算方程 ................................................................... 5 实验3 matlab绘图 ........................................................................... 8 实验4 角度弧度互换 ............................................

matlab上机指导练习

MATLAB

上机指导书

matlab上机指导练习

上机一 MATLAB编程环境

一、上机目的

1.熟悉MATLAB的操作环境及基本操作方法 2.熟悉MATLAB的通用参数设置

3.熟悉 MATLAB的搜索路径及设置方法 4.熟悉MATLAB帮助信息的查阅方法

二、上机内容和结果

1．利用菜单设置MATLAB的Command Window中字体大小，并更改输出格式。

示例：

matlab上机指导练习

结果：

2.在硬盘上创建以自己名字命名的文件夹，将当前路径修改为此文件夹 示例：

结果：

4．完成下列操作：

（1）在MATLAB命令窗口输入以下命令 >> x=0:pi/10:2*pi; >> y=sin(x);

（2）在工作空间选择变量y,在在工作空间窗口选择绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形。

matlab上机指导练习

结果：

5.利用帮助学习save,load命令的用法，将在工作区中变量全部保存在mypath.mat中，清空工作区，重载变量x,y查看变量信息，并把它们保存在mypath.mat中

结果：

matlab上机指导练习

6. 计算y=1.3^3*sin(pi/3)*sqrt(26) (1)结果用f

河北科技大学教案用纸

实验一 熟悉Matlab工作环境

1、熟悉Matlab的5个基本窗口 思考题：

（1）变量如何声明，变量名须遵守什么规则、是否区分大小写。

答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则：

? 变量名必须以字母开头，只能由字母、数字或下划线组成。 ? 变量名区分大小写。 ? 变量名不能超过63个字符。 ? 关键字不能作为变量名。

? 最好不要用特殊常量作为变量名。 （2）试说明分号、逗号、冒号的用法。

分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。

逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。 冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。 （3）linspace（）称为“线性等分”函数，说明它的用法。 LINSPACE Linearly spaced vector. 线性等分函数 LINSPACE(X1, X2) generates a row vector of 100 linearly equally spa

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后

1、画出z?sin(x2?y2)x?y22所表示的三维曲面，x,y的取值范围是[?8,8]。

要求采用不同的方法生成坐标网格，并注意如何回避分母不为零的问题。 方法1：

clear;x=-8:0.5:8; y=x'; X=ones(size(y))*x; Y=y*ones(size(x)); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; mesh(X,Y,Z); colormap(hot) xlabel('x'),ylabel('y'),zlabel('z') 方法2： x=-8:0.5:8;

[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; surf(X,Y,Z) surf(x,y,Z) 方法3： x=-8:0.5:8; y=x;

for i=1:length(y)

for j=1:length(x)

r=sqrt(x(j)^2+y(i)^2)+eps; z(i,j)=sin(r)/r; end end

surf(x,y,z);

2、通过M脚本文件，画出下列分段函数所表示的曲面

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输

第四次 上机作业

1、 从键盘输入一个4位整数，按照如下规则加密后输出。加密规则：每位数字

都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 Clear X=ones(1,4);

X（1）=input(’输入第一位：‘)； X（2）=input(’输入第二位：‘)； X（3）=input(’输入第三位：‘)； X（4）=input(’输入第四位：‘)； X=rem(7+x,10);

Y=1000.*x(3)+100.*x(4)+10.*x(1)+x(2)

2、 分别用if 和switch语句实现以下计算，其中a、b、c的值从键盘输入。

??ax2?bx?c,?y??asincb?x,?c?lnb?,x?0.5?x?1.51.5?x?3.5 3.5?x?5.5a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输入c：’)； If(x>=0.5&&x<=1.5) y=a.*x^2+b.*x+c Elseif(x>=1.5&&x<=3.5) y=a.*(sin(b))^c+x

Elseif(x>=3.5&&x<=5.5) y=log(abs(b+c./x)) end

a=input(‘请输入a：’)； b=input(‘请输入b：’)； c=input(‘请输

MATLAB上机作业1

对以下问题,编写M文件：

(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。 function f=qipaofa(x) for j=9:-1:1 for i=1:j

if(x(i)>x(i+1))

t=x(i);x(i)=x(i+1);x(i+1)=t; end end end f=x

x=round(10*rand(1,10)) qipaofa(x);

(2) 有一个4×5矩阵，编程求出其最大值及其所处的位置。 function f=zuidazhi(x) a=1;b=1;c=x(1,1); for i=1:4 for j=1:5

if x(i,j)>c

a=i;b=j;c=x(i,j); end end end

f=[c,a,b]

x=rand(4,5) zuidazhi(x)

20(3) 编程求?n!。

n?1function f=qiuhe(x) b=0; for i=1:x a=prod(1:i); b=b+a; end f=b

qiuhe(20)

(4)一球从100米高度自由落下,每次落地后

三、假设已知矩阵A，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给B矩阵，用A?magic(8)命令生成A矩阵，用上述命令检验一下结果是不是正确。

>> A=magic(8)

>> B=A(2:2:end,:)

五、选择合适的步距绘制出下面的图形。

（1）sin(1/t)，其中t?(?1,1)；（2）sin(tant)?tan(sint)，其中t?(??,?)。 1．>> t=[-1:0.0001:1];y=sin(1./t);plot(t,y) Warning: Divide by zero. >>

2．>> t=[-pi:0.001:pi];y=sin(tan(t))-tan(sin(t));plot(t,y) >>

七、试求出如下极限。 （1）lim(3?9)；（2）limx??xx1xxyxy?1?1x?0y?0；（3）lim1?cos(x2?y2)(x?y)e22x2?y2x?0y?0。

（1）>>symsx;f=(3^x+9^x)^(1/x);limit(f,x,inf) ans =9

（2）>>syms x y;f=x*y/(sqrt(x*y+1)-1);limit(limit(f,x,0),y,0) ans =2