MATLAB实验题答案

合集下载

matlab上机实验答案-整理版

matlab上机实验答案-整理版

第一次实验答案

1. 设

要求以0.01秒为间隔,求出y 的151个点,并求出其导数的值和曲线。 clc

clear

x=0:0.01:1.5;

y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)

y1=diff(y)

subplot(2,1,1)

plot(x,y)

subplot(2,1,2)

plot(x(1:150),y1)

2绘制极坐标系下曲线(a,b,n 自定数据)

clc

clear

a=10;

b=pi/2;

n=5;

theta=0:pi/100:2*pi;

rho=a*cos(b+n*theta);

polar(theta,rho)

3. 列出求下列空间曲面交线的程序

clc

clear

x=[-5:0.5:5];

[X,Y]=meshgrid(x);

z1=X.^2-2*Y.^2;

z2=X.*2-Y.*3;

xlabel('x')

ylabel('y')

zlabel('z')

surf(X,Y,z1)

hold on

surf(X,Y,z2)

k=find(abs(z1-z2)<0.5);

x1=X(k)

y1=Y(k)

z3=x1.^2-2*y1.^2

hold on

plot3(x1,y1,z3,'*')

⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2

212y x z -=y x z 322-=

4、设 ⎥⎦⎤⎢⎣

⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点,画出以x 为横坐标,y 为纵坐标的曲线,要求有图形标注。

MATLAB 课后实验答案

MATLAB 课后实验答案

实验一 MATLAB 运算基础

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

2sin851z e =+ (2) 221ln(1)2z x x =++,其中2120.455i x +⎡⎤=⎢⎥-⎣⎦

(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022

a a e e a z a a --+=++=--

(4) 2242011

122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0:0.5:2.5 解: M 文件:

z1=2*sin(85*pi/180)/(1+exp(2))

x=[2 1+2*i;-.45 5];

z2=1/2*log(x+sqrt(1+x^2))

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)

t=0:0.5:2.5;

z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 运算结果:

z1=2*sin(85*pi/180)/(1+exp(2))

x=[2 1+2*i;-.45 5];

z2=1/2*log(x+sqrt(1+x^2))

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)

MATLAB)课后实验答案

MATLAB)课后实验答案

实验一 MATLAB 运算基础

1、 先求下列表达式得值,然后显示MATLAB 工作空间得使用情况并保存全部变量。 (1) 012

2sin 851z e =+

(2) 21

ln(2z x =+,其中2120.45

5i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022

a a e e a z a a --+=++=--L (4) 224201

1

122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩

,其中t =0:0、5:2、5 解:

4、 完成下列操作:

(1) 求[100,999]之间能被21整除得数得个数。

(2) 建立一个字符串向量,删除其中得大写字母。

解:(1) 结果:

(2)、 建立一个字符串向量 例如:

ch='ABC123d4e56Fg9';则要求结果就是:

实验二 MATLAB 矩阵分析与处理

1、 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦

,其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵与对角阵,试通过数值计算验证22E R RS A O

S +⎡⎤=⎢

⎥⎣⎦。

解: M 文件如下;

5、 下面就是一个线性方程组:

1231112340.951

110.673450.521

11456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣

⎦⎢⎥⎢⎥⎣⎦

(1) 求方程得解。

(2) 将方程右边向量元素b 3改为0、53再求解,并比较b 3得变化与解得相对变化。

MATLAB)课后实验答案

MATLAB)课后实验答案

实验一 MATLAB运算基础 1. 先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0122sin851ze (2) 221ln(1)2zxx,其中2120.455ix (3) 0.30.330.3sin(0.3)ln,3.0,2.9,,2.9,3.022aaeeazaaL (4) 2242011122123ttzttttt,其中t=0:0.5:2.5 解: M文件: z1=2*sin(85*pi/180)/(1+exp(2)) x=[2 1+2*i;-.45 5]; z2=1/2*log(x+sqrt(1+x^2)) 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) t=0:0.5:2.5; z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 4. 完成下列操作: (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量,删除其中的大写字母。 解:(1) 结果: m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 (2). 建立一个字符串向量 例如: ch='ABC123d4e56Fg9';则要求结果是: ch='ABC123d4e56Fg9'; k=find(ch>='A'&ch<='Z'); ch(k)=[] 实验二 MATLAB矩阵分析与处理 1. 设有分块矩阵33322322ERAOS,其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵,试通过数值计算验证22ERRSAOS。 解: M文件如下; 5. 下面是一个线性方程组: 1231112340.951110.673450.52111456xxx ch = 123d4e56g9 (1) 求方程的解。 (2) 将方程右边向量元素b3改为0.53再求解,并比较b3的变化和解的相对变化。 (3) 计算系数矩阵A的条件数并分析结论。 解: M文件如下: 实验三 选择结构程序设计 1. 求分段函数的值。 2226035605231xxxxyxxxxxxx且且及其他 用if语句实现,分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。 解:M文件如下: 2. 输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。其中90分~100分为A,80分~89分为B,79分~79分为C,60分~69分为D,60分以下为E。 要求: (1) 分别用if语句和switch语句实现。 (2) 输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。 解:M文件如下 3. 硅谷公司员工的工资计算方法如下: (1) 工作时数超过120小时者,超过部分加发15%。 (2) 工作时数低于60小时者,扣发700元。 (3) 其余按每小时84元计发。 试编程按输入的工号和该号员工的工时数,计算应发工资。 解:M文件下 实验四 循环结构程序设计 1. 根据2222211116123nL,求π的近似值。当n分别取100、1000、10000时,结果是多少? 要求:分别用循环结构和向量运算(使用sum函数)来实现。 解:M文件如下: 运行结果如下: 2. 根据11113521ynL,求: (1) y<3时的最大n值。 (2) 与(1)的n值对应的y值。 解:M—文件如下: 3. 考虑以下迭代公式: 1nnaxbx 其中a、b为正的学数。 (1) 编写程序求迭代的结果,迭代的终止条件为|xn+1-xn|≤10-5,迭代初值x0=1.0,迭代次数不超过500次。 (2) 如果迭代过程收敛于r,那么r的准确值是242bba,当(a,b)的值取(1,1)、(8,3)、

MATLAB程序实验参考答案

MATLAB程序实验参考答案

MATLAB程序设计与应用(第二版)实验参考答案%实验一MATLAB运算基础

%第一题

%(1)

z1=2*sin(85*pi/180)/(1+exp(2))

%(2)

x=[2,1+2i;-0.45,5];

z2=0.5*log(x+sqrt(1+x.^2))

%(3)

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)

%(4)

t=0:0.5:2.5;

z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)

%第二题

A=[12 34 -4;34 7 87;3 65 7];

B=[1 3 -1;2 0 3;3 -2 7];

A+6*B

A-B+eye(size(A))

A*B

A.*B

A^3

A.^3

A/B

B\A

[A,B]

[A([1,3],:);B^2]

%第三题

A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]

C=A*B

F=size(C)

D=C(F(1)-2:F(1),F(2)-1:F(2))

whos

%第四题

%(1):

A=100:999;

B=rem(A,21);

C=length(find(B==0))

%(2):

A='lsdhKSDLKklsdkl';

MATLAB)课后实验答案

MATLAB)课后实验答案

实验一 MATLAB 运算基础

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

2sin851z e =+ (2) 221

ln(1)2z x x =++,其中2120.45

5i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022

a a e e a z a a --+=++=-- (4) 2242011

122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩

,其中t =0:0.5:2.5 解: M 文件:

z1=2*sin(85*pi/180)/(1+exp(2))

x=[2 1+2*i;-.45 5];

z2=1/2*log(x+sqrt(1+x^2))

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) t=0:0.5:2.5; z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)

运算结果:

z1=2*sin(85*pi/180)/(1+exp(2))

x=[2 1+2*i;-.45 5];

z2=1/2*log(x+sqrt(1+x^2))

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)

MATLAB全部实验及答案

MATLAB全部实验及答案

MATLAB全部实验及答案

实验一、MATLAB基本操作

实验内容及步骤

1、命令窗口的简单使用

(1)简单矩阵的输入

(2)求[12+2×(7-4)]÷32的算术运算结果

2、有关向量、矩阵或数组的一些运算

(1)设A=15;B=20;求C=A+B与c=a+b?

(2)设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与

A.*B?

A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘(3)设a=10,b=20;求i=a/b=0.5与j=a\b=2?

(4)设a=[1 -2 3;4 5 -4;5 -6 7]

请设计出程序,分别找出小于0的矩阵元素及其位置(单下标、全

下标的形式),并将其单下标转换成全下标。

clear,clc

a=[1 -2 3;4 5 -4;5 -6 7];

[x,y]=find(a<0);

c=[];

for i=1:length(x)

c(i,1)=a(x(i),y(i));

c(i,2)=x(i);

c(i,3)=y(i);

c(i,4)=(y(i)-1)*size(a,2)+x(i);

end

c

(5)在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?前面那

个是虚数矩阵,后面那个出错

(6)请写出完成下列计算的指令:

a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?

a^2= 22 16 16

25 26 23

26 24 28

matlab实验报告答案

matlab实验报告答案

matlab实验报告答案

篇一:参考答案Matlab实验报告

实验一Matlab基础知识

一、实验目的:

1. 熟悉启动和退出Matlab的方法。

2. 熟悉Matlab命令窗口的组成。

3. 掌握建立矩阵的方法。

4. 掌握Matlab各种表达式的书写规则以及常用函数的使

用。

二、实验内容:

1. 求[100,999]之间能被21整除的数的个数。(rem)

2. 建立一个字符串向量,删除其中的大写字母。(find)

3. 输入矩阵,并找出其中大于或等于5的元素。(find)

4. 不采用循环的形式求出和式i?1

三、实验步骤:

? 求[100,199]之间能被21整除的数的个数。(rem)

1. 开始→程序→Matlab

2. 输入命令:

?m=100:999;

?p=rem(m,21);

?q=sum(p==0)

ans=43

? 建立一个字符串向量,删除其中的大写字母。(find)

1. 输入命令:

?263i的数值解。(sum)

?k=input( ’,’s’);

Eie48458DHUEI4778

?f=find(k =’A’ k =’Z’);

f=9 10 11 12 13

?k(f)=[ ]

K=eie484584778

? 输入矩阵,并找出其中大于或等于5的元素。(find)

1. 输入命令:

?h=[4 8 10;3 空间。

2. 熟悉M文件与M函数的编写与应用。

3. 熟悉Matlab的控制语句。

4. 掌握if,switch,for等语句的使用。

二、实验内容:

1. 根据y=1+1/3+1/5+??+1/(2n-1),编程求:y 5时最

Matlab上机实验答案

Matlab上机实验答案

Column 61
0.3388 + 0.0000i
t 2 0 t 1 2 (4) z4 t 1 1 t 2 ,其中 t=0:0.5:2.5 2 t 2t 1 2 t 3
>> t=0:0.5:2.5; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^2-2.*t+1 ) z4 =
>> x=[2 1+2i;-0.45 5]; >> z2=1/2*log(x+sqrt(1+x^2))
z2 =
0.7114 - 0.0253i 0.2139 + 0.9343i
0.8968 + 0.3658i 1.1541 - 0.0044i
e0.3a e0.3a 0.3 a sin(a 0.3) ln , a 3.0, 2.9, , 2.9, 3.0 (3) z3 2 2
0.2000 0.1667 0.1429 0.1250 0.1111
1 1 1 1 1
1 2 3 4 5
1 3 6 10 15
1 4 10 20 35
1 5 15 35 70

MATLAB实验答案

MATLAB实验答案

2,

(1)A=[1 2;3 4 ];

B=[5 5;7 8 ];

A^2*B

(2) A=[1 2 3;4 5 6;7 8 9 ];

B=[1 0 0;0 2 0;0 0 3 ];

A\B,A/B

(3) A=[5+1i,2-1i,1;6*1i,4,9-1i ]; A1=A.',A2=A'

(4) A=[1 2 3;4 5 6;7 8 9 ];

B=A([1,2],[3]),C=A(2:end, : )

(5) M=magic(4),M(:,4)=[]

3,

(1) p1=[1 0 2 4];PS1=poly2str(p1,'x')

r=roots(p1)

(2) p=poly(A)

polyval(p,20)

poly2str(p,’A’)

4,

(1) t=0:.1:2*pi;

y=cos(t);

plot(t,y),grid

(2) t=0:.1:2*pi;

y1=cos(t-0.25);

y2=sin(t-0.25);

plot(t,y1,t,y2)

grid

t=0:.01:4*pi;

y=10*sin(t);

plot(t,y,'r:+')

grid

axis([0,10,-15,15])

xlabel('x','FontSize',16)

ylabel('y','FontSize',16)

title('正弦函数')

6,

(1)

sum=0;i=1;

while(sum<2000)

sum=sum+i;

i=i+1;

end

n=i-2

{ n=0;s=0;

while s<2000

n=n+1;

s=s+n;

end

n=n-1

}

(2)

解法1:function [ s i ] = mitifang1(n)

MATLAB-课后实验答案

MATLAB-课后实验答案

实验一 MATLAB 运算基础

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

2sin851z e =+

(2) 21ln(2z x =,其中2120.455i x +⎡⎤=⎢⎥-⎣⎦

(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022

a a e e a z a a --+=++=--

(4) 2242011

122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩,其中t =0:: 解:

2. 已知:

1234413134787,2033657327A B --⎡⎤⎡

⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦

求下列表达式的值:

`

(1) A+6*B 和A-B+I (其中I 为单位矩阵)

(2) A*B 和A.*B

(3) A^3和A.^3

(4) A/B 及B\A

(5) [A,B]和[A([1,3],:);B^2] 解:

运算结果:

3. 设有矩阵A 和B

12

3

4

53

166789101769,11

121314150

23416171819209

7021

2223242541311A B ⎡⎤⎡⎤

⎢⎥⎢⎥-⎢⎥⎢

⎥⎢⎥⎢⎥==-⎢⎥⎢

⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦

(1) 求它们的乘积C 。

(2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。

解:. 运算结果:

4. 完成下列操作:

(1) 求[100,999]之间能被21整除的数的个数。

(2) 建立一个字符串向量,删除其中的大写字母。解:(1) 结果:

MATLAB)课后实验答案-精简版

MATLAB)课后实验答案-精简版

实验一 MATLAB 运算基础

1. 先求下列表达式的值,然后显示MA TLAB 工作空间的使用情况并保存全部变量。 (1) 0

122sin 851z e

=+

(2) 21ln(2

z x =+

,其中2120.45

5i x +⎡⎤

=⎢

⎥-⎣⎦

(3) 0.30.330.3sin(0.3)ln

,

3.0, 2.9,,2.9,3.02

2

a

a

e

e a z a a --+=

++=--

(4) 2

24201

1

122123

t t z t t t t t ⎧≤<⎪

=-≤<⎨⎪-+≤<⎩

,其中t =0:0.5:2.5

2. 已知:

123441

3134

787,2

033

65

73

2

7A B --⎡⎤

⎡⎤

⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦

⎣⎦

求下列表达式的值:

(1) A+6*B 和A-B+I (其中I 为单位矩阵) (2) A*B 和A.*B

(3) A^3和A.^3

(4) A/B及B\A

(5) [A,B]和[A([1,3],:);B^2]

3. 设有矩阵A 和B

123453

0166789101769,1112

1314150

234

16171819209702122

23

24

254

13

11A B ⎡⎤

⎡⎤

⎢⎥⎢⎥-⎢

⎥⎢⎥⎢⎥⎢⎥==-⎢

⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦

⎣⎦

(1) 求它们的乘积C 。

(2) 将矩阵C 的右下角3×

2子矩阵赋给D 。 (3) 查看MA TLAB 工作空间的使用情况。 4. 完成下列操作:

(1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量,删除其中的大写字母。 (2). 建立一个字符串向量 例如:

Matlab实验习题集答案

Matlab实验习题集答案

1:用以上两种形式计算36sin 5e ++π算术运算结果。

>> 5^6+sin(pi)+exp(3)

ans =

1.5645e+004

>> x=5^6+sin(pi)+exp(3)

x = 1.5645e+004

2:已知矩阵

⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=2211,2121B A ,对它们做简单的关系与逻辑运算C=(A<B)&(A= =B)

>> A=[1 2;1 2];

>> B=[1 1;2 2];

>> C=(A<B)&(A==B) C =

0 0

0 0

3:对数7sin 5+=a 用五位定点、十五位定点以及有理数形式表示出来。

>> a=5+sin(7);

format short,a

a =

5.6570

>> a=5+sin(7);

>> format long,a

a =

5.656986598718789

4:直接输入创建矩阵

⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=98760154321A

>> A=[1 2 3;4 15 60;7 8 9]

A =

1 2 3

4 1

5 60

7 8 9

5:输入矩阵111111111⎛⎫ ⎪ ⎪ ⎪⎝

⎭。 %利用MATLAB 命令直接输入矩阵

OneMatrix=ones(3,3,1)

OneMatrix =

1 1 1 1 1 1 1 1 1

6:输入矩阵

00000 00000⎛⎫ ⎪⎝⎭

>> OneMatrix=ones(2,5,1);

ZeroMatrix=zeros(size(OneMatrix))

MATLAB数学实验答案(全)

MATLAB数学实验答案(全)

第一次练习

教学要求:熟练掌握Matlab 软件的基本命令和操作,会作二维、三维几何图形,能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令

vpa(x,n) 显示x 的n 位有效数字,教材102页

fplot(‘f(x)’,[a,b]) 函数作图命令,画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位(1-9班)或4位(10班以上) 1.1 计算30sin lim

x mx mx x →-与3

sin lim

x mx mx

x →∞- syms x

limit((902*x-sin(902*x))/x^3) ans =

366935404/3

limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos

1000

x

mx

y e =,求''y syms x

diff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =

(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算

22

11

x y e

dxdy +⎰⎰

dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.1394

1.4 计算4

2

2

4x dx m x +⎰ syms x

int(x^4/(902^2+4*x^2))//不定积分 ans =

(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//高阶导数

参考答案Matlab实验报告

参考答案Matlab实验报告

实验一 Matlab基础知识

一、实验目的:

1.熟悉启动和退出Matlab的方法。

2.熟悉Matlab命令窗口的组成。

3.掌握建立矩阵的方法。

4.掌握Matlab各种表达式的书写规则以及常用函数的使

用。

二、实验内容:

1.求[100,999]之间能被21整除的数的个数。(rem)

2.建立一个字符串向量,删除其中的大写字母。(find)

3.输入矩阵,并找出其中大于或等于5的元素。(find)

4.不采用循环的形式求出和式

63

1

2i

i=

的数值解。(sum)

三、实验步骤:

●求[100,199]之间能被21整除的数的个数。(rem)

1.开始→程序→Matlab

2.输入命令:

»m=100:999;

»p=rem(m,21);

»q=sum(p==0)

ans=43

●建立一个字符串向量,删除其中的大写字母。(find)

1.输入命令:

»k=input('’,’s’);

Eie48458DHUEI4778

»f=find(k>=’A’&k<=’Z’);

f=9 10 11 12 13

»k(f)=[ ]

K=eie

●输入矩阵,并找出其中大于或等于5的元素。(find)

1.输入命令:

»h=[4 8 10;3 6 9; 5 7 3];

»[i,j]=find(h>=5)

i=3 j=1

1 2

2 2

3 2

1 3

2 3

●不采用循环的形式求出和式的数值解。(sum)

1.输入命令:

»w=1:63;

»q=s um(2.^w)

q=1.8447e+019

- 1 -

实验二 Matlab 基本程序

一、 实验目的:

1. 熟悉Matlab 的环境与工作空间。

Matlab实验题目与答案

Matlab实验题目与答案

(温馨提示:实验课结束后,请将所有作业(题目、代码、结果)利用word 整

理成一个完整的实验报告,加上封面,打印,纸质档于18周周一交)

第一次上机作业

目的:

1. 掌握MATLAB 各种表达式的书写规则

2. 运行课堂上讲过的例子,熟悉矩阵、表达式的基本操作和运算。

作业:

1. 熟悉matlab 集成环境界面。回答以下问题,并操作相关的指令:

(1) 分别写出清除命令窗口和清除变量的指令。

答: clc 和clear

(2)在命令行输入命令后,matlab 的搜索过程是怎样的?

答: (1)检查该命令对象是不是一个变量。

(2)检查该命令对象是不是一个内部函数。

(3)检查该命令对象是否为当前目录下的程序文件。

(4)检查该命令对象是否为MATLAB 搜索路径中其他目录下的M 文件。

(3)什么是matlab 的当前工作目录?写出两种设置当前工作目录的方法? 答: 就是matlab 当前文件读取和存储的默认路径

(1)在当前目录窗口中更改

(2)在MATLAB 桌面工具栏中更改

(3)使用cd 命令:cd c:\mydir---将c :\mydir 设置为当前目录

(4)什么是matlab 的搜索路径?写出两种设置搜索路径的方法?

答: 指Matlab 运行文件时进行搜索的目录。

(1)用path 命令设置:

(2)用Set Path 对话框设置

(5)help 命令和doc 命令有什么作用,它们有什么区别?

答: help 命令:最基本的帮助命令,查询信息直接显示在命令窗口。

doc 命令:在帮助窗口中显示HTML 帮助文档,显示函数的详细用法及 例子,比help 命令更详细。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

result5 =

( 1 ) a = 1 : 2 : 5

a =

1 3 5

( 2 ) b = [ a' , a' , a' ;a ]

b =

1 1 1

3 3 3

5 5 5

1 3 5

( 3 ) c = a + b ( 2 , : )

c =

4 6 8

2、下列运算是否合法,为什么如合法,

结果是多少

>> result2=a*b

Error using *

Inner matrix dimensions must agree. >> result3=a+b

result3 =

3 6 2

58 11

>> result4=b*d

result4 =

31 22 22

40 49 13

31 22 22

40 49 13

-5 -8 7

>> result6=a.*b result6 =

2 8 -3

415 30

>> result7=a./b result7 =

>> result8=

Attempt to reference field of non-structure array.

>> result9=a.\b result9 =

>> result10=a92

result10 =

1 4 9

16 25 36

>> resultl 1=29a

result11 =

2 4 8

16 32 64

>>

result5=[b;c']*d 3、用MATLAB求解下面的的方程组。

1、求以下变量的值,并在MATLAB^验证。

1 2 x1

3 2 x2

11 5 x3

2 1

3 x4

>> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0]

>> B=B'

>> x=inv(A)*B

>> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5]

>> x2=inv(A1)*B2

7 2 1 2

9 15 3 2

2 2 11 5

1 3

2 13

(1)求矩阵A的秩(rank)

(2)求矩阵 A 的行列式(determinant)

(3)求矩阵 A 的逆(inverse)

(4)求矩阵 A 的特征值及特征向量

(eigenvalue and eigenvector)

>> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3)

>> a=det(A3) >> [V,D]=eig(A3)

10

n 10

查看y 的值)

m1=0;

for m=-10:10 m仁m1+2^m;

end

m1

m1 =

6、求分段函数的值。

用if 语句实现,算出下列表中x 对应的y 值。

x=input('enter x=');

if x<0

y=x A2+x-6;

elseif x>=0&&x<5

y=xA2-5*x+6;

else

y=xA2-x-1;

end

y

7、分别用if 和switch 语句实现,将百分

制成绩转换为成绩等级A、B、C、D、E。

其中90~1 00分为A,80~89 分为B,70~79

分为C,60~69 分为D,60 分以下为E。

对超出百分制范围的成绩,给出错误提示

信息。

if 结构程序:

x=input('please enter score=');

if x>=90&&x<=100

9

2

10

disp('A')

7 2

9 15

(1)

2 2

1 3 4 7 1 0

A

4、已知

2n 2 10 29

elseif x<80&&x>=70 otherwise

disp('C')

disp('error')

elseif x<70&&x>=60

end

disp('D')

8、思考题

elseif x<60&&x>=0

设计程序,完成成两位数的加、减、乘、除

disp('E') 四则运算,即产生两个两位随机整数,再输入一个运算符号,做相应的运算,

else 并显示相应的结果。

disp('error')

x=input(' 请输入运算符')

end

a=num2str(floor(rand(1)*90+10));

switch 结构程序:

a

x=input('please enter sco re=');

b=num2str(floor(rand(1)*90+10));

switch fix(x/10)

b

case{10,9}

if x=='+'

if x>100

y=a+b;

disp('error')

elseif x=='-'

else

y=a-b;

elseif x=='*'

disp('A')

y=a*b;

end

elseif x=='/'

case{8}

y=a/b;

disp('B')

else

case{7}

disp('error')

disp('C')

end

case{6}

y

disp('D') 9、启动MATLAB后,点击File|New|M-File ,

启动MATLAB的程序编辑及调试器

case{0,1,2,3,4,5} ( Editor/Debugger ),编辑以下程序,点

击File|Save 保存程序,注意文件名

相关文档
最新文档