matlab平时作业共享

第一章习题1

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

ABCDefgh

abcd-2xyz_33chana变量ABCDefgh

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

Clear为清除工作空间变量指令，包括内存中的变量和数据；

Clc为清叮嘱窗口指令，而保留工作空间内容；

Clf为清除当前图形指令。

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

的数据显示精度相同。"(2)MATLAB指令窗中显示的数值有效位

数不超过7位。”

解：(1)错误；

(2)错误，该显示数值有效位可以在菜单中设置或通过

vpa(X,n)进行设置。

-123一

&想要在MATLAB中产生二维数组5=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=[l23;456;789];

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

第三章习题3

1.要求在闭区间。24]上产生具有1。个等距采样点的一维数组。试用两

种不同的指令实现。

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

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

2.由指令rand(statd,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(15,2);D=zeros(l,15);

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

forii=l:3

forjj=l:5

ifA(ii,jj)>0.5

C(x,:)=[ii,jjJ;

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]^,绘制)，=1-Hos2/曲线。要求分别实行“标量

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

解：(1)标量循环运算法

plot(t,y)

fort=0:0.01:10

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

Plot(t,y);

grid;

holdon;

end如图:

(2)数组运算法

t=0:0.01:10;

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

plot(t,y)

如图：

6.先运行clear,formatlong,rand('state',l),A=rand(3,3),然后依据A

写出两个矩阵：一个对角阵B,其相应元素由A的对角元素构成；另一个

矩阵C,其对角元素全为0,而其余元素与对应的A阵元素相同。

A={0.952782149656620.5980.83681960067634

0.7000.8470.592

0.953877473592230.442818842235130.260}

B={0.9527821496566200

00.8470

000.260}

C={00.5980.83681960067634

0.70000.592

0.953877473592230.44281884223513

0}

第四章习题4

3.求函数/(工)="6”的数值积分s=并请接受符号计算尝试

复算。

解：

(1)数值积分：

clear

symsx

x=0:0.1:pi;

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

>>trapz(x,y)

ans=

5.0954

⑵符号计算：

symsx

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

Warning:Explicitintegralcouldnotbefound.

>Inat58

Isym=

5.761450876

4.用quad求取1的数值积分，并保证积分的确定精度为

10，

解:Y=quad('exp(-abs(x)).*abs(sin(x))',-5*pi,1.7*pi,le-9)

Y=

1.779

5.求函数/(0=(sin5r)2e)06/2-1.5rcos2r+1.即+0.5|在区间［-5,51中的最小值

点。

解：

clear

symsx

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

yd=diff(y,x);

a=solve(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=

-.6164860955e-2

yd_a=

.le-30

y_a=

.899099

y_m_5=

1.88547

y_p_5=

16.2715

最小值点即为a点.

第五章习题5

1.已知椭圆的长、短轴〃=4*=2（见图p5-l）,用“小红点线”画椭圆

x=acost

y=Z?sinl

图p5-l

解：如图:

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

P5-3所示的三城市上半年每月生产总值的累计直方图。（提示：运用

指令bar）

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

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

A170120180