matlab数学实验复习题有答案Word格式.docx
《matlab数学实验复习题有答案Word格式.docx》由会员分享,可在线阅读,更多相关《matlab数学实验复习题有答案Word格式.docx(17页珍藏版)》请在冰豆网上搜索。
![matlab数学实验复习题有答案Word格式.docx](https://file1.bdocx.com/fileroot1/2022-11/17/01651bbe-a4b3-453d-aa01-d6e57c1d606e/01651bbe-a4b3-453d-aa01-d6e57c1d606e1.gif)
10
A-3=
A.
A2
>
=
16
2536
tril
(A
)=45
6
49
6481
78
0(
triu
(A,
-1
)
dia
ig(A):
=0
5(
0A(:
2),=
5A(3,:
)=
0!
10、normcdf(1,1,2)=0.5%正态分布mu=1sigma=2,x=1处的概率
11、unifpdf([5,7],2,6)=【0.25;
0】
11、命令formatshort的作用保留小数点后四位而formatlong:
保留小数点后14位
12、formatrat的作用是最接近的有理数
12、interp1(x0,y0,x)的作用是求以x0,y0为节点数组,x为插值点数组的分段线性插值
13、13、[a,b,c,d]=fzero(fun,x0)中参数的涵义是a是变号点的近似值,b是对应,的函数值,c是停止运行的原因(c=1即为找到该点,c=0就是没有找到)d是一个结构变量,fun是求解方程的函数M文件,x0是零点或变号点附近的值。
14、龙格-库塔方法可用如下MATLAB^令求解微分方程
[t,x]=ode45(f,[a,b],x0),中参数的涵义是fun是求解方程的函数M文件,
[a.b]是输入向量即自变量的围a为初值,x0为函数的初值,t为输出指定的[a,b],x为函数值
15、写出下列命令的功能:
axisequal纵、横坐标轴采用等长刻度
text(1,2,‘y=sin(x)'
)在x=1,y=2处加上字符串y=sin(x);
holdon把新的plot产生的图形画在原来的图形上。
title('
y=sin(x)'
)在图形正上方加上字符串y=sin(x)
16、Matlab中自定义函数M文件的第一行必须以function开头;
17、二种数值积分的库函数名为:
quad;
quadl
18、unifrnd(1,2,3,4)的功能是:
随机生成3行4列均匀分布,每个元素服从(1,2)的矩阵_
19、binornd(20,0.3,3,4)的功能是随机生成3行4列服从(20,0.3)的二项分布的矩阵—
20、eig(A)的功能是矩阵_A的特征值
21、设x是一向量,则hist(x)的功能是作出将X十等分的直方图
22、interp1([1,2,3],[3,4,5],2.5)
Ans=4.5
23、建立一阶微分方程组
x(t)y(t)
3xy的函数M文件。
(做不出来)
4x3y
二、写出运行结果:
10
1、>
eye(3,4)=01
00
2、>
size([1,2,3])=1
3、设b=round(unifrnd(-5,5,1,4)),则=352-5>
[x,m]=min(b);
x=-5;
m=4
[x,n]=sort(b)
-5235
4312
mean(b)=1.25,median(b)=2.5,range(b)=10
4、向量b如上题,则
any(b),all(b<
2),all(b<
6)
Ans=101
5、>
[56;
78]>
[78;
56]=
12
6、若B34,则
7、>
diag(diag(B))=
&
>
[4:
-2:
1]*[-1,6]=-412
9、>
acos(0.5),atan
(1)ans=
1.6598
ans=
0.7448
10、>
norm([1,2,3])
Ans=3.3941
11、>
length([1,3,-1])=3
12、>
x=0:
0.4:
2;
plot(x,2*x,'
k*'
13、>
zeros(3,1);
ans=
11
14、
ones(3)=1
11,vander([2,3,5])=
25
16、
floor(1:
0.3:
3)=
11222
18
、
subplot(2,2,1);
fplot('
sin'
[0,2*pi]);
subplot(2,2,2);
plot([1,2,-1]);
x=linspace(0,6*pi);
subplot(2,2,3);
plot3(cos(x),sin(x),x);
subplot(2,2,4);
polar(x,5*sin(4*x/3));
19、>
t=linespace(0,2,11)
0.00.20.40.60.81.01.21.41.61.82.0
20、>
[a,b]=binostat(15,0.2)a=3b=2.4
y1=binopdf(5,10,0.7)=0.1029,y2=binocdf(5,10,0.7)=0.1503
21、>
log10([1,10,100])=[01^]
22、>
p=1;
fork=2:
3:
9p=p*k;
end;
pp=8
23、>
s=0;
9s=s+k;
ss=15
24、finline('
xA3exp(2*x)3'
);
feval(f,1)
Ans=3.8647
25、>
a1=norminv(0.6,3,4)a1=4.0134
26、>
unifinv(0.4,1,5),unifpdf(0.4,1,5),unifpdf(2,1,5)
Ans=2.600.25
27、>
A=[01-1;
210;
1-11];
01-1
1-11
A([1,3],:
01-1
A([3,1],:
)=1-11
01-1
A(2,:
)=210
-2*A(1,:
)=0-22
28、>
quad(‘sin(x)'
0,pi/2)=1.0000
29、>
trapz([3,4,6],[1,2,3])=6.5000
30、>
int('
x-sin(x)'
0,1)
Ans=
cos
(1)-1/2
31、>
round(3:
5),ceil(3:
5);
floor(3:
5)
334455
333445
limit(1+1/(3*x)Ax,inf)=1_
diff(sin(3*x)+xA3,2)=6*x-9*sin(3*x)
taylor(exp(3*x),5,1):
命令输入:
y=taylor(exp(3*x),x,1,'
Order'
5)
exp(3)+3*exp(3)*(x-1)+(9*exp(3)*(x-1)八2)/2+(9*exp(3)*(x-1F3)/2+(27*exp(3)*(x-1)八4)/8
a1=mod(15,4),b1=rem(15,4)=3,3_>
a2=mod(-15,-4),b2=rem(-15,-4)=-3,-3
a3=mod(15,-4),b3=rem(15,-4)=-1,
a4=mod(-15,4),b4=rem(-15,4)=1,-3
34、>
x=binornd(20,0.4,2,4)
87108
107912
sign(x),
1111
y=-poissrnd(8,2,4)
-16-108-7
-7-8-6-9
sign(y)
-1-1-1-1
35、>
[a1,b1]=binostat(20,0.4)a1=8b1=4.8
[a2,b2]=poisstat(8)ans=8,8
[a3,b3]=chi2stat(15)ans=[1530]
36、运行M文件:
chi2fig
n=5;
a=0.9;
xa=chi2inv(a,n);
x=0:
0.1:
15;
y=chi2pdf(x,n);
plot(x,y,'
b'
holdon;
xf=0:
xa;
yf=chi2pdf(xf,n);
fill([xf,xa],[yf,0],'
g'
text(xa*1.01,0.005,num2str(xa));
text(2.5,0.05,'
alpha=0.9'
'
fontsize'
20);
text(9,0.09,'
X~{\chi}A2(4)'
'
16);
x=[linspace(0,2*pi,100)];
y=sin(x).*cos(x);
y仁spline(xO,yO,x);
[x;
y;
y1]'
38、>
quadl('
exp(2*x).*log(3*x)'
1,3)
ans=
398.6352
39、x0=0:
2*pi/6:
2*pi;
y0=sin(x0).*cos(x0);
plot(x,y,'
k'
x,y1,'
b-'
注:
此处省略100组数
40、>
A=round(unifrnd(0,100,3,3));
[L,U]=lu(A)
L=
0.9897
0.4699
1.0000
0.1649
97.000080.000092.0000
035.804126.8247
00-89.6568
41、a=sparse([133],[235],[123],4,5);
s=full(a)
s=
、编程