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

  • 2015-12-18
  • 674
  • 0
  • 0

tumblr_n0kqo3XUUW1smz89ro1_1280

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

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’ ) %图形标注)’, ‘y5=x.^(+3.0)’, ‘y6=x.^(+4.0)’ ) %图形标
d3zd3t
x=0.1:0.1:10;
y1=log(x)/log(1/3);
y2=log(x)/log(1/2);
y3=log(x)/log(2);
y4=log(x)/log(3);
plot(x,y1,x,y2,x,y3,x,y4) % 洋红, 差号
Legend (‘y1=log1/3(x)’, ‘y2=log1/2(x)’, ‘y3=log2(x) ‘, ‘y4=log3(x)’ ) %图形标注
d3zd4t
x=-4:0.1:4;
subplot(1,3,1)
plot(x,x.^2.*sin(x)) ,grid on,title(‘用plot绘图y=x^2*sin(x)’) % 用plot绘图,网格打开,加上图名
subplot(1,3,2)
fplot(‘x.^2.*sin(x)’,[-4,4]) ,grid on,title(‘用fplot绘图y=x^2*sin(x)’ ) % 用fplot绘图,网格打开,加上图名 subplot(1,3,3)
ezplot(‘x.^2.*sin(x)’,[-4,4]) ,grid on,title(‘用ezplot绘图y=x^2*sin(x)’ ) % 用ezplot绘图,网格打开,加上图名
d3zd5t
% 用plot绘图
subplot(1,2,1)
x=0:0.2:5;
y= 3.*cos(x).*exp(sin(x));
plot(x,y,’yd–‘),grid on,title(‘用plot绘图y=3*cos(x)*exp(sin(x)) 注意:数据点!’)
% 用fplot绘图
subplot(1,2,2)
fplot(‘3*cos(x)*exp(sin(x))’,[0,5],’yd–‘) ,grid on,title(‘用fplot绘图y=3*cos(x)*exp(sin(x)) 注意:数据点!’) % 用fplot绘图
d3zd6t
x1=-pi:0.1:pi;
y1=x1.*sin(x1);
plot(x1,y1, ‘mo-.’)
hold on
x2=pi:0.1:4*pi;
y2=x2.*tan(1./x2).*sin((x2).^3);
plot(x2,y2, ‘rx-‘)
x3=-pi:0.1:4*pi;
y3=(x3).^3;
plot(x3,y3, ‘y*–‘)
x4=1:0.1:8;
y4=tan(x4)+sin(x4);
plot(x4,y4, ‘gs:’)
legend(‘x*sinx’,’x*tan(1/x)*sinx^3)’,’x^3′,’tanx+sinx’)
d3zd6t
subplot(2,2,1)
x1=-pi:0.1:pi;
y1=x1.*sin(x1);
plot(x1,y1, ‘mo-.’)
xlabel(‘横坐标:x,[-pi,pi]’)
ylabel(‘纵坐标:x*sinx’)
title(‘图1 函数y=x*sinx的曲线’)
subplot(2,2,2)
x2=pi:0.1:4*pi;
y2=x2.*tan(1./x2).*sin((x2).^3);
plot(x2,y2, ‘rx-‘)
xlabel(‘横坐标:x,[pi,4pi]’)
ylabel(‘纵坐标:x*tan(1/x)*sinx^3)’)
title(‘图2 函数y=x*tan(1/x)*sinx^3)的曲线’)
subplot(2,2,3)
x3=-pi:0.1:4*pi;
y3=(x3).^3;
plot(x3,y3, ‘y*–‘)
xlabel(‘横坐标:x,[-pi,4pi]’)
ylabel(‘纵坐标:x^3’)
title(‘图3 函数y=x^3的曲线’)
subplot(2,2,4)
x4=1:0.1:8;
y4=tan(x4)+sin(x4);
plot(x4,y4, ‘gs:’)
xlabel(‘横坐标:x,[1,8]’)
ylabel(‘纵坐标:tanx+sinx’)
title(‘图4 函数y=tanx+sinx的曲线’)
d3zd7t 绘制摆线
r=10 ;
t=0:0.1:4*pi;
x=r*(t-sin(t));
y=r*(1-cos(t));
plot(x,y)
第四章
d4zd1t 同p88例2
s=sym([‘(1+2*x-y)^2’]);
p=expand(s)
p1=collect(p, ‘y’)
p1=collect(p, ‘y’)
d4zd2t 参考p83例5
a=[4 -19 -4 1];
b=[-3 1];
[q,r]=deconv(a,b)
q=poly2str(q, ‘x’)
r=poly2str(r, ‘x’)
d4zd3t
p=sym(‘120-46*x-19*x^2+4*x^3+x^4′)
q=sym(’20*x+9*x^2+x^3’)
zdgy=gcd(p,q)
zxgb=lcm(p,q)
d4zd4t
%参考p99例题1.
p=[1 0 0 0 -1 2 -3];
r=roots(p)
%参考p99例题2.
s1=sym([‘x^6-x^2+2*x-3’]);
solve(s1, ‘x’)
d4zd5t
%画出方程对应的函数在[0,2]内的图形:
[x,y]=fplot(‘sin(x)-log(x+0.1)’,[0,2]);
plot(x,y)
grid on
%方程在1附近的根:
x=fzero(‘sin(x)-log(x+0.1)’,1)
%方程的根:
s1=sym([‘sin(x)-log(x+0.1)’]);
solve(s1, ‘x’)
d4zd7t
1
% syms x;
syms p;
p=input(‘输入第1个多项式: ‘,’s’)
syms q;
q=input(‘输入第2个多项式: ‘,’s’)
% p=sym(‘120-46*x-19*x^2+4*x^3+x^4′)
% q=sym(’20*x+9*x^2+x^3’)
syms zdgy;
syms zxgb;
zdgy=gcd(p,q)
zxgb=lcm(p,q)
2
p=sym(‘120-46*x-19*x^2+4*x^3+x^4′)
q=sym(’20*x+9*x^2+x^3′)
zdgy=gcd(p,q)
zxgb=lcm(p,q)
第五章
d5zd1t
a=rand(8,8);
a’ %求矩阵a的转置
a^(-1) %求矩阵a的逆
rank(a) %求矩阵a的秩
d5zd2t
A=[2 7 9 6 1 0;3 5 0 7 8 7;5 5 1 9 2 1;1 4 7 4 2 0;6 0 5 3 2 0]
B=[1 7 0 1 3 6 0;3 8 4 5 3 4 5;2 6 5 2 1 8 4;2 5 5 4 4 5 2;5 1 9 6 8 3 1;7 1 6 0 7 2 9]
rank(A*B) %求矩阵(A*B)的秩
rank(A) %求矩阵(A)的秩
rank(B) %求矩阵(A)的秩
rank(B*A) %求矩阵(B*A)的秩(无法计算B*A)
d5zd5t 参考p114-116 例题1、2、3
A=[2 2 -2;2 5 -4;-2 -4 5]
p=poly(A) % 求矩阵A的特征多项式
poly2str(p, ‘x’) % 以惯用形式表示求出的A的特征多项式p
[p,d]=eig(A)
%返回矩阵A的: 特征向量 矩阵p,
%返回矩阵A的: 特征值 矩阵d
%满足Av=vd
format short %5位定点表示
p*p’ %结果表明:p为“所求”得到的【正交的相似变换矩阵】
p’*A*p %能够将矩阵A对角化,即满足p’*A*p =以A的3个
%特征值1、1、10为对角元素的对角矩阵
d5zd7t 参考p118 例题5
A=[2 0 0;0 3 2;0 2 3]
[p,d]=eig(A)
%返回矩阵A的: 特征向量 矩阵p,
%返回矩阵A的: 特征值 矩阵d
%满足Av=vd
p*p’ %结果表明:p为“所求”得到的【正交的相似变换矩阵】
p’*A*p %能够将矩阵A对角化,即满足p’*A*p =以A的3个
%特征值1、2、5为对角元素的对角矩阵
syms x1 x2 x3
X=[x1; x2; x3];Y=p*X
d5zd8t 1参考p108 例题5
A=[2 5;1 3]
B=[4 -6;2 1]
X=B/A
d5zd8t 2 参考p108 例题5
A=[2 5;1 3]
B=[4 -6;2 1]
X=B/A
d5zd11t 1 参考p111 例题10
A=[2 1 -1 1;4 2 -2 1;2 1 -1 -1]
b=[1;2;1]
B=([A,b])
rref(B)
d5zd11t 2 参考p111 例题9
A=[2 1 -1 1;3 -2 1 -3;1 4 -3 5]
b=[1;1;-2]
B=([A,b])
rref(B)
d5zd16t
A=[2 -2 0;-2 1 -2;0 -2 0]
B=[2 1 1;1 2 1;1 1 2]
C=A+B
AT=eig(A)
BT=eig(B)
CT=eig(C)
X=AT(1,1)+AT(2,1)+AT(3,1)+BT(1,1)+BT(2,1)+BT(3,1)
Y=CT(1,1)+CT(2,1)+CT(3,1)

 

评论

还没有任何评论,你来说两句吧

发表评论

心若止水,方能惊起万丈澜。