matlab数学实验习题全部答案(胡良剑).doc

PAGE\*MERGEFORMAT33

数学实验答案
%Page20,ex1
(5)等于[exp(1),exp(2);exp(3),exp(4)]
(7)3=1*3,8=2*4
(8)a为各列最小值，b为最小值所在的行号
(10)1>=4,false,2>=3,false,3>=2,ture,4>=1,ture
(11)答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)
(12)答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)
%Page20,ex2
(1)a,b,c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b
(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码
%Page20,ex3
>>r=2;p=0.5;n=12;
>>T=log(r)/n/log(1+0.01*p)
T=
11.5813
%Page20,ex4
>>x=-2:0.05:2;f=x.^4-2.^x;
>>[fmin,min_index]=min(f)
fmin=
-1.3907%最小值
min_index=
54%最小值点编址
>>x(min_index)
ans=
0.6500%最小值点
>>[f1,x1_index]=min(abs(f))%求近似根--绝对值最小的点
f1=
0.0328
x1_index=
24
>>x(x1_index)
ans=
-0.8500
>>x(x1_index)=[];f=x.^4-2.^x;%删去绝对值最小的点以求函数绝对值次小的点
>>[f2,x2_index]=min(abs(f))%求另一近似根--函数绝对值次小的点
f2=
0.0630
x2_index=
65

>>x(x2_index)
ans=
1.2500
%Page20,ex5
>>z=magic(10)
z=
929918156774515840
9880714167355576441
4818820225456637047
8587192136062697128
869325296168755234
17247683904249263365
2358289914830323966
7961395972931384572
10129496783537444653
111810077843643502759
>>sum(z)
ans=
505505505505505505505505505505
>>sum(diag(z))
ans=
505
>>z(:,2)/sqrt(3)
ans=
57.1577
46.1880
46.7654
50.2295
53.6936
13.8564
2.8868
3.4641
6.9282
10.3923
>>z(8,:)=z(8,:)+z(3,:)
z=
929918156774515840
9880714167355576441
4818820225456637047
8587192136062697128
869325296168755234
17247683904249263365
2358289914830323966

83871011151198387101115119
10129496783537444653
111810077843643502759

%Page40ex1
先在编辑器窗口写下列M函数，保存为eg2_1.m
function[xbar,s]=ex2_1(x)
n=length(x);
xbar=sum(x)/n;
s=sqrt((sum(x.^2)-n*xbar^2)/(n-1));
例如
>>x=[81706551766690876177];
>>[xbar,s]=ex2_1(x)
xbar=
72.4000
s=
12.1124
%Page40ex2
s=log(1);n=0;
whiles<=100
n=n+1;
s=s+log(1+n);
end
m=n
计算结果m=37
%Page40ex3
clear;
F(1)=1;F(2)=1;k=2;x=0;
e=1e-8;a=(1+sqrt(5))/2;
whileabs(x-a)>e
k=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);
end
a,x,k
计算至k=21可满足精度
%Page40ex4
clear;tic;s=0;
fori=1:1000000
s=s+sqrt(3)/2^i;
end
s,toc
tic;s=0;i=1;
whilei<=1000000
s=s+sqrt(3)/2^i;i=i+1;

end
s,t