MATLAB程序设计与应用(刘卫国编)课后实验答案

实验MATLAB运算基础1、先求下列 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 达式得值，然后显示MATLAE工作空间得使用情况并保存全部变量(1)(2)，其中(3)Z30.3a0.3aeesin(a20.3)In0.3a23.0,2.9,川,2.9,3.0(4)，其中t=0:0、5:2、5解：M文件：zl=2*sin(85兴pi/180)/(1+exp(2))x=[21+2兴i;->455];z2=1/2*log(x+sqrt(1+x人2))a=—3、0:0、1:3、0;z3=(exp(0、3、*a)-exp(—0、3、*a))、/2、*sin(a+0、3)+1og((0、3+a)、/2)t=0:0、5:2、5;z4=(t〉=0&t<1)、*(t、八2)+(t>=1&tv2)、*(t、人2—1)+(t>=2&t<3)、*(t、八2-2*t+1)运算结果：z1=2*sin(85*pi/180)/(1+cxp(2))x=[21+2*i;—、455];z2=1/2兴log（x+sqrt（1+xA2））a=—3、0:0、1:3、0;z3=（exp（0、3、兴a）-exp（—0、3、*a））、/2、兴sin（a+0、3）+log（（0、3+a）、/2）t=0:0、5:2>5;z4=（t>=0&tv1）、*（t、A2）+（t>=1&t<2）、兴（t、A2—1）+（t>=2&t〈3）、兴（t、A2—2兴t+1）z1=0、2375z2=0、7114—0、0253i0、8968+0、3658i0、2139+0、9343i1、1541—0、0044iz3=Columns1through40、7388+3、1416i0、7696+3、1416i0、7871+3、1416i0、7913+3、1416iColumns5through80、7822+3、1416i0、7602+3、1416i0、7254+3、1416i0、6784+3、1416iColumns9through120、6196+3、1416i0、5496+3、1416i0、4688+3、1416i0、3780+3、1416iColumns13through160、2775+3、1416i0、1680+3、1416i0、0497+3、1416i—0、0771+3、1416iColumns17through20—0、2124+3、1416i—0、3566+3、1416i—0、5104+3、1416i—0、6752+3、1416iColumns21through24—0、8536+3、1416i—1、0497+3、1416i-1、2701+3、1416i-1、5271+3、1416iColumns25through28—1、8436+3、1416i-2、2727+3、1416i—2、9837+3、1416i—37、0245Columns29through32—3、0017—2、3085—1、8971-］、5978Columns33through36—1、3575—1、1531-0、9723—0、8083Columns37through40-0、6567—0、5151—0、3819-0、2561Columns41through44-0、1374—0>02550、07920、1766Columns45through480、26630、34780、42060、4841Columns49through520、53790、58150、61450、6366Columns53through560、64740、64700、63510、6119Columns57through600、57770、53270、47740、4126Column610、3388z4=00、250001、25001、00002、25002、已知：求下列表达式得值：(1)A+6*E与A-B+I(其中I为单位矩阵)(2)A兴E与A、*BAA3与A、人3A/E及B\A(5):A,B]与:A([1,3]B人2]解：M文件：A=[1234-4;34787;3657];B=[13—1;203;3—27];A+6、*BA-B+eye(3)A*BA、*BAA3A、A3A/BBA[A,B]:A([1,3],:);BA2]运算结果：A=[1234-4;34787;3657];B=[13-1;203;3—27];A+6、兴BA—B+eye(3)A*BA、兴BAA3A、A3A/BB\A:A,B]:A([1,3]，:);BA2]ans=1852-10467105215349ans=1231—3328840671ans=684462309—72596154-5241ans=1210246802619-13049ans=372262338244860424737014918860076678688454142118820ans=172839304—6427274625343ans=16、4000-13、60007、600035、8000-76、200050、200067、0000-134、000068、0000ans=109、4000—131、2000322、8000—53、000085、0000—171、0000—61、600089、8000—186、2000ans=1234-413—13478720336573—27ans=1234-436574511101920-540393043436585033、设有矩阵A与B求它们得乘积C.将矩阵C得右下角3X2子矩阵赋给D查瞧MATLAB工作空间得使用情况.解：、运算结果:E=(reshape(1:1:25,5,5))';F=[31016;17-69;023-4;970;41311];C=E*FH=C(3:5,2:3)C=9315077258335237423520397588705557753890717H=5203977055578907174、完成下列操作：(1)求［100,999］之间能被21整除得数得个数。(2)建立一个字符串向量，删除其中得大写字母解:(1)结果：m=100:999;n=find(mod(m,21)==0);length(n)ans=43(2)、建立一个字符串向量例如：ch='ABC123d4e56Fg9'则要求结果就是:ch二'ABC123d4e56Fg9‘；k=find(ch>='A'&h<='Z');ch(k)=口Ch=123d4e56g9实验二MATLAB矩阵分析与处理1、设有分块矩阵，其中E、R、0、S分别为单位矩阵、随机矩阵、零矩阵与对角阵，试通过数值计算验证。解：M文件如下；B圧<>-SH」00ah32OOOo21i1y0IIo42711OoooooOyooo_kyoooooooo_kyooooo00,53830,4oo—kyo—kyooo—kyoooooyo767o2oyo002996_k—ky32OOOOy96—k9oy§7001、00000、15642、88570001、0000000004、0000ans=0000000000000000000000000由ans，所以2、产生5阶希尔伯特矩阵H与5阶帕斯卡矩阵P,且求其行列式得值Hh与Hp以及它们得条件数Th与Tp，判断哪个矩阵性能更好。为什么？解：M文件如下：输出结果：H=1、00000、50000、33330、25000、20000、50000、33330、25000、20000、16670、33330、25000、20000、16670、14290、25000、20000、16670、14290、12500、20000、16670、14290、12500、1111P=111111234513610151410203515153570Hh=3、7493e—012Hp=1Th=4、7661e+005Tp=8、5175e+003因为它们得条件数Th>>Tp，所以pascal矩阵性能更好3、建立一个5X5矩阵，求它得行列式值、迹、秩与范数解：M文件如下：输出结果为A=17241815235714164613202210121921311182529d=507000065cl=6、8500c2=5、4618cinf=6、8500已知求A得特征值及特征向量，并分析其数学意义。解：M文件如图：输出结果为V=0、71300、28030、2733-0、6084—0、78670、87250、34870、55010、4050D=25、3169000—10、518200016、8351数学意义：V得3个列向量就是A得特征向量，D得主对角线上3个就是A得特征值,特别得,V得3个列向量分别就是D得3个特征值得特征向量•5、下面就是一个线性方程组：求方程得解.将方程右边向量元素b3改为0、53再求解，并比较b3得变化与解得相对变化。计算系数矩阵A得条件数并分析结论。解：M文件如下：输出结果：X=1、20000、60000、6000X2=1、20000、60000、6000C=1、3533e+003由结果,X与X2得值一样，这表示b得微小变化对方程解也影响较小，而A得条件数算得较小，所以数值稳定性较好，A就是较好得矩阵。6、建立A矩阵，试比较sqrtm(A)与$口rt(A),分析它们得区别。解：M文件如下：EditorUiitilled2*运行结果有:166TOC\o"1-5"\h\z20598b1=3、8891-0、3、29170、3855b2=4、00004、47213、0000b=16、000020、00009、00001812511023、2、14362、07602、44952、23612、82846、00005、00008、000021030、36981、73054、24263、46412、236118、000012、00005、0000分析结果知：sqrtm(A)就是类似A得数值平方根(这可由b1衣b1=A得结果瞧出)而sqrt(A)则就是对A中得每个元素开根号，两则区别就在于此。实验三选择结构程序 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 、实验目得1、掌握建立与执行M文件得方法。2、掌握利用if语句实现选择结构得方法•3、掌握利用switch语句实现多分支选择结构得方法4、掌握try语句得使用。、实验 内容 财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容 1、求分段函数得值。用if语句实现，分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时得y值.解:M文件如下：运算结果有：f(—5)y=14>f(—3)y=11>f(i)y=2>f(2)y=1>>f(2、5)y=-0、2500>>f(3)y=5>>f(5)y=192、输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分〜100分为A,80分~89分为B,79分~79分为C,60分~69分为D,60分以下为E.要求：(1)分别用if语句与switch语句实现。(2)输入百分制成绩后要判断该成绩得合理性，对不合理得成绩应输出出错信息解:M文件如下试算结果:score=88grade=Bscore=123错误:输入得成绩不就是百分制成绩3、硅谷公司员工得工资计算方法如下：工作时数超过120小时者，超过部分加发15%。工作时数低于60小时者，扣发700元。其余按每小时84元计发。试编程按输入得工号与该号员工得工时数，计算应发工资。解：M文件下4、设计程序,完成两位数得加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应得运算，并显示相应得结果。解：M文件如下;运算结果例：a=38b=3输入一个运算符：Ac=falsea=92b=0输入一个运算符：+c=51325、建立5X6矩阵,要求输出矩阵第n行元素。当n值超过矩阵得行数时，自动转为输出矩阵最后一行元素，并给出出错信息。解：M文件如下：Fil*EditTtxtGoCtllTo&LsDesktopWindowHtlp;aj必肴启巧4广口◎勻"・•■*耳Tx憑窗iI0eiou121314IE17elcAfinput「输人一亍5行网I]矩阵A);n=inputC儒人—正團由戸‘);ifn<50disp(Atrij:)),elsedi?p(A(5,:))；endscriptLn7Cel18OVR运算结果如下：输入一个5行6列矩阵234567]输入一正整数n=4A=[155;234576;222223;1123973;112输入一个5行6列矩阵3A=[175；2334576;222223;1123973;234输入一正整数567]n=6ansErrorusing==>manyinpudisptarguments、lasterr实验四循环结构程序设计、实验目得1、掌握利用f。r语句实现循环结构得方法.2、掌握利用wh1le语句实现循环结构得方法。3、熟悉利用向量运算来代替循环操作得方法。、实验内容1、根据，求n得近似值。当n分别取100、1000、10000时,结果就是多少？要求:分别用循环结构与向量运算（使用sum函数）来实现。解：M文件如下：运行结果如下：K〉＞%循环结构计算p1值y=0;n=input('n=');fori=1ny==y+1/i/i；endpi=sqrt(6*y)n=100Pi=3、1321n=100DPi=3、1406n=10000Pi=3、1415%向量方法计算Pi值n=input('n=：i=1>/(1:n)、A2；s=sum(i);Pi=sqrt(6*s)n=100Pi=3、1321n=1000Pi=3、1406n=10000Pi=3、14152、根据,求：(1)y<3时得最大n值.57⑵与（1）得n值对应得y值。解:M—文件如下：运行结果如下：K〉>y=0;n=0;whiley〈3n=n+1;y=y+1/(2*n-1);endynify>3n=n—1;endny=3、0033n=n=3、考虑以下迭代公式:其中a、b为正得学数。(1)编写程序求迭代得结果，迭代得终止条件为丨Xn+1—Xn|<10—5,迭代初值X0=1、0,迭代次数不超过500次.(2)如果迭代过程收敛于r,那么r得准确值就是，当(a,b)得值取(1,1)、(8,3)、(10,0、1)时，分别对迭代结果与准确值进行比较。解：M文件如下：运算结果如下;请输入正数a=1请输入正数b=1X=0、6180r=0、6180—4、7016r=0、6180-1、6180s=-0、0000—2、2361请输入正数a=8请输入正数b=3X=1、7016r=1、7016—1、6180r=1、7016—4、7016s=0.0-6、4031请输入正数a=10请输入正数b=0、1x=3、1127r13、1127—4、7016r=3、1127—3、2127s=-0、0000—6、32544、已知求f1〜f100中:(1)最大值、最小值、各数之与。(2)正数、零、负数得个数。解:M-文件以下就是运算结果：max(f)=437763282635min(f)=—899412113528sum(f)=-742745601951cl=49c2=2c3=495^若两个连续自然数得乘积减i就是素数，则称这两个边疆自然数就是亲密数对:该素数就是亲密素数•例如，2X3-1=5，由于5就是素数，所以2与3就是亲密数,5就是亲密素数。求［2,50］区间内：(1)亲密数对得对数•(2)与上述亲密数对对应得所有亲密素数之与。解：M文件:运算结果为：j=29s=23615实验五函数文件一、实验目得1、理解函数文件得概念.2、掌握定义与调用MATLAB函数得方法。二、实验内容1、定义一个函数文件，求给定复数得指数、对数、正弦与余弦，并在命令文件中调用该函数文件。解：M文件如下:函数fushu、M文件：function[e,l,s,c]=fushu(z)%fushu复数得指数，对数，正弦，余弦得计算%e复数得指数函数值%1复数得对数函数值%s复数得正弦函数值%c复数得余弦函数值e=exp(z)；l=1og(z)；s=sin(z);c=cos(z)；命令文件M:z=input('请输入一个复数z=/);[a,b,c,d]=fushu(z)运算结果如下：z=input('请输入一个复数z=');[a,b,c,d]=fushu(z)请输入一个复数z=1+ia=1、4687+2、2874ib=0、3466+0、7854ic=1、2985+0、6350id=0、8337—0、9889i2、一物理系统可用下列方程组来表示：从键盘输入m1、m2与B得值，求a1、a2、N1与N2得值.其中g取9、8,输入B时以角度为单位。要求：定义一个求解线性方程组AX=B得函数文件，然后在命令文件中调用该函数文件。解：M文件函数fc、M文件：functionX=fc(A,B)%fcfc就是求解线性方程得函数%AA就是未知矩阵得系数矩阵X=A\B;命令M文件：c1c;ml=input(/输入m1=');m2=inputC输入m2=/);theta=inputC输入theta=/);x=theta*pi/180;g=9、8;A=[m1*cos(x)—m1-sin(x)0m1*sin(x)0cos(x)00m2-sin(x)000—cos(x)1];B=[0;m1*g;0;m2衣g];X=fc(A,B)运算结果：输入m1=1输入m2=1输入theta=30X=7、84003、39486、789615、68003、一个自然数就是素数，且它得数字位置经过任意对换后仍为素数。例如13就是绝对素数。试求所有两位绝对素数。要求：定义一个判断素数得函数文件。解：M文件：函数prime、m文件function[p]=prime(p)%输入p得范围，找出其中得素数m=p(length(p));fori=2:sqrt(m)n=find(rem(p,i)==0&p〜=i);p(n)=[:;%将p中能被i整除,而却不等于i得元素，即下标为n得元素剔除，其余得即为素数endp;命令文件：clc；p=10:99;p=prime(p);%找出10到99内得所有素数p=l0*rem(p,10)+(p—rem(p,10))/10;%将卩素数矩阵每个元素个位十位调换顺序p=prime(p)%再对对换后得素数矩阵找出所有得素数运算结果:p=1131711373173797794、设，编写一个MATLAE函数文件fx、m,使得调用f(x)时，x可用矩阵代入，得出得f(x)为同阶矩阵•解：函数fx、m文件：functionf=fx(x)%fxfx求算x矩阵下得f(x)得函数值A=0、1+(x-2)、人2;B=0、01+(x-3)、A4;f=1、/A+1、/B;命令文件：c1c;x=input('输入矩阵x=/);f=fx(x)运算结果:>>x=input('输入矩阵x=/);f=fx(x)输入矩阵x=[72;125:f=0、043710、99010、01010、17245、已知当f(n)=n+10ln(n2+5)时，求y得值.当f(n)=1X2+2X3+3X4+、、、+nX(n+1)时，求y得值.解:⑴函数f、m文件：functionf=f(x)f=x+10*1og(xA2+5);命令文件：clc;nl=input(/n1=');n2=input('n2=z)；n3=input('n3=');y1=f(n1);y2=f(n2);y3=f(n3);y=yl/(y2+y3)运算结果如下：n1=40n2=30n3=20y=0、6390(2)、—函数g、m文件functions=g(n)fori=1:ng(i)=i*(i+1);ends=sum(g);命令文件：clc;nl=input(zn1=/);n2=input('n2=');n3=input('n3=');yl=g(nl)；y2=g(n2);y3=g(n3);y=y1/(y2+y3)运算结果如下：n1=40n2=30n3=20y=1、7662实验六高层绘图操作一、实验目得1、掌握绘制二维图形得常用函数。2、掌握绘制三维图形得常用函数。3、掌握绘制图形得辅助操作。、实验内容1、设，在x=0〜2n区间取101点，绘制函数得曲线。解：M文件如下：clc；x=linspace(0,2*pi,101)；y=(0、5+3*sin(x)、/(1+x、人2));plot(x,y)运行结果有：2、已知y仁x2,y2=cos(2x),y3=y1Xy2，完成下列操作:在同一坐标系下用不同得颜色与线型绘制三条曲线。以子图形式绘制三条曲线。分别用条形图、阶梯图、杆图与填充图绘制三条曲线•解：(1)M文件：clc;x=-pi:pi/100:pi;y1=x、A2;y2=cos(2*x);y3=y1、*y2;plot(x,y1,'b-'x,y2,'r:',x,y3，‘k—)运行结果:(2)M文件:clc;x=—pi:pi/100:pi;yl=x、人2；y2=cos(2*x);y3=y1、*y2；subplot(1,3,1)plot(x,y1,'b-');tit1e('y1=x人2');subplot(1,3,2)Plot(x,y2,'r:');title('y2=cos(2x)');subpot(1,3,3);plot(x,y3，‘k——')；title('y3=y1*y2/)；、运行结果:(3)M文件：clc;x=-pi:pi/100:pi；yl=x、人2；y2=cos(2*x);y3=y1、*y2;subplot(2,2,1）；plot(x,y1,'b-_/,x,y2,'r:',x,y3,'k—；subplot(2,2,2)bar(x,y1,'b')title('y1=xA2');subp1ot(2,2,3)；bar(x,y2,'rz);title('y2=cos(2x)')；subplot(2,2,4);bar(x,y3，/k');tit1e('y3=y1*y2')；由上面得M文件，只要依次将"bar"改为"stairs”、“stem”、“fill"，再适当更改区间取得点数，运行程序即可，即有下面得结果：3、已知在一50);P1ot(x,y)运行结果:Fil«EditlimInsertTool'sDtsktopWindowHelv□ESQfe]\|x<7◎42諂*□E|■El由图可瞧出，函数在零点不连续•4、绘制极坐标曲线p=asin(b+n0)，并分析参数a、b、n对曲线形状得影响解：M文件如下：clc;theta=O:pi/100:2*pi;a=input('输入a=');b=input(/输入b=/);n=input('输入n=');rho=a*sin(b+neta);polar(theta,rho,'m')采用控制变量法得办法，固定两个参数，变动第三个参数观察输出图象得变化。lu«llillaw如xrlXiiU1皿33drvUU2ddjbV■OkX-IQ20□k•、乂£3量/・<2D目■口tf=tRF2:15■••••••<•••••1':"•Z4HYi«beertloheyI«l|«i:广b^100■DmUlZ・!>;•uh丫一r«ubiila的［釦Y-Z0山的々\宀佝2JdXib•、•、少①O0|iO»F.bUi«Viw叶Tm1«]・・XqMy■judik\\CSk^-2OBBE运行结果:分析结果：由这8个图知道，当a,n固定时，图形得形状也就固定了，b只影响图形得旋转得角度；当a,b固定时,n只影响图形得扇形数，特别地，当n就是奇数时，扇叶数就就是n,当就是偶数时，扇叶数则就是2n个；当b,n固定时，a影响得就是图形大小，特别地，当a就是整数时，图形半径大小就就是a。5、绘制函数得曲线图与等高线。其中x得21个值均匀分布[—5,5]范围，y得31个值均匀分布在]0,10],要求使用subplot(2,1,1)与subplot(2,1,2)将产生得曲面图与等高线图画在同一个窗口上。解：M文件：clc；x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x)、*cos(y)、*exp(-sqrt(x、人2+y、A2)/4);subplot(2,1,1);surf(x,y,z);titleC曲面图')；subplot(2,1,2);surfc(x,y,z);title('等高线图/);|、1口1x|I.ViawDw&ktcipndowHwip□臼H◎1h\-X'1|□目■□曲面图等高銭图0-56、绘制曲面图形，并进行插值着色处理。解:M文件：clc;s=0:pi/100:pi/2;t=0:pi/100:3*pi/2;[s，口=meshgrid(s,t);x=cos(s卜*cos(t);y=cos(s)、*sin(t);z=sin(s);subplot(2,2,1);mesh(x,y,z);title('未着色得图形’);subplot(2,2,2);surf(x,y,z);title('shadingfaceted(缺省)');subplot(2,2,3);surf(x,y,z)；shadingflat;title(h'adingflat');subplot(2,2,4);surf(x,y,z);shadinginterp;title(zshadinginterp';)运行结果有:运行结果:实验七低层绘图操作、实验内容1、建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有得菜单项，而且在按下鼠标器得左键之后显示出LeftButtonPressed字样.解:M文件如下：clc;hf=figure('color',[100],、、、'WindowButtonDownFen'/disp(''LtButtonPressed、/')';运行结果：J>1蓉・1应・EZxuiartTodsH-eIaX叶lK.imJh>wUwlp左击鼠标后:2、先利用默认属性绘制曲线y=x2e2x，然后通过图形句柄操作来改变曲线得颜色、线型与线宽，并利用文件对象给曲线添加文字标注^解：M文件：uduj总 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 方差。最大元素与最小元素。(3)大于0、5得随机数个数占总数得百分比。解：M文件：clc;x=fand(1,30000);mu=mean(x)%求这30000个均匀分布随机数得平均值sig=std(x)%求其标准差(Tiy=length(find(x>0、5));%找出大于0、5数得个数5得所占百分比P=y/30000%大于0、运行结果：mu=0、499488553231043sig=0、288599933559786P:0、49942、将100个学生5门功课得成绩存入矩阵P中，进行如下处理：分别求每门课得最高分、最低分及相应学生序号。分别求每门课得平均分与标准方差。5门课总分得最高分、最低分及相应学生序号。将5门课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh.提示:上机调试时，为避免输入学生成绩得麻烦，可用取值范围在]45,95]之间得随机矩阵来表示学生成绩。解:M文件：clc;t=45+50*rand(100,5);P=fix(t);%生成100个学生5门功课成绩[x,l]=max(P)%x为每门课最高分行向量，l为相应学生序号[y,k:=min(P)%y为每门课最低分行向列，k为相应学生序号mu=mean(P)%每门课得平均值行向量sig=std(P)%每门课得标准差行向量s=sum(P,2)%5门课总分得列向量[X,m]=max(s)%5门课总分得最高分X与相应学生序号m[Y,n]=min(s)%5门课总分得最低分Y与相应学生序号n[zcj,xsxh]=sort(s)%zcj为5门课总分从大到小排序，相应学生序号XsXh运行结果:3、某气象观测得某日6:00〜18:00之间每隔2h得室内外温度（0C）如实验表1所示。实验表1室内外温度观测结果（°C）时间h6?8?101??2?14?1618??室内温度t118、020?0?22、025、030?02?8、024、0室外温度t2?15、01?9、024?028?0?34、032、0?30、0试用三次样条插值分别求出该日室内外6:温度(叱)。解：M文件：clc;h=6:2:18;t1=[18、020、022、025、030、028、024、0];t2=[15、019、024、028、034、032、030、0];T1=interp1(h,t1,'spline')%室内得3次样条插值温度T2=interp1(h,t2,'spline')%室外得3次样条插值温度运行结果：T1=Columns1through340、070344、113048、1705Columns4through632654、288564、588360、4512Column752、2444T2=Columns1through334、028442、090252、2444Columns4through660、451272、940868、7503Column764、58834、已知1gx在[1,101］区间10个整数采样点得函数值如实验表2所示。实验表2lgx在10个采样点得函数值X1112131415161718191101lgx01、04141、32221、49141、61281、70761、78531、85131、90851、95102、0043试求1gx得5次拟合多项式p(x),并绘制出lgx与p(x)在［1,101］区间得函数曲线。解：M文件：x=1：10:101;y=lg10(x)；P=po1yfit(x,y,5)y1=poly/val(P,x)；p1ot(x,y,/:o'x,y1,—*/运行结果:Warning:Polynomia1isbadlyconditioned>AddpointswithdistinctXva1ues,reducethedegreeofthepolynomial,ortrycenteringandsca1ingasdescribedinHELPPOLYFIT、>Inpo1yfitat80P=0、0000—0、00000、0001—0、00580、1537—0、1(这里出现警告就是提示不必用5价函数就已经可以完美拟合了，就是可以降价拟合。)在［1,101］得区间函数图像5、有3个多项式P1(x)=x4+2x3+4x2+5,P2(x)=x+2,P3(x)=x2+2x+3，试进行下列操作：(1)求P(X)=P1(x)+P2(x)P3(x)。(2)求P(x)得根。当x取矩阵A得每一元素时，求P(x)得值.其中：当以矩阵A为自变量时，求P(x)得值。其中A得值与第(3)题相同。解:M文件：clc;clear;p1=[1,2,4,0,5];p2=：1,2]；P3=[1,2,3]；p2=[0,0,0,p2];p3=[0,0,p3];p4=conv(p2,p3);%p4就是p2与p3得乘积后得多项式np4=length(p4);np1=length(p1);p=[zeros(1,np4—np1)p1]+p4%求戸(x)=p1(x)+p2(x)x=roots(p)%求p(x)得根A=[—11、2—1、4;0、7523、5;052、5];y=polyval(p,A)%x取矩阵A得每一元素时得p(x)值运行结果:p=3871100001X=—1、3840+1、8317i—1、3840-1、8317i—0、1160+1、4400i—0、1160-1、4400iy=1、0e+003*0、01000、03820、01250、02230、09700、41220、01101、24600、1644实验九数值微积分与方程数值求解、实验目得1、掌握求数值导数与数值积分得方法2、掌握代数方程数值求解得方法。3、掌握常微分方程数值求解得方法。、实验内容1、求函数在指定点得数值导数。解：M文件：clc;clear;x=1；i=1；f=inline('de([xx人whilex<=3、01g(i)=f(x)；i=i+1;x=x+0、01;xA3;12*x3*xA2;026*x])');%以0、01得步长增加，可再缩小步长提高精度g；t=1:0、01:3、01;dx=diff(g)/0、01;f仁dx(1)f2=dx(101)f3=dx(length(g)—1)%差分法近似求导%x=1得数值倒数%x=2得数值倒数%x=3得数值倒数运行结果:fl6、0602f224、1202f354、18022、用数值方法求定积分。得近似值。(2)解:M文件：clc;clear；f=inline(/sqrt(cos(t、人2)+4*sin(2*t)、人2+1)/);I仁quad(f,0,2*pi)g=in1ine('(1+x)、/(1+x、A2)');I2=quad(g,0,2*pi运行结果：=10、4285=0、99973、分别用3种不同得数值方法解线性方程组。解:M文件：clc;c1ear;A=:65—25;9-14-1;342—2;3-902];b=[-413111]'；x=A\by=inv(A)*b[L,U:=lu(A);z=(L\b)运行结果:x=0、6667-1、00001、50C0—0、0000y=c、6667一1、000C1、50C0—-0、0000z二—0、6667-1、0C001、5000—0、00004、求非齐次线性方程组得通解。解：M文件function[x,y]=line—solution(A,b)[m,n]=size(A);y=[];ifnorm(b)>0%非齐次方程组ifrank(A)==rank([A,b])ifrank(A)==ndisp(/有唯一解x')；x=A\b;elsedisp('有无穷个解，特解x,基础解系y');x=A\b;y=null(A,'endelsedisp('无解/);x=[]；endelse%齐次方程组disp(/有零解X，);x=zeros(n,1);ifrank(A)〈ndisp(/有无穷个解，基础解系y');y=null(A,'));endendclc;clear;formatratA=[2731;3522;9417];b=[642]'[x,y]=line—solution(A,b)运行结果：有无穷个解，特解X,基础解系yWarning:Rankdeficient,rank=2,tol=8、6112e—015、>Inline_solutionat11x=—2/1110/II00y=1/11—9/115/111/111001所以原方程组得通解就是:，其中为任意常数。5、求代数方程得数值解。3x+sinx—ex=0在xo=1、5附近得根。在给定得初值x0=1,y0=1,zo=1下，求方程组得数值解。?解:M文件:functiong=f(x)g=3*x+sin(x)—exp(x);clc;c1ear;fzero('f'3)结果就是:ans=1289/682⑵、M文件：functionF=fun(X)x=X(1);y=x(2)；z=X(3);F(1)=sin(x)+yA2+1og(z)—7;F(2)=3*x+2-zA3+1;F(3)=x+y+z-5;X=fsolve('myfun'[,1,1,1]',(otimset('Display'off))运行结果：X=909/10731735/7281106/6256、求函数在指定区间得极值。(1)在(0,1)内得最小值。(2)在]0,0]附近得最小值点与最小值。解：M文件：functionf=g(u)x=u(1);y=u(2);f=2*x、人3+4*x、*y人3—10*x、*y+y、a2;clc;clear;formatlongf=inline('(xA3+cos(x)+x*log(x))/exp(x)f;[x,fmin1]=fminbnd(f,0,1)[U,fmin2]=fminsearch(fgf[0,0])运行结果x=0、522288340666172fmin1=0、397363464998461U=1、0810、833488282765738fmin2=—3、32447、求微分方程得数值解.?解：M文件：functionxdot=sys(x,y)xdot=:y(2);(5*y(2)-y(1))/x:;clc;clear;x0=1、0e—9;xf=20;[x,y]=ode45('sys/,[x0,xf],[00]);[x,y]运行结果：xyyans=0、0000000、5000001、0000001、5000002、0000002、5000003、0000003、5000004、0000004、5000005、0000005、5000006、0000006、5000007、0000007、5000008、0000008、5000009、0000009、50000010、00000010、50000011、00000011、50000012、00000012、50000013、00000013、50000014、00000014、50000015、00000015、50000016、00000016、50000017、00000017、50000018、00000018、50000019、00000019、50000020、000000&求微分方程组得数值解，并绘制解得曲线。解：令y1=x,y2=y,y3=z;这样方程变为：，自变量就是tM文件：functionxdot=sys(x,y)xdot=[y(2)*y(3);—y(l)*y(3);—0、51*y(1)*y(2)];clc;clear;t0=0;tf=8;[x,y]=ode23('sy',[t0,tf],:0,1,1])plot(x,y)运行结果：x=00、00010、00050、00250、01250、06250、16320、30330、48290、71620、98491、26101、56781、95502、32872、70243、01533、29213、48893、64523、75383、86243、99414、16454、38354、65374、92655、22455、5861TOC\o"1-5"\h\z6、03026、34286、65556、93717、15417、32387、45027、57657、70427、87068、0000y=01、00001、00000、00011、00001、00000、00051、00001、00000、00251、00001、00000、01250、999!91、00000、06240、99800、99900、16210、98680、99330、29650、95500、97730、45630、88980、94540、63500、77220、89120、79440、60690、82330、90690、42030、76170、97800、20660、71550、9975-0、06440、70160、9450—0、32580、73770、8127-0、58170、81410、6303—0、77550、89270、4130-0、90980、95520、2324-0、97160、98580、0795-0、99580、9980—0、0289—0、99860、9994-0、1367—0、98960、9949—0、2640-0、96340、9817—0、4187—0、90690、9538—0、5935-0、80340、9053-0、7644-0、64270、8373—0、8859—0、46090、7738-0、9656—0、25420、7235-0、9985—0、00140、7003—0、94950、30920、7343—0、84950、52470、7943-0、68580、72560、8712—0、48020、87510、9387—0、28880、95540、9778—0、12540、99020、99530、00020、99810、99930、12570、99010、99530、24940、96630、98330、40160、91360、95720、51000、85780、9305图形:I.lluiI3£JlmW工±w*L工卫*1匚£1kJkkwv企血gvtIJulp・Uc213J|^.-”，?7壬＞乞山"口OtE)].口实验十符号计算基础与符号微积分、实验目得1、掌握定义符号对象得方法•2、掌握符号表达式得运算法则以及符号矩阵运算。3、掌握求