matlab上机作业
- 格式:doc
- 大小:133.84 KB
- 文档页数:8
第二次 上机作业
1、 求下列矩阵的主对角线元素、上三角阵、下三角阵、秩、范数、条件数和迹。
(1)⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡--=901511250324153211A (2)⎥⎦⎤⎢⎣⎡-=2149.824343.0B 1.
A=[1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9]
D=diag(A)
C=triu(A)
B=tril(A)
E=rank(A)
F=trace(A)
a1=norm(A,1)
a2=norm(A,inf)
a3=norm(A,inf)
c1=cond(A)
c1=cond(A,1)
c2=cond(A,2)
c3=cond(A,inf)
2.
B=[0.43,43,2;-8.9,4,21]
D=diag(B)
C=triu(B)
B=tril(B)
E=rank(B)
F=trace(B)
a1=norm(B,1)
a2=norm(B,inf)
a3=norm(B,inf)
c1=cond(B)
c1=cond(B,1)
c2=cond(B,2)
c3=cond(B,inf)
2、 求矩阵A 的特征值和相应的特征向量。
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=225.05.025.0115.011A
A=[1,1,0.5;1,1,0.25;0.5,0.25,2]
[V ,D]=eig(A)
3、 下面是一个线性方程组:
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.06/15/14/15/14/13/14/13/12/1321x x x (1) 求方程的解。
(2) 将方程右边向量元素3b 改为0.53,再求解,并比较3b 的变化和解的相
对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]
B=[0.95,0.67,0.52]
X=inv(A)*b
c1=cond(A,1)
c2=cond(A,2)
c3=cond(A,inf)
4、 利用Matlab 提供的randn 函数生成符合正态分布的10×5随机矩阵A,进行
如下操作:
(1)A 各列元素的均值和标准方差
(2)A 的最大元素和最小元素
(3)求A 每行元素的和以及全部元素之和
(4)分别对A 的每列元素按升序、每行按降序排列
X=randn(10,5)
M=mean(X)
D=std(X)
m=max(X)
n=min(X)
P=sum(X,2)
sum(p)
Ans
5、按要求对指定函数进行插值
(1)按下表用三次样条方法插值计算0~900范围内整数点的正弦值和0~
750范围内整数点的正切值。
表1 特殊角的正弦与正切值表
α(度)0 15 30 45 60 75 90 α
sin0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000 α
tan0 0.2679 0.5774 1.0000 1.7320 3.7320
(2)用三次多项式方法插值计算1~100之间整数的平方根
表2 1~100内特殊值的平方根
N 1 4 9 16 25 36 49 64 81 100 N 1 2 3 4 5 6 7 8 9 10 1.
X=0:15*pi/180:pi/2;
sin(X)
tan(X)
Format long
Interp1(x,sin(x),’spline’)
Columns 1 through 4
Columns 5 through 7
2.
N=[1 4 9 16 25 36 49 64 81 100]
n=sqrt(N)
format long
i nterp1(N,n,’cubic’)
6、将10个学生5门功课的成绩存入矩阵P中,进行如下处理:
(1)分别求每门功课的最高分、最低分及相应的学号。
(2)分别求每门功课的平均值和标准方差。
(3) 5门课总分的最高分、最低分及相应的学号。
(4)将5门课总分按从大到小的顺序存入abc中,相应的学生序号存入mmw。
提示:为了避免输入学生成绩的麻烦,可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。
X=45+(95-45)*rand(5,10)
M1=max(x,[],2)
[M1,l1]=max(x,[],2)
M2=min(x,[],2)
[M2,l2]=min(x,[],2)
X1=mean(x)
X2=std(x)
N1=sun(X,2)
N2=sun(X,1)
[n2,l3]=max(x,[],1)
7、某气象站观测某日6:00~8:00之间每隔2h内温度(0C)如实验表1所示
实验表1 室内外温度观测结果(0C)
试用三次样条插值分别求出该日室内外6:30~17:30之间每隔2h各个点的近似温度(0C)。
h=6:2:18;
t1=[18,20,22,25,30,28,24]
t2=[15,19,24,28,34,32,30]
format long
interp1(h,t1,’splime’)
interp1(h,t2,’splime’)
8、已知一组实验数据如下表所示: