matlab习题与答案,Matlab编程与应用习题和一些参考答案.doc

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

Matlab上机实验一、二1 安装Matlab软件。2 验证所学内容和教材上的例子。3.求下列联立方程的解 a3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10; b4;4;9;4; cabc 5.2226 4.4570 1.4718 1.59944设求C1A*B;C2A*B;C3A.*B,并求上述所有方阵的逆阵。 A1 4 8 13;-3 6 -5 -9;2 -7 -12 -8; B5 4 3 -2;6 -2 3 -8;-1 3 -9 7; C1A*BC1 19 -82 30 12 27 3 -38 54 29 C2A*BC2 -15 16 -24 36 63 -17 93 -105 22 6 117 -60 19 46 84 -10 C3A.*BC3 5 16 24 -26 -18 -12 -15 72 -2 -21 108 -56 invC1ans 0.0062 0.0400 -0.0106 -0.0046 0.0169 0.0030 0.0168 0.0209 0.0150 invC2Warning Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND 8.997019e-019.ans 1.0e015 * -0.9553 -0.2391 -0.1997 0.2700 0.9667 0.2420 0.2021 -0.2732 -0.4473 -0.1120 -0.0935 0.1264 -1.1259 -0.2818 -0.2353 0.3182 invC3 Error using invMatrix must be square.5设 ，把x02间分为101点，画出以x为横坐标，y为纵坐标的曲线。 xlinspace0,2*pi,101; ycosx*0.51x.23*sinx; plotx,y,r6产生86阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 (mean var)arandn8,6meanavarakmeanak1meankiones8,6i1i*k1i2a-i1i3i2.*i2gmeani3g2meang或者ureshapea,1,48;p1meanup2varu7.设xrcost3t,yrsint3,分别令r2,3,4，画出参数t010区间生成的xy曲线。 tlinspace0,10; r12; x1r1*cost3*t; y1r1*sint3; r23; x2r2*cost3*t; y2r2*sint3; r34; x3r3*cost3*t; y3r3*sint3; plotx1,y1,r,x2,y2,b,x3,y3,m8.设fxx5- 4x4 3x2- 2x 61 在x-2,8之间取100个点，画出曲线，看它有几个过零点。(提示用polyval 函数) xlinspace2,8,100; ypolyval1 0 -4 3 -2 6,x; plotx,y,b,x,0,y2 用roots函数求此多项式的根。t1 0 -4 3 -2 6prootst9.建立一个字符串向量，然后对该向量做如下处理1 取第15个字符组成的子字符串。2 将字符串倒过来重新排列。3 统计字符串中小写字母的个数。10.利用帮助查找limit函数的用法，并自己编写，验证几个函数极限的例子。Matlab上机实验二1.验证所授课程的课件。2.设xsint, ysinnta，1 若a1，令n 1,2,3,4,在四个子图中分别画出其曲线。2 若n2,取a0，/3，/2，及，在四个子图中分别画出其曲线。1a1; xsint;y1sin1*ta;y2sin2*ta;y3sin3*ta;y4sin4*ta; subplot2,2,1;plotx,y1;subplot2,2,2;plotx,y2;subplot2,2,3;plotx,y3;subplot2,2,4;plotx,y42n2;xsint;y1sin2*t0;y2sin2*tpi/3;y3sin2*tpi/2;y4sin2*tpi;subplot2,2,1;plotx,y1;subplot2,2,2;plotx,y2;subplot2,2,3;plotx,y3;subplot2,2,4;plotx,y4注本题好像题目本身就有问题，因为“t” 取值不明，所以运行不了，画不出图来。转化为符号式来画图也是不可行的。(纯粹个人意见，可能我想错了方向，会做的同学请上传一下正确的做法)3绘制在定义域x-2,2,y-2,2内的曲面。(利用meshgrid)x-22;yx;X,Y meshgridx,y;ZX2*exp-X2Y2;meshX,Y,Z 4求代数方程3x54x47x32x29x120的所有根。(利用roots函数)p3 4 7 2 9 12;rootsp5把1开五次方，并求其全部五个根。提示解x5-10p1 0 0 0 0 -1;rootsp6画出在上的曲线。x1-30.01-1;y1-x1.2-4.*x1-3/2; x2-10.011;y2-x2.21; x310.013;y3-x3.24.*x3-3/2;plotx1,y1;hold on; plotx2,y2;hold on; plotx3,y3 Matlab上机实验三1.设,求xsymx; yx*exp2*x;diffy,202.设,求xsymx; ysinx4cosx4;diffy,103.设,求xsymx; ysymy;zy4*sinx4x4*cosy4;intz,x4xsymx;int1/x21*x2x,x5xsymx;intsqrt9-x23/x6,x6xsymx;intxx3,x,1,27.求级数的和 symn;symsumn2/2n,n,1,inf8.利用函数int计算二重不定积分xsymx; ysymy;zxy*exp-x*y;aintz,x;inta,yMatlab上机实验四1. 编写一个M函数文件计算下列级数的和x和n的值由指令给出。test.mfunction ftest x the xn; n the nn; j1; f20;for i1n; j1;for k1i; jj*k; f1xi/j;endf2f2f1;endff21;2编制m文件，等待键盘输入，输入密码123，密码正确，显示输入密码正确，程序结束；否则提示，重新输入。方法一function fcheckpasswordPlease the passwordn;if password123 disp输入密码正确else check;end方法二function pwn0;k0;n请输入密码整数n;whilek0if n123 disp输入正确 k1;else k0; n密码错误，请重新输入密码整数n;endend3编制m文件，输入n(正整数)，显示所有小于n的质数。第一m文件(用于判断某数是否为素数)function leappanduannleap1;for i2sqrtn1 if modn,i0 leap0; break; endend第二个m文件(用于读进n个数并调用panduan.m判断读数是否为素数，并把是素数的显示出来)function zhishunfor m2n if panduanm1 dispm; endend4编写一个M函数文件fun_esx，计算如下函数，其中参数可以为标量，也可以为向量。在Matlab指令窗口键入如下命令检验该函数 fun_es3function fun_esxy0.5.*expx./3-x.2.*sinx5下面程序实现什么功能for n100200if remn,770 continueendbreakendn 求100到200之间，第一个77的整倍数6编写m函数文件求半径为r的圆的面积和周长。function fcirclerSpi*r2L2*pi*r7数组199中素数(即质数)元素的和第一m文件(用于判断某数是否为素数)function leappanduannleap1;for i2sqrtn1 if modn,i0 leap0; break; endend第二个m文件(用于读进199并调用panduan.m判断读数是否为素数，并把是素数的加起来，最后输出)function totalnk0;for m1n if panduanm1 kkm; endendk8从键盘输入若干个数，当输入0时结束输入，求这些数的平均值以及平方和。 求n个数的平均值及平方和function pingfang1 l0; total10; total20; a an; whilea0 ll1; total1total1a; total2total2a2; a; end avgtotal1/l total2 9输入x,y的值，并将它们的值互换后输出。function fchangex,yx xn;y yn;ax;xy;ya;x,y10利用rand产生10个随机数，利用for循环对其进行排序(从大到小)。disp排列前arand10,1for i110for j110-iif aj1,1aj,1baj,1;aj,1aj1,1;aj1,1b;endendenddisp排列后aMatlab上机实验四1. 编写一个M函数文件计算下列级数的和x和n的值由指令给出。function yjiechengn 定义 jiecheng函数s1;for i1n ss*i;end ys;function zy1x请输入x值;n请输入n值;m1;i0n;for mmx.i/jiechengi; 调用上面的公式end m2编制m文件，等待键盘输入，输入密码123，密码正确，显示输入密码正确，程序结束；否则提示，重新输入。function fzy2t keykey is 123 if t123 goodelse Your are fool, again zy2; 从新输入；end3编制m文件，输入n(正整数)，显示所有小于n的质数。function fmyzy31nksqrtn;flag1;for i2k; if modn,i0 flag0; break; endendfflag;判断一个数字是否是素数；function myzy32nfor i2nif myzy31i1调用myzy31，判断是否是素数，如果是就输出； i else i; i后面有“；”就是不输出了； endend4编写一个M函数文件fun_esx，计算如下函数，其中参数可以为标量，也可以为向量。在Matlab指令窗口键入如下命令检验该函数 fun_es3function ffun_esnf0.5.*expn./3-n.*n.*sinn;use the dot not only adapt to numbers but also matrix 5下面程序实现什么功能for n100200if remn,770 continueendbreakendn实现输出在100到200之间第一个被77整除的数字，就是154；6编写m函数文件求半径为r的圆的面积和周长。function fmyzy6r r dispThe aear is dispr*r*pidispThe perimeter isdisp2*pi*r7数组199中非素数(即质数)元素的和在这个函数中令n99即可；function myzy32nsum0；function myzy32nfor i2nif myzy31i1调用myzy31，判断是否是素数，如果是就输出； i sumsumi;求和 else i; i后面有“；”就是不输出了； endendsum求和输出；for i2nif myzy31i1调用myzy31，判断是否是素数，如果是就输出； i yy1; else i; i后面有“；”就是不输出了； endend8从键盘输入若干个数，当输入0时结束输入，求这些数的平均值以及平方和。function myzy8sum0;sum20;n2;while n0 k the numberend with 0 sumksum; sum2k.2sum2; nk;标志位传递；endsumsum2 输出；9输入x,y的值，并将它们的值互换后输出。function myzy8xPlease the x yPlease the y tx; 保护数据;xy;yt;disp“after change ”xy 输出;10利用rand产生10个随机数，利用for循环对其进行排序(从大到小)。function myzy10 输入10即可；k the number of randn numbers xrandn1,ktemp0;for j1k-1; for ik-jk-1; if xixi1 txi; xixi1; xi1t; else end endendx all by 穿了袜子才变身