中南大学操作题集.docVIP

操作题集 操作题1 《Matlab与科学计算》P28 例2.1 要求计算水在温度为0℃，20℃，40℃，60℃，80℃时的黏度，已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2) 其中μ0为0℃水的黏度，值为1.758*10-3。 要求分别以长短格式显示计算结果。 muw0=1.785e-3; a=0.03368; b=0.000221; t=0:20:80; muw=muw0./(1+a*t+b*t.^2) muw = 0.0018 0.0010 0.0007 0.0005 0.0003 操作题2 《Matlab与科学计算》P31 例2.2 数值数组【1：5】和字符串转换示例。 a=[1:5]; b=num2str(a); a*2 ans = 2 4 6 8 10 b*2 ans = Columns 1 through 9 98 64 64 100 64 64 102 64 64 Columns 10 through 13 104 64 64 106 str2num(b)*2 ans = 2 4 6 8 10 操作题3 《Matlab与科学计算》P44 例2.9 比较用左除和右除法分别求解恰定方程的解。 rand(seed,12); a=rand(100)+1.e8; x=ones(100,1); b=a*x; cond(a) ans = 5.0482e+011 tic;x1=b/a;t1=toc t1 = 0.0940 er1=norm(x-x1) er1 = 139.8328 re1=norm(a*x1-b)/norm(b) re1 = 4.3095e-009 tic;x1=a\b;t1=toc t1 = 0.0160 tic;x1=a\b;t1=toc t1 = 0 er2=norm(x-x1) er2 = 1.6120e-004 re1=norm(a*x1-b)/norm(b) re1 = 4.5617e-016 操作题4 《Matlab与科学计算》 P48 例2.14 计算矩阵magic(3)的指数，并比较不同函数的结果。 b=magic(3); expm(b) ans = 1.0e+006 * 1.0898 1.0896 1.0897 1.0896 1.0897 1.0897 1.0896 1.0897 1.0897 expm1(b) ans = 1.0e+003 * 2.9800 0.0017 0.4024 0.0191 0.1474 1.0956 0.0536 8.1021 0.0064 操作题5 《Matlab与科学计算》 P50 例2.18 计算下面矩阵[-149 -50 -154；537 180 546；-27 -9 -25]的特征值条件数。 a=[-149 -50 -154;537 180 546;-27 -9 -25] a = -149 -50 -154 537 180 546 -27 -9 -25 [V,D,s]=condeig(a) V = 0.3162 -0.4041 -0.1391 -0.9487 0.9091 0.9740 -0.0000 0.1010 -0.1789 D = 1.0000 0 0 0 2.0000 0 0 0 3.0000 s = 603.6390 395.2366 219.2920 操作题6《Matlab与科学计算》P62 例2.29 对矩阵a = 1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20 进行三角抽取使之分别变成 0 0 0 0 1 1 1 0 1 0 0 0