对的对数。（2）与上述亲密数对对应的所有亲密素数之和。
f实验五
二、实验内容
4设
f
xx21201
1x34001
，编写一个
MATLAB
函数文件fxm，使得调用fx时x可用矩阵代入，得出的fx为
同阶矩阵。
f5已知
y
f40
f30f20
1当f
10l
25时，求y的值。
2当时f
122334
1，求y的值。
1
2
实验六
1
设
y

05
3si
1x
x
2

cos
x
，在x02区间取101点，绘
制函数的曲线。
f4．绘制极坐标曲线
asi
b
并分析参数a、
b、
对曲线形状的影响。
f以上五张截图分别是a1b1
1、2、3、4、7时的情况，不
难发现，当
为奇数时画出的图有奇数个环，而当
为偶数时
画出的图有该偶数的两倍个环。参数a控制极坐标的半径，参
数b可对图进行角度旋转。
6绘制曲面图形，并进行插值着色处理
xcosScost

ycos
S
si

t
zsi
S
0s0t3
2
2
实验七2利用曲面对象绘制曲面
fvxt10e001xsi
2000t02x，先利用默认属性绘制曲线，然后通过图形句柄操作来改变曲线的颜色、线型和线宽，并利用文字对象给曲线添加文字标注。
实验八1利用MATLAB提供的ra
d函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：（1）均值和标准方差。（2）最大元素和最小元素。（3）大于05的随机数个数占总数的百分比。
f1
2
（3）
2某气象观测站测得某日6001800之间每隔2h的室内外温度（℃）如实验表1所示。实验表1室内外温度观测结果（℃）
时间h681012141618室内温度180200220250300280240
t1室外温度150190240280340320300
t2
试用三次样条插值分别求出该日室内外6301730之间每隔2h各点的近似温度（℃）。
f5有3个多项式，P1xx42x34x25，P2xx2，
时进行下列操作：
（1）求PxP1xP2xP3x。
（2）求Px的根。
（3）当x取矩阵A的每一元素时，求Px的值。其
中：
11214
A0752
35

0525
（4）当以矩阵A为自变量时，求Px的值。其中A的
值与第（3）题相同。
1
2
f2
3
实验九
1求函数在指定点的数值导数。
xx2x3fx12xx2，x123
026x
2用数值方法求定积分。
（1）
2
I10
cost24si
2t21dt的近似值。
3分别用3种不同的数值方法解线性方程组。
f6x5y2z5z4
9xy4zu13

3x

4y

2z

2u
1
3x9y2u11
直接解法：
LU分解：
通解法：
4求非齐次线性方程组的通解。
32xx1157xx2223xx332xx4464

9x1

4x2

x3

7x4

2
f5求代数方程的数值解。
（2）
在给定的初值x01，
y0

1
，
z
0

1
下，求方程组的数
值解。
f
r