Matlab程序设计与应用第二版刘卫国.docx
《Matlab程序设计与应用第二版刘卫国.docx》由会员分享，可在线阅读，更多相关《Matlab程序设计与应用第二版刘卫国.docx（22页珍藏版）》请在冰豆网上搜索。
Matlab程序设计与应用第二版刘卫国
实验一：
T1：
%%第一小题
z1=2*sin（85*pi/180）/（1+exp
（2））
%%第二小题
x=[2,1+2i;-0.45,5];
z2
%%第三小题
a=-3.0:
0.1:
3.0;
0.3）+log（（0.3+a）/2）
%%第四题
t=0:
0.5:
2.5
z4=（t>=0&t<1）.*（t.^2）+（t>=1&t<2）.*（t
.^2-1）+（t>=2&t<3）.*（t.^2-2*t+1）
T2：
A=[12,34,-4;34,7,87;3,65,7]
B=[1,3,-1;2,0,3;3,-2,7]
disp（'A+6*B='）;
disp（A+6*B）;
disp（'A-B+I='）;disp（A-B+eye（3））;
disp（'A*B='）;
disp（A*B）;
disp（'A.*B='）;
disp（A.*B）;
disp（'A^3='）;
disp（A^3）;
disp（'A.^3='）;
disp（A.^3）;
disp（'A/B='）;
disp（A/B）;
disp（'B\A='）;
disp（B\A）;
disp（'[A,B]='）;
disp（[A,B]）;
disp（'[A（[1,3],:
）;B^2]='）;
disp（[A（[1,3],:
）;B^2]）;
T3:
z=1:
25;
A=reshape（z,5,5）';
B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C（3:
5,2:
3）
T4-1:
a=100:
999;
b=find（rem（a,21）==0）;
c=length（b）
T4-2:
a=input（'请输入一个字符串：
','s'）;
b=find（a>='A'&a<='Z'）;
a（b）=[];
disp（a）;
实验二：
T1：
E=eye（3）,R=rand（3,2）,O=zeros（2,3）,S=d
iag（[1,2]）;
A=[E,R;O,S]
disp（'A^2='）;
disp（A^2）;
disp（'[E,R+RS;O,S^2]'）;
B=[E,R+R*S;O,S^2]
T2：
H=hilb（5）
P=pascal（5）Hh=det（H）
Hp=det（P）
Th=cond（H）
Tp=cond（P）
a=abs（Th-1）;
b=abs（Tp-1）;
ifa>b
disp（'帕萨卡矩阵P性能更好'）;elseifadisp（'希尔伯特矩阵H性能更好'）;else
disp（'两个矩阵性能相同'）;
end
T3:
a=1:
25;
A=reshape（a,5,5）
disp（'行列式的值：
'）;
disp（det（A））;
disp（'矩阵的秩：
'）;
disp（rank（A））;
disp（'矩阵的迹：
'）;
disp（trace（A））;disp（'矩阵的范数：
'）;
disp（norm（A））;
T4:
A=[-29,6,18;20,5,12;-8,8,5]
[V,D]=eig（A）
T5:
B=[0.95,0.67,0.52]'
X1=A\B
B（3）=0.53
X2=A\B
disp（'系数矩阵A的条件数：
'）;
disp（cond（A））;
T6:
a=1:
25;
A=reshape（a,5,5）
disp（'B1=sqrtm（A）'）;
B1=sqrtm（A）
disp（'B2=sqrt（A）'）;
B2=sqrt（A）
disp（'B1*B1'）;
B1*B1disp（'B2.*B2'）;
B2.*B2
实验三：
T1:
x=-5.0:
2:
5.0
fori=1:
length（x）;
if（x（i）<0）&（x（i）~=-3）
y（i）=x（i）^2+x（i）-6;
elseif
（x（i）>=0）&（x（i）<5）&（x（i）~=2）&（x（i）~=3）;
y（i）=x（i）^2-5*x（i）+6;
else
y（i）=x（i）^2-x（i）-1;
end
end
y
T2-if:
s=input（'pleaseenterthescore:
'）;while（s<0||s>100）
disp（'thescoreisnotreasonable'）;s=input（'pleaseenterthescore:
'）;
end
if（s>=90&s<=100）;
disp（'A'）;
elseif（s>=80&s<90）;
disp（'B'）;
elseif（s>=70&s<80）;
disp（'C'）;
elseif（s>=60&s<70）;
disp（'D'）;
else
disp（'E'）;
end
T2-switch:
s=input（'pleaseenterthescore:
'）;
while（s<0||s>100）
disp（'thescoreisnotreasonable'）;
s=input（'pleaseenterthescore:
'）;
end
switchfix（s/10）
case{9,10}
disp（'A'）;
case{8}disp（'B'）;
case{7}
disp（'C'）;
case{6}
disp（'D'）;
case{0,1,2,3,4,5}
disp（'E'）;
end
T3:
t=input（'请输入工时time='）;
ift>120
w=120*84+（t-120）*1.15*84;elseift<60
w=t*84-700;
else
w=84*t;
end
disp（'应发工资为：
'）;
disp（w）;
T4:
a=10+floor（rand
（1）*89）
b=10+floor（rand
（1）*89）s=input（'请输入+或-或*或/','s'）;
while（s~='+'&&s~='-'&&s~='*'&&s~
='/'）
disp（'输入的符号错误，请重新输入'）;
s=input（'请输入+或-或*或/','s'）;
end
switchs
case{'+'}
c=a+b;
case{'-'}
c=a-b;
case{'*'}
c=a*b;
case{'/'}
c=a/b;
end
c
T5:
A=rand（5,6）
n=input（'请输入n:
'）;
while（n<1）
disp（'输入的n有误，请重新输入'）;n=input（'请输入n:
'）;
end
ifn>5
B=A（5,:
）
else
B=A（n,:
）
End
实验四：
T1:
n=input（'请输入n:
'）;
m=0;
fori=1:
n;
m=m+1/i^2;
end
pii=sqrt（6*m）
T1-2:
n=input（'请输入n:
'）;
m=1./（1:
n）.^2;
s=sum（m）;
pii=sqrt（6*s）
T2:
n=1000;y=0;
fori=1:
n;
y=y+1/（2*i-1）;
ify>=3
disp（'最大n值为：
'）;
i-1
disp（'对应的y值为：
'）;
y-1/（2*i-1）
break;
end
end
T3:
formatlong
a=input（'请输入a:
'）;
b=input（'请输入b:
'）;
x
（1）=a/（1.0+b）;i=1;x
（2）=a/（x（i）+b）;
while（abs（x（i+1）-x（i））>0.00001&i<500）
i=i+1;
x（i+1）=a/（x（i）+b）;
end
r1=（-b+sqrt（b^2+4*a））/2;
r2=（-b-sqrt（b^2+4*a））/2;
disp（['x',num2str（i+1）,'=',num2str（x（i+1）,8）]）;
disp（['x',num2str（i）,'=',num2str（x（i）,8）]）;disp（['r1=',num2str（r1,8）,',r2=',num2str（r2,8）]）;
formatshort
T4:
f
（1）=1;f
（2）=0;f（3）=1;i=4;
while（i<=100）
f（i）=f（i-1）-2*f（i-2）+f（i-3）;
i=i+1;
end
fmax=max（f）
fmin=min（f）
fsum=sum（f）
fp=length（find（f>0））
fn=length（find（f<0））
f0=length（find（f==0））
T5:
j=0;s=0;
fori=3:
50
m=i*（i-1）-1;
ifisprime（m）==1
m
s=s+m;
j=j+1;
end
end
disp（['亲密数的对数：
',num2str（j）]）;
disp（['所有亲密素数之和：
',num2str（s）]）;
实验五：
fn1:
functionf=fn1（n）
f=n+10*log（n.^2+5）;
fn2:
functionf=fn2（n）
f=0;
fori=1:
n
f=f+i*（i+1）;
end
fx:
functionf=fx（x）
f=1./（（x-2）.^2+0.1）+1./（（x-3）.^4+0.01）;
T1:
x=input（'请输入一个复数：
'）;fe=exp（x）
flg=log（x）
fsin=sin（x）
fcos=cos（x）
T2:
m1=input（'请输入m1:
'）;
m2=input（'请输入m2:
'）;
t=input（'请输入theta:
'）;
A=[m1*cos（t）,-m1,-sin（t）,0;m1*sin（t）,0,...
cos（t）,0;0,m2,-sin（t）,0;0,0,-cos（t）,1]B=[0,m1*9.8,0,m2*9.8]'
disp（'[a1,a2,N1,N2]='）;
disp（[A\B]'）;
T3:
j=1;n=1;k=1;
fori=10:
99%挑选出10~99中所有的素数
ifisprime（i）==1
m（j）=i;
j=j+1;
end
end
fort=1:
length（m）;
s（n）=10*rem（m（t）,10）+fix（m（t）/10）;%
挑选出的素数进行位置交换
ifisprime（s（n））==1%再挑选出交
换绝对素数
r（k）=m（t）;
k=k+1;
end
n=n+1;
end
disp（['所有两位绝对素数：
',num2str（r）]）
T5:
y1=fn1（40）/（fn1（30）+fn1（20））%函数
调用
y2=fn2（40）/（fn2（30）+fn2（20））
实验六：
T1:
x=linspace（0,2*pi,101）;
y=（0.5+3*sin（x）./（1+x.^2））.*cos（x）;
plot（x,y）;T2:
x=linspace（0,2*pi,100）;
y1=x.^2;
y2=cos（2*x）;
y3=y1.*y2;
subplot（2,2,1）;
plot（x,y1,'r:
',x,y2,'y-',x,y3,'b-.'）;subplot（2,2,2）;
plot（x,y1,'r:
'）;
subplot（2,2,3）;
plot（x,y2,'y-'）;
subplot（2,2,4）;
plot（x,y3,'b-.'）;
T2-3:
x=linspace（0,2*pi,10）;
y1=x.^2;
y2=cos（2*x）;
y3=y1.*y2;
subplot（3,4,1）;
bar（x,y1,'r'）;axis（[0,7,0,40]）;
subplot（3,4,2）;
stairs（x,y1,'y'）;axis（[0,7,0,40]）;
subplot（3,4,3）;
stem（x,y1,'b'）;axis（[0,7,0,40]）;
subplot（3,4,4）;
fill（x,y1,'g'）;axis（[0,7,0,40]）;
subplot（3,4,5）;
bar（x,y2,'r'）;axis（[-1,7,-1,1]）;
subplot（3,4,6）;
stairs（x,y2,'y'）;axis（[-1,7,-1,1]）;
subplot（3,4,7）;
stem（x,y2,'b'）;axis（[-1,7,-1,1]）;
subplot（3,4,8）;
fill（x,y2,'g'）;axis（[-1,7,-1,1]）;
subplot（3,4,9）;
bar（x,y3,'r'）;axis（[0,7,-30,50]）;
subplot（3,4,10）;
stairs（x,y3,'y'）;axis（[0,7,-30,50]）;
subplot（3,4,11）;
stem（x,y3,'b'）;axis（[0,7,-30,50]）;
subplot（3,4,12）;
fill（x,y3,'g'）;axis（[0,7,-30,50]）;
T3:
x=linspace（-5,5,100）;
y=[];forx0=x;
ifx0<=0;
y=[y,（x0+sqrt（pi））/exp
（2）];else
end
end
subplot（1,2,1）;
plot（x,y）;
fori=1:
length（x）
ifx（i）<=0
y（i）=（x（i）+sqrt（pi））/exp
（2）;else
end
end
subplot（1,2,2）;
plot（x,y）;
T4:
a=input（'请输入a:
'）;
b=input（'请输入b:
'）;
n=input（'请输入n:
'）;
theta=0:
0.001:
2*pi;
y=a.*sin（b+n.*theta）;
polar（theta,y）;
T5:
x=linspace（-5,5,21）;
y=linspace（0,10,31）;
[x,y]=meshgrid（x,y）;
z=cos（x）.*cos（y）.*exp（-sqrt（x.^2+y.^2）/
4）;
subplot（1,2,1）;
surf（x,y,z）;
subplot（1,2,2）;
contour3（x,y,z）;
T6:
s=0:
pi/2;
t=0:
3*pi/2;
[s,t]=meshgrid（s,t）;
x=cos（s）.*cos（t）;
y=cos（s）.*sin（t）;z=sin（s）;
surf（x,y,z）;
shadinginterp;
实验八：
T1:
A=rand（1,30000）;
disp（'均值为：
'）;
mean（A）
disp（'标准方差为：
'）;
std（A）
disp（'最大元素为：
'）;
max（A）
disp（'最小元素为：
'）;
min（A）
k=find（A>0.5）;
b=length（k）;
disp（'大于0.5的随机数个数占总数的百分比：
'）;
sprintf（'%f%%',b/300）
T2:
P=45+50*rand（100,5）
disp（'最高分及序号：
'）;
[Y,U]=max（P,[],1）
disp（'最低分及序号：
'）;
[Y1,U1]=min（P,[],1）
disp（'每门课平均分：
'）;
mean（P,1）
disp（'每门课的标准方差：
'）;
std（P,0,1）
s=sum（P,2）;
disp（'总分最高分：
'）;
[Y2,U2]=max（s,[],1）
disp（'总分最低分：
'）;
[Y3,U3]=min（s,[],1）
[zcj,xsxh]=sort（s,1,'descend'）%desce
nd降序
T3:
h=6:
2:
18;
t1=[18.0,20.0,22.0,25.0,30.0,28.0,24.0];
t2=[15.0,19.0,24.0,28.0,34.0,32.0,30.0];
t=6.5:
2:
17.5;
T1=interp1（h,t1,t,'spline'）
T2=interp1（h,t2,t,'spline'）
T4:
x=1:
10:
101;
y=log10（x）;
p=polyfit（x,y,5）
x1=1:
0.01:
101;
y1=log10（x1）;
y2=polyval（p,x1）;
plot（x1,y1,'r-',x1,y2,'b-.'）;
T5:
P1=[1,2,4,0,5];
P2=[0,0,0,1,2];
P3=[0,0,1,2,3];
P4=conv（P2,P3）;
sp1=length（P1）;
sp2=length（P4）;
P=[zeros（1,sp2-sp1）,P1]+P4
x=roots（P）
A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];Ax=polyval（P,A）
AA=polyvalm（P,A）
实验九：
T1:
f=inline（'det（[xx^2x^3;12*x3*x^2;02
6*x]）'）;
i=1;x=1;
whilex<=3.01
g（i）=f（x）;
i=i+1;
x=x+0.01;
end
dx=diff（g）/0.01;
dx
（1）
dx（101）
dx（length（g）-1）
T2:
f1=inline（'sqrt（cos（t.^2）+4*sin（2*t）.^2
+1）'）
f2=inline（'log（1+x）./（1+x.^2）'）
I1=quad（f1,0,2*pi）
I2=quad（f2,0,2*pi）
T3:
A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2];
B=[-4,13,1,11]';
x1=A\B
[L,U]=lu（A）;
x2=U\（L\B）[Q,R]=qr（A）;
x3=R\（Q\B）
fun:
functionF=fun（X）;
x=X
（1）;
y=X
（2）;
z=X（3）;
F
（1）=sin（x）+y^2+log（z）-7;
F
（2）=3*x+2^y-z^3+1;
F（3）=x+y+z-5;
T5:
f=inline（'3*x+sin（x）-exp（x）'）;
fzero（f,1.5）
X=fsolve（'fun',[1
1]',optimset（'Display','off'））
实验十：
T1:
x=sym（'6'）;
y=sym（'5'）;
z=（x+1）/（sqrt（3+x）-sqrt（y））
T2:
1
symsxy;
t=sym（'5135'）;
factor（x^4-y^4）
factor（t）
T3:
symsbeta1beta2x;
simplify（sin（beta1）*cos（beta2）-cos（beta
1）*sin（beta2））
simplify（（4*x^2+8*x+3）/（2*x+1））
T4:
symsabcdefghk;
p1=[0,1,0;1,0,0;0,0,1];
p2=[1,0,0;0,1,0;1,0,1];
A=[a,b,c;d,e,f;g,h,k];
B=p1*p2*A
B1=inv（B）
B1*B;
tril（B）
det（B）
T5:
symsxtay
f1=（x*（exp（sin（x））+1）-2*（exp（tan（x））-1））/sin（x）^3;
limit（f1,x,0）%
（1）
f2=（sqrt（pi）-sqrt（acos（x）））/sqrt（x+1）;limit（f2,x,-1,'right'）%
（2）
f3=（1-cos（2*x））/x;
diff（f3,x,1）%（3）
diff（f3,x,2）
A=[a^x,t^3;t*cos（x）,log（x）];
diff（A,x,1）&（4）
diff（A,t,2）
diff（diff（A,x,1）,t,1）
f5=（x^2-2*x）*exp（-x^2-y^2-x*y）;yx=-diff（f5,x,1）/diff（f5,y,1）%（5）fxy=diff（diff（f5,x,1）,y,1）
x=sym（'0'）;y=sym（'1'）;
eval（fxy）
T6:
x=sym（'x'）;
f1=1/（1+x^4+x^8）;
int（f1,x）
f2=1/（asin（x））^2/sqrt（1-x^2）;int（f2,x）
f3=（x^2+1）/（x^4+1）;
int（f3,x,0,inf）
f4=exp（x）*（1+exp（x））^2;
int（f4,x,0,log
（2））
实验十一：
T1：
symsnx;
S1=symsum（1/（2*n-1）,n,1,10）
S2=symsum（n^2*x^（n-1）,n,1,inf）
T2：
x=sym（'x'）;
f=log（x）;
taylor（f,x,6,1）
T3：
x1=solve（'log（1+x）-5/（1+sin（x））=2'）
x2=solve（'x^2+9*sqrt（x+1）-1=0'）
x3=solve（'3*x*exp（x）+5*sin（x）-78.5=0'）
[x4,y4]=solve（'sqrt（x^2+y^2）-100=0','3
*x+5*y-8=0'）
T4：
dsolve（'D2y+4*Dy+29*y=0','y（0）=0','Dy（0）=15','x'）
T5：
[x,y,z]=dsolve（'Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z','t'）欢迎下载，谢谢观看！
资料仅供参考学习
升级会员 