MATLAB 题库

《MATLAB编程技术》题库

基础知识

1、多条命令可以放在一行，中间用逗号或者分号分开，逗号会显示语句的结果（回

显），分号不回显。如果一行语句太长，在结尾处加三个点…续行；注释语句不能续行，注释语句另起一行时前面加%

2、通常，在MATLAB内部，每个数据元素都是用双精度来存储的，format命令

只影响数据的输出格式，而不影响数据的计算和存储。

已知圆周率的值为：3.14159265358979：

>> clear all; format short; pi

ans =

3.1416

>> format long; pi

ans =

3.14159265358979

3、A(:) 将A的每一列元素堆叠起来，成为一个列向量，这也是MATLAB变量的

内部存储格式。

4、一些函数的功能：

linspace(a,b,n)与冒号表达式a:(b-a)/(n-1):b等价

>> linspace(0,10,6)

ans =

0 2 4 6 8 10

>> A = 1:18; B = reshape(A,3,6); B(:,[2,4])=[]

B =

1 7 13 16

2 8 14 17

3 9 15 18

>> A = [1:3; 2:4]; A>=2

ans =

0 1 1

1 1 1

>> A = [1:3; 2:4]; B=[3,2,1; 4,3,2]; A>B ans =

0 0 1

0 0 1

>> disp('It''s a book')

It's a book

>> x=3; y=4; m = '[x*y, x+y]'; eval(m) ans =

12 7

>> A=5; A=A==10

A =

>> x=-5; y=3; [rem(x,y),mod(x,y)]

ans =

-2 1

>> x=-18.98; [fix(x),floor(x),ceil(x),round(x)] ans =

-18 -19 -18 -19

>> x = 1:3; y=4:6; 2.^[x; y]

ans =

2 4 8

16 32 64

>> x = 1:3; y=4:6; x.^y

ans =

1 3

2 729

pwd 返回当前路径

clc 清除当前的命令行窗体显示内容

dir 显示当前路径下的文件

length(A),给出A的行数和列数中的较大者。ndims(A),给出A的维数

基本编程

一、判断1000以内有多少个素数，显示这些素数，并给出所用的时间。

clear all; close all; clc; tic

p = primes(1000); m = size(p,2);

disp(['There are ',num2str(m),' primes less than 1000.They are as follows:']); disp(p);

toc;

二、给出解下列方程或者方程组的程序：

1、??

?

??=++=++=++321222321z y x z y x z y x ；

solve('x+y+z=1','x+y^2+z^3=2','x^2+y^2+z^2=3')

2、0322=++c bx a ax 其中a 是未知数； solve('a*x^2 + a^2*b*x + c^3 = 0','a')

3、??

?

??=++=++=-+36542321z y x z y x z y x （此线性方程组要求用两种方法）。

solve('x+y-z=1','x+2*y+3*z=2','4*x+5*y+6*z=3')

A = [1,1,-1; 1,2,3; 4,5,6];

B = [1; 2; 3]; S=A\B;

x = S(1); y=S(2); z=S(3);

三、画图：

1、绘制如右图所示的13角形；设置线的属性为“蓝色虚线 线宽为4”。 figure(1);

plot(cos(0:pi*10/13:10*pi),sin(0:pi*10/13:10*pi),'b--','LineWidth', 4 );

2、绘制一条三维曲线??

?

??<<+=+=π100)10/()cos()10/()sin(z z z y z z x 。

figure(2);

z=0:pi/50:10*pi; x=sin(z)./(z+10); y=cos(z)./(z+10);

plot3(x,y,z,'LineWidth',3);

3、 画出2

222)

sin(),(y x y x y x f ++=

在区域??

???<<1515

y x 上的图形（Z>0.35部分设置为0.35）。

figure(3);

[x,y] = meshgrid(-15:0.5:15,-15:0.5:15); z=cos(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2)+eps; z(find(abs(z)>0.35))=0.35; mesh(z); xlabel('x'); ylabel('y');

4、绘制如右图所示的13角形，设置线的属性为“红色虚线 线宽为5”；

figure(1);

plot(cos(0:pi*12/13:12*pi),sin(0:pi*12/13:12*pi),'r--','LineWidth',5);

四、用两种方法求下列积分：

?20)(ln dx

x

Jifen = quad('log(x)',0,2)

x = 0+eps:0.0001:2;

y = log(x);

Jifen= 2*sum(y)/20000 ?20dx

e x

Jifen = quad('exp(x)',0,2)

x = 0:0.001:2;

y = exp(x);

Jifen = 2*sum(y)/2001

五、对下列数据进行曲线拟合（5次多项式）并绘制图形显示拟合效果。

x= 0:2:14;

y= [1.2, 1.8, 2.3, 5, 4, 3.3, 2.0, 1.0];

P = polyfit(x, y, 5);

C = 0:0.2:14;

D = polyval(P, C);

plot(x, y, 'r-', C, D, '*G-');

六、符号运算编程

求F = e x的逆函数

求F = ln(x)的逆函数

F(x) = e x；G(x)=sin(x); 给出复合函数F(G(x))。

F = cos(2x2)+5; 令x=7; 求符号表达式F的值。

F = ax3+by2+cx+d; 其中y是自变量，求符号表达式F微分。

F = cos(x)+ln(y); 其中y是自变量，求符号表达式F在[2,3]上的积分。

F = log10(5x); 画出符号表达式F在[2,3]上的图形。

F = cos(x)+ln(y); 其中y是自变量，求符号表达式F微分。

F = sym('exp(x)');

finverse(F)

F = sym('log(x)');

finverse(F)

F = sym('exp(x)');

G = sym('sin(x)');

compose(F,G)

F = sym('cos(2*x^2)+5');

x=sym('7');

eval(F)

F = sym('a*x^3+b*y^2+c*x+d');

diff(F,'y')

F = sym('cos(x)+log(y)');

int(F,'y',2,3)

F = sym('log10(5*x)');

ezplot(F,[2,3])

F = sym('cos(x)+log(y)');

diff(F,'y')

七、矩阵的运算

已知 A = [1,0; 2,2]; B = [1,0; 0,2]; 试求：

A*B A/B A./B A\B A.^B 的值。

1 0 1 0 1 NaN 1 0 1 1

2 4 2 1 Inf 1 -1 1 1 4

已知 A = [1,0; 0,2]; B = [1,0; 2,2]; 试求：

A*B A/B A./B A\B A.^B的值。

1 0 1 0 1 NaN 1 0 1 1

4 4 -2 1 0 1 1 1 0 4

综合编程（15）

☆八、综合编程：TXT文件的读写操作。

>> [a1,a2,a3,a4]=textread('test1.txt','%s%s%s%s','headerlines',4)

idin=fopen('test.txt'); % 打开test2.txt文件

fidout=fopen('mkmatlab.txt','w'); % 创建MKMATLAB.txt文件

while ~feof(fidin) % 判断是否为文件末尾

tline=fgetl(fidin); % 从文件读行

if double(tline(1))>=48&&double(tline(1))<=57 % 判断首字符是否是数值fprintf(fidout,'%s\n\n',tline); % 如果是数字行，把此行数据写入文件MKMATLAB.txt

continue % 如果是非数字继续下一次循环

end

end

fclose(fidout);

MK=importdata('MKMATLAB.txt'); % 将生成的MKMATLAB.txt文件导入工作空间，变量名为MK，

☆九、综合编程：Excel读写。

clear all

close all

clc

[Nm,Nm,c] = xlsread('C:\名单.xls');

[Nm,Nm,C] = xlsread('C:\成绩库.xls');

[n,Nm] = size(c);

[N,Nm] = size(C);

for i= 2:n

c{i,3} = [];

for j = 1:N

if c{i,1} == C{j,2}

c{i,3} = [c{i,3},C{j,4}];

end

end

c{i,3} = max(c{i,3});

end

xlswrite('c:\成绩单.xls',c);

☆十、综合编程：矩阵、元胞、结构体数组的访问和操作，mat文件的使用，动态字段方法。

论述（15）

十一、简述M脚本文件和M函数文件的主要区别。

（1）M脚本文件名不必与函数名相同；（2）M函数文件调用时需要提供函数参数，而M脚本文件只需要提供文件名；

（3）M脚本文件创建的每个变量都会保存到Matlab工作区，而M函数文件只有输出结果才保存到Matlab工作区中；

（4）文件格式上的一些区别

十二、谈谈在MATLAB 编程过程中的一些问题：设置断点、如何单步调试、如何优化提高程序运行速度。

点击在要设置断点的行左侧的-，将其变成圆圈后即设置了断点，也可按F12设置断点。条件断点在debug菜单下的set/modify conditional breakpoint进行设置。

设置断点后按F5运行程序，会在断点处停止运行，按F10可以单步运行调试，以上内容也可在Debug菜单下和工具条（第一条）中找到。

----------------------------------------------- The End-------------------------------------- 综合编程题和论述题，一定有。老师要求实验报告考试当场上交。