数学实验练习题杨振华(MATLAB)

Matlab数学实验一——matlab初体验一、实验目的及意义[1] 熟悉MATLAB软件的用户环境；[2] 了解MATLAB软件的一般目的命令；[3] 掌握MATLAB数组操作与运算函数；通过该实验的学习，使学生能熟悉matlab的基础应用，初步应用MATLAB软件解决一些简单问题。

二、实验内容1．认识matlab的界面和基本操作2．了解matlab的数据输出方式（format）3. MATLAB软件的数组（矩阵）操作及运算练习；三、实验任务根据实验内容和步骤，完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会）完成如下题目，并按照实验报告格式和要求填写实验报告1．在command window中分别输入如下值，看它们的值等于多少，并用matlab的help中查询这些缺省预定义变量的含义，用中文写出它们的意义。

i j eps inf nan pi realmax realmin2．分别输入一个分数、整数、小数等，（如：a=1/9）,观察显示结果，并使用format 函数控制数据的显示格式，如：分别输入format short、format long、format short e、format long g、format bank、format hex 等，然后再在命令窗口中输入a，显示a的值的不同形式，并理解这些格式的含义。

3．测试函数clear、clc的含义及所带参数的含义（利用matlab的help功能）。

4. 写出在命令窗口中的计算步骤和运行结果。

（1）计算1.2210(ln log)81eππ+-；>> (log(pi)+log(pi)/log(10)-exp(1.2))^2/81 >>ans =0.0348(2)>> x=2;y=4;>> z=x^2+exp(x+y)-y*log(x)-3z =401.6562（3）输入变量135.3,25a b⎡⎤==⎢⎥⎣⎦,在工作空间中使用who，whos，并用save命令将变量存入”D:\exe01.mat”文件。

数学实验第一次测验题及参考答案一、写出下列MATLAB命令的运算结果.1. a=[1 2 3 4].\ [2 1 3 2], b=[1 2 3 4].* [2 1 3 2]a =2.0000 0.5000 1.0000 0.5000b =2 2 9 82. A=[1 2; 3 4]; a=A(2, :) , b=A(: , 2), c=A(:)'a =3 4b =24c =1 32 43. a=linspace(0, 10, 5), b=0:2.5:10a =0 2.5000 5.0000 7.5000 10.0000b =0 2.5000 5.0000 7.5000 10.00004. A=[1,2,3; 4,5,6]; b=A(2, 3), c=A([1 2], [1,3]) , d=size(A)b =6c =1 34 6d =2 35. x=7.1; [fix(x), floor(x) ceil(x) round(x)]ans = 7 7 8 7x=-7.1; [fix(x), floor(x) ceil(x) round(x)]ans = -7 -8 -7 -76. syms x y, z='x^2+x*y-y'; zx=diff(z, x), zy=diff(z, y)zx =2*x+yzy =x-17. syms t, x=[exp(t), cos(t), sin(t)]; dx=diff(x, t)dx =[ exp(t), -sin(t), cos(t)]8. syms x y, f=x*y; y=x^2; f1=subs(f, 'y', y)f =x^39. a=[1 2 3]; b=[4 5 6]; x=dot(a, b), y=cross(a, b)x = 32y = -3 6 -310. v=[1, 2, 3, 4]; A=diag(v); det(A)ans =24二、写出下列MATLAB 命令的实验目的. 1. f=inline('sum(1./(1:n).^2)', 'n'); f(10) 求∑=nk k121当10=n 时的值.2. syms x, limit(sin(x)/x, x, 0) 求极限xxx sin lim0→.3. x=0: 0.1: 5; y=zeros(1, length(x)); plot(x, y, 'r')画出函数0=y ，50≤≤x 的图形.4. syms x, f=3*x^2-x+2; c=roots([3, -1, 2]) 求多项式方程0232=+-x x 的解.syms x, x=-3:0.1:3; f=inline('x.^2 -4*x-10'); c1=fzero(f, [-3,0]) 求函数1042--x x 在区间]0,3[-的零点.5. f='x/(1+x.^2)'; [xmin, ymin]=fminbnd(f, -10,10) 求函数21x xf +=在区间]10,10[-的最小值.6. syms x, int('x^2*sin(x)', x, 0,1) 计算定积分⎰12sin xdx x 的值.syms x, quad('sin(x)./x', 0.001, 1) 求定积分⎰10sin dx x x的近似值.7. ezplot3('cos(t)', 'sin(t)', 't', [0, 8*pi])作出空间曲线⎪⎩⎪⎨⎧===t z t y t x sin cos 当]8,0[π∈t 时的图形.8. [x,y]=meshgrid(-2:0.1:2, -2:0.1:2); z=x.^2-y.^2; contour(x, y, z, 20) 作函数22y x z -=的等高线.9. s= solve( 'x^2 + x*y + y = 3', 'x^2 - 4*y + 3 = 0'), [s.x, s.y]求方程组 ⎩⎨⎧=+-=++034322y x y xy x 的解.10. syms x y; int( int(x*y^2, y, 0, x), x, 0, 1) 计算重积分⎰⎰1002xdx dy xy 的值.dblquad(inline('sin(exp(x.*y))'), 0, 1, 0, 1) 求重积分⎰⎰101)sin(dx dy exy的近似值.三、为下列实验目的写MATLAB 指令.1. 作出星形线t x 3cos 2=,t y 3sin 2=(π20≤≤t )的图形.ezplot('2*cos(t)^3','2*sin(t)^3',[0,2*pi])作出摆线)sin (2t t x -=, )cos 1(2t y -=(π40≤≤t )的图形.ezplot('2*(t-sin(t))','2*(1-cos(t))',[0,4*pi])2. 用极坐标命令, 作出五叶玫瑰线θρ5sin 4=的图形.theta=0:0.1:2*pi; rho=4*sin(5*theta); polar(theta,rho)用隐函数命令作出椭圆322+=+xy y x 的图形和双曲线3322+=+xy y x 的图形. ezplot('x^2+y^2-x*y-3',[-3,3,-3,3])ezplot('x^2+y^2-3*x*y-3',[-10,10,-10,10]) 3. 求!100的近似值. prod(1:100) 求∑=101!1n n 的近似值. s=0;for n=1:10s=s+prod(1./(1:n)); end s4. 求极限 1512lim 33++∞→n n n .syms nlimit((2*n^3+1)/(5*n^3+1),n,inf)5. 求函数x x x f 3cos sin )(=的一阶导数. 并求)1(f '. syms xdaoshu=diff('sin(x)*cos(3*x)') x=1;daoshuzhi=eval(daoshu)求由方程0122222=++++-y x y xy x 确定的隐函数的导数. 并求)1(f '. syms x yz=2*x^2-2*x*y+y^2+x+2*y+1; daoshu=-diff(z,x)/diff(z,y) x=1;daoshuzhi=eval(daoshu)6. 求由参数方程⎩⎨⎧==ty tx 33sin cos 确定的函数的导数. syms t x=cos(t)^3; y=sin(t)^3;daoshu=diff(y,t)/diff(x,t) daoshu=simple(daoshu)求由参数方程⎪⎩⎪⎨⎧+=+=3231616t t y t t x 确定的函数的导数. syms tx=6*t/(1+t^3); y=6*t^2/(1+t^3);daoshu=diff(y,t)/diff(x,t)7. 求⎰)(cos 02)(x dx x w dxd .diff(int('w(x)',0,(cos(x))^2))8. 求⎰-ππ2222sin xdx ex.syms xint('exp(2*x)*sin(2*x)^2',x,-pi,2*pi)9. 绘制曲面22y x z +=的图形 x=-2:0.1:2; y=-2:0.1:2;[x,y]=meshgrid(x,y); z=x.^2+y.^2; surf(x,y,z)作出空间曲线t t x cos =,t t y sin =,t z 2=(π60≤≤t )的图形. ezplot3('t*cos(t)','t*sin(t)','2*t',[0,6*pi])10. 作双曲抛物面4122y x z -=的图形，其中66≤≤-x ，1414≤≤-y . x=-6:0.1:6;y=-14:0.1:14;[x,y]=meshgrid(x,y); z=x.^2-y.^2./4; mesh(x,y,z)11. 设y xy z )1(+=,求x z ∂∂，yz ∂∂. syms x yz='(1+x*y)^y'; diff(z,x)，diff(z,y)设),(y xy f z =,求22x z ∂∂,22y z ∂∂,yx z∂∂∂2.syms x y z='f(x*y,y)'; zxx=diff(z,x,2) zyy=diff(z,y,2)zxy=diff(diff(z,x),y)12. 求级数∑∞=12k kk的和.syms k;s1=symsum(k/(2^k),k,1,inf) 求级数∑∞=-12)12(1k k 的和. syms k;s2=symsum(1/((2*k-1)^2),k,1,inf)13. 求函数)1ln()1(x x ++的6阶麦克劳林多项式. clear; syms xser3=taylor((1+x)*log(1+x),7)四、为下列实验目的编写MATLAB 程序.1. 作出分段函数⎪⎩⎪⎨⎧>≤<--≤+=1,11,11,23)(2x x x x x x f 的图形.y=[];for x=-3:0.1:3; if x<=-1y=[y, 3+2*x]; elseif x<=1y=[y, 1]; elsey=[y, x^2]; end endx=-3:0.1:3; plot(x,y)作出分段函数⎩⎨⎧>≤+=0,0,sin 1)(x e x x x h x 的图形.y=[ ];for x=-4:0.1:4 if x<=0y=[y, 1+sin(x)]; elsey=[y, exp(x)]; end endx=-4:0.1:4; plot(x, y)2. 设D 为由2=+y x ,y x =,2=y 所围成的有界区域.(1) 作出区域D 的图形;(2) 根据图形确定积分限. 先对x 后对y 积分计算⎰⎰Ddxdy xy2.x=0:0.01:2; y1=2-x; y2=x.^2; y3=2;plot(x, y1,'r-', x, y2,'b-', x, y3,'k-') syms x y;int(int(x*y^2, x, 2-y, sqrt(y)), y, 1,2) 3. 计算⎰⎰+-Dy x dxdy e )(22, 其中D 为122≤+y x .提示：用极坐标变换 22y x r +=，⎰⎰⎰⎰-+-=πθ201)(222rdrd e dxdy er Dy xtheta=0:0.1:2*pi; x=cos(theta); y=sin(theta);plot(x,y,'r-')clear;syms r theta; f=exp(-(r^2))*r ;int(int(f,r,0,1),theta,0,2*pi)4. 求曲线x x x g 2sin )(=)0(π≤≤x 与x 轴所围成的图形分别绕x 轴, y 轴旋转所成的旋转体体积，并写出作该旋转体的图形的指令. (1) 在图形绕x 轴旋转时，体积dx x g v ⎰=ππ02)(.syms xv=int('pi*(x*sin(x)^2)^2', x, 0, pi)v= 1/8*pi^4-15/64*pi^2 即)86415(22ππ+-=v .(2) 图形绕y 轴旋转时，体积 dx x xg v ⎰=ππ0)(2syms xquad('2*pi*x.^2*sin(x).^2', 0, pi) 则得到体积的近似值为 ans = 27.5349(3) 绕x 轴旋转所得旋转体图形的作图指令 r=0:0.1:pi; t=-pi:0.1:pi; [r, t]=meshgrid(r,t); x=r;z=r.*sin(t).*sin(r).^2; y=r.*cos(t).*sin(r).^2; surf(x,y,z)title('绕x 轴旋转');xlabel('x 轴'); ylabel('y 轴'); zlabel('z 轴')(4) 绕y 轴旋转所得旋转体图形的作图指令 r=0:0.1:pi; t=-pi:0.1:pi; [r, t]=meshgrid(r,t); x=r.*cos(t); z=r.*sin(t); y=r.*sin(r).^2; surf(x,y,z)title('绕y 轴旋转');xlabel('x 轴'); ylabel('y 轴'); zlabel('z 轴')5. 作出抛物柱面2y x =和平面1=+z x 相交的图形.u=-2:0.1:2; v=-2:0.1:2; [u,v]=meshgrid(u,v); x=u.^2; y=u; z=v;mesh(x,y,z) hold onx1=-2:0.1:2; y1=-2:0.1:2; [x1,y1]=meshgrid(x1,y1); z1=1-x1;mesh(x1,y1,z1) hold off6. 作出球面22222=++z y x 和柱面1)1(22=+-y x 相交的图形.解：提示：柱面1)1(22=+-y x 的参数方程：⎪⎪⎪⎩⎪⎪⎪⎨⎧=⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛==⎪⎭⎫ ⎝⎛=+=z z t t t y t t x 2cos 2sin 2sin 2cos 2)cos(12t=0:0.1:pi; r=0:0.1:2*pi;[r,t]=meshgrid(r,t); x=2*sin(t).*cos(r); y=2*sin(t).*sin(r); z=2*cos(t); mesh(x,y,z) hold onu=-pi/2:0.1:pi/2; v=-3:0.1:3;[u,v]=meshgrid(u,v); x1=2*cos(u).^2; y1=sin(2*u); z1=v;mesh(x1,y1,z1) hold off 7. 求⎰Lds z y x f ),,(，其中y x z y x f 10301),,(2++=,路径L 为: t x =,2t y =,23t z =，20≤≤t .解：首先把曲线积分化为定积分. 因为dt z y x ds t t t 222++=，clear; syms t;x= t; y=t^2; z= 3*t^2; f=sqrt(1+30*x^2+10*y);f1=f*sqrt(diff(x,t)^2+diff(y,t)^2+diff(z,t)^2); s=int(f1,t,0,2)8. 计算曲面积分⎰⎰∑++ds zx yz xy )(，其中∑为锥面22y x z +=被柱面 x y x 222=+所截得的有限部分.解：因为面积微元dxdy z z ds y x 221++=，根据曲面积分化作二重积分的计算公式, 并采用极坐标计算重积分. 注意到投影曲线x y x 222=+的极坐标方程为t r cos 2=,22ππ≤≤-t .clear;syms x y z r t f=x*y+y*z+z*x; z= sqrt(x^2+y^2); f=subs(f, 'z ', z);mj=sqrt(1+diff(z,x)^2+diff(z,y)^2); x= r*cos(t); y = r*sin(t); f=eval(f); mj=eval(mj); f1=f*mj*r;int(int(f1, r, 0, 2*cos(t)), t,-pi/2,pi/2) 输出为64/15*2^(1/2).五．其它练习1. 作函数2470722234++-+=x x x x y 及其二阶导函数在区间]7,8[-上的图形，并求函数的凹凸区间和拐点.syms xddf=diff(x^4+2*x^3-72*x^2+70*x+24,x,2) ddf =12*x^2+12*x-144>>ezplot('x^4+2*x^3-72*x^2+70*x+24',[-8,7])>>ezplot('12*x^2+12*x-144',[-8,7])>> x=-8:0.1:8;y1=x.^4+2*x.^3-72*x.^2+70.*x+24;y2=12*x.^2+12.*x-144;y3=zeros(1,length(x));plot(x,y1,'r',x,y2,'b',x,y3,'g')ddf=inline('12*x^2+12*x-144');c1=fzero(ddf,[-8,0])c2=fzero(ddf,[0,7])c1 =-4c2 =3>> ddf=inline('12*x^2+12*x-144')>> zhi1=ddf(-5),zhi2=ddf(0),zhi3=ddf(4)zhi1 =96zhi2 =-144zhi3 =96>> f=inline('x^4+2*x^3-72*x^2+70*x+24','x'); zhi4=f(-4),zhi5=f(3)zhi4 =-1280zhi5 =-279函数的下凸（向上凹）区间为[-8,-4],[3,+∞]，下凹（向上凸）区间为[-4,3];拐点(-4,-1280), (3,-279).2. 设⎪⎭⎫ ⎝⎛=-πx e x f x cos )(162，45sin )(23+=x x g . 作它们在区间],0[π上的图形. 并求方程)()(x g x f =在该区间内的近似根.x=0:0.1:pi;y1=exp(-x.^2/16).*cos(x./pi);y2=sin(x.^(3/2))+5/4; plot(x,y1,'r',x,y2,'b')clear; syms xf=inline('exp(-x^2/16)*cos(x/pi)-(sin(x^(3/2))+5/4)');x0=fzero(f,2.5)x0 =2.5411>> x0=fzero(f,3)x0 =2.9746 为所求近似根.3. 求曲面14),(22++=y x y x k 在点)2164,21,41(处的切平面方程，并把曲面和它的切平面作在同一图形里.syms x y zF='4/(x^2+y^2+1)-z';f=diff(F,x);g=diff(F,y);h=diff(F,z);x=1/4;y=1/2;z=64/21;a=eval(f);b=eval(g);c=eval(h);x=-1:0.1:1;y=-1:0.1:1;[x,y]=meshgrid(x,y);z1=a*(x-1/4)+b*(y-1/2)+64/21;z2=4*(x.^2+y.^2+1).^(-1);mesh(x,y,z1)hold onmesh(x,y,z2)4. 求函数22y x z +=在条件0122=-+++y x y x 下的极值.syms x y rg=x^2+y^2;h=x^2+y^2+x+y-1;la=g+r*h;lx=diff(la,x)ly=diff(la,y)lr=diff(la,r)输出lx = 2*x+r*(2*x+1)ly = 2*y+r*(2*y+1)lr = x^2+y^2+x+y-1输入solve('2*x+r*(2*x+1)=0', '2*y+r*(2*y+1)=0', 'x^2+y^2+x+y-1=0','x,y,r')得到输出ans =r: [2x1 sym]x: [2x1 sym]y: [2x1 sym]再分别输入rxy得到r =[ -1+1/3*3^(1/2)][ -1-1/3*3^(1/2)]x =[ 1/2*3^(1/2)-1/2][ -1/2-1/2*3^(1/2)]y =[ 1/2*3^(1/2)-1/2][ -1/2-1/2*3^(1/2)]即有解：1(33r =--,1(12x =-,1(12y =--;1(33r =-+,1(12x =-+,1(12y =-+. 因此有两个极值可疑点. 再输入x = 1/2*3^(1/2)-1/2;y = 1/2*3^(1/2)-1/2;f1=eval(g)x = -1/2-1/2*3^(1/2);y = -1/2-1/2*3^(1/2);f2=eval(g)得到输出0.26793.7321即得到两个可能是条件极值的函数值{2+3,2-3}. 但是否真的取到条件极值呢？可利用等高线作图来判断：输入[x,y]=meshgrid(-2:0.1:2,-2:0.1:2);z=x.^2+y.^2;contour(x,y,z,30)hold onezplot('x^2+y^2+x+y-1')5. 求∑∞=+-021)3(4n nn n x 的收敛域与和函数. clear;syms n x;a1=4^(2*n)*(x-3)^n/(n+1);a2=subs(a1, n, n+1);p=limit(a2/a1, n, inf)输出为p =16*x-48注意这里对a2和a1都没有加绝对值. 因此上式的绝对值小于1时，幂级数收敛，大于1时发散. 为了求出收敛区间的端点，输入x=solve('abs(16*x-48)=1')输出为x =[ 49/16][ 47/16] 由此可知16491647<<x 时收敛，1647<x 或1649>x 时发散. 为了判断端点的敛散性，输入 simplify(subs(a1, 'x', 49/16))得到x 为右端点时幂级数的一般项为ans =1/(n+1), 因此当1649＝x 时发散. 再输入 simplify(subs(a1, 'x', 47/16))输出结果为ans =1/(n+1)*(-1)^n, 因此当4716x ＝时, 级数收敛. 也可以在收敛域内求得这个级数的和函数. 输入clear;syms n xs4=symsum(4^(2*n)*(x-3)^n/(n+1), n, 0, inf )输出为 s4 =-1/(16*x-48)*log(49-16*x)。

MATLAB数学实验答案（全）第⼀次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作⼆维、三维⼏何图形，能够⽤Matlab 软件解决微积分、线性代数与解析⼏何中的计算问题。

补充命令vpa(x,n) 显⽰x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形在下⾯的题⽬中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin lim x mx mxx →∞-syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x y edxdy +??dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +? syms xint(x^4/(902^2+4*x^2))//不定积分 ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//⾼阶导数syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最⾼次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=⽤循环语句编程给出该数列的前20项（要求将结果⽤向量的形式给出）。

MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。

要求输出结果精确至小数点后两位。

答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。

要求输出结果精确至小数点后两位。

答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。

matlab数学实验考试题及答案一、选择题（每题2分，共10分）1. MATLAB中用于生成0到1之间均匀分布的随机数的函数是？A. randB. randiC. randnD. randperm答案：A2. 下列哪个命令可以计算矩阵的行列式？A. detB. rankC. eigD. inv答案：A3. MATLAB中用于求解线性方程组的命令是？A. solveB. linsolveC. fsolveD. ode45答案：A4. 在MATLAB中，如何创建一个3x3的单位矩阵？A. eye(3)B. ones(3)C. zeros(3)D. identity(3)答案：A5. MATLAB中用于绘制二维图形的函数是？A. plotB. surfC. meshD. contour答案：A二、填空题（每题3分，共15分）1. MATLAB中，使用________函数可以计算矩阵的迹。

答案：trace2. 若要在MATLAB中创建一个从1到10的向量，可以使用________函数。

答案：1:103. MATLAB中，使用________函数可以计算矩阵的特征值。

答案：eig4. 若要在MATLAB中绘制一个正弦波，可以使用________函数。

答案：sin5. MATLAB中，使用________函数可以计算矩阵的逆。

答案：inv三、简答题（每题10分，共20分）1. 描述MATLAB中如何使用循环结构来计算并打印1到100之间所有奇数的和。

答案：可以使用for循环结构，初始化一个变量sum为0，然后遍历1到100之间的每个数，使用模运算符判断是否为奇数，如果是，则将其加到sum上，最后打印sum的值。

2. 简述MATLAB中如何使用条件语句来检查一个数是否为素数，并打印出所有小于100的素数。

答案：可以使用for循环遍历2到99之间的每个数，对于每个数，使用一个while循环检查它是否有除1和它本身之外的因数，如果没有，则使用if语句判断该数是否为素数，如果是，则打印该数。

数学实验练习2.1画出下列常见曲线的图形。

（其中a=1,b=2,c=3）1、立方抛物线3xy=解：x=-5:0.1:0;y=(-x).^(1/3);y=-y;x=0:0.1:5;y=[y,x.^(1/3)];x=[-5:0.1:0,0:0.1:5];plot(x,y)2、高斯曲线2x e=y-解：fplot('exp(-x.^2)',[-5,5])3、笛卡儿曲线)3(13,1333222axy y x t at y t at x =++=+=解：ezplot('x.^3+y.^3-3*x*y',[-5,5])xyx.3+y.3-3 x y = 0或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)4、蔓叶线)(1,1322322xa x y t at y t at x -=+=+=解：ezplot('y.^2-x.^3/(1-x)',[-5,5])xyy.2-x.3/(1-x) = 0或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi;x=t-sin(t); y=2*(1-cos(t)); plot(x,y)6、星形线)(sin ,cos 32323233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3;plot(x,y)或ezplot('x.^(2/3)+y.^(2/3)-1',[-1,1])xyx.2/3+y.2/3-1 = 07、螺旋线ct z t b y t a x ===,sin ,cos 解：t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z) grid on8、阿基米德螺线θa r = 解：x =0:0.1:2*pi; r=x; polar(x,r)902701809、对数螺线θa e r = 解：x =0:0.1:2*pi; r=exp(x); polar(x,r)90270180010、双纽线))()((2cos 22222222y x a y x a r -=+=θ 解：x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-(x.^2-y.^2)',[-1,1]) grid onxy(x.2+y.2).2-(x.2-y.2) = 011、双纽线)2)((2sin 222222xy a y x a r =+=θ 解：x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-2*x*y',[-1,1]) grid onxy(x.2+y.2).2-2 x y = 012、心形线)cos 1(θ+=a r 解：x =0:0.1:2*pi; r=1+cos(x); polar(x,r)90270练习2.21、求出下列极限值。

matlab数学实验复习题(有答案)复习题1、写出3个常用的绘图函数命令2、inv （A ）表示A 的逆矩阵；3、在命令窗口健入clc4、在命令窗口健入clear 5、在命令窗口健入6、x=-1：0.2：17、det （A ）表示计算A 的行列式的值；8、三种插值方法：拉格朗日多项式插值，分段线性插值，三次样条插值。

9、若A=123456789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，则fliplr （A ）=321654987⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦A-3=210123456--⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦A ．^2=149162536496481⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦tril （A ）=100450789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦ triu （A ，-1）=123456089⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦diag （A ）=100050009⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦A(:,2),=258A(3,:)=369 10、normcdf （1，1，2）=0.5%正态分布mu=1，sigma=2，x=1处的概率[t,x]=ode45(@f,[a,b],x0),中参数的涵义是@fun 是求解方程的函数M 文件，[a,b]是输入向量即自变量的范围a 为初值，x0为函数的初值，t function 开头；1721、设x ）的功能是作出将X 十等分的直方图22、interp1([1,2,3],[3,4,5],2.5)Ans=4.523、建立一阶微分方程组⎩⎨⎧+='-='yx t y y x t x 34)(3)(2的函数M 文件。

(做不出来)二、写出运行结果：1、>>eye(3，4)=1000010000102、>>size([1,2,3])=1；33、设b=round （unifrnd （-5，5，1，4）），则=3 5 2 -5 >>[x,m]=min(b)；x=-5；m=4,[x,n]=sort(b)-5 2 3 5 4 3 1 2mean(b)=1.25，median （b ）=2.5，range （b ）=104、向量b 如上题，则>>any(b),all(b<2),all(b<6)Ans=1 0 15、>>[5 6;7 8]>[7 8;5 6]=00116、若1234B ⎡⎤=⎢⎥⎣⎦，则 7、>>diag(diag(B))=10048、>>[4:-2:1].*[-1,6]=-4 129、>>acos(0.5),atan(1)ans=1.047197551196598ans=0.78539816339744810、>>norm([1,2,3])Ans=3.74165738677394111、>>length（[1，3，-1]）=312、>>x=0:0.4:2;plot(x,2*x,’k*’)13、>>zeros(3,1);ans=14、>>ones(3)=111111111，vander([2,3,5])=421931255116、>>floor（1：0.3：3）=1 1 1 12 2 218、>>subplot(2,2,1); fplot('sin',[0,2*pi]);subplot(2,2,2);plot([1,2,-1]);>>x=linspace(0,6*pi);subplot(2,2,3);plot3(cos(x),sin(x),x);>>subplot(2,2,4);polar(x,5*sin(4*x/3));19、>>t=linespace（0，2，11）0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.020、>>[a,b]=binostat(15,0.2)a=3 b=2.4>>y1=binopdf(5,10,0.7)=0.1029,y2=binocdf(5,10,0.7)=0.15031 1 1 11 1 1 1>>y=-poissrnd(8,2,4)-16 -10 8 -7-7 -8 -6 -9>>sign(y)-1 -1 -1 -1-1 -1 -1 -135、>>[a1,b1]=binostat(20,0.4) a1=8 b1=4.8 >>[a2,b2]=poisstat(8)ans=8,8>>[a3,b3]=chi2stat(15)ans=[15 30]36、运行M文件：chi2fign=5;a=0.9;xa=chi2inv(a,n);x=0:0.1:15;y=chi2pdf(x,n);plot(x,y,'b');hold on;xf=0:0.1:xa;yf=chi2pdf(xf,n);fill([xf,xa],[yf,0],'g');text(xa*1.01,0.005,num2str(xa));text(2.5,0.05,'alpha=0.9','fontsize',20); text(9,0.09,'X~{\chi}^2(4)','fontsize',16);37、>>t=linspace(0,2*pi);>>polar(t,3*t,’g*’)38、>>quadl(’exp(2*x).*log(3*x)’,1,3)ans =398.635239、x0=0:2*pi/6:2*pi;y0=sin(x0).*cos(x0);x=[linspace(0,2*pi,100)];y=sin(x).*cos(x);y1=spline(x0,y0,x); [x;y;y1]'plot(x,y,'k',x,y1,'b-')注：此处省略100组数据40、>>A=round(unifrnd(0,100,3,3));>>[L,U]=lu(A)L =0.9897 0.4699 1.00000.1649 1.0000 01.0000 0 0U =97.0000 80.0000 92.00000 35.8041 26.82470 0 -89.656841、a=sparse([1 3 3],[2 3 5],[1 2 3],4,5);s=full(a)s =0 1 0 0 00 0 0 0 00 0 2 0 30 0 0 0 0三、编程1、 分别用矩形公式、梯形公式、辛普森公式、Gauss-Lobatto 公式及随机模拟方法计算数值积分/230sin 2x e xdx π⎰，并与符号运算计算的结果进行比较。