MATLAB程序设计方案与应用第二版实验参考答案Word下载.docx
《MATLAB程序设计方案与应用第二版实验参考答案Word下载.docx》由会员分享，可在线阅读，更多相关《MATLAB程序设计方案与应用第二版实验参考答案Word下载.docx（33页珍藏版）》请在冰豆网上搜索。
F=size（C）
D=C（F
（1）-2:
F
（1）,F
（2）-1:
F
（2））
whos
%第四题
%
（1）:
A=100:
999。
B=rem（A,21）。
C=length（find（B==0））
%
（2）:
A='
lsdhKSDLKklsdkl'
。
k=find（A>
='
A'
&
A<
Z'
A（k）=[]
%实验二MATLAB矩阵分析与处理
E=eye（3）。
R=rand（3,2）。
O=zeros（2,3）。
S=diag（[2,3]）。
A=[E,R。
O,S]。
A^2
B=[E,（R+R*S）。
O,S^2]
H=hilb（5）
P=pascal（5）
Hh=det（H）
Hp=det（P）
Th=cond（H）
Tp=cond（P）
%第三题:
A=fix（10*rand（5））
H=det（A）
Trace=trace（A）
Rank=rank（A）
Norm=norm（A）
%第四题:
A=[-29,6,18。
20,5,12。
-8,8,5]
[V,D]=eig（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
B=[0.95,0.67,0.53]'
x=inv（A）*B
%（3）:
cond（A）
%第五题方法二:
A=hilb（4）
A（:
1）=[]
A（4,:
）=[]
B=[0.95,0.67,0.52]'
X=inv（A）*B
B1=[0.95,0.67,0.53]'
X1=inv（A）*B1
N=cond（B）
N1=cond（B1）
Na=cond（A）%矩阵A为病态矩阵
%第六题
A=[1,4,9。
16,25,36。
49,64,81]
B=sqrtm（A）
C=sqrt（A）%sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算
%实验三选择程序结构设计
%第一题程序一
x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0]。
y=[]。
%建立存放所有y值的矩阵
forx0=x
ifx0<
0&
x0~=-3
y=[y,x0*x0+x0-6]。
elseifx0>
x0<
5&
x0~=2&
x0~=3
y=[y,x0*x0-5*x0+6]。
else
y=[y,x0*x0-x0-1]。
end
end
x%输出所有x
y%输出所有y
%第一题程序二
x=[-5,-3,1,2,2.5,3,5]。
fora=1:
7
ifx（a）<
x（a）~=-3
y=[y,（x（a））^2+x（a）-6]。
elseifx（a）>
x（a）<
x（a）~=2&
x（a）~=3
y=[y,（x（a））^2-5*x（a）+6]。
else
y=[y,x（a）*x（a）-x（a）-1]。
%第二题程序一:
x=input（'
请输入一个百分制成绩：
'
ifx>
100|x<
disp（'
您输入的成绩不是百分制成绩，请重新输入。
ifx<
=100&
x>
=90
elseifx<
=89&
=80
B'
=79&
=70
C'
=69&
60
D'
E'
%第二题程序二：
s=input（'
请输入一个成绩（0分到100分之间）：
%s用于存放成绩
while1%判断输入成绩的合理性
ifs<
0|s>
100
disp（'
输入的成绩需在0到100之间，请重新输入：
）
s=input（'
break。
switchfix（s/10）%对成绩做出等级判断
case{9,10}
case8
case7
case6
otherwise
n=input（'
请输入员工工号：
h=input（'
该员工工作时数是：
ifh>
120
x=（h-120）*84*（1+0.15）+120*84。
elseifh<
x=h*84-700。
x=h*84。
disp（[num2str（n）,'
号员工'
'
的应发工资为'
num2str（x）]）。
%第四题（还可以用switch语句实现）
a=fix（10+（99-10）*rand（1,2））%产生两个随机整数
x=a
（1）。
y=a
（2）。
t=input（'
请输入运算符号：
s'
ift=='
+'
z=x+y。
elseift=='
-'
z=x-y。
*'
z=x*y。
/'
z=x/y。
disp（[num2str（x）,t,num2str（y）,'
num2str（z）]）%输出运算结果
%第五题
a=rand（5,6）%产生5x6的随机矩阵
请输入您要输出矩阵的第几行：
ifn>
5
超出了矩阵的行数，矩阵的最后一行为：
a（5,:
disp（['
矩阵的第'
num2str（n）,'
行为：
]）
a（n,:
%实验四循环结构程序设计
s=0。
n=?
fori=1:
n
s=s+1/i/i。
PI=sqrt（6*s）
pi
a=1:
n。
b=1./a.^2。
PI=sqrt（6*sum（b））
y=0。
n=1。
while（y<
y=y+1/（2*n-1）。
n=n+1。
y=y-1/（2*（n-1）-1）
n=n-2
a=input（'
a=?
b=input（'
b=?
Xn=1。
Xn1=a/（b+Xn）。
n=0。
whileabs（Xn1-Xn）>
1e-5
Xn=Xn1。
Xn1=a/（b+Xn）。
n=n+1。
ifn==500
Xn1
r1=（-b+sqrt（b*b+4*a））/2
r2=（-b-sqrt（b*b+4*a））/2
ifi==1
f（i）=1。
elseifi==2
f（i）=0。
elseifi==3
f（i）=f（i-1）-2*f（i-2）+f（i-3）。
max（f）
min（f）
sum（f）
length（find（f>
0））
length（find（f==0））
length（find（f<
%第五题：
fori=2:
49
b=i*（i+1）-1。
m=fix（sqrt（b））。
forj=2:
m
ifrem（b,j）==0
break
ifj==m
s=s+b。
s
%实验五函数文件
functiony=mat1（x）%建立函数文件mat1.m
y=[exp（x）,log（x）,sin（x）,cos（x）]。
%在命令窗口调用上述函数文件：
y=mat1（1+i）
%第二题程序一
function[a,b,N,M]=shiyanwu2（m,n,t）
A=[m*cos（t*pi/180）,-m,-sin（t*pi/180）,0。
m*sin（t*pi/180）,0,cos（t*pi/180）,0。
0,n,-sin（t*pi/180）,0。
0,0,-cos（t*pi/180）,1]。
B=[0,9.8*m,0,9.8*n]。
C=inv（A）*B'
a=C
（1）。
b=C
（2）。
N=C（3）。
M=C（4）。
%在命令窗口调用该函数文件:
m1=input（'
m1='
m2=input（'
m2='
theta=input（'
theta='
[a1,a2,N1,N2]=shiyanwu2（m1,m2,theta）
%第二题程序二
functionX=mat2（m1,m2,t）
g=9.8。
A=[m1*cos（t*pi/180）,-m1,-sin（t*pi/180）,0。
m1*sin（t*pi/180）,0,cos（t*pi/180）,0。
0,m2,-sin（t*pi/180）,0。
B=[0。
m1*g。
0。
m2*g]。
X=inv（A）*B。
%在命令窗口调用该函数文件：
X=mat2（1,1,60）
functionflag=mat3（x）
flag=1。
sqrt（x）
ifrem（x,i）==0
flag=0。
fori=10:
99
j=10*rem（i,10）+fix（i/10）。
ifmat3（i）&
mat3（j）
disp（i）
functiony=fx（x）
y=1./（（x-2）.^2+0.1）+1./（（x-3）.^4+0.01）。
y=fx
（2）
a=[1,2。
3,4]。
y=fx（a）
functionf1=mat5（n）
f1=n+10*log（n*n+5）。
%在命令窗口中调用该函数文件：
y=mat5（40）/（mat5（30）+mat5（20））
%
（2）方法一
functionf2=mat6（n）
f2=0。
f2=f2+i*（i+1）。
%在命令窗口中调用该函数文件如：
y=mat6（40）/（mat6（30）+mat6（20））
%
（2）方法二
functionf2=mat7（n）
i=1:
m=i.*（i+1）。
f2=sum（m）。
y=mat7（40）/（mat7（30）+mat7（20））
%实验六高层绘图操作
%第一题：
x=linspace（0,2*pi,101）。
y=（0.5+3*sin（x）./（1+x.^2））.*cos（x）。
plot（x,y）
%第二题：
x=linspace（-2*pi,2*pi,100）。
y1=x.^2。
y2=cos（2*x）。
y3=y1.*y2。
plot（x,y1,'
b-'
x,y2,'
r:
x,y3,'
y--'
text（4,16,'
\leftarrowy1=x^2'
text（6*pi/4,-1,'
\downarrowy2=cos（2*x）'
text（-1.5*pi,-2.25*pi*pi,'
\uparrowy3=y1*y2'
subplot（1,3,1）。
%分区
plot（x,y1）。
title（'
y1=x^2'
%设置标题
subplot（1,3,2）。
plot（x,y2）。
y2=cos（2*x）'
subplot（1,3,3）。
plot（x,y3）。
y3=x^2*cos（2*x）'
x=linspace（-2*pi,2*pi,20）。
subplot（2,2,1）。
bar（x,y1）。
y1=x^2的条形图'
subplot（2,2,2）。
stairs（x,y1）。
y1=x^2的阶梯图'
subplot（2,2,3）。
stem（x,y1）。
y1=x^2的杆图'
subplot（2,2,4）。
fill（x,y1,'
r'
%如果少了'
则会出错
y1=x^2的填充图'
%其他的函数照样做。
x=-5:
0.01:
5。
%起始设y为空向量
ifx0<
=0%不能写成x0=<
y=[y,（x0+sqrt（pi））/exp
（2）]。
%将x对应的函数值放到y中
y=[y,0.5*log（x0+sqrt（1+x0^2））]。
%第四题：
a='
b='
n='
t=-2*pi:
2*pi。
r=a*sin（b+n*t）。
polar（t,r）
x=linspace（-5,5,21）。
y=linspace（0,10,31）。
[x,y]=meshgrid（x,y）。
%在[-5,5]*[0,10]的范围内生成网格坐标
z=cos（x）.*cos（y）.*exp（-sqrt（x.^2+y.^2）/4）。
subplot（2,1,1）。
surf（x,y,z）。
subplot（2,1,2）。
contour3（x,y,z,50）。
%其中50为高度的等级数，越大越密
ezsurf（'
cos（s）*cos（t）'
cos（s）*sin（t）'
sin（s）'
[0,0.5*pi,0,1.5*pi]）。
%利用ezsurf隐函数
shadinginterp%进行插值着色处理
%实验七低层绘图操作
h=figure（'
MenuBar'
figure'
color'
WindowButtonDownFcn'
LeftButtonPressed'
）'
x=-2:
2。
y=x.^2.*exp（2*x）。
h=line（x,y）。
set（h,'
linestyle'
:
linewidth'
2）
text（1,exp
（2）,'
y=x^2*exp（2*x）'
0.00001:
0.001。
[t,x]=meshgrid（t）。
v=10*exp（-0.01*x）.*sin（2000*pi*t-0.2*x+pi）。
axes（'
view'
[-37.5,30]）。
h=surface（t,x,v）。
v=10*exp（-0.01*x）.*sin（2000*pi*t-0.2*x+pi）'
xlabel（Ct'
）,ylabel（'
x'
）,zlabel（'
v'
x=0:
y1=sin（x）。
y2=cos（x）。
y3=tan（x）。
y4=cot（x）。
plot（x,y4）。
cylinder（5）。
light（'
Position'
[0,1,1]）。
materialshiny
%实验八数据处理与多项式运算
A=rand（1,30000）。
b=mean（A）
std（A,0,2）
max（A）
min（A）
30000
ifA（i）>
0.5
p=n/30000
A=45+51*rand（100,5）。
[Y,U]=max（A）
[a,b]=min（A）
m=mean（A）
s=std（A）
sum（A,2）
[Y,U]=max（ans）
[a,b]=min（ans）
[zcj,xsxh]=sort（ans）
h=6:
2:
18。
x=6.5:
17.5。
t1=[18,20,22,25,30,28,24]。
t2=[15,19,24,28,34,32,30]。
T1=spline（h,t1,x）
T2=spline（h,t2,x）
x=1:
101。
y1=log10（x）。
p=polyfit（x,y1,5）
y2=polyval（p,x）。
p1=[1,2,4,0,5]。
p2=[1,2]。
p3=[1,2,3]。
p=p1+[0,conv（p2,p3）]%为使两向量大小相同，所以补0
A=roots（p）
A=[-1,1.2,-1.4。
0.75,2,3.5。
0,5,2.5]。
polyval（p,A）
polyvalm（p,A）
实验十
程序：
x=sym（'
6'
y=sym（'
5'
z=（x+1）/（sqrt（3+x）-sqrt（y））
1、分解因式
（1）
symsxy。
A=x^4-y^4。
factor（A）
（2）
factor（sym（'
5135'
））
3、化简表达式
symsbeta1beta2
y=sin（beta1）*cos（beta2）-cos（beta1）*sin（beta2）
simple（y）
symsx
y=（4*x^2+8*x+3）/（2*x+1）
5、用符号方法求下列极限或导数
f=（x*（exp（sin（x））+1）-2*（exp（tan（x））-1））/（sin（x））
limit（f）
y=（sqrt（pi）-sqrt（acos（x）））/（sqrt（x+1））。
limit（f,x,-1,'
right'
（3）
y=（1-cos（2*x））/x。
y1=diff（y）
y2=diff（y,x,2）
6、用符号方法求下列积分
f=1/（1+x^4+x^8）
int（f）
f=1/（（（asin（x））^2）*sqrt（1-x^2））
f=（x^2+1）/（x^4+1）
int（f,x,0,inf）
（4）
f=exp（x）*（1+exp（x））^2
y=int（f,x,0,log
（2））
double（y）
0
0
版权申明
本文部分内容，包括文字、图片、以及设计等在网上搜集整理。
版权为潘宏亮个人所有
Thisarticleincludessomeparts,includingtext,pictures,anddesign.CopyrightisPanHongliang'
spersonalownership.
用户可将本文的内容或服务用于个人学习、研究或欣赏，以及其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关法律的规定，不得侵犯本网站及相关权利人的合法权利。
除此以外，将本文任何内容或服务用于其他用途时，须征得本人及相关权利人的书面许可，并支付报酬。
Usersmayusethecontentsorservicesofthisarticleforpersonalstudy,researchorappreciation,andothernon-commercialornon-profitpurposes,butatthesametime,theyshallabidebythe
升级会员 