2023年matlab平时作业共享

第一章习题1

3.请指出如下5个变量名中，哪些是合法的?abed-2,xyz_3,ABCDefgh

»abcd_2xyz_33»chan»a变量ABCDefgh

6.指令clear,elf,clc各有什么用处？

Clear为清除工作空间变量指令,涉及内存中的变量和数据；

Clc为清命令窗口指令，而保存工作空间内容；

Clf为清除当前图形指令。

7.以下两种说法对吗？(1)“MATLAB的数值表达精度与其指令窗中的数据显示精度

相同。"(2)MATLAB指令窗中显示的数值有效位数不超过7位。”

解：(1)错误；

(2)错误，该显示数值有效位可以在菜单中设立或通过vpa(X,n)进行设立。

'123'

8.想要在MATLAB中产生二维数组S=456,下面哪些指令能实现目的？

789

S=[l,2,3；4,5,6;7,8;9]

S=[l23;456;789]

S=[l,2,3；4,5,6;7,8,9]

解:S=[123;456；789]；

S=[1,2,3；4,5,6;7,8,9]。

第三章习题3

1.规定在闭区间[0,2万]上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。

解：a=linspace(0,2*pi,10)

a=0：2*pi/9:2*pi

2.由指令rand('state',0),A=rand(3,5)生成二维数组A,试求该数组中所有大于0.5的

元素的位置，分别求出它们的“全下标”和“单下标”。

解：rand('state',0),A=rand(3,5);

L(A>0.5)=NAN

程序如下：

%注释:最终求出的该数组中所有大于0.5的元素的“全下标”信息存储于C中

%%%%“单下标”信息存储于D中

clear;

ii=0；jj=0;x=l;

C=zeros(15z2);D=zeros(1,15);

rand(^tate*,0),A=rand(3,5)；

forii=l:3

forjj=l:5

ifA(iizjj)>0.5

C(x,:)=[ii,jj];

D(x)=(ii-l)*5+jj;

x=x+l;

end

end

end

C,D

结果如下：

C=

11

15

22

24

25

31

32

33

34

00

00

00

00

00

00

D=

157910111213140

00000

5.在时间区间［0,10］中，绘制y=l-eScos2♦曲线。规定分别采用“标量循环运

算法”和“数组运算法”编写两段程序绘图。

解：（1）标量循环运算法

plot(t,y)

fort=0:0.01:10

y=l—exp(—0.5*t)*cos(2*t);

plot(t,y);

grid;

holdon;

(2)数组运算法

t=0:0.01:10；

y=I-exp(-0.5*t).*cos(2*t);

plot(t,y)

如图:

6.先运营clear,format1ong,rand('state',1),A=rand(3,3),然后根据A写出

两个矩阵：一个对角阵B,其相应元素由A的对角元素构成;另一个矩阵C,其对角元素全为

0,而其余元素与相应的A阵元素相同。

A=(0.620.5980.34

0.000.070.592

0)

B={0.6200

00.070

000.60}

C={00.5980.34

0.0000.592

0.230.130)

第四章习题4

.3

3.求函数/(%)=6仙=的数值积分s=「/(为公，并请采用符号计算尝试复算。

J0

解：

(1)数值积分：

clear

symsx

x=0:0.l:pi;

y=exp(sin(x).A3);

»trapz(x,y)

ans=

5.0954

(2)符号计算：

symsx

Isym=vpa(int(exp(sin(x)A3),x,0,pi))

Warning:Exp1icitintegralcou1dnotbefbund.

>Iat58

Isym=

5.

4.用quad求取J；：e-W|sinx\dx的数值积分，并保证积分的绝对精度为104。

解：Y=quadCexp(-abs(x)).*abs(sin(x))',-5*pi,1.7*pi,le-9)

Y=

1.79

5.求函数八。=($皿5。26°。城一1.5北0$2£+1.81+0.5|在区间［—5,5］中的最小值点。

解：

clear

symsx

y=((sin(5*x))八2)*(exp(0.06*xA2))-1.5*x*cos(2*x)+1.8*abs(x+0.5);

yd=diff(y,x);。。。

a=so1ve(yd)。。

yd_a=vpa(subs(yd,x,xs0),6)。

y_a=vpa(subs(y,x,xs0),6>

y_m_5=vpa(subs(y,x,-5),6)

y_P_5=vpa(subs(y,x,5),6)

a=

e-2

yd_a=

.1e-30

y_a=

.899099

y_m_5=

1.88547

y_p-5=

16.2715

最小值点即为a点.

第五章习题5

1.已知椭圆的长、短轴。=4/=2(见图p5-l),用“小红点线”画椭圆《~.

y-bsint

图p5-l

解:如图:

3.A,B,C三个城市上半年每月的国民生产总值见表p5.1。试画出如图p5—3所示的三城市

上半年每月生产总值的累计直方图。（提醒:使用指令bar）

表p5.1各城市生产总值数据（单位:亿元）

城市1月2月3月4月5月6月

A170120180200190220

B120100110180