Matlab编程与应用习题和一些参考答案

1、.word 可编辑.专业.专注4y 7z 12w7y 4z 2wx 8z 5w 96x 5y 2z 10w a=3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10; b=4;4;9;4; c=ab28,求 3= A*B ;C2=A *B;C3=A.*并,求上述所7有方阵的逆阵。 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; 3= A*B C2=A*B C3=A.*B in v(C1) in v(C2) in v(C3)3sin x5.设 y cosx 0.52,把 x=02 n

2、间分为 101 点，画出以 x 为横坐标,y为纵坐标的(1 x )曲线。 x=li nspace(0,2*pi,101); y=cos(x)*(0.5+(1+x.A2)3*si n(x); plot(x,y,r)Matlab上机实验一3.求下列联立方程的解3x5x1213.word 可编辑.专业.专注6.产生 86 阶的正态分布随机数矩阵 R1,求其各列的平均值和均方差。并求该矩阵全体数的 平均值和均方差。(mean var)a=ra ndn( 8,6)mean(a)var(a)k=mea n(a)k1=mea n(k)i=o nes(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=

3、mea n(i3)g2=mea n(g)或者u=reshape(a,1,48);p1=mea n(u)p2=var(u)7.设 x=rcost+3t,y=rsint+3, 分别令 r=2,3,4，画出参数 t=010 区间生成的 xy 曲线。 t=li nspace(0,10); r1=2;.word 可编辑.专业.专注 x1=(r1*cos(t)+3*t); y1=r1*si n( t)+3; r2=3; x2=(r2*cos(t)+3*t); y2=r2*si n( t)+3; r3=4; x3=(r3*cos(t)+3*t); y3=r3*si n( t)+3; plot(x1,y1,r

4、,x2,y2,b,x3,y3,m)8.设 f(x)=x5- 4x4+3x2- 2x+ 6(1) 在 x=-2,8之间取 100 个点，画出曲线，看它有几个过零点。(提示：用 polyval 函 数) x=li nspace(2,8,100); y=polyval(1 0 -4 3 -2 6,x); plot(x,y,b,x,0,y)(2) 用 roots 函数求此多项式的根。t=1 0 -4 3 -2 6p=roots(t)9建立一个字符串向量，然后对该向量做如下处理：(1) 取第 15 个字符组成的子字符串。(2) 将字符串倒过来重新排列。.word 可编辑.专业.专注(3) 统计字符串中小

5、写字母的个数 a=abcdef123456; subch=a(1:5)subch =abcde revch=a(e nd:-1:1)revch =654321fedcba k=fi nd(a=a &a a(k)=a(k)-(a-A); char(a) len gth(a) len gth(k)10.利用帮助查找 limit 函数的用法，并自己编写，验证几个函数极限的例子limit(si n(x)/x)returns 1limit(x-2)/(xA2-4),2)returns1/4limit(1+2*t/x)A(3*x),x,i nf)returnsexp(6*t)limit(1/x,x

6、,0,right)returnsinflimit(1/x,x,0,left)returns-inflimit(si n(x+h)-si n(x)/h,h,0)returnscos(x)v = (1 + a/x)Ax, exp(-x);.word 可编辑.专业.专注Matlab 上机实验二1. 验证所授课程的课件。2.设 x=sint, y=sin(nt+a) ,(1) 若 a=1 ,令 n =1,2,3,4,在四个子图中分别画出其曲线。(2) 若 n=2,取 a=0,n/3，冗/2，及 n,在四个子图中分别画出其曲线。(1)a=1;x=si n( t);y1=si n(1*t+a);y2=si

7、 n(2*t+a);y3=si n(3*t+a);y4=si n(4*t+a);subplot(2,2,1);plot(x,y1);subplot(2,2,2);plot(x,y2);subplot(2,2,3);plot(x,y3);subplot(2,2,4);plot(x,y4)n=2;x=si n( t);y1=si n(2*t+0);y2=si n(2*t+pi/3);y3=si n(2*t+pi/2);y4=si n(2*t+pi);subplot(2,2,1);plot(x,y1);subplot(2,2,2);plot(x,y2);subplot(2,2,3);plot(x,y

8、3);subplot(2,2,4); plot(x,y4)2 2、3.绘制z x2e(x y)在定义域 x=-2,2,y=-2,2内的曲面。(利用 meshgrid )x=-2:2;y=x;X,Y= meshgrid(x,y);Z=XA2*exp(-(XA2+YA2);mesh(X,Y,Z)4.求代数方程 3X5+4X4+7X3+2X2+9X+12=0 的所有根。(利用 roots 函数)p=3 4 7 2 9 12;roots(p)5.把 1 开五次方,并求其全部五个根。（提示：解 x5-1=0）p=1 0 0 0 0 -1;roots(p).word 可编辑.专业.专注6,y(x24x 3

9、)/2x21123 x 1x 1,画出在1 x 33,3上的曲线.word 可编辑.x1=-3:0.01:-1;y 1= (-x1.A2-4.*x1-3)/2; x2=-1:0.01:1;y2=-x242+1;x3=1:0.01:3;y3=(-x3.A2+4.*x3-3)/2;plot(x1,y1);hold on; plot(x2,y2);hold on; plot(x3,y3)Matlab 上机实验三1. 设】 - :求 syms x y=x*exp(2*x);Diff(y,10)2. 设 y =总曲.十 閃/龙,求 Jsyms x y=si n( x)A4+cos(x)A4; yAl04

10、 _ 444-3. 设盅=y 肌总扯十 K佃y,求柚 syms x y z=yA4*si n(x)A4+xA4*cos(yF4; diff( z,x)14.2-2-dx(x 1)(x x) syms x y y=1/(xA2+1)*(xA2+x); in t(y).word 可编辑. in t(sqrt(9-xA2)A3)/xA6)专业.专注.word 可编辑.专业.专注2dx131X X syms x in t(1/(x+xA3),1,2) simple(a ns) syms n symsum( n+2)/2A n,n-1,i nf) simple(a ns)8.利用函数 int 计算二重不

11、定积分 syms x y in t(i nt(x+y)*exp(-x*y)9.设方程(1)x3cos(a) 0，( 2)x3cos(x) 0及(3)x3cos(ax) 0，用符号运算工具箱函数分别求 x 的解。用 a=0.5 代入，求 x 的数值解，并与用 roots 函数所求的结果进行比较。 syms x1 x2 x3 x1=solve(x1A3+cos(0.5)=0); x2=solve(x2A3+cos(x2)=0); x3=solve(x3A3+cos(0.5*x3)=0); x1,x2,x3 p=1 0 0 cos(0.5); roots(p)7求级数的和n 12n（利用帮助，查找

12、symsum ）(x y)exydxdysolve.word 可编辑.专业.专注10.求微分方程的解：y-y=x , y(2)=4 , y(1)=7 , y(1)=8 syms x y dsolve(D3y-D2y=x,D2y(2)=4,Dy(1)=7,y(1)=8)11.设x zsin3z, y zcos3z,要求在 z=010 区间内画出 x,y,z 三维曲线。 z=0:0.01:10; x=z. *si n( 3*z); y=z.*cos(3*z); plot3(x,y,z)2 2x1 2e八,画出定义域 x=-2,2,y=-2,2 x,y=meshgrid(-2:0.1:2); z=x

13、.A2.*exp(-x.A2-y.A2); mesh(x,y,z)13.设 z=0.05x-0.05y+0.1; 画出 z 的曲面(平面)图。 X,Y=meshgrid(1:10:200,1:10:200); Z=0.05*X-0.05*Y+0.1; surf(X,Y,Z)For 循环: s=0;1 .编写程序计算：(1) 使用一个 for/end 循环，但不能使用数组运算符 sum。(2) 使用数组运算符和 sum ,但不能使用任何 for/end 循环答：(1): 1100 相加12.设z内的曲面图。.word 可编辑.专业.专注for i=1:100s=s+i;ends(2): 1100

14、 相加Sum 函数: i=1:100; s=sum(i)2.编制 m 文件，等待键盘输入，输入密码 123 ,密码正确，显示输入密码正确，程序结束; 否则提示，重新输入，直到正确为止。答:function mypasswordk=0;n=input（等待键盘输入n）;while（k=0）if n=123disp（输入密码正确）k=1;elsen=input（密码错误，重新输入n）;k=0;endend3. 编制 m 文件，输入 n (正整数)，显示所有小于 n 的质数。.word 可编辑.专业.专注答：function zhishum=input(输入正整数 m=n);p=2:m-1;for

15、i=2:sqrt(m)n=fin d(rem(p,i)=0&p=i);P( n)=；endp4.编写一个 M 函数文件 fun_es(x)，计算如下函数：y 0.5ex/3x2sin(x)，其中参数可以为 标量，也可以为向量。在 Matlab 指令窗口键入如下命令检验该函数：fun_es(3)答: function fun_es(x)y=0.5*exp(x/3)-x.A2.*si n(x)5.利用 rand 产生 10 个随机数，利用 for 循环对其进行排序(从大到小)。答: function pailien=ran d(1,10)for j=1:10for i=j+1:10if n

16、(j) n(i)a=n(j);n (j)=n(i);.word 可编辑.专业.专注n (i)=a;endendenddisp(重新排列后:)n6.编写 m 函数文件求半径为 r 的圆的面积和周长答：function yua n(r)s=pi*rA2c=2*pi*r7.数组1 : 99中非素数（即质数）元素的和答：fun cti on feizhishuhem=99;p=2:m;for i=2:sqrt(m).word 可编辑.专业.专注n=fin d(rem(p,i)=0&p=i);P(n )=；ends=sum(1:99)-sum(p)8.从键盘输入若干个数(每个数用回车结束)，当输

17、入 0 时结束输入，求这些数的平均值以及平方和。(结束标志 0 不计算在内)答：function abck=0;s=0;i=0;y=0;n=in put();while(k=0)if n=0k=0;i=i+1;s=s+n;y=y+nT;n=i nput();elsek=1;a=s/i.word 可编辑.专业.专注yendend9.输入 x,y 的值，并将它们的值互换后输出答：function jiaohuanx=i nput(x=n);y=i nput(y=n);z=x;x=y;y=z;xy10. 下面程序实现什么功能？for n=100:200if rem( n, 77)=0con ti nueendbreak.word 可编辑.专业.专注end答：在 100 到 200 之间输出第一个能被 77 整除的整数x=i nput(x=n);n=in put( n=n);y=