中南大学matlab课后习题

ans = 01 00 01
e&d
ans = 11 01 11
e|d
ans = 11 11 11
~e|~d
ans = 00 10 00
4,产生 5 阶希尔伯特矩阵 H 和 5 阶帕斯卡矩阵 P，且求其行列式的值 Hh 和 Hp 以及他们的条件数 Th 和 Tp，判断哪个矩阵性能更好，为什么 format rat H=hilb(5)
V= 1
-1227/1438 290/593
D= -3595/142 0 0
361/1013 -1
551/788
0 -3755/357
0
493/1574 1
331/713
0 0 4697/279
1/5 1/6 1/7 1/8 1/9
1 5 15 35 70
第三章
1． 从键盘输入一个 3 位整数，将它反向输出。如输入 639，输出为 936.
plot(x,y)
200
150
100
50
0
-50
-100
-150
-200 -10 -8 -6 -4 -2 0 2 4 6 8 10
(2) x=-10::10; ezplot('x^2+2*y^2-64',[-8,8]); grid on;
x2+2 y2-64 = 0 8
6
4
2
y
0
-2
-4
-6
-8
-8
a=
a=
y=
(2)当 f(n)=1*2+2*3+3*4+ +n*(n+1)时，y 的值是多少。 function b=n1b(n) b=0; for i=1:n b=b+i*(i+1); end y=n1b(40)/(n1b(30)+n1b(20)) y=
第四章 1.绘制下列曲线。 （1）y=x-
x=-10::10; y=x-x.^3/6;
clear
m=input('请输入一个三位数：');
m1=fix(m/100);
m2=rem(fix(m/10),10);
m3=rem(m,10);
m=m1+m2*10+m3*100;
disp(m)；
请输入一个三位数：453
354
2．输入一个百分制成绩，要求输出成绩等级 A、B、C、D、E。其中 90~100 分为 A，80~89 分为 B，
a^2-b+eye(3)
ans = -18 -217 17 22 533 109 21 867 526
（2）a*b
ans = 14 14 16 -10 51 21 125 328 180
a.*b
ans = -8 15 4
0 35 24 -9 122 0 b*a
ans = -11 0 -15 7 228 53 3 -1 28
（3）a/b
ans =
b\a ans =
（4）[a,b]
ans = -1 5 -4 8 0782 3 61 7 -3
[a([1,3],:);b^2]
3 -1 53 20
ans = -1 5 -4 3 61 7 73 37 1 17 37 13 -20 1 9
3.已知 a，完成下列操作。 a=[23,10,,0;41,-45,65,5;32,5,0,32;6,,54,]; （1） 输出 a 在[10,25]范围内的全部元素。 k=find(a>10&a<25)
120 150 180
90 8 60
6
4
30
2
0
210
330
240
300
270
4.绘制下列三维图形。
（1）
,0
.
t=0:pi/50:2*pi;
x=exp(-t./20).*cos(t); y=exp(-t./20).*sin(t); z=t; plot3(x,y,z); grid on;
ห้องสมุดไป่ตู้
8
6
4
2
0 1
switch fix(c/10)
%用 switch 语句实现
case{9,10}
m='A';
case{8}
m='B';
case{7}
m='C';
case{6}
m='D';
casenum2cell(0:5)
m='E';
otherwise
m='输入出错';
end
disp(m);
请输入一个百分制分数：98
a(k)
k= 1
ans = 23
(2)取出 a 前 3 行构成矩阵 b，前两列构成矩阵 c，右下角 3*2 子矩阵构成矩阵 d，b 与 c 的乘积构 成矩阵 e。 b=a(1:3,:)
b= 0
0 c=a(:,1:2)
c=
d=a(2:4,3:4)
d=
0
e=b*c
e= +003 *
(3)分别求表达式额 e<d,e&d,e|d 和~e|~d 的值。 e<d
70~79 分为 C,60~69 分为 D，60 分以下为 E。要求：
（1）分别用 if 语句和 switch 语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
c=input('请输入一个百分制分数：');
if c>=90&c<=100
%用 if 语句实现
disp('A');
第二章 1·求下列表达式的值。 （1）w=sqrt(2)*(1+*10^-6)
w=
（2）a=;b=5;c=; x=(2*pi*a+(c+b)/(pi+a*b*c)-exp(2))/(tan(b+c)+a); x
x=
（3）a=;b=; y=2*pi*a^2*[(1-pi/4)*b-4)*a]; y
y=
Hp = 1
Th=cond(H)
Th = 476607
Tp=cond(P)
Tp = 178868/21
P 矩阵性能更好，因为 Tp 更接近 1. 5.已知 A，求 A 的特征值及特征向量，并分析其数学意义。 A=[-29,6,18;20,5,12;-8,8,5]
A=
-29
6
18
20
5
12
-8
8
5
[V,D]=eig(A,'nobalance')
第五章 1， 利用 matlab 提供的 randn 函数生成符合正态分布的 10*5 随机矩阵 A，进行以下操作： （1） A 各列元素的均值和标准方差。 （2） A 的最大元素和最小元素。 （3） 求 A 每行元素的和以及全部元素的和。 （4） 分别对 A 的每列元素按升序，每行元素按降序排列。
-6
-4
-2
0
2
4
6
8
x
2,设
, -π≤t≤π，
在同一图形窗口采用子图的形式绘制条形图，阶梯图，杆图和对数坐标图等不同图形，并对 不同图形加标注说明。 t=-pi:pi/10:pi; y=1./(1+exp(-t));
subplot(2,2,1); bar(t,y); title('条形图(t,y)'); axis([-pi,pi,0,1]); subplot(2,2,2); stairs(t,y,'b'); title('阶梯图(t,y)'); axis([-pi,pi,0,1]); subplot(2,2,3); stem(t,y,'k'); title('杆图(t,y)'); axis([-pi,pi,0,1]); subplot(2,2,4); loglog(t,y,'y'); title('对数坐标图(t,y)');
0.5
0 -0.5
-1 -1
0 -0.5
1 0.5
（2）z=5,|x| 5,|y| [x,y]=meshgrid(-5:5); z=zeros(11)+5; mesh(x,y,z); shading interp;
5.要求应用插值着色处理。
5.播放一个直径不断变化的球体。 [x,y,z]=sphere(20); surf(x,y,z); axis off; shading interp; m=moviein(20); for i=1:20 axis([-i,i,-i,i,-i,i]) m(:,i)=getframe; end movie(m,4);
（4）t=[2,1-3i;5,]; z=1/2*exp(2*t)*log(t+sqrt(1+t^2)); z
z= +004 * ---
2，已知 a，b，求下列表达式的值。 a=[-1,5,-4;0,7,8;3,61,7];b=[8,3,-1;2,5,3;-3,2,0]; （1）a+6*b
ans = 47 23 -10 12 37 26 -15 73 7
Warning: Negative data ignored > In title at 42
In title at 23
条 形 图 (t,y)
1
1
阶 梯 图 (t,y)
0.5
0.5
0
-2
0
2
杆 图 (t,y) 1
0
-2
0
2