当前位置：文档库 › 数学实验答案-1

数学实验答案-1

1.（1）

[1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2.

A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1]

X=(B+2*A)/2

A=[-4 -2 0 2 4;-3 -1 1 3 5]

abs(A)>3

A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5]

det(A),eig(A),rank(A),inv(A)

求计算机高手用matlab解决。

>> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5]

求|A|

>> abs(A)

ans =

2 3 2 4

1 2 3 2

3 2 3 4

0 4 2 5

求r（A）

>> rank(A)

ans =

求A-1

>> A-1

ans =

-3 2 1 3

0 -3 2 1

2 1 2 3

-1 3 -3 4

求特征值、特征向量

>> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵V

V =

0.7335 0.7335 -0.3804 - 0.0312i -0.3804 + 0.0312i -0.0024 + 0.5329i -0.0024 - 0.5329i -0.3907 - 0.0001i -0.3907 + 0.0001i -0.3166 - 0.0283i -0.3166 + 0.0283i -0.8280 -0.8280 -0.0556 - 0.2718i -0.0556 + 0.2718i 0.0301 - 0.1235i 0.0301 + 0.1235i

D =

-3.1766 + 0.6201i 0 0 0 0 -3.1766 - 0.6201i 0 0 0 0 5.1766 + 0.7101i 0 0 0 0 5.1766 - 0.7101i

将A的第2行与第3列联成一行赋给b

>> b=[A(2,:),A(:,3)']

b =

1 -

3 2 2 3 3 -2

1．

a=round(unifrnd(1,100))

i=7;

while i>=0

i=i-1;

b=input('请输入一个介于0到100的数字： '); if b==a

disp('You won!');

break;

else if b>a

disp('High');

else if b

disp('Low');

end

结果

a =

请输入一个介于0到100的数字： 50

Low

请输入一个介于0到100的数字： 75

Low

请输入一个介于0到100的数字： 85

请输入一个介于0到100的数字： 82

You won!

2．

clear all；clc;

n=input('请输入数字 n=');

n1=floor(n/100); %取出百位数字n1

n2=mod(floor(n/10),10); %取出十位数字n2

n3=mod(n,10) ; %取出个位数字n3

if n1^3+n2^3+n3^3==n

fprintf('%d是“水仙花数”', n) % 注意输出格式前须有%符号else

fprintf('%d不是“水仙花”', n) % 注意输出格式前须有%符号end

结果

请输入数字 n=234

234不是“水仙花数”>>

3．

price=input('请输入商品价格');

switch fix(price/100)

case {0,1} %价格小于200

rate=0;

case {2,3,4} %价格大于等于200但小于500

rate=3/100;

case num2cell(5:9) %价格大于等于500但小于1000 rate=5/100;

case num2cell(10:24) %价格大于等于1000但小于2500 rate=8/100;

case num2cell(25:49) %价格大于等于2500但小于5000 rate=10/100;

otherwise %价格大于等于5000

rate=14/100;

end

price=price*(1-rate) %输出商品实际销售价格结果

请输入商品价格250

price =

242.5000

Function f=myfun(x)

x=input;

s=pi*x*x

l=pi*x^2

4、

Function y=circle(r)

s=pi*x*x

l=pi*x^2

syms r

s=pi*r*r

l=2*pi*r

5. function fibonacci(n,m) f(1)=1;f(2)=1;

for i=3:max(n,m)

f(i)=f(i-1)+f(i-2); end

fprintf('第%d项',m)

x=f(m)

fprintf('前%d项',n)

s=f(1:n)

COMMAND WINDOW输入：fibonacci(20,50)

1.绘制])4,0[)(3sin(3

π∈=x x e y x 的图像，要求用蓝色的星号画图；并且画出器官包络线3x e y ±=的图像，用红色的点划线画图。

2.用fplot 和ezplot 命令绘出函数)21sin(32t e

y t +=-在区间]10,0[上的图像。 3.在同一图像窗口画三个子图要求使用指令gtext,axis,legend,title,xlabel,和ylabel ：（3）]8,1[,sin 1

∈=x x e y x

x=0:pi/25:4*pi;

y1=exp(x/3).*sin(3*x);y2=exp(x/3);y3=-exp(x/3);

plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.')

t=1:0.1:10

y=exp(-2*t/3).*sin(1+2*t);

plot(t,y);

figure

fplot('exp(-2*t/3).*sin(1+2*t)',[1,10])

ezplot('exp(-2*t/3).*sin(1+2*t)',[1,10])

x=1:1/50:8;

y=exp(1./x).*sin(x);

subplot(1,3,3),plot(x,y,'b-');

legend('y=exp(1/x)sinx'),grid on ,

title('y=exp(1/x)sinx'),

xlabel('x 轴'),ylabel('y 轴')

gtext('真棒'),axis([x1 x8 y1 y10])

x=0:pi/50:2*pi;

y1=sin(x);y2=cos(x);y3=sin(2*x);

plot(x,y1,'k*--',x,y2,'rs-',x,y3,'bo--'),grid on

title('曲线y1=sin(x),y2=cos(x)与y3=sin(2*x)')

xlabel('x 轴'),ylabel('y 轴')

gtext('y1=sin(x)'),gtext('y2=cos(x)'),gtext('y3=sin(2*x)')

legend('y1=sin(x)','y2=cos(x)','y3=sin(2*x)')

5.绘制圆锥螺线的图像并加各种标注，圆锥螺线的参数方程为：

????

?????=≤≤==t

z t t t t y t

t x 2)

200(,6

cos 6cos π

6.在同一图形窗口画半径为1的球面，柱面122=+y x 以及极坐标]2,0[,4sin 21πρ∈=

t t 7.用mesh 与surf 命令绘制三维曲面223y x z +=的图像，并使用不同的着色效果及光照效

果。

8.绘制由函数14

1692

22=++z y x 形成的立体图，并通过改变观测点获得该图形在各个坐标平面上的投影。

9.画三维曲面)2,2(52

2≤≤---=y x y x z 与平面3=z 的交线。

t=1:pi/50:20*pi;

x=t.*cos(pi/6.*t);

y=t.*sin(pi/6.*t);

z=2*t;

plot3(x,y,z);

grid on ;

title('圆锥螺线')

xlabel('x 轴');ylabel('y 轴');zlabel('z 轴');

axis square

6. v=[-2 2 -2 2 -2 2];

subplot(1,3,1),sphere

title('以半径为1的球面'),xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴'),axis(v)

subplot(1,3,2),cylinder

title('柱面'),xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')

subplot(1,3,3),t=0:pi/100:2*pi;

polar(t,(1/2)*sin(4*t))

title('p=(1/2)*sin(4t)')

[X,Y]=meshgrid(-8:0.5:8);

Z=X.^2+3*Y.^2;

subplot(1,2,1);mesh(X,Y,Z);shading interp

subplot(1,2,2);surf(X,Y,Z);shading flat

[xx,yy,zz]=sphere(40);

x=xx*3;

y=yy*4;

z=zz*2;

surf(x,y,z)

axis equal

《数学实验》试题答案

北京交通大学海滨学院考试试题课程名称：数学实验2010-2011第一学期出题教师：数学组适用专业: 09机械, 物流, 土木, 自动化班级：学号：姓名：选做题目序号： 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。解: 由题意每月的成兔与幼兔的数量如下表所示： 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序： x=zeros(1,24); x(1)=1;x(2)=1； for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用

已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解：根据求积分的过程，我们先对区间[0,1]进行n 等分，然后针对函数x e 取和，取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ，其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点，则当10n =时，1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算，当10n =0时，100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?，当10n =000时，10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图，Matlab 命令如下： x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10)；%n=10 symsum(exp(k/100)/100,k,1,100)；%n=100 symsum(exp(k/10000)/10000,k,1,10000)；%n=10000

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“”练习题要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;

>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =

0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1

6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5

8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =

数学实验答案-1

1.（1） [1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2. A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1] X=(B+2*A)/2 3. A=[-4 -2 0 2 4;-3 -1 1 3 5] abs(A)>3 % 4. A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5] det(A),eig(A),rank(A),inv(A) 求计算机高手用matlab解决。 >> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5] 求|A| >> abs(A) ans = （ 2 3 2 4 1 2 3 2 3 2 3 4 0 4 2 5 求r（A） >> rank(A) ans =

4 求A-1 《 >> A-1 ans = -3 2 1 3 0 -3 2 1 2 1 2 3 -1 3 -3 4 求特征值、特征向量 >> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵V , V = - + + - - + - + - + - + D = { + 0 0 0 0 - 0 0 0 0 + 0 0 0 0 - 将A的第2行与第3列联成一行赋给b >> b=[A(2,:),A(:,3)'] b = 《 1 - 2 3 2 2 3 3 -2

1． a=round(unifrnd(1,100)) i=7; while i>=0 i=i-1; b=input('请输入一个介于0到100的数字：'); if b==a ￥ disp('You won!'); break; else if b>a disp('High'); else if b