内容发布更新时间 : 2024/6/8 10:48:18星期一 下面是文章的全部内容请认真阅读。

练习1 基础练习

一、矩阵及数组操作：

1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）。

解eye(3) eye(15,8) ones(3) ones(15,8) zeros(3) zeros(15,8) X1=(1+(1-(-1))*rand(3)) X2=(1+(1-(-1))*rand(15，8)) Y1=1+sqrt(4)*randn(3) Y2=1+sqrt(4)*randn(15，8)

2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。

解：a=fix((10-0+1)*rand(10)+0) x=5 i=find(a>=x) n=length(i)

3．在给定的矩阵中删除含有整行内容全为0的行，删除整列内容全为0的列。

解：A=[1 2 2 3 4;0 0 0 0 0;4 5 6 2 3;0 0 0 0 0;]；

A([2,4],:)=[] B=A'； B(:,[2,4])=[]

二、绘图：

4．在同一图形窗口画出下列两条曲线图像： y1=2x+5； y2=x^2-3x+1， 并且用legend标注。 解：x=0:0.01:10;

y1=2*x+5; y2=x.^2-3*x+1;

1

plot(x,y1,x,y2) legend('y1', 'y2')

80706050403020100-10 0y1y2 12345678910

5．画出下列函数的曲面及等高线： z=x^2+y^2+sin(xy)．

解：[x,y]=meshgrid(0:0.25:4*pi); z=x.^2+y.^2+sin(x.*y)； contour3(x,y,z); meshc(x,y,z)

4003002001000151010500515

2

三、程序设计：

6．编写程序计算（x在[-3，3]，间隔0.01）

解：x=input('??ê?è?xμ??μ:'); if x>=-3&x<-1 y=(-x.^2-4*x-3)/2; elseif x>=-1&x<1 y=-x.^2+1; elseif x>=1&x<=3 y=(-x.^2+4*x-3)/2; else y='error' end y

7．有一列分数序列：

求前15项的和。

解：a=1; b=2; sum=0; for k=1:15 c=b/a; sum=sum+c; t=b; b=a+b; a=t; end sum

8．用至少三种方法编写函数实现求任意整数n的阶乘。

解：法一：s=factorial(n)

法二：gamma(n) 求出的是（n-1）！ 法三：n=input('please input n:');

x=1:n;

3