MatLab数学矩阵练习册答案.docx

MatLab数学矩阵练习册答案.docx

《MatLab数学矩阵练习册答案.docx》由会员分享，可在线阅读，更多相关《MatLab数学矩阵练习册答案.docx（12页珍藏版）》请在冰豆网上搜索。

MatLab数学矩阵练习册答案

MatLab练习册

（1）参考答案

1,（5*2+1.3-0.8）*10^2/25答案：

42

2，y=

答案：

2.2618

sin（10*pi）*exp（-0.3+4^2）+log（23）/log（4）

3，z=10*sin（pi/3）*...

cos（pi/3）答案：

4.3301

4，x=sin（223/3），y=x^2，z=y*10；

求x+2y-5z答案：

-37.6025

5，

计算

答案：

-2.1379+1.6552i

z1=2+7*i

z2=2*i

z3=5*exp（2*pi*i）

z=（z1*z2）/（z2+z3）

6,建立起始值=3,增量值=5.5,终止值=44的一维数组x答案：

（填写程序语句）

x=（3:

5.5:

44）

7,建立等差一维数组x：

首项为0,末项为

项数为15答案：

（填写程序语句）

x=linspace（0,pi,15）

8,计算行列式的值

答案：

-8

A=[3,2,3;4,2,6;7,8,1]

det（A）

9,矩阵

矩阵

;分别求出AxB及A与B中对应元素之间的乘积.

答案：

A=[3,2,3;4,2,6;7,8,1]

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

A*B

161616

262626

262626

A.*B

323

8412

21243

10,写出计算Sin（30o）的程序语句.

答案：

（填写程序语句）sin（pi*30/180）或sin（pi/6）

11,

计算f

（2）,f

（1）f

（2）+f2（3）答案：

x=2

f2=（x^3-2*x^2+x-6.3）/（x^2+0.05*x-3.14）

f2=-4.4792

x=1

f1=（x^3-2*x^2+x-6.3）/（x^2+0.05*x-3.14）

x=3

f3=（x^3-2*x^2+x-6.3）/（x^2+0.05*x-3.14）

f1*f2+f3^2

-12.6023

MatLab练习册

（2）参考答案

1,在[0,4pi]画sin（x）,cos（x）（在同一个图象中）;其中cos（x）图象用红色小圆圈画.并在函数图上标注“y=sin（x）”,“y=cos（x）”,x轴,y轴,标题为“正弦余弦函数图象”.

答案：

x=linspace（0,4*pi,100）;

y=sin（x）;

plot（x,y）;

gtext（'y=sin（x）'）;%图形注解,注意要用鼠标定位

holdon;

y=cos（x）;

plot（x,y,'ro'）;

gtext（'y=cos（x）'）;%图形注解

xlabel（'x轴'）;%x轴注解

ylabel（'y轴'）;%y轴注解

title（'正弦余弦函数图象'）;%图形标题

2,画出

所表示的三维曲面。

的取值范围是

答案：

x=-8:

0.5:

8;%定义自变量x的一维刻度向量

y=x;%定义自变量y的一维刻度向量

[X,Y]=meshgrid（x,y）;

R=sqrt（X.^2+Y.^2）;%计算中间变量

Z=sin（R）./R;%计算相应的函数值

mesh（X,Y,Z）;%或用surf（X,Y,Z）;画图

3,绘制椭圆

表面图,x在[-2,2],y在[-2,2]范围.

答案：

x=-2:

0.3:

2;

y=-2:

0.3:

2;

[X,Y]=meshgrid（x,y）;

Z=sqrt（4-X.^2/9-Y.^2/4）;

surf（X,Y,Z）;

4,绘制环面:

答案：

u=linspace（0,2*pi,30）;

v=u;

[U,V]=meshgrid（u,v）;

X=（1+cos（U））.*cos（V）;

Y=（1+cos（U））.*sin（V）;

Z=sin（U）;

mesh（X,Y,Z）;

surf（X,Y,Z）;

MatLab练习册（3）

学生信息：

系别：

班级：

姓名：

1,利用二分法编程求方程

在[0,3]内的根.

答案：

1（填写程序语句）

a=0;

b=3;

x=（a+b）/2;

whileabs（a-b）>10^（-6）

iff（a）*f（x）<0

b=x;

elseiff（b）*f（x）<0

a=x;

end

x=（a+b）/2;

end

x

%%%%%%%%%%%%%%%%%%%%%%%%%下面是定义的函数，存为f.m

functiony=f（x）

y=x.^3-3*x.^2-x+3;

2,利用二分法编程求方程tan（x）-x-0.5=0在[0,1.5]内的根.

答案：

0.975（填写程序语句）

a=0;

b=1.5;

x=（a+b）/2;

whileabs（a-b）>10^（-6）

iff2（a）*f2（x）<0

b=x;

elseiff2（b）*f2（x）<0

a=x;

end

x=（a+b）/2;

end

x

%%%%%%%%%%%%%%%%%%%%%%%%%下面是定义的函数，存为f2.m

functiony=f2（x）

y=tan（x）-x-0.5;

3,利用

公式求

的近似值,直到最后一项小于

答案：

3.14（填写程序语句）

pi=0;n=1.0;s=1;t=1;

whileabs（t）>10.^（-6）

pi=pi+t;

n=n+2;

s=-s;

t=s/n;

end

4.0*pi

4,列出所有水仙花数,水仙花数是一个三位数,其各位数字立方和等于该数本身.例如:

答案：

（填写程序语句）

fori=1:

9

forj=0:

9

fork=0:

9

if（100*i+10*j+k）==（i.^3+j.^3+k.^3）;

disp（100*i+10*j+k）;

end

end

end

end

5,方程

的根.

答案：

1.076（填写程序语句）

6,求积分

答案：

0.6559（填写程序语句）

MatLab练习册（5）参考答案

1,求下面的优化问题:

min-5x1+4x2+2x3

6x1-x2+x3<=8

x1+2x2+4x3<=10

3>=x1>=-1;

2>=x2>=0;

x3>=0;

答案：

（填写程序语句）

c=[-5,4,2,];

a=[6-11;124];

b=[810];

xlb=[-100];

xub=[32inf];

x0=[0,0,0];

x=lp（c,a,b,xlb,xub,x0）;

disp（['最优值为:

',num2str（c*x）]）

x=

1.3333

0

0

最优值为:

-6.6667

2,某车间生产A和B两种产品.为了生产A和B,所需的原料分别为2个和3个,而所需要的工时分别为4个和2个单位,现在可以应用的原料为100个单位,工时为120个单位,每生产一台A和B分别可获得利润6元和4元,应当安排生产A,B各多少台,才能获得最大的利润?

答案：

（填写程序语句）

c=[-6,-4];

a=[2,3;4,2];

b=[100,120];

xlb=[0,0];

xub=[inf,inf];

x=lp（c,a,b,xlb,xub,xo）

disp（['最优值为:

',num2str（-c*x）]）

x=

20.0000

20.0000

最优值为:

200

3,min

x1+x2=0;

1.5+x1x2-x1-x2<=0;

-x1x2-10<=10

答案：

（填写程序语句）

function[f,g]=fg（x）

f=exp（x

（1））*（4*x

（1）^2+2*x

（2）^2+4*x

（1）*x

（2）+2*x

（2）+1）;

g

（1）=x

（1）+x

（2）;

g

（2）=1.5+x

（1）*x

（2）-x

（1）-x

（2）;

g（3）=-x

（1）*x

（2）-20;

x0=[-1,1];

options（13）=1;

[x,options]=constr（'fg',x0,options）

options（8）

x=

-1.22471.2247

options（8）

ans=

1.8951

4,对以下数据分别作二次,三次多项式拟合,并画出图形.

x=1:

16;

y=[4,6.4,8,8.4,9.28,9.5,9.7,9.86,10,10.2,10.32,10.42,10.5,10.55,10.58,10.6];

答案：

（填写程序语句）

二次多项式拟合

x=1:

1:

16;

y=[4,6.4,8,8.4,9.28,9.5,9.7,9.86,10,10.2,10.32,10.42,10.5,10.55,10.58,10.6];

a=polyfit（x,y,2）

a=

-0.04451.07114.3252

ezplot（'-0.0445*x^2+1.0711*x+4.3252'）

三次多项式拟合

x=1:

1:

16;

y=[4,6.4,8,8.4,9.28,9.5,9.7,9.86,10,10.2,10.32,10.42,10.5,10.55,10.58,10.6];

a=polyfit（x,y,3）

a=

0.0060-0.19632.13462.5952

ezplot（'0.0060*x^3-0.1963*x^2+2.1346*x+2.5952'）