MATLAB程序设计与应用第二版刘卫国实验部分Word文档下载推荐.docx
《MATLAB程序设计与应用第二版刘卫国实验部分Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《MATLAB程序设计与应用第二版刘卫国实验部分Word文档下载推荐.docx（15页珍藏版）》请在冰豆网上搜索。
A=[1,2,3,4,5;
6,7,8,9,10;
11,12,13,14,15;
16,17,18,19,20;
21,22,23,24,25];
B=[3,0,16;
17,-6,9;
0,23,-4;
9,7,0;
4,13,11];
C=A*B
D=C（3:
end,2:
end）
Namevalueclass
A<
5×
5double>
double
B<
3double>
C<
D[520397;
705557;
…]double
4.
X=100:
999;
length（find（~rem（X,21）））
y='
ahHGJjskjJjsdsaSa'
;
M=y（find（y<
='
z'
&
y>
a'
））
实验二
实验三
1x=input（'
请输入x的值'
）
ifx<
0&
x~=-3
y=x^2+x-6;
elseifx>
x<
5&
x~=2&
x~=3
y=x^2-5*x+6;
else
y=x^2-x-1;
end
y
2
x=input（'
请输入成绩'
ifx>
=90&
=100
disp（'
等级为A'
=80&
=89
等级为B'
=70&
=79
等级为C'
=60&
=69
等级为D'
=59
等级为E'
成绩有误'
switch（x）
case{90,100}
case{80,89}
case{70,79}
case{60,69}
case{0,59}
otherwise
成绩错误'
3
Num=input（'
输入员工的工号:
'
）;
Time=input（'
该工号员工的工作时数:
salary=input（'
该员工底薪:
ifTime>
=120
salary=salary*（1+0.15）;
elseifTime<
=60
salary=salary-700;
elseifTime>
60&
Time<
120
salary=Time*84;
salary=salary
实验四
1：
程序：
n=input（'
输入n:
fora=1:
n;
b（a）=1/（a^2）;
pi=sqrt（6*sum（b））
2：
.
sum=0;
forn=1:
200;
b（n）=1/（2*n-1）;
sum=sum+b（n）;
if（sum>
=3）
break;
end
n-1
f
（1）=1;
max=0;
f
（2）=0;
min=0;
f（3）=1;
c=0;
zhengshu=0;
fushu=0;
zero=0;
forn=4:
100
f（n）=f（n-1）-2*f（n-2）+f（n-3）;
form=1:
if（f（m）>
max）
max=f（m）;
if（f（m）<
min）
min=f（m）;
c=sum（f）;
if（f>
0）
zhengshu=zhengshu+1;
if（f（m）<
fushu=fushu+1;
if（f（m）==0）
zero=zero+1;
if（f（m）>
disp（'
最大值'
max
最小值'
min
和'
c
正数的个数'
zhengshu
负数的个数'
fushu
零的个数'
zero
实验五
1，函数文件：
文件名：
Yunsuan
function[z,d,s,c]=yunsuan（a,b）
z=exp（a+b*i）;
d=log（a+b*i）;
s=sin（a+b*i）;
c=cos（a+b*i）;
disp（['
复数'
num2str（a）,'
+'
num2str（b）,'
i的指数为:
num2str（z）,'
，对数为：
num2str（d）,'
，正弦为：
num2str（s）,'
，余弦为：
num2str（c）]）
函数调用：
clear
请输入复数实部的值'
y=input（'
请输入复数虚部的值'
[x,y]=yunsuan（x,y）;
2函数文件：
zhu
function[p,q,w]=zhu（p,q,w）
A=[p*cos（w）,-p,-sin（w）,0;
p*sin（w）,0,cos（w）,0;
0,q,-sin（w）,0;
0,0,-cos（w）,1];
B=[0;
9.8*p;
0;
9.8*q];
X=A\B
m1=input（'
m1='
m2=input（'
m2='
th=input（'
theta='
theta=pi*th/180;
[m1,m2,theta]=zhu（m1,m2,theta）;
4函数文件：
hanshu
function[y]=hanshu（x）
y=1/（（x-2）.^2-0.1）+1./（（x-3）.^4+0.01）;
实验六
1、
、
x=0:
0.02*pi:
2*pi;
y=（0.5+3*sin（x）./（1+x.*x））.*cos（x）
plot（x,y）
2、
（1）
x=（0:
2*pi）'
y1=x.*x;
y2=cos（2*x）;
y3=y1.*y2;
plot（x,y1,'
y-'
）,holdon
plot（x,y2,'
c-'
plot（x,y3,'
k-'
）,holdoff
（2）
subplot（2,2,1）
stairs（x,y1）;
r-'
subplot（2,2,2）
stairs（x,y2）;
b-'
subplot（2,2,3）
stairs（x,y3）;
（3）
（1）绘制条形图
0.2*pi:
bar（x,y1）
subplot（2,2,2）
bar（x,y3）
（2）绘制梯形图
0.05*pi:
（4）杆图：
stem（x,y1）;
stem（x,y2）;
stem（x,y3）;
（5）填充图
程序：
x=（0:
fill（x,y1,'
r'
fill（x,y2,'
b'
fill（x,y3,'
g'
3、：
x=linspace（-5,5,1000）;
y=[];
forx0=x
ifx0<
=0
y=[y,（x0+sqrt（pi））./exp
（2）];
else
y=[y,0.5*log（x0+sqrt（1+x0.*x0））];
4、
a=input（'
请输入a的值'
b=input（'
请输入b的值'
请输入n的值'
theta=0:
0.01:
rho=a.*sin（b+n.*theta）
polar（theta,rho）
5、
x=-linspace（-5,5,21）;
y=linspace（0,10,31）;
[x,y]=meshgrid（x,y）
z=cos（x）.*cos（y）.*exp（-0.25*sqrt（x.^2+y.^2））
subplot（1,2,1）;
title（'
surf（x,y,z）'
surf（x,y,z）;
subplot（1,2,2）;
surfc（x,y,z）'
surfc（x,y,z）;
实验八
x=rand（1,30000）;
EX=mean（x）%返回向量x的算术平均值
S=std（x）%返回向量x的标准方差
MAX=max（x）%返回向量x的最大元素
MIN=min（x）%返回向量x的最小元素
p=sum（x>
0.5）/30000%满足条件的百分比
p=100*rand（100,5）;
[MAXOBJECT,U]=max（p）%求每门课的最高分,U记录最大元素的行号
[MAXOBJECT,V]=min（p）%求每门课的最低分,V记录最小元素的行号
EX=mean（p）%每门课的平均分
S=std（p）%每门课的标准方差
SUM_1=sum（p,2）;
%求每名学生的总成绩
[MAXSUM,W]=max（SUM）%求五门课总分的最高分,W记录最高分学生的序号
[MINSUM,Y]=min（SUM）%求五门课总分的最低分，Y记录最低分学生的序号
[zcj,xsxh]=sort（SUM_1）;
zcj（end:
-1:
1）;
xsxh（end:
实验九
2、程序：
I1=quad（'
sqrt（cos（t.^2）+4*sin（2*t）.^2+1）'
0,2*pi）
I2=quad（'
（log（1+x））./（1+x.*x）'
0,1）
3、
A=[6,5,-2,5;
9,-1,4,-1;
3,4,2,-2;
3,-9,0,2]
b=[-4;
13;
1;
11]
x=inv（A）*b
实验十
x=sym（'
6'
y=sym（'
5'
z=（x+1）/（sqrt（3+x）-sqrt（y））
2、分解因式
symsxy;
A=x^4-y^4;
factor（A）
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'
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））
（3）
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）
升级会员 