Matlab与数学实验(第二版)(张志刚 刘丽梅 版) 习题答案

Matlab与数学实验（第二版）（张志刚刘丽梅版）习题答案

（1,3,4,5章）

第一章

d1zxt1

用format的不同格式显示2*Pi，并分析格式之间的异同。

a=2*pi ;

disp('***(1) 5位定点表示2*pi：')

format short , a % 5位定点表

disp('***(2) 15位定点表示2*pi：')

format long , a % 15位定点表

disp('***(3) 5位浮点表示2*pi：')

format short e , a % 5位浮点表示

disp('***(4) 15位浮点表示2*pi：')

format long e , a % 15位浮点表示

disp('***(5) 系统选择5位定点和5位浮点中更好的表示2*pi：')

format short g , a % 系统选择5位定点和5位浮点中更好的表示

disp('***(6) 系统选择15位定点和15位浮点中更好的表示2*pi：')

format long g , a % 系统选择15位定点和15位浮点中更好的表

disp('***(7) 近似的有理数的表示2*pi：')

format rat , a % 近似的有理数的表

disp('***(8) 十六进制的表示：')

format hex , a % 十六进制的表

disp('***(9) 用圆角分(美制)定点表示2*pi：')

format bank , a % 用圆角分(美制)定点表示

d1zxt2利用公式求Pi的值。

sum=0 ;

n=21;

for i = 1:4:n % 循环条件

sum= sum+(1/i) ; % 循环体

end

diff=0 ;

for j = 3:4:(n-2) % 循环条件

diff= diff+(1/j) ; % 循环体

end

pai=4*(sum-diff)

d1zxt3 编程计算1!+3!+...+25！的阶乘。

% 方法1：利用“while循环”来计算1!+3!+...+25！的值。

sum=0;

i=1;

while i<26

prd=1;

j=1;

while j<=i % 求数i的阶乘

prd=prd*j;

j=j+1;

end

sum=sum+prd;

i=i+2;

end

disp('1!+3!+...+25！的和为：')

sum

% 方法2：利用“阶乘函数”来计算1!+3!+...+25！的值。

sum=0;

for i=1:2:25

sum=sum+factorial(i); % 求阶乘函数factorial(i) end

disp('1!+3!+...+25！的和为：')

sum

d1zxt4 编M文件计算n的阶乘。

n = 6 % 可以任意输入n的值？

prd=1;

j=1 ;

while j<=n % 求数n的阶乘

prd=prd*j;

j=j+1;

end

disp('n的阶乘n！为：')

prd

d1zxt4 JCHS 编阶乘计算n的阶乘。

function d1zxt4JCHS(n)

n

disp('编阶乘计算n的阶乘为：')

prd=1;

j=1 ;

while j<=n % 求数n的阶乘

prd=prd*j;

j=j+1;

end

d1zxt5_1计算表达式的值。

exp(123)+1234^34*log2(3)/cos(21/180*pi)

value=exp(123)+1234^34*log2(3)/cos(21/180*pi)

d1zxt5_2 计算函数的值。

function [y]=d1zxt5_2(x)

x

disp('此时对应x 的y 值为：')

y=tan(-x^2)*acos(x);

第二章

d3zl1

x=0:0.1:4;

y1=x.^(-1.0 )

y2=x.^(-0.5)

y3=x.^(+0.5)

y4=x.^(+2.0 )

y5=x.^(+3.0)

y6=x.^(+4.0)

plot(x,y1,x,y2,x,y3,x,y4,x,y5,x,y6) % 洋红，差号

Legend ('y1=x.^(-1.0 )', 'y2=x.^(-0.5) ', 'y3=x.^(+0.5) ', 'y4=x.^(+2.0 )', 'y5=x.^(+3.0)', 'y6=x.^(+4.0)' ) %图形标注

d3zd2t

x=-10:0.1:10;

y1=(1/3).^x

y2=(1/2).^x

y3=(2).^x

y4=(3).^x

plot(x,y1,x,y2,x,y3,x,y4) % 洋红，差号

Legend ('y1=(1/3)^x', 'y2=(1/2)^x ', 'y3=(2)^x ', 'yy4=(3)^x' ) %图形标注