matlab语言习题集(3)

>>h=x1-x2+x3+8*x4-1; >>i=7*x1+x2-2*x3+2*x4-5; >>[x1,x2,x3,x4]=solve(f,g,h,i) x1 =

13/17 x2 =

320/187 x3 =

212/187 x4 =

19/187 或

>>a=[2 -3 1 2;1 3 0 1;1 -1 1 8;7 1 -2 2]; >>b=[8 6 1 5]'; >>xx=inv(a)*b

5．绘制函数曲线，要求写出程序代码（12分）。

（1）在区间[0:2π]均匀的取50个点，构成向量；

（2）在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求y1曲线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号；

（3）为坐标轴做相应的标注，分别是“时间”和“函数值”； （4）给图行加标题，为“三角函数波形图”。 答：

>>t=linspace(0,2*pi,50);

>>y1=sin(2*t-0.3);y2=3*cos(t+0.5); >>plot(t,y1,'r-.o',t,y2,'b--*');

>>xlabel('时间');ylabel('函数值'); >>title('三角函数波形图')

6．分析下面程序运行后s1,n1,s2,n2,m的值。

s1=0;s2=0;n1=0;n2=0;

x=[1,-4,-8,3,10,-9,7,-3,10, 8,-5,-2,2,0 ]; m=length(x); for i=1:m if x(i)<0

s1=s1+x(i); n1=n1+1; else

s2=s2+x(i); n2=n2+1; end end

s1,n1,s2,n2,m

请回答s1,n1,s2,n2,m的结果。 答：

s1 = -31，n1 = 6，s2 = 41，n2 = 8，m = 14 7．A = [1 2; 3 4]; B = [1 2; 1 4]; 写出下面几条指令的运行结果：C = A*B;D = A+B;E = A.*B。 答：

C =

3 10 7 22 D =

2 4 4 8 E =

1 4 3 16

8．(1) 写出关系运算符中的等号、不等号、小于号、大于号、小于等于号和大于等于号； (2) 写出逻辑操作中逻辑“与”、逻辑“或”及逻辑“非”的符号；

(3) 并用语句表达“如果a等于 b且c等于0就让d 等于3，否则如果a大于b且c=1让d等于 0，其余情况让d等于3”。 答：

(1) ==, ~=, <, >, <=, >=; (2) &, |, ~; (3)

if a==b&c==0 d=3;

else if a>b&c= =1 d=0; else d=3; end 9．设A=?aij?10*10和B=?bij?10*10是两个10行10列的矩阵（数组），试说明命令A*B，A/B，A\\B，A.*(B.^A)，A./B，A.\\B和sin(A).*exp(B)的涵义。 答：

矩阵A乘以矩阵B，矩阵A乘以矩阵B的逆，矩阵A的逆阵乘以矩阵B，矩阵B元素A次幂乘以矩阵A的对应元素，矩阵A元素除以矩阵B的对应元素，矩阵B元素除以矩阵A的对应元素，矩阵A元素的正弦与矩阵B对应元素的指数相乘。 10．设A为14*17 数组，B为一个行数大于3的数组，请给出：

(1) 删除A的第4、8、12三列的命令； (2) 删除B的倒数第3行的命令。 答：

(1) >>A(:,[4 8 12])=[] (2) >>B(end-2,:)=[]

11．请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5]。 答：

>>repmat(1:5,100,1)

12．已知变量A='ilovematlab';B='matlab'; 请找出：

(1) B在A中的位置；

(2) 把B放在A后面，形成C='ilovematlabmatlab'。 答：

>>strfind(A,B):(strfind(A,B)+length(B)-1) >>C=[A,B]

13．你知道几种数据矩阵的建立方式？请简述并举例。（6分） 答：

①不规则矩阵可直接在命令窗口输入>> X=[1 3 5;7 9 11;3 6 8] X =

1 3 5 7 9 11 3 6 8 ②规则矩阵的创建

增量矩阵创建方式，一般格式为X=（First:Increment:Last） >> X=(0:2:10) X =

0 2 4 6 8 10

14．在命令窗口中输入“a＝[2 2 2 2]”，请问：

（1）可以用“a”直接表示多项式吗？若能，请写出该多项式的数学表达式； （2）如果输入命令poly(a)，则表示什么意思？

（3）写出将第(2)问的多项式转换成自变量为“x”的符号形式的完整命令。 答：

（1）能，2x3?2x2?2x?2。

（2）求取矩阵特征多项式系数向量。 （3）

>> poly2sym(a,'x') ans =

2*x^3 + 2*x^2 + 2*x + 2 15．依次写出下列指令的运行结果。（6分）

>>a=zeros(3,3) >>a(3,4)=5^2 >>a(6)=100 >>a(:,3) >>a(:,:) >>size(a) 答：

a =

0 0 0 0 0 0 0 0 0

a =

0 0 0 0 0 0 0 0 0 0 0 25

a =

0 0 0 0 0 0 0 0 0 100 0 25

ans = 0 0 0

ans =

0 0 0 0 0 0 0 0 0 100 0 25

ans =

3 4

16．写出下列指令的运行结果。（5分）

>>A=1:2:7.5 >>B=1:4 >>C=[B;A] >>D=fliplr(C) >>E=sort(D) 答：

A =

1 3 5 7

B =

1 2 3 4

C =

1 2 3 4 1 3 5 7

D =

4 3 2 1 7 5 3 1

E =

4 3 2 1 7 5 3 1

17．写出下列指令的运行结果。（3分）

>>X=[2,0;2,4];Y=[-1 -4;1 0];A=2; >>X.*Y >>X*Y >>X.*A 答：

ans =

-2 0 2 0

ans =

-2 -8 2 -8

ans =

4 0 4 8

18．写出下列指令的运行结果。（4分）

>>a=eye(4); >>b=pascal(4) >>tril(b,0) >>a

b =

1 1 1 1 1 2 3 4 1 3 6 10 1 4 10 20

ans =

1 0 0 0 1 2 0 0 1 3 6 0 1 4 10 20

ans =

0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 19．写出下列指令的运行结果。（3分）

>>a=[5 2 -6;6 7 8; 12 10 8;20 40 30]; >>median(a) ans =

9.0000 8.5000 8.0000

>>cumsum(a) ans =

5 2 -6 11 9 2 23 19 10 43 59 40

>>diff(a,2) ans =

5 -2 -14 2 27 22

20．现有下列脚本程序，请分析该“for”循环执行情况并写出该程序运行结果。（4分）

for m=1:3 for n=1:3

H(m,n)=m^2+n end end

答：

生成三行三列的数组。

m=1,n=1,2,3时，依次生成H =[ 2 3 4] m=2,n=1,2,3时，依次生成H =[ 5 6 7] m=3,n=1,2,3时，依次生成H =[ 10 11 12] 最后结果 H =

2 3 4 5 6 7 10 11 12 21．设有程序：

>>A=[2,0,1,0;1,0,1,5;0, 6 ,2, 3]; >>B=[1,2,-1;1,3,1]; >>C=A(1:3,[1,2]);

>>D=size(C)- length(B(2,:)); >>F= C.*B?;

将运行结果写出来。

D =

0 -1 ； F= 2 0 2 0 0 6

。 22．设有程序：

>>A=[2,3,2,1;1,2,2,1; 1,2,3,-2]; >>A1=sum(A>A(2,1)); >>A(2,:)=[]; >>A2=A.^2-4;

将运行结果写出来。

A1= 1 3 3 0 ； A2=

0 5 0 -3

-3 0 5 0

。 23．设有程序：

>>P=[1,2,5,4;2,2,1,3]; >>max(P,[ ],1); >>mean(P,2);

将运行结果写出来。 max(P,[ ],1)=

2 2 5 4

；

mean(P,2)=

3 2

。 24．设有程序：

x=0; n=1;