《马特来不实验及答案》由会员分享，可在线阅读，更多相关《马特来不实验及答案（6页珍藏版）》请在金锄头文库上搜索。
1、实验一、二1、上机运行教材 2.1 节中的例子，直至理解为止。2、自定义一个变量，然后分别用 8 种不同的数字显示格式显示查看。3、自定义一个矩阵，分别对其进行超维调用和超维赋值，比较其结果，并分析产生这种结果的原因。4、自定义一个矩阵，给矩阵的全行赋值；取出矩阵中某一个或几个元素，然后组成一个新的矩阵。5、分别用不同的方法给复数赋值。6、自定义一个复数矩阵，分别实现共轭转置、共轭、转置等运算。7、自己举例，完成对 eye, ones, zeros, magic, linspace, logspace 等函数的实验。答案：2 format3 A = rand(5);A(5,6)A(5,6)=1
2、;4 A(1,:) = 6; c=A(1,3,2,4);%取出 1、3 行与 2、4 列交点上的元素5 %共两种方法：逐个元素赋值，实部和虚部分别赋值。实验三1、 在0,2 内产生均匀分布的 200 个点，形成 1200 向量赋给变量 X。2、 在0,10内每间隔 3 产生向量 Y。3、 产生 6*6 阶的随机数矩阵 R，分别求其逆矩阵，矩阵的行列式。4、 设 x=x1;x2;x3;x4，A*x =B，求下列方程组的解 x：12344123475896108x5、 a)列出 2*2 阶的单位矩阵 I，4*4 阶的魔方矩阵 M 和 4*2 阶的全幺矩阵A，2*4 全零矩阵 B。b)将这些矩阵拼接
3、为 6*6 阶的矩阵 C，并将 C 输出至 txt 文本文档mtest.txt中（在实验报告中提交 txt 文档截屏图像） 。c)取出 C 的第 2、4、6 行，组成 3*6 阶的矩阵，取出 C 的第 2、4、6 列，组成 6*3 阶的矩阵。6、 已知矩阵：a=1 2 3;3 4 5;2 4 6 b=3 6 2;3 7 9;6 2 8 请运算 a.*b；a*b；a./b ；a/b；a.b ；ab ；a2 ；a.2 ；查看输出结果并进行比较，分析普通乘除运算和元素群运算的区别。7、 设 ，89132576431A 782934165B求 C1=AB， C2=AB， C3=A.*B 并求它们的逆阵
4、。答案：1 X = linspace(0,2*pi,200);2 Y = 0:3:10;3 R = rand(6);v = inv(R);d = det(R);%R = randn(6);%R = randi(-16,16,4,6)4 A = 3,4,-7,-12;5,-7,4,2;1,8,0,-5;-6,5,-2,10;B = 4;-3;9;-8;X = inv(A)*B;%X = AB5a I = eye(2)M = magic(4)A = ones(4,2)B = zeros(2,4)b: C = I,B;A,Msave test.txt -ascii Cc: c= C(2,4,6,:)
5、 ;%取行值c= C(:,2,4,6) %取列值实验四8、 用 for 循环语句实现 1100 求和。9、 编制一个函数文件，利用 while-end 循环实现 1100 求和。10、 用 for 循环语句实现一个求 n 阶乘的函数文件。11、 绘制参数方程 y=sin(t),z=cos(t)在 t=0,7区间的三维曲线，并加入网格线（提示：利用grid on 命令） 。12、 用 hold on 命令在同一个窗口绘制曲线 y=sin(t)，y1=sin(t+0.25) y2=sin(t+0.5)，其中t=0 10，要求利用不用的颜色和线形将三条曲线区分开来。答案：1sum = 0 ;for
6、i = 1:100sum = sum + i;end4t = 0:.01:7;y = sin(t);z = cos(t);plot3(t,y,z)5t=0:.01:10;y=sin(t);plot(t,y,+g);hold ony1=sin(t+0.25);plot(t,y1,-r);hold ony2=sin(t+0.5);plot(t,y2,.b);实验五1、给定两个多项式 a(x)=6x5+2x4-8x3-3x2+7x-4，b(x)=x 3-7x2+4x+2，分别求a*b，a+b 和 a/b。2. 给定矩阵 a，分别求其三角分解、正交分解和奇异值分解。a=3,4,7;2,9,12;6,8
7、,9;3、产生 42 阶的正态随机数矩阵 R1，求其各列的平均值和标准差，并求 全体的平均值和标准差；对 R1 的各列进行排序（按照 降序排列） ，排列后的结果输出至 sort.txt4. 设方程的根为 x=-3,-5,-8,-9，求它所对应的 x 多项式的系数。5、求代数方程 的所有根543237910xx6.y=x5-4x4+3x2-2x+6 a) 当 x=-2,8之间取值时（取 100 个点） ，求函数的值，并画出曲线。 7、给定一个信号：y=6sin(2t-1)*cos(4t+2)，t 的取值范围为0,2 ，取值间隔为 0.01，试画出其幅频和相频特性。答案：1a = 6 2 -8 -
8、3 7 -4;b = 1 -7 4 2;c = conv(a,b)a+0 0 bf,g = deconv(c,a)2l,u=lu(a);%三角分解q,r=qr(a);%正交分解u,s,v=svd(a);%奇异值分解3 R1=randn(4,2);m=mean(R1);s=std(R1);ma=mean(R1(:);sa=std(R1(:);sort(R1,descend);4x=-3,-5,-8,-9;a=poly(x);5a=3 4 7 2 9 12;x=roots(a)6 a=1 -4 0 3 -2 6;x=linspace(-2,8);y=polyval(a,x);fplot(x5-4*
9、x4+3*x2-2*x+6,-2,8)y =x.5-4*x.4+3*x.2-2*x+6;plot(x,y)ezplot(x5-4*x4+3*x2-2*x+6,-2,8)b=roots(a)7t=0:.01:2*pi;y=6*sin(2*t-1).cos(4*t+2);am=abs(y);ang=angle(y);subplot(1,2,1);plot(t,am);title(幅频特性)subplot(1,2,2);plot(t,ang);title(相频特性)实验六2、设 ，在 x1,4区间内：321()sin5cosfxx1）画出它的曲线；2）求出它的过零点的值，并判断是否存在极大值、极小值
10、，并分别求出它们的值。2、设 ：542()36fxx1）求 x=-2,8之间函数的值（ 取 100 个点） ，画出曲线，看它有几个过零点（提示：利用 polyval 函数） ；2）用 roots 函数求此多项式的根。3、将上题写成一个函数文件 f1.m，用 fzero 函数求它的过零点；若在该函数中加一项 xsin(x)，求过零点。4、已知线型常微分方程： ， 8695yyuu在输入信号 u(t)为单位阶跃信号时的解析解，并绘制出其曲线。5、用 eval 命令执行字符串：s=y=magic(3)。6、编写程序，要求利用 for 循环及 eval 语句实现：yn=magic(n)(n=3,4,5
11、)。7、已知三次曲线的方程为 y=ax3+cx，请采用编制函数文件、由主程序调用的方法，探讨参数 a 和 c 对其图形的影响。 （提示：仿照例 5-1-1（2）的求解方法） 。答案：1. （1） editfunction y=quxian(x)y=x.3-2*x2*sin(x)+5*x*cos(x)+1./x;fplot(quxian,1,4)(2) a1=fzero(quxian,1)a 1=1.5117 a2= fzero(quxian,2.5)a2=2.6095 m=fminbnd(quxian,1,4)m =2.1589 quxian(2.1589)ans =-3.21872(1) a
12、=1,-4,3,-2,6 x=linspace(-2,8) t=polyval(a,x) plot(x,t)(2) roots(a)ans =3.0000 1.6956 -0.3478 + 1.0289i3function y=(x)y=x.5-4*x4+3*x2-2*x+6;m=fzero(f1,-1)m=-1.2607n=fzero(f1,1)n=1.3479z=fzero(f1,4)z=3.79994 a=1 0 8 6 0 0;b=6 -9 5;r p k=residue(b,a); t=0:10;ys=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)+r(3)*exp(p(3)*t)+r(4)*exp(p(4)*t)+r(5)*exp(p(5)*t); plot(t,ys)5. s=y=magic(3)s =y=magic(3) eval(s)y =8 1 63 5 74 9 26.for n=3:5s=y=magic(n);eval(s)end
