matlab程序题复习

1.编写m 文件要求利用matlab 数值绘图功能画出衰减振荡曲线t e

y t 3sin 3-

=及其它的包络线3

0t e y -

=。t 的取值范围是]4,0[π，数值采样间隔点为pi/50，t e

y t 3sin 3

-

=在图中用红色是

线段绘出，3

0t e

y -=用蓝色虚线绘出，且两函数作在同一张图上。

1.t=0:pi/50:4*pi;y0=exp(-t/3);

y=exp(-t/3).*sin(3*t);figure;

plot(t,y,'-r')hold on;

plot(t,y0,':b');plot(t,-y0,':b');hold off;

title(‘y=exp(-t/3)sin(3t)’);xlabel(‘x’)ylabel(‘y’)

2.利用matlab 数值绘图功能，画出2

2

22)

sin(y

x y x z ++=

所表示的三维曲面。y x ,的取值范围是

]8,8[-，要求去除当x,y 均等于0时，函数z 的不连续点，数值坐标采样间隔为0.5.

2.clear;x=-8:0.5:8;

y=x';

［X Y］=meshgrid(x,y);R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;surf(X,Y,Z);(mesh(X,Y,Z))colormap(cool)

xlabel('x'),ylabel('y'),zlabel('z')

Title(‘z=sin(sqrt(x^2+y^2)/sqrt(x^2+y^2));

3.利用matlab 数值绘图功能，绘制t 在区间［-2pi 2pi］时函数y=sin(t)/t 的图形，采样间

隔为pi/10，要求利用逻辑运算求出x=0时函数近似极限，并修补图形缺口。

3.t=-2*pi:pi/10:2*pi;

y=sin(t)./t;

tt=t+(t==0)*eps;yy=sin(tt)./tt;

subplot(1,2,1),plot(t,y),axis([-7,7,-0.5,1.2]),xlabel('t'),ylabel('y')；

subplot(1,2,2),plot(tt,yy),axis([-7,7,-0.5,1.2])

4.编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。

4.function c=collatz(n)

c=n;

while n>1

if rem(n,2)==0

n=n/2;

else

n=3*n+1;

end

c=[c n];

end

5.有一周期为4 的正弦波上叠加了方差为0.1的正态分布的随机噪声的信号，用循环结构编制一个三点线性滑动平均的程序。（提示：①用0.1*randn(1,n)产生方差为0.1的正态分布的随机噪声；②三点线性滑动平均就是依次取每三个相邻数的平均值作为新的数据，如x1(2)=(x(1)+x(2)+x(3))/3，x1(3)=(x(2)+x(3)+x(4))/3……）。

5.t=0:pi/50:4*pi;

n=length(t);

y=sin(t)+0.1*randn(1,n);

ya(1)=y(1);

for i=2:n-1

ya(i)=sum(y(i-1:i+1))/3;

end

ya(n)=y(n);

plot(t,y,'c',t,ya,'r','linewidth',2)

6.有一组测量数据如下表所示，数据具有y=x2的变化趋势，用最小二乘法求解y。并要求图示之。

x1 1.52 2.53 3.54 4.55

y-1.4 2.73 5.98.412.216.618.826.2

6.>>x=[11.522.533.544.55]'

>>y=[-1.42.735.98.412.216.618.826.2]'

>>e=[ones(size(x))x.^2]

>>c=e\y

>>x1=[1:0.1:5]';

>>y1=[ones(size(x1)),x1.^2]*c;

>>plot(x,y,'ro')

hold on;

plot(x1,y1,'k');

xlabel(x);

ylabel(y);

title('拟合曲线');legend;

7.2

2

y x

xe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口

中绘制出三维线图、网线图、表面图和带渲染效果的表面图。7>>[x,y]=meshgrid([-2:.2:2]);

>>z=x.*exp(-x.^2-y.^2);>>mesh(x,y,z)>>subplot(2,2,1),plot3(x,y,z)>>title('plot3(x,y,z)')>>subplot(2,2,2),mesh(x,y,z)>>title('mesh (x,y,z)')>>subplot(2,2,3),surf(x,y,z)

>>title('surf (x,y,z)')

>>subplot(2,2,4),

surf(x,y,z),shading interp

>>title('surf (x,y,z),shading interp')

8.利用matlab 绘图功能，分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x),要求作在同一张图上，并标注标题。8.解：

x=0:pi/10:2*pi;y=2*sin(x);subplot(2,2,1);bar(x,y,'g');

title('bar(x,y,''g'')');axis([07–22]);subplot(2,2,2);stairs(x,y,'b');title('stairs(x,y,''b'')');axis([07–22]);

subplot(2,2,3);stem(x,y,'k');

title('stem(x,y,''k'')');axis([07–22]);subplot(2,2,4);fill(x,y,'y');

title('fill(x,y,''y'')');axis([07–22]);

必考9.使用MATLAB 画一个圆心在原点、半径等于10的圆，并在圆周上依逆时钟方向取任意四点A 、B 、C 、D ，将线段AB 、AC 、AD 、BC 、BD 、CD 用直线画出。计算线段AB 、AC 、AD 、BC 、BD 、CD 的长度。

clear all;t=linspace(0,2*pi,50);r=10;x=r*cos(t);

y=r*sin(t);

a=r*[cos(0.50*pi),sin(0.50*pi)];

b=r*[cos(0.90*pi),sin(0.90*pi)];

c=r*[cos(1.25*pi),sin(1.25*pi)];

d=r*[cos(1.80*pi),sin(1.80*pi)];

plot(x,y,'b',a(1),a(2),'.k',b(1),b(2),'.k',c(1),c(2),'.k',d(1), d(2),'.k');

axis image

text(a(1),a(2),'A');

text(b(1),b(2),'B');

text(c(1),c(2),'C');

text(d(1),d(2),'D');

line([a(1),b(1)],[a(2),b(2)],'color','r');t=(a+b)/2;text(t(1), t(2),'AB');

line([b(1),c(1)],[b(2),c(2)],'color','r');t=(b+c)/2;text(t(1), t(2),'BC');

line([c(1),d(1)],[c(2),d(2)],'color','r');t=(c+d)/2;text(t(1), t(2),'CD');

line([d(1),a(1)],[d(2),a(2)],'color','r');t=(d+a)/2;text(t(1), t(2),'DA');

line([a(1),c(1)],[a(2),c(2)],'color','r');t=(a+c)/2;text(t(1), t(2),'AC');

line([b(1),d(1)],[b(2),d(2)],'color','r');t=(b+d)/2;text(t(1), t(2),'BD');

ab=sqrt((a(1)-b(1)).^2+(a(2)-b(2)).^2);

fprintf('ab=%f\n',ab);

bc=sqrt((b(1)-c(1)).^2+(b(2)-c(2)).^2);

fprintf('bc=%f\n',bc);

cd=sqrt((c(1)-d(1)).^2+(c(2)-d(2)).^2);

fprintf('cd=%f\n',cd);

ad=sqrt((a(1)-d(1)).^2+(a(2)-d(2)).^2);

fprintf('ad=%f\n',ad);

ac=sqrt((a(1)-c(1)).^2+(a(2)-c(2)).^2);

fprintf('ac=%f\n',ac);

bd=sqrt((b(1)-d(1)).^2+(b(2)-d(2)).^2);

fprintf('bd=%f\n',bd);

必考10.试写一函数regPolygon(n)，其功能为画出一个圆心在(0,0)、半径为1的圆，并在圆内画出一个内接正n边形，其中一顶点位于(0,1)。例如regPolygon(8)可以画出如下之正八边型：

10.>>function regpoly(n)

vertices=[1];

for i=1:n

step=2*pi/n;

vertices=[vertices,exp(i*step*sqrt(-1))];

end

plot(vertices,'-o');

axis image

%画外接圆

hold on

theta=linspace(0,2*pi);

plot(cos(theta),sin(theta),'-r');

hold off

axis image

11.请用surf指令来画出下列函数的曲面图：z=x*exp(-x2-y2)其中x在[-2,2]间共等切分为21点，y在[-1,1]间共等切分为21点，所以此曲面共有21*21=441个点。请以曲面的斜率来设定曲面的颜色。请以曲面的曲率来设定曲面的颜色。

11.x=linspace(-2,2,21);%在x轴[-2,2]之间取21点

y=linspace(-1,1,21);%在y轴[-1,1]之间取21点

[xx,yy]=meshgrid(x,y);%xx和yy都是21×21的矩阵

zz=xx.*exp(-xx.^2-yy.^2);%计算函数值，zz也是21×21的矩阵subplot(1,3,1)

surf(xx,yy,zz);axis image

subplot(1,3,2)

surf(xx,yy,zz,gradient(zz));axis image

subplot(1,3,3)

surf(xx,yy,zz,del2(zz));axis image

12.试写一函数regStar(n)，其功能为画出一个圆心在(0,0)、半径为1的圆，并在圆内

画出一个内接正n星形，其中一顶点位于(0,1)。

12.解：

function regStar(n)

vertices=[1];

for i=1:n

step=2*pi*floor(n/2)/n;

vertices=[vertices,exp(i*step*sqrt(-1))];

end

plot(vertices,'-o');

%画外接圆

hold on

theta=linspace(0,2*pi);

plot(cos(theta),sin(theta),'-r');

hold off

axis image

13.利用matlab 函数功能，编写function 文件要求利用matlab 数值绘图功能画出衰减振荡曲线t sin e

y t

55

-=及其它的包络线5

0t e y -

=。t 的取值范围是]4,0[π，数值采样间隔点为pi/50，

t sin e

y t 55

-

=在图中用红色是线段绘出，5

0t e y -=用蓝色虚线绘出，且两函数作在同一张图上,

并说明其调用方式。

38.function

y=test(t);

y0=exp(-t/3);

y=exp(-t/3).*sin(3*t);figure;

plot(t,y,'-r')hold on;

plot(t,y0,':b');plot(t,-y0,':b');hold off;

调用t=0:pi/50:4*pi;y=test(t)

14.有一正弦衰减数据y=cos(x).*exp(-x/5)，其中x=0:pi/5:4*pi，利用matlab 三次样条插值函数进行插值，要求每两个数据点间插入3个数，并作出拟合后曲线。

14.>>x0=0:pi/5:4*pi;>>y0=cos(x0).*exp(-x0/5);>>x=0:pi/20:4*pi;>>y=spline(x0,y0,x);>>plot(x0,y0,'or',x,y,'b')title（‘曲线拟合’）；xlabel（‘x’）；ylabel（‘y’）；

15.利用matlab 数值绘图功能，画出2

2

223333y

x )

y x sin(f ++=

所表示的三维曲面。y x ,的取值

范围是],[66-，要求去除当x,y 均等于0时，函数z 的不连续点，数值坐标采样间隔为0.2。15.clear all;x=-6:0.2:6;

y=x';

［X Y］=meshgrid(x,y);R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;

surf(X,Y,Z);(mesh(X,Y,Z))colormap(cool)

xlabel('x'),ylabel('y'),zlabel('z')

16.利用matlab 指令计算表达式(

)

2

23

230y x e y x f ---=的梯度并利用数值绘图功能绘出梯度

图。（x,y 的取值范围［-33］采样间隔0.1）

16.>>

v =-3:0.1:3;

>>[x,y]=meshgrid(v);

>>z=30*(x.^2-y.^3).*exp(-x.^2-y.^2);>>[px,py]=gradient(z,.1,.1);>>contour(x,y,z)>>hold on

>>quiver(x,y,px,py)>>hold off

17.有一周期为3π的正弦波上叠加了方差为0.2的正态分布的随机噪声的信号，利用matlab 循环结构编制一个三点线性滑动平均的程序去除噪声。（提示：①用0.2*randn(1,n)产生方差为0.2的正态分布的随机噪声；②三点线性滑动平均就是依次取每三个相邻数的平均值作为新的数据，如x1(2)=(x(1)+x(2)+x(3))/3，x1(3)=(x(2)+x(3)+x(4))/3……）。

17.t=0:pi/50:3*pi;n=length(t);

y=sin(t)+0.2*randn(1,n);ya(1)=y(1);for i=2:n-1

ya(i)=sum(y(i-1:i+1))/3;end

ya(n)=y(n);

plot(t,y,'c',t,ya,'r','linewidth',2)

18.分段函数else x /x /,x /)x cos()/cos(y π

ππππππ23

8372305<<<<<

?

??=,x 取样区间［04π］，区间内采样500个点，要求利用逻辑关系运算对函数y=sin(x)逐段解析计算，并绘图。（给出绘图指令即可）。

18.t=linspace(0,4*pi,500);y=cos(t);

z1=((t2*pi)).*y;

w=(t>pi/3&t<2*pi/3)+(t>7*pi/3&t<8*pi/3);w_n=~w;

z2=w*cos(pi/5)+w_n.*z1;

plot(t,z2,'-b'),axis([010-11])

19.有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi，用三次样条法进行插值。

19.>>x0=0:pi/5:4*pi;>>y0=sin(x0).*exp(-x0/10);>>x=0:pi/20:4*pi;>>y=spline(x0,y0,x);>>plot(x0,y0,'or',x,y,'b')

20.计算表达式()

2

2

e 1053y x y x z ---=的梯度并利用matlab 数值绘图功能绘图。（x,y 的取值范

围［-22］采样间隔0.2）

20.>>

v =-2:0.2:2;

>>[x,y]=meshgrid(v);

>>z=10*(x.^3-y.^5).*exp(-x.^2-y.^2);>>[px,py]=gradient(z,.2,.2);>>contour(x,y,z)>>hold on

>>quiver(x,y,px,py)>>hold off

21.分段函数sin(/3)/32,7/38/302sin()x x y x x else πππππππ<<<

=<

,x 取样区间［03π］，区间内采样300个点，要求利用逻辑关系运算对函数y=sin(x)逐段解析计算，并绘图。（给出绘图指令即可）。

21.t=linspace(0,3*pi,500);y=sin(t);

z1=((t2*pi)).*y;

w=(t>pi/3&t<2*pi/3)+(t>7*pi/3&t<8*pi/3);w_n=~w;

z2=w*sin(pi/3)+w_n.*z1;

plot(t,z2,'-b'),axis([010-11])

22.有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。并在结果图中添加标题-at e =y ，并用箭头线标识出各曲线a 的取值并添加标题-at e =y 和图例框。

22.>>t=0:0.5:10;

>>y1=exp(-0.1*t);>>y2=exp(-0.2*t);

>>y3=exp(-0.5*t);

>>plot(t,y1,'-ob',t,y2,':*r',t,y3,'-.^g')

>>title('\ity\rm=e^{-\itat}')

>>title('\ity\rm=e^{-\itat}','FontSize',12)

>>text(t(6),y1(6),'\leftarrow\ita\rm=0.1','FontSize',11) >>text(t(6),y2(6),'\leftarrow\ita\rm=0.2','FontSize',11) >>text(t(6),y3(6),'\leftarrow\ita\rm=0.5','FontSize',11) >>title('\ity\rm=e^{-\itat}','FontSize',12)

>>legend('a=0.1','a=0.2','a=0.5')

概念格matlab程序

clear all U=形式背景 [m,n]=size(U); waiyan={}; waiyan{1}=1:m; temp1=[]; temp={}; for i=1:n temp1(1)=1; for j=2:m if(U(j,i)==1) temp1=[temp1,j];%得到形式背景 end end temp{i}=temp1; temp1=[]; end a=intersect(cell2mat(temp(1)),cell2mat(waiyan(1))); b{1}=a; waiyan=[waiyan,b]; b=[]; for i=2:n m=size(waiyan); for j=1:m(2) a=intersect(cell2mat(temp(i)),cell2mat(waiyan(j)));%求子集 if(~isziji(a,waiyan));%若原来不存在该概念 waiyan=[waiyan,a];%则加入该概念 end end end disp('概念格构造为：'); for i=1:max(size(waiyan)) t=num2str(cell2mat(waiyan(i))); disp(strcat('{',t,'}')); end ----------------------------------------------------------------------------- function logic=isziji(a,waiyan) m=size(waiyan); m1=max(size(a)); for i=1:max(m) m2=max(size(cell2mat(waiyan(i))));

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程序大全答案

频率特性类题目 1 一个系统的开环传递函数为 ，试绘制其当K=5、30时系统的开环频率特性Nyquist 图，并判断系统的稳定性。 w=linspace(0.5,5,1000)*pi; sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600) figure(1) nyquist(sys1,w) title('system nyquist charts with k=5') figure(2) nyquist(sys2,w) title('system nyquist charts with k=30') 由图可知K=5时，开环Nyquist 曲线没有包围（-1，j0）点，所以系统稳定。 K=30时，开环Nyquist 曲线包围（-1，j0）点，所以系统不稳定。 2系统开环传递函数为 ，建立其零极点增益模型， 然后分别绘制当K=5、K=30时系统的开环频率特性Bode 图，并判断系统的稳定性。 sys1=zpk([],[0 -10 -2],100) sys2=zpk([],[0 -10 -2],600) figure(1) [Gm1,Pm1,Wcg1,Wcp1]=margin(sys1) bode(sys1) title('system bode charts with k=5'),grid figure(2) [Gm2,Pm2,Wcg2,Wcp2]=margin(sys2) bode(sys2) title('system bode charts with k=30'),grid 因为K=5时，Wcg>Wcp,所以系统稳定。 K=10时，Wcg

《Matlab语言》期末考试试题样例

一、填空题：（每空2分）大约20空 （注：写出完成要求所需的命令语句或题目给定语句执行后输出的 结果） 1. Matlab 中用来表示不确定值的变量是 NaN ；默认变量i 表示 虚数单位 。 2. Matlab 中符号 ； （分号） 可以使命令行不显示运行的结果，表示该行为注释行的符号为 % （百分号） 。 3. Matlab 将工作空间所有变量保存在默认目录下的命令是 Save 。(命令复习) 4. 骰子会掷出1到6的随机数字，用函数rand 和取整函数(Important)模拟骰子某次投掷出的数字的命令为 ceil(6*rand) ，模拟两 个骰子某次投掷的数字组合之和情况的命令为 。 5. 生成一个含有10个元素的行向量x ，其元素在[-1,1]之间均匀分布的的命令为 x=2*rand(1,10)-1 ，使用向量x 的后5个元素作为对角线 元素生成对角阵A 的命令为 A=diag(x(end-4:end)) ，将A 的主对 角线下方的元素（不含主对角线）全部填充为2014的命令为 B=ones(5); 2014*B-triu(B)+A 6. 为函数22 110at t x +=分布定义内联函数命令是 ，定义匿名函数的命令是 7. 在循环结构中使用 break 语句可以终止当前循环。 8. Matlab 用于绘制统计直方图的命令是 hist ；用于绘制y 轴为对数坐标的单对数坐标图的命令为 semilog 。

二、解答题：大约七道题目 1. 在[0,π2]区间上绘制如下图所示的sin(x)与sin(x 2)的曲线，并如图为图像添加标题以及为坐标轴添加标注。 2.绘制上面曲面的图形（提示：该曲面由上下两部分组成：） 3.创建命令式的M 文件计算圆形的面积S ，提示用户输入圆形的半径r ，并输出圆形的面积。（注意输入输出都要对用户添加一定的提示语） %Function Circle 4.请利用matlab 的符号工具箱做下面的运算，并写出相应的命令。 (1) 求极限x e x L x x -+=→10)1(lim ； (2)求积分：dx x x ?∞+02sin ； Syms x; f=((1+x)^(1/x)-exp(1)/x; limit(f,x,0) 2222, [2,2], [1,1]z x y x y =+∈-∈-

Matlab 期末考试题库(共12套卷)

第一套 华东交通大学2015-2016学年第一学期Matlab期末考试 一、填空题（4*6=24分） 1、在MATLAB命令窗口中的“>>”标志为MATLAB的提示符，“│”标志为提示符。 2、符号表达式sin(2*a+t)+m 中独立的符号变量为__ _。 3、在通常情况下，左除x=a\b是的解，右除x=b/a是的解，一般情况下，。 4、为了使两个plot的图形在同一个坐标显示，可以使用_____命令进行图形保持；可以使用_ __命令为图形添加网格。 5、倘若要是对x进行赋值，从5到25，间隔是0.1,为；倘若要是对x进行赋值，从2到20，中间间隔100个点,为 6、A=[1,2;3,4];B=[1,0;0,1]; A*B= _ , A.*B=_ _； 二、选择题（4*6=24分） 7、如果x=1: 2 : 8,则x(1)和x(4)分别是( ) A．1，8 B．1, 7 C．2, 8 D．2, 7 8、运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为( ) c=input('请输入一个字符','s'); if c>='A' & c<='Z' disp(setstr(abs(c)+abs('a')-abs('A'))); elseif c>='a'& c<='z' disp(setstr(abs(c)- abs('a')+abs('A'))); elseif c>='0'& c<='9' disp(abs(c)-abs('0')); else disp(c); end A. 9 B. 8 C. 7 D. 6 9、MATLAB表达式2*2^3^2的结果是() A．128 B．4096 C. 262144 D．256 10、在循环结构中跳出循环，执行循环后面代码的命令为( ) (A) return (B) break (C) continue (D) keyboard 11、在图形指定位置加标注命令是（） A. title(x,y,?y=sin(x)?); B. xlabel(x,y,?y=sin(x)?); C. text(x,y,?y=sin(x)?); D. legend(x,y,?y=sin(x)?); 12、下列哪个变量的定义是不合法的（）。 (A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz 三、程序题（52分） 13（10分）、在同一窗口绘制6个子图，已知t=0.1:0.1:2*pi，每个子图分别绘制（其中y1=sin(t),y2=cos(2t),y3=sin(t)+cos(t),y4=3t, y5=2ln(t),y6=et）； 14（10分）、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price<200 没有折扣 200≤price<500 3%折扣 500≤price<1000 5%折扣 1000≤price<2500 8%折扣 2500≤price<5000 10%折扣

matlab考试题及答案

%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (iA(i)) y=A(i); end end disp('最小元素是：') disp(y) % 7、某系统由四个典型环节构成，试根据四阶龙格—库塔公式，求输出量y的阶跃响应，取仿真时间Tf=10s， %步长h=。

MATLAB简单程序大全

MATLAB简单程序大全 求特征值特征向量 A=[2 3 4;1 5 9;8 5 2] det(A) A' rank(A) inv(A) rref(A) eig(A)%求特征值和特征向量 卫星运行问题 h=200,H=51000,R=6378; a=(h+H+2*R)/2; c=(H-h)/2; b=(a^2-c^2)^(1/2); e=c/a; f=sqrt(1-exp(2).*cos(t)^2); l=int(f,t,0,pi/2) L=4*a.*l 动态玫瑰线 n=3;N=10000; theta=2*pi*(0:N)/N; r=cos(n*theta); x=r.*cos(theta); y=r.*sin(theta); comet(x,y) 二重积分 syms x y f=x^2*sin(y); int(int(f,x,0,1),y,0,pi) ezmesh(f,[0,1,0,pi]) 函数画图 syms x;f=exp(-0.2*x)*sin(0.5*x); ezplot(f,[0,8*pi])

玫瑰线 theta=0:0.01:2*pi; r=cos(3*theta); polar(theta,r,'r') 求x^2+y^2=1和x^2+z^2=1所围成的体积 syms x y z R r=1; Z=sqrt(1-x^2); y0=Z; V=8*int(int(Z,y,0,y0),x,0,1) 求导数及图像 f='1/(5+4*cos(x))'; subplot(1,2,1);ezplot(f) f1=diff(f) subplot(1,2,2);ezplot(f1) 绕x轴旋转 t=(0:20)*pi/10; r=exp(-.2*t).*sin(.5*t); theta=t; x=t'*ones(size(t)); y=r'*cos(theta); z=r'*sin(theta); mesh(x,y,z) colormap([0 0 0]) 某年是否闰年 year=input('input year:='); n1=year/4; n2=year/100; n3=year/400; if n1==fix(n1)&n2~=fix(n2) disp('是闰年') elseif n1==fix(n1)&n3==fix(n3) disp('是闰年') else

MATLAB 期末复习试题

6. 用两种方法制出下列行向量。 并将A B 连接成为两种新的矩阵。 A= 1 2 3 4 5 6 7 8 9 B=-5 -2 1 4 7 10 13 16 19 8， 创建3行3列的零矩阵，并用赋值的方式使其成为9行9 9， 将习题5中创建的矩阵，将其第3行到第4行，第2列到第5列取出，为一2行4列的子矩阵，并将此子矩阵中的第2行 第3 10，求 11. 用 3. 将上题中的结果分别赋值给A,B,C,D,E ，并将ABCDE 组成一个5行3列的矩阵。找出该矩阵中所有绝对值大于9的元素，将 对应元素重新赋值为0；再找出该矩阵中所有绝对值小于1的元素，并将对应元素重新赋值为1. 4. 向量操作时MATLAB 的主要部分，使用已给出的向量ok 来做下面练习。 ok = [0 1 2 3 4 5 0 0 -1 -2 -3 -4 -5 0 ] 创建一个新的向量oka ，使其为ok 的转置； 创建一个新的向量okb ，使其为ok 中得前半部分与后半部分对换的结果。（即okb 中的元素为ok 的后半部分“0 -1 -2 -3 -4 -5 -” 在前，前半部分“0 1 2 3 4 5 0”在后； 创建一个新向量okc ，使其包含ok 中所有小于0的元素，且元素的排序与ok 创建一个新向量okd ，使其包含ok 中所有等于0或者等于1的元素的下标值 写出生成下图所示波形的MATLAB 程序。图中虚线为余弦波，x 范围为 处被削顶。 ]3,3[],3,3[ππππ

假设a 是这样一组数组：a(n)=a(n-1)2-a(n-2)2+a(n-3)2，且a(1)=0；a(2)=1；a(3)=2，试分别用循环指令for 以及while 来寻求该数组中第一个大于10000的元素。 分别用四个图形表示连续调制波形y=sin(t)sin(9t)，t 的范围为0~11。其中，第一个图形，在t 的范围内均匀取12个点画出y,图 形用红色的点显示；第二个图形，在t 范围内均匀取101个点，图形用蓝色的点显示；第三个图形，把第一个图形的图用红色线条显示；第四图形，把第二个图形的图用蓝色线条显示。 建立一个函数文件jie.m ，其作用为计算关于n 再建立一个m 主程序文件，调用已建立的jie.m ，计算s=2+22+23+…+220 。 function y=jie(n) y=2^n; return ////// clear clc s=0; for n=1:20 y=jie(n); s=s+y; end

Matlab考试题库答案

填空题 1、标点符号; %—用来表示该行为注释行。 可以使命令行不显示运算结果， 2、x 为0 ～4pi ，步长为 的向量，使用命令 x=0:*pi:4*pi 创建。 3、输入矩阵A=[1323?57569 ] ，使用全下标方式用A(2,2)取出元素“-5 ”，使用单下标方式用A(5)取出元素“-5 ”。 4、符号表达式sin(2*a+t)+m 中独立的符号变量为 t 。 5、M 脚本文件和M 函数文件的主要区别是M 脚本文件没有函数定义和M 函数文件有函数定义_______。 6. 设x 是一维数组，x 的倒数第3个元素表示为x(_end-2_) 设y 为二维数组，要删除y 的第34行和48列，可使用命令y(34,:)=[] ； y(:,48)=[]； 7. 将变量x 以Ascii 文本格式存储到文件，应使用命令 save _x ； 8. 在while 表达式, 语句体, End 循环语句中，表达式的值 非零 时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 9.要从键盘读入一个字符串并赋值给变量x ，且给出提示“Who is she ”，应使用命令x=input(‘Who is she’,’s’) ； 10．设A=(a ij )m?n 和B=(b ij )m?n 和C=(c ij )m?n 均为m*n 矩阵，且存在于WorkSpace 中，要产生矩阵D=(a ij ?c ij b ij c ij ?)m?n ，可用命令D=(A-C)/B.^C ， 计算|(A ′B )?1| 可用命令det(inv(A’*B) 11. 在MATLAB 命令窗口中的“>>”标志为MATLAB 的 命令行 提示符，“│”标志为 输入 提示符。 12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。 A+B ； A.*B ； A==B ； ans= [2,2,6;5,10,0;7,9,11] ans= [1,0,9;4,25,0;0,8,18] ans= [1,0,1;0,1,1;0,0,0] 13.已知A 是矩阵，求A 的对角矩阵函数是 diag ； 求A 的下三角矩阵函数是 tril 。 的程序文件和Simulink 模型文件的扩展名分别是.m 、.mdl 最基本的绘图函数为plot() 16. A=[1,2,3;4,5,6]; A(:,[1,3])=[];A= [2；5] 17. fix= -1 , round= -2 . 18.已知A=[1 3;4 6];C=[A,eye(sizeA);A ,ones(sizeA)] 则Ｃ＝ [1 3 1 0；4 6 0 1；1 3 1 1；4 6 1 1]

Matlab程序代码

Matlab程序代码： clc; clear; N=20; T=0.1 t=0:T:N m=length(t) syms x1 x2 x3 fx=[0;x1+x2^2;x1-x2] gx=[exp(x2);exp(x2);0] hx=x3; R=10*eye(1) Q=[10 0 0;0 1 0;0 0 1] A=[0 1 0;0 0 1;0 0 0] B=[0;0;1] SS=B*inv(R)*B' [p1,p2,lamp,perr,wellposed,P]=aresolv(A,Q,SS) z1=hx z2=[diff(hx,x1) diff(hx,x2) diff(hx,x3)]*fx z3=[diff(z2,x1), diff(z2,x2), diff(z2,x3)]*fx ax=[diff(z3,x1), diff(z3,x2), diff(z3,x3)]*fx bx=[diff(z3,x1), diff(z3,x2), diff(z3,x3)]*gx z=[z1;z2;z3] k=inv(R)*B'*P %diff(z)=A*z+B*v=(A-B*K)*Z %x(0)=[1;0;0] abk=A-B*k x1(1)=1 x2(1)=0 x3(1)=1 z1(1)=x3(1) z2(1)=x1(1)-x2(1) z3(1)=-(x1+x2^2) for i=2:m z1(i)=z1(i-1)+T*(abk(1,1)*z1(i-1)+abk(1,2)*z2(i-1)+abk(1,3)*z3(i-1)) z2(i)=z2(i-1)+T*(abk(2,1)*z1(i-1)+abk(2,2)*z2(i-1)+abk(2,3)*z3(i-1)) z3(i)=z3(i-1)+T*(abk(3,1)*z1(i-1)+abk(3,2)*z2(i-1)+abk(3,3)*z3(i-1))

MATLAB期末上机试题带答案

MATLAB 期末上机考试试题带答案版 姓名： 学号： 成绩： 1.请实现下图： 50 100 150 200 250 -1 -0.8-0.6-0.4-0.200.20.40.60.81x 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.请实现下图： 2 4 6 8 10 12 14 -1-0.8-0.6-0.4-0.200.20.40.60.81x y '自己的姓名' sin(x)sin(2*x) x=linspace(0,4*pi,1000); y1=sin(x);

matlab基础练习题及答案

第1章 MATLAB 基础 MATLAB 操作桌面有几个窗口如何使某个窗口脱离桌面成为独立窗口又如何将脱离出去的窗口重新放置到桌面上 答：在MATLAB 操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop 菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 如何启动M 文件编辑/调试器 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动。在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器。 存储在工作空间中的数组能编辑吗如何操作 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 命令历史窗口除了可以观察前面键入的命令外，还有什么用途 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 在MATLAB 中有几种获得帮助的途径 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 在MATLAB 中如何建立矩阵??? ???194375，并将其赋予变量a 答：在Command Window 窗口输入操作: >> a=[5 7 3;4 9 1] 有几种建立矩阵的方法各有什么优点 答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷； （2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； （3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； （4）通过数据文件建立，该方法可以调用由其他软件产生数据。 在进行算术运算时，数组运算和矩阵运算各有什么要求 答：进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 数组运算和矩阵运算的运算符有什么区别 答：在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 计算矩阵??????????897473535与??? ?? ?????638976242之和。 答： 求??? ???+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。 答： 计算??????=572396a 与?? ? ???=864142b 的数组乘积。 答：

MATLAB复习题.doc

MATLAB期末考试题目 ------电气二班 一．选择题。 1．下列哪个变量的定义是不合法的（A） A. abcd-3 B. xyz_3 C. abcdef D. x3yz 2．在同一个图形窗口中画出三行两列的图形并指定当前可编辑图形为第3个图 的正确的 MATLAB命令是： A.subplot(3,2,1)B subplot(3,3,) C.subplot(3,2,3 ) D.subplot(1,3,2) 3.x=[1234],y=[1111],x.*y的结果是：(B) A．1234 B.出错信息 C.4 3 2 1 D.1111 4．角度x=[30 45 60]，计算其正弦函数的运算为（D） A. sin(deg2rad(x)) B. sin(x) C. sin(x) D. sin(deg2rad(x)) 5．i=2;a=2i;b=2*i;c=2*sqrt(-1);执行程序后，a,b,c的值分别是多少？（C） A.a=4,b=2,c=2.0000i B.a=4,b=2.0000i,c=2.0000i C.a=2.0000i,b=4,c=2.0000i D.a=2.0000i,b=2.0000i,c=2.0000i 6. legend()命令是(A) A．在图中标注“图例” B.．在图中加上网格线C．在图中加上标题 D．打印图形 7.命令plot(x,y,′m:*′),其中字符串“m:o”中,第一个字符“m”(A) A. 表示曲线为m符号构成的线 B. 表示曲线颜色为紫色 C. 表示曲线x坐标的最大值为m变量的值 D. 表示曲线y坐标的最大值为m变量的值 8. 下列可作为matlab合法变量名的是(D) A、合计 B、123 C、@h D、xyz_2a 9. 使用语句t=0:7生成的是()个元素向量(A) A、8 B、7 C、6 D、5 10. 输入字符串或选项时，要用(D)括住 A、() B、【】 C、{} D、‘’ 11 . 清空 Matlab 工作空间内所有变量的指令是（A）。 A.clc B. cls C.clear D. clf

MATLAB期末复习联练习题

M ATL A B(R2014a)仿真程序设计 一、选择题 1、下列哪个变量的定义是不合法（） A）a b c d-4 B）xyz_4a C）a bc def D）xy4z 2、在循环结构中跳出循环，执行循环后面的代码的命令为（） A）re tur n B）conti n u e C）break D）key b o ard 3、清空MATL AB工作空间内所有变量的指令是（） A）c l c B）c ls C）cl ear D）c l f 4、清空MATL AB中图形窗的指令是（） A）c l c B）c ls C）cl ear D）c l f 5、下列哪条指令是求矩阵的行列式的值（） A）inv B）d ia g C）e i g D）de t 6、用ro un d函数四舍五入对数组[1.49 6.39 3.93 8.51]取整，结果为（）A）[1 6 3 8] B）[1 6 4 8]C）[1 6 4 9]D）[2 7 4 9] 7、已知a=2:2:8,b=2:5。下列运算表达式中错误的为（） A）a*b B）a.*b C）a’*b D）a-b 8、角度x=[15 30 45 60]，计算其正弦函数的运算为（） A）S IN(de g2rad(x))B）si n(de g2rad(x))C）si n(x) D）SIN(x) 9、在M ATL A B中，下列变量最大数是（） A）-p i B）rea l m i n C）rea lm ax D）eps 10、在MATL AB中，下列变量比0大的最小数是（） A）p i B）rea l m i n C）real max D）i 11、下列表示可以作为M AT L A B的变量的是（） A）a b c d-2x B）x yz2016_2a C）@hy D）3h o ur 12、指出下列错误的指令（） A) sy ms a b B)sym s a, b C)sy m s(‘a',‘b')D)syms(‘a','b','p o si ti ve') 13、下列哪条指令是求逆矩阵的行列式的值（）。 A) de t B)d ia g C) e i g D) i nv 14、清除命令窗口（co m m a n d w i nd o ws）空间的命令是（） A）c l c B）c lear C）c l f D）de le te 15、已知函数文件如下，则fa c to r(4)=（） f u n c ti on f=fa c tor(n) i f n<=1 f=1;

Matlab源程序代码

正弦波的源程序： (一)，用到的函数 1,f2t函数 function x=f2t(X) global dt df t f T N %x=f2t(X) %x为时域的取样值矢量 %X为x的傅氏变换 %X与x长度相同并为2的整幂 %本函数需要一个全局变量dt(时域取样间隔) X=[X(N/2+1:N),X(1:N/2)]; x=ifft(X)/dt; end 2,t2f函数。 function X=t2f(x) global dt df N t f T %X=t2f(x) %x为时域的取样值矢量 %X为x的傅氏变换 %X与x长度相同,并为2的整幂。 %本函数需要一个全局变量dt(时域取样间隔) H=fft(x); X=[H(N/2+1:N),H(1:N/2)]*dt; end (二),主程序。 1，%（1）绘出正弦信号波形及频谱 global dt df t f N close all k=input('取样点数＝2^k, k取10左右'); if isempty(k), k=10; end f0=input('f0=取1(kz)左右'); if isempty(f0), f0=1; end N=2^k; dt=0.01; %ms df=1/(N*dt); %KHz T=N*dt; %截短时间

Bs=N*df/2; %系统带宽 f=[-Bs+df/2:df:Bs]; %频域横坐标 t=[-T/2+dt/2:dt:T/2]; %时域横坐标 s=sin(2*pi*f0*t); %输入的正弦信号 S=t2f(s); %S是s的傅氏变换 a=f2t(S); %a是S的傅氏反变换 a=real(a); as=abs(S); subplot(2,1,1) %输出的频谱 plot(f,as,'b'); grid axis([-2*f0,+2*f0,min(as),max(as)]) xlabel('f (KHz)') ylabel('|S(f)| (V/KHz)') %figure(2) subplot(2,1,2) plot(t,a,'black') %输出信号波形画图grid axis([-2/f0,+2/f0,-1.5,1.5]) xlabel('t(ms)') ylabel('a(t)(V)') gtext('频谱图') 最佳基带系统的源程序： (一)，用到的函数 f2t函数和t2f函数。代码>> (二),主程序 globaldt t f df N T close all clear Eb_N0 Pe k=input('取样点数＝2^k, k取13左右'); if isempty(k), k=13; end z=input('每个信号取样点数＝2^z, z

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) 产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏） 答案： a=10*rand(1,10)-5; b=sort(a,'descend') 1.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5] 2. 已知变量：A=’ilovematlab’；B=’matlab’, 请找出： （A）B在A中的位置。 （B）把B放在A后面，形成C=‘ilovematlabmatlab’ 3. 请修改下面的程序，让他们没有for循环语句！ A=[1 2 3; 4 5 6; 7 8 9]; [r c]=size(A); for i=1:1:r for j=1:1:c if (A(i,j)>8 | A(i,j)<2) A(i,j)=0; end end end 4. 请把变量A=[1 2 3; 4 5 6; 7 8 9]写到文件里，写完后文件看起来是这样的 1 2 3 4 5 6 7 8 9 5.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。 6.编写M文件，从Yahoo网站批量读取至在2008年9月份的每日收盘价（提示：使用字符串函数）。 7. 将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。 8.已知资产每日回报率为，标准差为，资产现在价值为亿，求5％水平下资产的10天在险价值(Var)。 =[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB中最简单的方法计算b，注意最简单哦。 1、求下列联立方程的解

matlab期末考试试卷及答案206

MATLAB期末考试题目 一．选择题。 1．下列哪个变量的定义是不合法的（A） A. abcd-3 B. xyz_3 C. abcdef D. x3yz 2．在同一个图形窗口中画出三行两列的图形并指定当前可编辑图形为第3个图 的正确的 MATLAB命令是： A.subplot(3,2,1)B subplot(3,3,) C.subplot(3,2,3 ) D.subplot(1,3,2) 3.x=[1234],y=[1111],x.*y的结果是：(B) A．1234 B.出错信息 C.4 3 2 1 D.1111 4．角度x=[30 45 60]，计算其正弦函数的运算为（D） A. sin(deg2rad(x)) B. sin(x) C. sin(x) D. sin(deg2rad(x)) 5．i=2;a=2i;b=2*i;c=2*sqrt(-1);执行程序后，a,b,c的值分别是多少？（C） A.a=4,b=2,c=2.0000i B.a=4,b=2.0000i,c=2.0000i C.a=2.0000i,b=4,c=2.0000i D.a=2.0000i,b=2.0000i,c=2.0000i 6. legend()命令是(A) A．在图中标注“图例” B.．在图中加上网格线C．在图中加上标题 D．打印图形 7.命令plot(x,y,′m:*′),其中字符串“m:o”中,第一个字符“m”(A) A. 表示曲线为m符号构成的线 B. 表示曲线颜色为紫色 C. 表示曲线x坐标的最大值为m变量的值 D. 表示曲线y坐标的最大值为m变量的值 8. 下列可作为matlab合法变量名的是(D) A、合计 B、123 C、@h D、xyz_2a 9. 使用语句t=0:7生成的是()个元素向量(A) A、8 B、7 C、6 D、5 10. 输入字符串或选项时，要用(D)括住 A、() B、【】 C、{} D、‘’ 11 . 清空 Matlab 工作空间内所有变量的指令是（A）。 A.clc B. cls C.clear D. clf