文档库 最新最全的文档下载
当前位置:文档库 › 2015下MATLAB语言上机练习内容

2015下MATLAB语言上机练习内容

2015下MATLAB语言上机练习内容
2015下MATLAB语言上机练习内容

MATLAB 语言上机实验内容

注意:每个实验奇数题做实验,偶数题做作业(做实验时一同带往实验室作为实验预习题提交)

实验1: MATLAB 基础与符号计算

1、指令>>A= [2 2 ; 2 2];

>>A .^2= 。 >>A^2= 。 A=[2 2 ;2 2 ];

A.^2 A^2 ans = 4 4 4 4 ans = 8 8 8 8

2、数组>>x=10:-2:2; y=1:5;

>>(x<5)&(y>2),ans= 。 >> y==2|x~=2, ans= 。

x=10:-2:2;y=1:5; (x<5)&(y>2) y==2|x~=2 ans =

0 0 0 1 1 ans =

1 1 1 1 0

3、(1)对于0>x ,求1

2011122+∞

=∑?

?

?

??+-+k k x x k 。(提示:理论结果为x ln )

syms x k

f=(2/(2*k+1))*((x-1)/(x+1))^(2*k+1);

symsum(f,k,0,inf) ans =

piecewise([abs(x + 1)*(abs(x - 1) - abs(x + 1)) < 0, 2*atanh((x - 1)/(x + 1))])

(2)求∑∞=0k k

x 的符号解,并进而用该符号解求∑∞

=-0

)31(k k ,∑∞=0)

1(k k

π,∑∞

=0

3k k

的准

确值。 (提示:注意subs 的使用)

s1=symsum(x^k,k,0,inf) s2=symsum((-1/3)^k,k,0,inf) s3=symsum((1/pi)^k,k,0,inf) s4=symsum(3^k,k,0,inf) s1 =

piecewise([1 <= x, Inf], [abs(x) < 1, -1/(x - 1)]) s2 = 3/4 s3 =

18014398509481984/12280237370259325 s4 = Inf

4、(1)通过符号计算求t t y sin )(=的导数

dt

dy

。 syms t

f=abs(sin(t)) s1=diff(f,t)

s2=limit(s1,t,0,'left') s3=limit(s1,t,pi/2) f =

abs(sin(t)) s1 =

abs(1,sin(t))*cos(t) s2 = -1 s3 = 0

(2)然后根据此结果,求

-

=0

t dt

dy 和

2

π

=t dt dy 。

5、指令>>A= [2 1 ; 1 3];

>>A .^2-A= 。 >>A^2+A= 。

A=[2 1;1 3]; A.^2-A ans = 2 0 0 6 A^2+A ans = 7 6 6 13

6、数组>>x=10:-2:2; y=5:-1:1;

>>(x>=4)&(y>=2),ans= 。 >> y==3|x~=3, ans= 。

x=10:-2:2;y=5:-1:1; (x>=4)&(y>=2) y==3|x~=3 ans =

1 1 1 1 1 ans =

1 1 1 1 0

7计算二重积分

??

+2

1

1

222

)(x dydx y x 。

syms x y

f=int(int((x^2+y^2),y,1,x^2),x,1,2)

f = 1006/105

8在]2,0[π区间,画出dt t

t

x y x

?

=

sin )(曲线,并计算)5.4(y 。 syms t x

y=int(sin(t)/t,t,0,x) y=sinint(x)

ezplot(y,[0,2*pi]) f=subs(y,x,4.5) y = sinint(x) y = sinint(x) f = 1.6541

1

2

34

5

6

00.20.40.60.811.21.41.61.8x

sinint(x)

9设系统的冲激响应为t

e t h 3)(-=,求该系统在输入t t u cos )(=,0≥t 作用下的输出。(提

示:运用卷积进行计算)

syms t tao s ht=exp(-3*t); ut=cos(t);

hu=subs(ut,t,tao)*subs(ht,t,t-tao); y1=int(hu,tao,0,t) y1 =

(3*cos(t))/10 - (3*exp(-3*t))/10 + sin(t)/10

10、求f={A (1-|t|/τ)|t|≦τ 的 Fourier 变换,并画出 A= 2,τ=2时的幅频谱。 0 |t|>τ

syms A t w;

syms tao positive;

ht=heaviside(t)-heaviside(t-tao);htt=heaviside(t+tao)-heaviside(t); ft=A*(1-t/tao)*ht+A*(1+t/tao)*htt; fw=fourier(ft,t,w)

fw2=subs(fw,[A tao],[2,2]) ezplot(fw2,[-10 10]) fw =

4*A/tao/w^2*sin(1/2*tao*w)^2 fw2 =

4/w^2*sin(w)^2

-10

-8

-6

-4

-2

02

4

6

8

10

00.511.522.533.54w

4/w 2 sin(w)2

11、方程求方程的解

初始条件为,,4)0(r 3,r(0),423'2'''===++-t e r r r

syms t

y=dsolve('D2r+3*Dr+2*r=4*exp(-2*t)','r(0)=3,Dr(0)=4','t') y =

14*exp(-t) - 11*exp(-2*t) - 4*t*exp(-2*t)

12、求

T k ke k f )(λ-=的Z 变换表达式。

.syms k z r T f=k*exp(-r*k*T) fZ=ztrans(f,k,z); f=k*exp(-r*k*T) pretty(fZ) f =

k*exp(-r*k*T) f =

k*exp(-r*k*T)

z exp(-r T) ---------------- 2 (z - exp(-r T))

13、求方程2,12

2

==+xy y x 的解。

[x,y]=solve('x^2+y^2=1','x*y=2','x','y') x =

((15^(1/2)*i)/2 + 1/2)^(1/2)/2 - ((15^(1/2)*i)/2 + 1/2)^(3/2)/2 - ((15^(1/2)*i)/2 + 1/2)^(1/2)/2 + ((15^(1/2)*i)/2 + 1/2)^(3/2)/2 (1/2 - (15^(1/2)*i)/2)^(1/2)/2 - (1/2 - (15^(1/2)*i)/2)^(3/2)/2 - (1/2 - (15^(1/2)*i)/2)^(1/2)/2 + (1/2 - (15^(1/2)*i)/2)^(3/2)/2 y =

((15^(1/2)*i)/2 + 1/2)^(1/2) -((15^(1/2)*i)/2 + 1/2)^(1/2) (1/2 - (15^(1/2)*i)/2)^(1/2) -(1/2 - (15^(1/2)*i)/2)^(1/2)

14求边值问题

1)0(,0)0(,34,43==+-=+=g f g f dx

dg g f dx df 的解。 syms x;

[f,g]=dsolve('Df=3*f+4*g,Dg=-4*f+3*g','f(0)=0,g(0)=1','x') f =

exp(3*x)*sin(4*x) g =

exp(3*x)*cos(4*x)

15求4

633

)(2

3++++=

s s s s s F 的Laplace 反变换。

syms s t

f=(s+3)/(s^3+3*s^2+6*s+4); ft=ilaplace(f,s,t)

ft =

(2*exp(-t))/3

-

(2*exp(-t)*(cos(3^(1/2)*t)

-

(3^(1/2)*sin(3^(1/2)*t))/2))/3

实验2: 数组(矩阵)及数值计算

1、指令>>A=ones(4,4); >>B=A(3,[2 3]); >>C=A(3,:)<1; 写出输出结果:

>>B=_________________________。 >>C=_________________________。

A=ones(4,4);

B=A(3,[2 3])

C=A(3,:)<1

B =

1 1

C =

0 0 0 0

2、指令>>B=magic(3);写出求矩阵B的特征值D和特征向量V的语句:

B=magic(3);

[V,D]=eig(B)

V =

0.0846 0.4928 0.7071 0.5000

0.0846 0.4928 -0.7071 0.5000

-0.7815 -0.3732 0 0.5000

0.6124 -0.6124 0 0.5000

D =

-0.0000 0 0 0

0 -0.0000 0 0

0 0 0 0

0 0 0 4.0000

3、>>A=rand(100,100);用向量化运算的方法,用1条语句实现求矩阵A中大于0.4的所有

元素个数:

> A=rand(100,100);

A1=sum(sum(A>0.4))

A1 =

5984

2,0[ 上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。4要求在闭区间]

0:2*pi/9:2*pi

ans =

0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832

x=linspace(0,2*pi,10)

x =

0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832

5由指令rand('state',0),A=rand(3,5)生成二维数组A,试求该数组中所有大于0.5的元

素的位置,分别求出它们的“全下标”和“单下标”。

rand('state',0)

A=rand(3,5)

si=find(A>0.5)

[ri,ci]=ind2sub(size(A),si)

[rj,cj]=find(A)

A =

0.9501 0.4860 0.4565 0.4447 0.9218

0.2311 0.8913 0.0185 0.6154 0.7382

0.6068 0.7621 0.8214 0.7919 0.1763

si =

1

3

5

6

9

11

12

13

14

ri =

1

3

2

3

2

3

1

2 ci =

1

1

2

2

3

4

4

5 5

rj =

1

2

3

1

2

3

1

2

3

1

2

3

1

2

3 cj =

1 1 1

2 2

3 3 3

4 4 4

5 5 5

6在时间区间 [0,10]中,绘制t e y t 2cos 15.0--=曲线。要求分别采取“标量循环运算法”和“数组运算法”编写两段程序绘图。

t=linspace(0,10,200); N=length(t); for k=1:N

y1(k)=1-exp(-0.5*t(k))*cos(2*t(k)); end

subplot(1,2,1),plot(t,y1),xlabel('t'),ylabel('y1'),grid on y2=1-exp(-0.5*t).*cos(2*t);

subplot(1,2,2),plot(t,y2),xlabel('t'),ylabel('y2'),grid on

7先运行指令x=-3*pi:pi/15:3*pi; y=x; [X,Y]=meshgrid(x,y); warning off; Z=sin(X).*sin(Y)./X./Y; 产生矩阵Z 。(1)请问矩阵Z 中有多少个“非数”数据?(2)用指令surf(X,Y,Z); shading interp 观察所绘的图形。(3)请写出绘制相应的“无裂缝”图形的全部指令。

clear,

x=-3*pi:pi/15:3*pi;

y=x;

[X,Y]=meshgrid(x,y); warning off;

Z=sin(X).*sin(Y)./X./Y;

A=sum(sum(isnan(Z)))

surf(X,Y,Z); shading interp A =

181

x=-3*pi:pi/15:3*pi;y=x; [X,Y]=meshgrid(x,y);

X=X+(X==0)*eps;

Y=Y+(Y==0)*eps;

Z=sin(X).*sin(Y)./X./Y;

surf(X,Y,Z);shading interp

8、指令>>A=3*ones(5,5);

>>B=A(3,[1 2]);

>>C=A(:,3)>1;

写出输出结果:

>>B=_________________________。

>>C=_________________________。

A=3*ones(5,5);

B=A(3,[1 2])

C=A(:,3)>1

B =

3 3

C =

1

1

1

1

1

9、指令>>B=magic(4), 按要求分别写出指令:

求矩阵B的行列式:_______________________。

求矩阵B的逆矩阵:_______________________。

求矩阵B的迹:。

求矩阵B的特征值D和特征向量V:。

B=magic(4);

A1=det(B) A2=inv(B) A3=trace(B) [V,D]=eig(B) A1 = 0 A2 =

1.0e+014 *

0.9382 2.8147 -2.8147 -0.9382 2.8147 8.4442 -8.4442 -2.8147 -2.8147 -8.4442 8.4442 2.8147 -0.9382 -2.8147 2.8147 0.9382 A3 = 34 V =

-0.5000 -0.8236 0.3764 -0.2236 -0.5000 0.4236 0.0236 -0.6708 -0.5000 0.0236 0.4236 0.6708 -0.5000 0.3764 -0.8236 0.2236 D =

34.0000 0 0 0 0 8.9443 0 0 0 0 -8.9443 0 0 0 0 0.0000

10采用数值计算方法,画出dt t

t

x y x

?

=

sin )(在]10 ,0[区间曲线,并计算)5.4(y 。 d=0.5; tt=0:d:10;

t=tt+(tt==0)*eps; y=sin(t)./t; s=d*trapz(y) ; ss=d*(cumtrapz(y)); plot(t,y,t,ss,'r'),hold on y4_5=ss(find(t==4.5))

yi=interp1(t,ss,4.5), plot(4.5,yi,'r+'); yy=quad('sin(t)./t',0,4.5); yy=quadl('sin(t)./t',0,4.5); warning off y4_5 =

1.654178716404084 yi =

1.654178716404084 >>

11用quad 求取

dx x e

x

sin 7.15?

--π

π

的数值积分,并保证积分的绝对精度为910-。

format long d=pi/8;

t=-5*pi:d:1.7*pi;

s1=quad('exp(-abs(x)).*abs(sin(x))',-5*pi,1.7*pi,1e-9)

s1 =

1.087849437547788

12求函数5.08.12cos 5.1)5(sin )(2

06.02++-=t t t e

t t f t

在区间]5,5[-中的最小值点。 t1=-5; t2=5;

ft=@(t)(sin(5*t))^2*exp(0.06*t^2)-1.5*t*cos(2*t)+1.8*abs(t+0.5); [t,fval,exitflag,output]=fminbnd(ft,t1,t2) t =

-1.284977032658338 fval =

-0.186048010633265

exitflag = 1

output =

iterations: 12 funcCount: 13

algorithm: 'golden section search, parabolic interpolation' message: [1x112 char]

13设0)

0(,1)0(,1)(2)(3)(2

2===+-dt

dy y t y dt t dy dt t y d ,用数值法和符号法求5.0)(=t t y 。 format long tspan=[0,1]; y0=[1;0];

dydt=@(t,y)[y(2);-2*y(1)+3*y(2)+1]; [tt,yy]=ode45(dydt,tspan,y0);

y0_5=interp1(tt,yy(:,1),0.5,'spline') y0_5 =

0.789580207901271

S = dsolve('D2y-3*Dy+2*y = 1','y(0) = 1','Dy(0) = 0') ys0_5=subs(S,0.5) S =

1/2-1/2*exp(2*t)+exp(t) ys0_5 =

0.7896

14求矩阵Ax b 的解,A 为4 阶魔方阵,b 是(4×1)的全 1 列向量。

A=magic(4); b=ones(4,1); x=A\b; xg=null(A) xg = 0.2236 0.6708 -0.6708 -0.2236

15求矩阵b Ax =的解,A 为4阶魔方阵,????

?

???????=4321b 。

clear all A=magic(4); b=(1:4)'; rref([A,b]) ans =

1 0 0 1 0 0 1 0 3 0 0 0 1 -3 0

0 0 0 0 1

16 求解二元函数方程组??

?=+=-0

)cos(0

)sin(y x y x 的解。(用符号法解)

S=solve('sin(x-y)=0','cos(x+y)=0','x','y') S.x, S.y S =

x: [2x1 sym] y: [2x1 sym] ans = -1/4*pi 1/4*pi ans = -1/4*pi 1/4*pi

>>

实验3: 图形绘制及m 文件编程

1. 已知椭圆的长、短轴2,4==b a ,用“小红点线”画椭圆??

?==t

b y t

a x sin cos 。

解:t=0:pi/20:2*pi; x=4*cos(t);

y=2*sin(t);

plot(x,y,'r.',’marksize ’,10) axis equal

2. 根据表达式θ

ρcos 1-=绘制如图p5-2的心脏线。(提示:采用极坐标绘线指令polar ) 解:theta=0:pi/50:2*pi;

rho=1-cos(theta); h=polar(theta,rho,'r-')

做此类题目注意把握图形的信息,看图标以及图形的题目(图标用legend,题目用title ,并注意里面需用单引号)

3. A,B,C 三个城市上半年每个月的国民生产总值见表p5.1。试画出如图p5-3所示的三城

市上半年每月生产总值的累计直方图。(提示:使用指令bar )

表p5.1 各城市生产总值数据(单位:亿元) 城市 1月 2月 3月 4月 5月 6月 A 170 120 180 200 190 220 B 120 100 110 180 170 180 C 70

50

80

100

95

120

解:x=(1:6);

y=[170,120,180,200,190,220;120,100,110,180,170,180;70,50,80,100,95,120]'; bar(x,y,'stacked'); colormap(cool); legend('a','b','c'); axis([0,7,0,600]) shg

做这类题注意数字的对列,即将数字以列的形式一一对

4. 用绿实线绘制)sin(t x =,)cos(t y =,t z =的三维曲线,曲线如图p5-5所示。

(提示:使用plot3指令)

-1

-0.5

0.5

1

-1

1

5

10

15

解:t=0:pi/20:4*pi;

x=sin(t); y=cos(t); z=t;

plot3(x,y,z,'g-') box on

做这题需注意的是图形窗的开与闭,若闭合用box on ,若开用box off

5. 在区域]3,3[,-∈y x ,绘制2

24y

x xe z --=的如图p5-6的三维(透视)网格曲面。(不得

使用ezmesh )

解:x=-3:3;

y=x;

[X,Y]=meshgrid(x,y); z=4.*X.*exp(-X.^2-Y.^2); mesh(X,Y,z)

hold on,mesh(X,Y,z),colormap(hot) hold off hidden off

axis equal,axis off

做这题需注意图形并茂,网格面用mesh,光滑面用surf,必要时注意图形的坐标范

围,用axis ,括号里用方括号

6. 在]

4,4[,ππ-∈y x 区间里,根据表达式y

x y x z ++=)

sin(,绘制如图p5-7所示的曲面。

解:x=-4*pi:pi/20:4*pi;

y=x;

[X,Y]=meshgrid(x,y);

z=sin(X+Y)./(X+Y+(X+Y==0)*eps); surf(X,Y,z); shading interp view(30,45)

若此类图形的定义域不能为零,注意机器零eps 的使用 7. 试用图解法回答:

(1)方程组???

??

=+=++0

))cos(sin(1.0)1(22y x y x y 有多少个实数解?

(2)求出离0

,0==y x 最近、且满足该方程组的一个近似解。 解:(1)ezplot('y/(1+x^2+y^2)-0.1',[-2*pi,2*pi,-pi/3,3.5*pi])

hold on

ezplot('sin(x+cos(y))',[-2*pi,2*pi,-pi/3,3.5*pi]) 可看到6个交点,及方程组有6个解。 (2)syms x y;

s=solve('y/(1+x^2+y^2)-0.1','sin(x+cos(y))')

8.一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。

解:syms x

8.for m=100:999 x=fix(m/100); y=rem(fix(m/10),10); z=rem(m,10);

if m==x*x*x+y*y*y+z*z*z disp(m) end end

9.Fibonacci 数组的元素满足Fibonacci 规则:1

2+++=k k k a a a ,),2,1( =k ;且121=

=a a 。现要求该数组中第一个大于10000的元素。

Matlab上机_测试题c答案讲课教案

<<<<<<精品资料》》》》》 Matlab 上机测试题答案 班级 姓名 学号 要求: (1)请将程序和运行结果粘贴到题目下面,写清楚哪部分是程序,哪部分是运行结果,运行结果包括图形,图形请适当缩小,然后将word 文档传到教育在线,注意保存格式是2003的doc ,文件名”自111_110110_张三” (2)在最开始前请将输入Matlab 命令窗口中输入”diary ”命令,随后在workspace 窗口中可以发现一个名为“diary ”的文件。在老师要求结束练习后,在命令窗口输入”diary off ”命令,然后将文件“diary ”上传到教育在线。 (3)注意题号不要弄乱。 --------------------------------------------------------------------- 1、 矩阵Y= ???? ? ???? ???3472123100451150425 ,给出元素1的全下标和单下标,并用函数练习全下标和单下标的转换,求出元素100的存储位置。取出子矩阵? ? ? ? ??21301,并求该矩阵的维数。 解:命令为: Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10) sub2ind([4 3],2,3) [i,j]=ind2sub([4 3],10) find(Y==100) sub2ind([4 3],3,2) B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y) --------------------------------------------------------------------- 2、已知矩阵A=[1 0 -1 ;2 4 1; -2 0 5],B=[0 -1 0;2 1 3;1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ,A/B 、A\B 解:命令为: A=[1 0 -1 ;2 4 1; -2 0 5] B=[0 -1 0;2 1 3;1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B --------------------------------------------------------------------- 3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵,并计算两者的乘积。 解:命令为: A=eye(3,4) B=8*ones(4) C=A*B --------------------------------------------------------------------- 4、创建矩阵a=????? ???????------7023021.5003.120498601 ,取出其前两列构成的矩阵b ,取出前两行构成矩阵c ,转置矩阵b 构成矩阵d ,计算a*b 、c

matlab上机练习(附答案)

1.以下两种说法对吗? (1)MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。 (2)MATLAB指令窗中显示的数据有效位数不超过七位。 2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同? DIARY filename causes a copy of all subsequent command window input and most of the resulting command window output to be appended to the named file. If no file is specified, the file 'diary' is used. DIARY OFF suspends it. DIARY ON turns it back on. DIARY, by itself, toggles the diary state. Use the functional form of DIARY, such as DIARY('file'), when the file name is stored in a string. 3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上?“位于搜索路 径上的目录”与“当前目录”在MATLAB中的功用相同吗? 4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组?用直 接键入法?用数组编辑器?用M文件编辑器? 5.运用数组算术运算符去掉下面程序里的for/end循环: x=11:15 for k=1:length(x) z(k)=x(k)^2+2.3*x(k)^0.5; end x=11:15 z1=x.^2+2.3*x.^0.5 6.不使用数组算术运算符,重写下面的程序代码: x=[2 1 4] z=1./(1+x.^2) x=2;k=1; while i<=4, z2(k)=1/(1+i^2); i=i+1;x=x+1; end 7.某公司销售电脑打印机的价格方案如下: ()如果顾客只买一台打印机,则一台的基本价格为$150。 ()如果顾客购买两台以上打印机,则第二台价格为$120。 ()第三台以后,每台$110。

matlab基础练习题(带答案)

Matlab基础练习题 常量、变量、表达式 1、MATLAB中,下面哪些变量名是合法的?() (A)_num(B)num_(C)num-(D)-num 2、在MATLAB中,要给出一个复数z的模,应该使用()函 数。 (A)mod(z) (B)abs(z) (C)double(z) (D)angle(z) 3、下面属于MATLAB的预定义特殊变量的是?() (A)eps(B)none (C)zero(D)exp 4、判断:在MATLAB的内存工作区中,存放一个英文字符'a'需要 占用1个字节,存放一个中文字符‘啊’需要占用2个字节。(错,都是2个字节) 5、判断:MATLAB中,i和j都是虚数单位,它们之间没有什么 区别。(对) 6、判断:MATLAB中,pi代表圆周率,它等于3.14。(错,后面 还有很多位小数) 7、在MATLAB中,若想计算的值,那么应该在MATLAB 的指令窗中输入的MATLAB指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。

8、在MATLAB中,a = 1,b = i,则a占_8__个字节,b占_16_个 字节,c占________字节。 9、在MATLAB中,inf的含义是__无穷大__,nan的含义是__非数 (结果不定)___。 数组 1、在MATLAB中,X是一个一维数值数组,现在要把数组X中的所 有元素按原来次序的逆序排列输出,应该使用下面的()指令。 (A)X[end:1] (B)X[end:-1:1](C)X(end:-1:1)(D)X(end:1) 2、在MATLAB中,A是一个字二维数组,要获取A的行数和列数, 应该使用的MATLAB的命令是()。 (A)class(A)(B)sizeof(A)(C)size(A)(D)isa(A) 3、在MATLAB中,用指令x=1:9生成数组x。现在要把x数组的 第二和第七个元素都赋值为0,应该在指令窗中输入()(A)x([2 7])=(0 0) (B)x([2,7])=[0,0] (C)x[(2,7)]=[0 0] (D)x[(2 7)]=(0 0) 4、在MATLAB中,依次执行以下指令:clear;A=ones(3,4); A(:)=[-6:5];这时,若在指令窗中输入指令b=A(:,2)',那么,MATLAB输出的结果应该是() (A)b=-3-2-1(B)b=-2 -1 0 1 (C)b=-5 -1 3 (D)b = -5-214

matlab基础练习题带答案

Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中,下面哪些变量名是合法的?( ) (A )_num (B )num_ (C )num- (D )-num 2、 在MA TLAB 中,要给出一个复数z 的模,应该使用( )函数。 (A )mod(z) (B )abs(z) (C )double(z) (D )angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是?( ) (A )eps (B )none (C )zero (D )exp 4、 判断:在MA TLAB 的内存工作区中,存放一个英文字符 'a' 需要占用1个字节,存放 一个中文字符‘啊’需要占用2个字节。( 错,都是2个字节 ) 5、 判断:MA TLAB 中,i 和j ( 对 ) 6、 判断:MA TLAB 中,pi 代表圆周率,它等于3.14。( 错,后面还有很多位小数 ) 7、 在MA TLAB 中,若想计算的5 1)3.0sin(21+= πy 值,那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中,a = 1,b = i ,则a 占_8__个字节,b 占_16_个字节,c 占________字 节。 9、 在MA TLAB 中,inf 的含义是__无穷大__,nan 的含义是__非数(结果不定)___。 数组 1、 在MA TLAB 中,X 是一个一维数值数组,现在要把数组X 中的所有元素按原来次序 的逆序排列输出,应该使用下面的( )指令。 (A )X[end:1] (B )X[end:-1:1] (C )X (end:-1:1) (D )X(end:1) 2、 在MA TLAB 中,A 是一个字二维数组,要获取A 的行数和列数,应该使用的MATLAB 的命令是( )。 (A )class(A) (B )sizeof(A) (C )size(A) (D )isa(A) 3、 在MATLAB 中,用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0,应该在指令窗中输入( ) (A )x([2 7])=(0 0) (B )x([2,7])=[0,0] (C )x[(2,7)]=[0 0] (D )x[(2 7)]=(0 0) 4、 在MA TLAB 中,依次执行以下指令:clear;A=ones(3,4); A(:)=[-6:5];这时, 若在指令窗中输入指令b=A(:,2)',那么,MATLAB 输出的结果应该是( ) (A )b = -3 -2 -1 (B )b = -2 -1 0 1 (C )b = -5 -1 3 (D )b = -5 -2 1 4 5、 在MA TLAB 中,A = 1:9,现在执行如下指令L1 = ~(A>5),则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。

MATLAB上机习题一

MATLAB上机习题一 请按以下步骤完成上机实验: 1)在FTP上下载“MATLAB上机习题一.doc”文件,所有习题列在该文件内; 2)在MATLAB中完成所有习题,并将屏幕截图粘贴到相应习题后面; 3)如果习题是问答题,请将答案写在题目后; 4)如果有的习题要求提供脚本文件,请将脚本文件内容拷贝到相应习题后; 5)将文件保存并重命名为“自己的学号-姓名”,例如“20110771-张三.doc”; 6)上传该文件到FTP的相关目录。 1. 运行MATLAB软件,拖放、关闭界面上的子窗口,并恢复到原始试图。 2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。 3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。 4. 通过使用帮助确定内置变量ispc 的含义。 6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。 5. 创建double类型的变量,并进行计算 1)a=87,b=190,计算a+b、a-b、a*b;

2)创建uint8类型的变量,数值与(1)中相同,进行相同的计算,观察计算结果与预想的是否一致。 6. 计算如下表达式: 1)()sin 60 2)3e 3)3cos 4π?? ??? 4)2 7562323336 +?-?? 7. 设u=2,v=3,计算: 1)4log uv v 2)() 2 2u e v v u +- 3 8. 计算如下表达式: 1)()()3542i i -+ 2)()sin 28i - 9. 判断下面语句的运算结果,并与MATLAB 计算结果做比较:

1)420< 2)420≤ 3)420== 4)4~20= 5)''''b B < 10. 设a=39,b=58,c=3,d=7,判断下面表达式的值与MATLAB 计算结果做比较: 1)a b > 2)a b < 3)&&a b b c >> 4)a d == 5)|a b c > 6)~~d 11. 判断下列变量名哪些是合法的MATLAB 变量: 1)fred 7)fred! 2)book_1 8)book-1 3)2ndplace 9)Second_Place 4)#1 10)No_1 5)vel_5 11)vel.5 6)tan 12)while 提示,可以使用isvarname 对上述变量名进行检验。 12. 编写脚本程序,命名为barbell.m ,完成如下计算: 如图1所示的杠铃,

MATLAB上机考试题(一)

(1)在MATLAB的命令窗口中执行_____命令,将命令窗口的显示内容清空。() A.clear B.clc C.echo off D.cd (2)在MATLAB的命令窗口中执行_____命令,使数据输出显示为十六进制表示。() A.format long B.format rat C.format hex D.format short e (3)下列变量名中_____是合法的。() A.x*y,a,1 B.x\y,a1234 C.end,1 bcx D.char_1,i,j (4)已知x=0:5,则x有_____个元素。() A.5 B.6 C.7 D.8 (5)一下运算符中哪个的优先级最高_____。() A./ B.^ C.~= D.& (6)使用检测函数isnumeric(10)的结果是_____。() A.1 B.0 C.false D.true (7)三维图形中默认视角是_____度。() A.方位角=0 俯仰角=90 B.方位角=90 俯仰角=0 C.方位角=37.5 仰俯角=30 D.方位角=0 仰俯角=180 (8)将符号表达式化简为因式分解因式分解因式分解因式分解形式,使用_____函数。() A.collect B.expand C.horner D.factor (9)运行以下命令,则_____描述是正确的。()>>syms a b c d >>A=[a b;c d] A.A占用的内存小于100B B.创建了5个符号变量 C.A占用的内存是a b c d的总和 D.不存在 (10)已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。 (11)运行命令bitor(8,7)的结果是_____。 (12)运行以下命令: >>x=0:10; >>y1=sin(x); >>y2=5*sin(x); >>y3=[10*sin(x );20*sin(x)]; >>plot(x,y1,x,y2,x,y3) 则在一个图形窗口中,可以看到_____条曲线。 (13)符号表达式“g=sym(sin(a*z)+cos(w*v))”中的自由符号变量是_____。 (14)运行以下命令: >>syms t >>f1=1/t >>limitf1_r=limit(f1,'t','0','right'); 则函数limitf1_r趋向0的右极限为_____。 15.在MATLAB的命令窗口中执行______命令,使数值5.3显示为5.300000000000000e+000 A. format long B. format long e C. format short D. format short e 16.下列变量名中______是合法的。A.char_1,i,j B.1_1, a.1 C.x\y,a1234 D.end,1bcx 17.已知x=0:9,则x有_____个元素。 A.12 B.11 C.10 D.9 18.产生对角线上为全1其余为0的2行3列矩阵的命令是______ A. ones(2,3) B. ones(3,2) C. eye(2,3) D. eye(3,2) 19.已知数组a= [1 2 3 4 5 6 7 8 9] ,则运行a(:,1)=[]命令后______ A. a变成行向量 B. a数组为2行2列 C. a 数组为3行2列 D. a数组中没有元素3 20.按含义选出各个函数名:表示4舍5入到整数的是____,表示向最接近0取整的是____,表示向最接近-∞取整的是____,表示向最接近∞取整的是_____ A. round(x) B. fix(x) C. floor(x) D. ceil(x) 21.已知a=0:5,b=1:6,下面的运算表达式出错的为______ A. a+b B. a./b C. a’*b D. a*b 22.已知s=’显示”hello”’,则s的元素个数是______ A. 12 B. 9 C. 7 D.18

MATLAB上机实验练习题答案

数学建模 MATLAB上机实验练习题 1、给出一个系数矩阵A[234;541;132],U=[123],求出线性方程组的一个精确解。 2、给出两组数据x=[00.30.81.11.62.3]’y=[0.820.720.630.600.550.50]’,我们可以简单的 认为这组数据在一条衰减的指数函数曲线上,y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式,并且在图形窗口画出这条曲线,已知的点用*表示。 3、解线性方程 4、通过测量得到一组数据: 5、已知一组测量值 6、从某一个过程中通过测量得到: 分别采用多项式和指数函数进行曲线拟合。 7、将一个窗口分成四个子窗口,分别用四种方法做出多峰函数的表面图(原始数据法,临近 插值法,双线性插值法,二重三次方插值法) 8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。分别使用不同的 颜色,线形和标识符。 9、下面的矩阵X表示三种产品五年内的销售额,用函数pie显示每种产品在五年内的销售额

占总销售额的比例,并分离第三种产品的切片。 X=19.322.151.6 34.270.382.4 61.482.990.8 50.554.959.1 29.436.347.0 10、对应时间矢量t,测得一组矢量y t00.30.8 1.1 1.6 2.3 y0.50.82 1.14 1.25 1.35 1.40 采用一个带有线性参数的指数函数进行拟合,y=a0+a1e-t+a2te-t,利用回归方法求出拟合函数,并画出拟合曲线,已知点用圆点表示。 11、请创建如图所示的结构数组(9分) 姓名编号指标 江明顺071023身高:176,体重:82 于越忠060134身高:168,体重:74 邓拓050839身高:182,体重:77 12、创建如图所示的元胞数组。(9分) 13、某钢材厂从1990年到2010年的产量如下表所示,请利用三次样条插值的方法计算1999年该钢材厂的产量,并画出曲线,已知数据用‘*’表示。要求写出达到题目要求的MATLAB 操作过程,不要求计算结果。 年份19901992199419961998200020022004200620082010 产量(万吨)75.99591.972105.711123.203131.669150.697179.323203.212226.505249.633256.344 14、在一次化学动力学实验中,在某温度下乙醇溶液中,两种化合物反应的产物浓度与反应时间关系的原始数据如下,请对这组数据进行三次多项式拟合,并画出拟合曲线,已知数据如下。 time=[2.55.07.510.013.017.020.030.040.050.060.070.0] res=[0.290.560.771.051.361.522.002.272.813.053.253.56]

Matlab 上机题及答案

1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 for m=100:999 m1=fix(m/100); %求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10); %求m的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2.从键盘输入若干个数,当输入0时结束输入,求这些数的平均值和它们之和。 sum=0; n=0; val=input('Enter a number (end in 0):'); while (val~=0) sum=sum+val; n=n+1; val=input('Enter a number (end in 0):'); end if (n > 0) sum mean=sum/n end 3. 若一个数等于它的各个真因子之和,则称该数为完数,如6=1+2+3,所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m); end end 4. 从键盘上输入数字星期,在屏幕上显示对应英文星期的单词。 function week n=input('input the number:'); if isempty(n) errror('please input !!')

end if n>7|n<1 error('n between 1 and 7') end switch n case 1 disp('Monday') case 2 disp('Tuesday') case 3 disp('Wednesday') case 4 disp('Thursday') case 5 disp('Friday') case 6 disp('Saturday') case 7 disp('Sunday') end 5. 某公司销售电脑打印机的价格方案如下: ()如果顾客只买一台打印机,则一台的基本价格为$150。 ()如果顾客购买两台以上打印机,则第二台价格为$120。 ()第三台以后,每台$110。 写一段程序分别计算出购买1--10台打印机所需的钱数。打印机台数可以在程序开始处指定,或通过input命令读入。运行程序,计算出购买10台打印机的总价格。 写出程序,生成分别购买1--10台打印机所需价格的图表(使用fprintf命令输出图表,不允许手算)。 x=input('请输入购买的打印机台数:'); for m=1:x if m<=1 y(m)=150*m; elseif m<=2 y(m)=150+120*(m-1); else y(m)=150+120+110*(m-2); y(1,m)=y(m); end end y(x) plot(1:m,y,'r*--')

MATLAB上机题汇总(基本题,抽3题)合肥学院

1设置matlab 的工作环境,将工作目录设置为d:\work ,添加搜索目录d:\example 2在matlab 的命令窗口里完成如下计算,其中t 的值分别取-1,0,1,表达式如下: 4/3)2 t y e π-= 3自行产生一个5行5列的数组,得到最中间的三行三行矩阵。 4用magic 产生一个5*5的矩阵,将这个矩阵的第二行与第三行互换位置 5求方程组的根 x 1+4x 2-3x 3=2 2x 1+5x 2-x 3=11 x 1+6x 2+x 3=12 6已知:一个多项式的系数向量是p=[1 -6-72 -27],求这个多项式的根。 7已经两个多项式的系数分别是:[1 2 3 4]和[1 4 9 16],请求这两个多项式的乘积,及商和余数。 8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式 9 A=[2 3 4; 1 5 7; 6 2 5] 用什么函数,保证第一列排序的时候,其他列跟着变化。最后的结果是: [1 5 7; 2 3 4; 6 2 5] 10补充题:电路分析 电路分析常常涉及到对方程组的求解,这些方程常常是利用描述电流进入和离开节点的电流方程,或者描述电路中网络回路上的电压的电压方程得出的。下图描述了3个回路电压的方程

方程式如下: 11121222132423432 532()0()()0()0V R I R I I R I I R I R I I R I I R I V -++-=??-++-=??-++=? 假设5个电阻值为已知,2个电压值也为已知,求3个电流值。 11自行产生一个5行5列的数组,用两种方法得到最中间的三行三行矩阵。 12根据a=reshape(-4:4,3,3)做一个矩阵,然后(1)取出所有大于0的元素构成一个向量(可推广到从一个矩阵里筛选出符合条件的元素组成一个向量)(2)将原矩阵中大于0的元素正常显示,而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示,符合条件的值不变)。(思考:大于0的正常显示,小于等于0的用-10来表示) 13建立如下一个元胞数组,现在要求计算第一个元胞第4行第2列加上第二个元胞+第三个元胞里的第二个元素+最后一个元胞的第二个元素。 a={pascal(4),'hello';17.3500,7:2:100} 14建立一个结构体的数组,包括3个人,字段有姓名,年龄,分数,其中分数由随机函数产生一个3行10列的数据表示了有10门课程,每门课程有三个阶段的分数。问题是: 问题1,如何找到第2个人的分数并显示出来 问题2,如何找到第2个人的每门课程3个阶段的平均分数并显示出来 问题3,全班同学(指这3个学生)的10门课程的每门课程的平均分如何计算出来?要求放到一个数组里。 问题4,找到这个班所有同学的姓名放到一个元胞数组里 15给定一个图像文件,格式是jpg ,通过inportdata 引入这个文件,查看引入后数据保存是一个数组,是100*100*3的一个数组,这是一个三维的数组,表明有100*100个点,每个点有RGB 三个方面决定的。现在,要求将这些点的中间部分30个象素宽度的正方形用黑点表示。 16手动构造一个长宽页是1000*1000*3的图像,每个图像的点有三个0-255之间的随机值构成。然后使用image 函数显示这个图像,再用imwrite 函数保存这个图像,使这个图像用图片浏览器也能够查看。 17给定一个矩阵 a = 1 2 3 4 5 6 编写一个M 函数,要求输入是a ,输出有三个:平均数,标准差,秩。程序运行

Matlab上机练习参考答案

Matlab 上机练习二 班级 学号 姓名 按要求完成题目,并写下指令和运行结果。(不需要画图) 1、 求??????+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。 >> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i]; >> x’ ans = - - - + + - - + + - 2、计算??????=572396a 与?? ????=864142b 的数组乘积。 >> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8]; >> a.*b

ans = 12 36 3 8 42 40 3、 对于B AX =,如果??????????=753467294A ,??????????=282637B ,求解X 。 >> A=[4 9 2;7 6 4;3 5 7]; >> B=[37 26 28]’; >> X=A\B X = 4、 ??????-=463521a ,??????-=263478b ,观察a 与b 之间的六种关系运算的结果。 >> a=[1 2 3;4 5 6]; >> b=[8 –7 4;3 6 2]; >> a>b ans = 0 1 0

1 0 1 >> a>=b ans = 0 1 0 1 0 1 >> a> a<=b ans = 1 0 1 0 1 0 >> a==b ans = 0 0 0 0 0 0

>> a~=b ans = 1 1 1 1 1 1 5、[]7.0 = - a,在进行逻辑运算时,a相当于什么样的逻辑量。 8 2.0 5- 相当于a=[1 1 0 1 1]。 6、角度[]60 x,求x的正弦、余弦、正切和余切。 = 30 45 >> x=[30 45 60]; >> x1=x/180*pi; >> sin(x1) ans = >> cos(x1) ans = >> tan(x1) ans =

MATLAB上机答案

一熟悉Matlab工作环境 1、熟悉Matlab的5个基本窗口 思考题: (1)变量如何声明,变量名须遵守什么规则、是否区分大小写。 答:变量一般不需事先对变量的数据类型进行声明,系统会依据变量被赋值的类型自动进行类型识别,也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则: 变量名必须以字母开头,只能由字母、数字或下划线组成。 变量名区分大小写。 变量名不能超过63个字符。 关键字不能作为变量名。 最好不要用特殊常量作为变量名。 (2)试说明分号、逗号、冒号的用法。 分号:分隔不想显示计算结果的各语句;矩阵行与行的分隔符。 逗号:分隔欲显示计算结果的各语句;变量分隔符;矩阵一行中各元素间的分隔符。 冒号:用于生成一维数值数组;表示一维数组的全部元素或多维数组某一维的全部元素。 (3)linspace()称为“线性等分”函数,说明它的用法。 LINSPACE Linearly spaced vector.线性等分函数 LINSPACE(X1,X2)generates a row vector of100linearly equally spaced points between X1and X2. 以X1为首元素,X2为末元素平均生成100个元素的行向量。 LINSPACE(X1,X2,N)generates N points between X1and X2. For N<2,LINSPACE returns X2. 以X1为首元素,X2为末元素平均生成n个元素的行向量。如果n<2,返回X2。 Class support for inputs X1,X2: float:double,single 数据类型:单精度、双精度浮点型。 (4)说明函数ones()、zeros()、eye()的用法。 ones()生成全1矩阵。 zeros()生成全0矩阵。 eye()生成单位矩阵。 2、Matlab的数值显示格式

MATLAB期末上机试题带答案

MATLAB 期末上机考试试题带答案版 姓名: 学号: 成绩: 1.请实现下图: x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图: x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,'k--') grid on xlabel('x') ylabel('y')

title('sin(x)') legend('y=sin(x)') y2=cos(x); subplot(2,2,2) plot(x,y2,'r--') grid on xlabel('x') ylabel('y') title('cos(x)') legend('y=cos(x)') y3=tan(x); subplot(2,2,3) plot(x,y3,'k-') grid on xlabel('x') ylabel('y') title('tan(x)') legend('y=tan(x)') y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel('x') ylabel('y') title('cot(x)') legend('y=cot(x)') 3.解方程组: a=[3 2 1;1 -1 3;2 4 -4];b=[7;6;-2] ; x=a\b 4.请实现下图:

x y x=linspace(0,4*pi,1000); y1=sin(x); y2=sin(2*x); plot(x,y1,'--',x,y2,'b*') grid on xlabel('x');ylabel('y');title('耿蒙蒙') legend('sin(x)','sin(2*x)') 5.请在x ,y 在(-2,2)内的z=xexp (-x 2-y 2) 绘制网格图 [x,y]=meshgrid(-2:0.1:2); z=x.*exp (-x.^2-y.^2); mesh(x,y,z) 6.请实现peaks 函数: -5 5 x Peaks y [x,y]=meshgrid(-3:1/8:3); z=peaks(x,y); mesh(x,y,z)

MATLAB上机题汇总(较难题,抽1题)合肥学院

1.根据下面的方程绘制图形,y=cos(x)+x ,x 取值-pi-pi ,然后用多项式拟合函数进行拟合。 2.要求编写一个函数文件,将一个给定图片的所有颜色都变深(变深多少自己斟酌),其中黑色为0,白色为255 3.根据下面的方程绘制图形,y=cos(2x)+sin(2x),然后用5次多项式进行拟合。 4.要求编写一个M 函数,完成下列函数的求解:y =sin (x );y =sin (x +pi );y =sin (x +pi /3);y =sin (x +pi/2)自变量x 的变化范围从-pi 到pi ,要求四个图像画在同一个图形窗体中。 5.M 文件的编写 计算分段函数 ?????≠+=+++=10101)1cos(2x x x x x x x y 6.函数文件的编写 编写函数文件,求半径为r 的圆的周长和面积。

7.题目见教材p85,其中的数据可以给出文本文件或者由下面的式子表达 a=[10*rand(24,1)+15,1*rand(24,1)+30,2*rand(24,1)+45],这是24行3列的一个矩阵,第一列表示温度,第二列表示湿度,第三列表示气压。要求如下: (1) 将温度数据取出 (2) 找出一天中最高温度 (3) 按温度进行从小到大(默认)排序 (4) 求出最高温度的时刻 (5) 按温度进行从大到小排序 8.已知Fibonacci (斐波那)数列由式子??=+=--4,3,21k a a a k k k 可以生成,其中初值为121==a a ,试编写出生成某项Fibonacci 数值的MATLAB 函数,并求20a 的值。 9.编写M 函数,判断给定的三边能够构成一个三角形?并求其周长 10.给定一个声音文件,如1.wav ,引入这段声音,将音量改为2倍之后,再保存成2.wav

2010专第一次上机作业: MatLab练习题-基本运算

一、书上命令练习: 二、MatLab 练习1 注意: 1,请按要求作题,写出程序并把结果附上,并保存文件. 2,仔细填写 “学生信息”. 3,把此word 文档的文件名更名为你的 “学号 姓名(*)”-*表示第几次作业.在下课前5分钟提交. 1, (5*2+1.3-0.8)*10^2/25 答案: 2,y =23log )10sin( 4) 43.0(2+?+-e π 答案: 3,z = 10sin(π/3) ... cos(π/3) 答案: 4,x = sin(223/3), y = x^2,z = y*10 ; 求x+2y-5z 答案: 5,i e z i z i z ?==+=π23215,2,72,计算3 22 1z z z z z += 答案: 6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x (填写程序语句) 答案: 7, 建立等差一维数组x :首项为0,末项为π,项数为15 (填写程序语句) 答案:

8, 计算行列式的值1 876243 23=A 答案: 9, 矩阵??????????=187624323A ,矩阵?? ?? ? ?????=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积的程序语句。 答案 10,写出计算 Sin(36o )的程序语句. 答案:

三、MatLab 练习2-矩阵及数组运算 1.输入矩阵。 2.输入一个与A同阶的随机矩阵B。要求元素为整数。 3.输入一个数组b=(1 3 5 7 9 11 ) 4.输入与A同阶的单位矩阵E 5.输入与B同阶的零矩阵Z 6.构建一个大矩阵 7.将 46-70这25个整数填入一个五行五列的矩阵数表G中,使其各行、各 列以及主对角线和次对角线的和等相等。 8.分别求d1=|A|、d2=|B|、d3=|D|、d4=|G| 9.对以上四个矩阵A、B、D、G中的可逆矩阵求逆。 10.(1) 计算d5=A+B (2) 计算d6=A*B 11.取A的1,3,5行做一个子阵a1. 12.取A的2,4,6行,4,5,6列上的交叉元素做一个子阵a2.

matlab练习题和答案

matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233,,1443678,i,,,,2357,,,,2335542,i,,,, A,1357B,,,2675342, i,,3239,,,,189543,,,,1894,, 再求出它们的乘积矩阵C,并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后,调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];

MATLAB期末上机试题带答案

MATLAB期末上机试题带答案

MATLAB 期末上机考试试题带答案版 姓名: 学号: 成绩: 1.请实现下图: x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图:

x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,'k--') grid on xlabel('x') ylabel('y') title('sin(x)') legend('y=sin(x)') y2=cos(x); subplot(2,2,2) plot(x,y2,'r--') grid on xlabel('x') ylabel('y') title('cos(x)') legend('y=cos(x)') y3=tan(x); subplot(2,2,3) plot(x,y3,'k-') grid on xlabel('x')

ylabel('y') title('tan(x)') legend('y=tan(x)') y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel('x') ylabel('y') title('cot(x)') legend('y=cot(x)') 3.解方程组: a=[3 2 1;1 -1 3;2 4 -4];b=[7;6;-2] ; x=a\b 4.请实现下图:

matlab上机练习及答案

第二、三次上机练习: 目的: 运行课本第四章及课堂上讲过的例子,掌握Matlab 的流程控制语句、函数及脚本文件的编程、调试方法。 作业: 1、 完成下列操作: 1) 求[100,999]之间能被21整除的个数。 2) 建立一个字符串向量(要求字符串向量中必须包含自己的姓名首字母,大小写均可),删除其中的大写字母 2. 编写脚本文件,实现 用magic(6)产生一矩阵,用for 循环指令求解其所有元素的和。 3. 定义一个函数文件,求∑=n i m i 1 ,要求在函数文件中包含能够通过help 查询到的说明;然后调用该函数文件求∑∑∑===++ 10 1501210011k k k k k k 的值。 4. 已知) 7.1cos(12 ln )7.1sin(++++=x x x y π,当x 取-3.0,-2.9,-2.8,…,2.8,2.9,3.0时, 1) 求各点的函数值; 2) 求这些数据的平均值; 5、求分段函数的值。 222 603 565231x x x x y x x x x x x x ?+-<≠-?=-+≤<≠≠??--?, 且, 0且及, 其它 用if 语句实现,分别输出x=-5.0,-3.0,1.0,2.0,2.5时的值。 6、输入一个百分制成绩,要求输出成绩等级A 、B 、C 、D 、E 。其中90分~100分为A ,80分~89分为B ,70分~79分为C ,60分~69分为D ,60分以下为E 。 要求:分别用if 语句和switch 语句实现。 7、根据2 22221 111...,6123n ππ=++++求的近似值。当n 分别取100、1000、10000时,结 果是多少?(要求:分别用循环结构和向量运算来实现) 8、已知n=1时,f 1=1;n=2时,f 2=0;n=3时,f 3=1;n>3时,f n =f n-1 -2f n-2+ f n —3; 求f 1~ f 100中,最大值、最小值以及各数之和。 9、某公司员工的工资计算方法如下: 1)工作时数超过120小时者,超过部分加发15%; 2)工作时数低于60小时者,扣发700元;

相关文档