MATLAB_程序设计与应用第二版_刘卫国_实验部分

实验一

1.

（1）

z1=2*sin(85*pi/180)/(1+exp(2))

（2）

>> x=[2,1+2i;-0.45,5];

>> z2=(1/2)*log(x+sqrt(1+x*x))

(3)

>> a=-3.0:0.1:3.0;

>> z3=((exp(0.3*a)-exp(-0.3*a))/2).*sin(a+0.3)+log((0.3+a)/2)

(4)

t=0:0.5:2.5;

>>

z4=(t.*t).*(t>=0&t<1)+(t.*t-1).*(t>=1&t<2)+(t.*t-2.*t+1).*(t>=2&t<3) 2.

(1)

>> A=[12,34,-4;34,7,87;3,65,7];

>> B=[1,3,-1;2,0,3;3,-2,7];

>> A+6*B

(2)

>> A*B

(3)

>> A^3

(4)

>> A/B

(5)

>> [A,B]

3.

(1)

>>

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

(2)

>> D=C(3:end,2:end)

(3)

Name value class

A <5×5double> double

B <5×3double> double

C <5×3double> double

D [520 397;705 557;…] double

4.

(1)

>> X=100:999;

>> length(find(~rem(X,21)))

(2)

>> y='ahHGJjskjJjsdsaSa';

>> M=y(find(y<='z'&y>='a'))

实验二

实验三1 x=input('请输入x的值')

if x<0&x~=-3

y=x^2+x-6;

elseif x>=0&x<5&x~=2&x~=3

y=x^2-5*x+6;

else

y=x^2-x-1;

end

y

2

（1）

x=input('请输入成绩')

if x>=90&x<=100

disp('等级为A')

elseif x>=80&x<=89

disp('等级为B')

elseif x>=70&x<=79

disp('等级为C')

elseif x>=60&x<=69

disp('等级为D')

elseif x>=0&x<=59

disp('等级为E')

else

disp('成绩有误')

end

(2)

x=input('请输入成绩')

switch(x)

case{90,100}

disp('等级为A')

case{80,89}

disp('等级为B')

case{70,79}

disp('等级为C')

case{60,69}

disp('等级为D')

case{0,59}

disp('等级为E')

otherwise

disp('成绩错误')

end

>>

3

Num=input('输入员工的工号:');

Time=input('该工号员工的工作时数:'); salary=input('该员工底薪:');

if Time>=120

salary=salary*(1+0.15);

elseif Time<=60

salary=salary-700;

elseif Time>60&Time<120

salary=Time*84;

end

salary=salary

实验四1：

程序：n=input('输入n:');

for a=1:n;

b(a)=1/(a^2);

end

pi=sqrt(6*sum(b))

2：.

程序： sum=0;

for n=1:200;

b(n)=1/(2*n-1);

sum=sum+b(n);

if(sum>=3)

break;

end

end

n-1

4.

程序：f(1)=1;max=0;

f(2)=0;min=0;

f(3)=1;c=0;

zhengshu=0;

fushu=0;

zero=0;

for n=4:100

f(n)=f(n-1)-2*f(n-2)+f(n-3);

end

for m=1:100

if(f(m)>max)

max=f(m);

end

if(f(m)

min=f(m);

end

c=sum(f);

if(f>0)

zhengshu=zhengshu+1;

end

if(f(m)<0)

fushu=fushu+1;

end

if(f(m)==0)

zero=zero+1;

end

if(f(m)>0)

zhengshu=zhengshu+1;

end

if(f(m)<0)

fushu=fushu+1;

end

if(f(m)==0)

zero=zero+1;

end

disp('最大值')

max

disp('最小值')

min

disp('和')

c

disp('正数的个数')

zhengshu

disp('负数的个数')

fushu

disp('零的个数')

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

x=input('请输入复数实部的值');

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

函数调用：

程序：

clear

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:0.02*pi:2*pi)';

y1=x.*x;

y2=cos(2*x);

y3=y1.*y2;

plot(x,y1,'r-'),hold on

plot(x,y2,'b:'),hold on

plot(x,y3,'gp'),hold off

（2）

程序：

x=(0:0.02*pi:2*pi)';

y1=x.*x;

y2=cos(2*x);

y3=y1.*y2;

subplot(2,2,1)

stairs(x,y1);

plot(x,y1,'r-'),hold on

subplot(2,2,2)

stairs(x,y2);

plot(x,y2,'b:'),hold on

subplot(2,2,3)

stairs(x,y3);

plot(x,y3,'gp'),hold off

（3）

（1）绘制条形图

程序：

x=(0:0.2*pi:2*pi)';

y1=x.*x;

y2=cos(2*x);

y3=y1.*y2;

subplot(2,2,1)

stairs(x,y1);

bar(x,y1)

subplot(2,2,2)

stairs(x,y2);

bar(x,y1)

subplot(2,2,3)

stairs(x,y3);

bar(x,y3)

（2）绘制梯形图

程序：

x=(0:0.05*pi:2*pi)';

y1=x.*x;

y2=cos(2*x);

y3=y1.*y2;

subplot(2,2,1)

stairs(x,y1);

subplot(2,2,2)

stairs(x,y2);

subplot(2,2,3)

stairs(x,y3);

（4）杆图：

程序：

x=(0:0.05*pi:2*pi)';

y1=x.*x;

y2=cos(2*x);

y3=y1.*y2;

subplot(2,2,1)

stem(x,y1);

subplot(2,2,2)

stem(x,y2);

subplot(2,2,3)

stem(x,y3);

（5）填充图

程序：

x=(0:0.05*pi:2*pi)';

y1=x.*x;

y2=cos(2*x);

y3=y1.*y2;

subplot(2,2,1)

fill(x,y1,'r');

subplot(2,2,2)

fill(x,y2,'b');

subplot(2,2,3)

fill(x,y3,'g');

3、：

程序：

x=linspace(-5,5,1000);

y=[];

for x0=x

if x0<=0

y=[y,(x0+sqrt(pi))./exp(2)];

else

y=[y,0.5*log(x0+sqrt(1+x0.*x0))];

end

end

plot(x,y)

4、

程序：

a=input('请输入a的值')

b=input('请输入b的值')

n=input('请输入n的值')

theta=0:0.01:2*pi;

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);

title('surfc(x,y,z)');

surfc(x,y,z);

实验八

1、

程序：

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 %满足条件的百分比

2、

程序：

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:-1:1);

实验九

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

实验十

1、

程序：

x=sym('6');

y=sym('5');

z=(x+1)/(sqrt(3+x)-sqrt(y))

2、分解因式

（1）

程序：

syms x y;

A=x^4-y^4;

factor(A)

（2）

程序：

factor(sym('5135'))

3、化简表达式

（1）

程序：

syms beta1 beta2

y=sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2)

simple(y)

（2）

程序：

syms x

y=(4*x^2+8*x+3)/(2*x+1)

simple(y)

5、用符号方法求下列极限或导数

（1）

程序：

syms x

f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(sin(x)) limit(f)

（2）

程序：

syms x

y=(sqrt(pi)-sqrt(acos(x)))/(sqrt(x+1));

limit(f,x,-1,'right')

（3）

程序：

syms x

y=(1-cos(2*x))/x;

y1=diff(y)

y2=diff(y,x,2)

6、用符号方法求下列积分

（1）

程序：

syms x

f=1/(1+x^4+x^8)

int(f)

（2）

程序：

syms x

f=1/(((asin(x))^2)*sqrt(1-x^2))

int(f)

（3）

程序：

syms x

f=(x^2+1)/(x^4+1)

int(f,x,0,inf)

（4）

程序：

syms x

f=exp(x)*(1+exp(x))^2

y=int(f,x,0,log(2))

double(y)