MATLAB上机体会

河北科技大学教案用纸

实验一 熟悉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

三、假设已知矩阵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

Matlab上机实验答案

实验一 MATLAB运算基础

1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。

2sin850(1) z1? 1?e2>> z1=2*sin(85*pi/180)/(1+exp(2)) z1 =

0.2375 (2) z2?1?2i??21 ln(x?1?x2)，其中x???5?2??0.45>> x=[2 1+2i;-0.45 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 =

0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i

e0.3a?e?0.3a0.3?az?sin(a?0.3)?ln,a??3.0,?2.9,(3) 322,2.9,3.0

>> a=-3.0:0.1:3.0;

>> z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) (>> z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) 可