1、C2=AB;C3=A.*B,并求它们的逆阵。程序 A=1,4,8,13;-3,6,-5,-9;2,-7,-12,-8; B=5,4,3,-2;6,-2,3,-8;-1,3,-9,7;C1= A*B,C2 = A*B, C3 = A.*Binv(C1),inv(C2),inv(C3)3. a. 列出22阶的单位矩阵I, 44阶魔方矩阵M和42阶的全幺矩阵A,全零矩阵Bb. 将这些矩阵拼接为66阶的矩阵C:c. 求出C的第2,4,6行,组成36阶的矩阵C1,及第2,4,6,裂,组成63阶的矩阵C2,d. 求D=C1C2及D1=C2C1.程序 I=eye(2),A=ones(4,2),B=zeros
2、(4,2),M=magic(4),C=I,A;B,M C1=C(2,4,6,:),C2=C(:,2,4,6,) D=C1*C2, D1=C2*C14设 把x=02间分为101点,画出以x为横坐标,y为纵坐标的曲线程序 x=linspace(0,2*pi,101)y=cos(x).*(0.5+3*sin(x)./(1+x.2);plot(x,y),grid5求代数方程3x5+4x4+7x3+2x2+9x+12=0的所有根。a=3,4,7,2,9,12;r=roots(a)r = -0.8612 + 1.4377i -0.8612 - 1.4377i 0.6737 + 1.0159i 0.6737
3、 - 1.0159i -0.9583 6把1开五次方,并求其全部五个根。(提示:解x5-1=0)a1=1,0,0,0,0,-1;r1=roots(a1)r1 = -0.8090 + 0.5878i -0.8090 - 0.5878i 0.3090 + 0.9511i 0.3090 - 0.9511i 1.0000 7. 设方程的根为x=-3,-5,-8,-9,求它们对应的x多项式的系数。 a=poly(-3,-5,-8,-9)a=1 25 223 831 1080即 a(x) = x4 + 25x3+ 223x2+ 831x+10808. 设微分方程求输入u(t)=(t)时的输出y(t).解:
4、 系统传递函数为输入信号的拉普拉斯变换为U(s)=1,故输出信号的拉普拉斯变换为:求拉普拉斯反变换的程序为: b=1;a=1,2,5,4,3;r,p=residue(b,a) t=0:0.1:10;y=r*exp(p*t);plot(t,y),grid9. 产生86阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求全体的平均值和均方差。程序 x=randn(8,6)m=mean(x) mm=mean(mean(x) s=std(x) ss=std(x(:)10. 产生46阶的均匀分布随机数矩阵R,要求其元素在1到16之间取整数值。并求此矩阵前四列组成的方阵的逆阵。 R=rand(4,
5、6)x=ceil(16*R)y=inv(x(:,1:4)11. 设x=r cos t+3t, y=r sint+3,分别令r=2,3,4,画出参数t=010区间生成的xy曲线。for r=2:4x=r*cos(t)+3*t;y=r*sin(t)+3;plot(x,y)hold onendhold offgrid12. 设x=sin t, y=sin(Nt+), a) 若=常数,令N = 1,2,3,4,在四个子图中分别画出其曲线b) 若N=2,取=0,/3,/2,及,四个子图中分别画出其曲线 程序a x=sin(t);y=sin(t+pi/2);subplot(2,2,1),plot(x,y)
6、y=sin(2*t+pi/2);subplot(2,2,2),plot(x,y)y=sin(3*t+pi/2);subplot(2,2,3),plot(x,y)y=sin(4*t+pi/2);subplot(2,2,4),plot(x,y)程序by=sin(2*t);y=sin(2*t+pi/3);y=sin(2*t+pi);13设f(x)=x5- 4x4 +3x2- 2x+ 6(1)x=-2,8之间函数的值(取100个点),画出曲线,看它有几个过零点。用polyval 函数)(2)用roots函数求此多项式的根。 a=1,-4,0,3,-2,6;x=linspace(-2,8,100);f=
7、polyval(a,x);plot(x,f),gridr=roots(a) 14.设, 要求在z=010区间内画出x,y,z三维曲线。程序 z=0:x=z.*sin(3*z);y=z.*cos(3*z);plot3(x,y,z),grid15设求定义域x=-2,2,y=-2,2内的z值(网格取0.1间隔), x=-2:2;y=-2: X=x*ones(1,length(y); Y=ones(length(x),1)*y;X,Y=meshgrid(-2:2,-2:2); z=X.*X.*exp(-X.*X-Y.2); mesh(z),grid16. 设 z1=0.05x-0.05y+0.1; 画
8、出z1 的曲面(平面)图,迭合在上题的图中。 z1=0.05*x-0.05*y+0.1; hold on, mesh(z1);20. 设(1) 画出它在x=0,4 区间内的曲线。求出它的过零点的值。(2) 求此曲线在x轴上方及下方的第一块所围的面积的大小。函数程序:function y=pbf32(x)y=x.3-2*x.2.*sin(x)+5*x.*cos(x)+1./x;主程序: fplot(pbf32,0,4) grid x1=fzero(,1.5)x1 = 1.5117 x2=fzero(,2.5)x2 = 2.6095 s1=quad8(,0,x1)s1 = Inf s2=quad8(,x1,x2)s2 = -2.279421 已知微分方程:,若y(0)=1,求它在x=0,5区间内的数值积分,并画出曲线。function Dy=pbf33(x,y) Dy=x.2./y-x.*cos(y);x,Iy=ode23(pbf33,0,5,1); plot(x,Iy), grid22 eval命令执行字符串s=y=magic(3). s=magic(3) M= eval(s)M = 8 1 6 3 5 7 4 9 2
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1