《2022年MATLAB程序设计方案与应用实验参考答案.docx》由会员分享，可在线阅读，更多相关《2022年MATLAB程序设计方案与应用实验参考答案.docx（34页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精品学习资源封面欢迎下载精品学习资源作者： PanHongliang仅供个人学习MATLAB 程序设计与应用（其次版）试验参考答案% 试验一 MATLAB 运算基础%第一题%（ 1 ） z1=2*sin85*pi/180/1+exp2欢迎下载精品学习资源%（ 2 ）x=2,1+2i；-0.45,5；z2=0.5*logx+sqrt1+x.2%3a=-3.0:0.1:3.0；z3=exp0.3*a-exp-0.3*a/2.*sina+0.3+log0.3+a/2%4t=0:0.5:2.5；z4=t.2.*t=0&t=1&t=2&t=A&A=Z；Ak=% 试验二 MATLAB矩阵分析与处理%第一题
2、E=eye3；R=rand3,2；欢迎下载精品学习资源O=zeros2,3；S=diag2,3；A=E,R；O,S；A2B=E,R+R*S；O,S2%其次题H=hilb5 P=pascal5Hh=detH Hp=detP Th=condH Tp=condP%第三题 : A=fix10*rand5 H=detATrace=traceA Rank=rankA Norm=normA%第四题 :A=-29,6,18；20,5,12；-8,8,5 V,D=eigA%数学意义略%第五题方法一 :%1:A=1/2,1/3,1/4；1/3,1/4,1/5；1/4,1/5,1/6；b=0.95,0.67,0.5
3、2；x=invA*b%2:B=0.95,0.67,0.53；x=invA*B%3:condA%第五题方法二 : A=hilb4 A:,1=A4,:=B=0.95,0.67,0.52；X=invA*BB1=0.95,0.67,0.53；X1=invA*B1 N=condB N1=condB1Na=condA %矩阵 A为病态矩阵%第六题欢迎下载精品学习资源A=1,4,9；16,25,36；49,64,81B=sqrtmAC=sqrtA %sqrtm函数是以矩阵为单位进行运算，sqrt函数是以矩阵中的元素进行运算%试验三选择程序结构设计%第一题程序一x=-5.0,-3.0,1.0,2.0,2.5,
4、3.0,5.0；y=；%建立存放全部y 值的矩阵for x0=xif x0=0&x05&x0=2&x0=3y=y,x0*x0-5*x0+6；elsey=y,x0*x0-x0-1；x y%第一题程序二%输出全部 x输出全部 yx=-5,-3,1,2,2.5,3,5；y=；for a=1:7if xa=0&xa100|x0disp您输入的成果不是百分制成果，请重新输入；elseif x=90 dispA；elseif x=80dispB；elseif x=70dispC；elseif x60欢迎下载精品学习资源dispD；elsedispE；end end%其次题程序二：s=input请输入一个成
5、果（0 分到 100分之间）： ； %s用于存放成果while 1%判定输入成果的合理性if s100disp输入的成果需在 0 到 100 之间，请重新输入： s=input请输入一个成果（0 分到 100 分之间）： ；elsebreak；end endswitch fixs/10%对成果做出等级判定case 9,10 dispAcase 8dispB case 7dispC case 6dispD otherwisedispEend%第三题n=input请输入员工工号：；h=input该员工工作时数是： ；if h120x=h-120*84*1+0.15+120*84；elseif h5
6、disp超出了矩阵的行数，矩阵的最终一行为： a5,:elsedisp矩阵的第 ,num2strn,行为： an,:end% 试验四 循环结构程序设计%第一题程序一s=0 ；n=inputn=.；for i=1:ns=s+1/i/i；end PI=sqrt6*s pi%第一题程序二n=inputn=.；a=1:n；b=1./a.2；PI=sqrt6*sumb pi%其次题y=0 ；n=1 ；whiley1e-5 Xn=Xn1；Xn1=a/b+Xn；n=n+1；if n=500break；end endn Xn1r1=-b+sqrtb*b+4*a/2 r2=-b-sqrtb*b+4*a/2%第四
7、题for i=1:100if i=1fi=1；elseif i=2fi=0；elseif i=3 fi=1；elsefi=fi-1-2*fi-2+fi-3；end end maxfminfsumf lengthfindf0 lengthfindf=0 lengthfindf0%第五题： s=0 ； n=0 ；for i=2:49b=i*i+1-1；m=fixsqrtb；for j=2:mif remb,j=0 break欢迎下载精品学习资源end endif j=mn=n+1 ；s=s+b ；endend ns% 试验五 函数文件%第一题function y=mat1x%建立函数文件mat1.
8、my=expx,logx,sinx,cosx；%在命令窗口调用上述函数文件：y=mat11+i%其次题程序一function a,b,N,M=shiyanwu2m,n,tA=m*cost*pi/180,-m,-sint*pi/180,0；m*sint*pi/180,0,cost*pi/180,0； 0,n,-sint*pi/180,0；0,0,-cost*pi/180,1；B=0,9.8*m,0,9.8*n；C=invA*B；a=C1； b=C2； N=C3； M=C4 ；%在命令窗口调用该函数文件:m1=inputm1=；m2=inputm2=；theta=inputtheta=；a1,a2
9、,N1,N2=shiyanwu2m1,m2,theta%其次题程序二function X=mat2m1,m2,t g=9.8；A=m1*cost*pi/180,-m1,-sint*pi/180,0；m1*sint*pi/180,0,cost*pi/180,0； 0,m2,-sint*pi/180,0； 0,0,- cost*pi/180,1；B=0 ；m1*g ；0 ；m2*g；X=invA*B；%在命令窗口调用该函数文件： X=mat21,1,60%第三题function flag=mat3x欢迎下载精品学习资源flag=1；for i=2:sqrtx if remx,i=0flag=0；b
10、reak；end end%在命令窗口调用该函数文件： for i=10:99j=10*remi,10+fixi/10；if mat3i&mat3j dispiend end%第四题function y=fxxy=1./x-2.2+0.1+1./x-3.4+0.01；%在命令窗口调用该函数文件： y=fx2a=1,2； 3,4；y=fxa%第五题%1function f1=mat5n f1=n+10*logn*n+5；%在命令窗口中调用该函数文件：y=mat540/mat530+mat520%2 方法一function f2=mat6n f2=0；for i=1:nf2=f2+i*i+1；end
11、%在命令窗口中调用该函数文件如： y=mat640/mat630+mat620%2 方法二function f2=mat7n i=1:n；m=i.*i+1；f2=summ；end%在命令窗口中调用该函数文件如： y=mat740/mat730+mat720% 试验六 高层绘图操作欢迎下载精品学习资源%第一题：x=linspace0,2*pi,101；y=0.5+3*sinx./1+x.2.*cosx；plotx,y%其次题：%（ 1 ）x=linspace-2*pi,2*pi,100；y1=x.2；y2=cos2*x；y3=y1.*y2；plotx,y1,b-,x,y2,r:,x,y3,y-t
12、ext4,16,leftarrow y1=x2；text6*pi/4,-1,downarrow y2=cos2*x；text-1.5*pi,-2.25*pi*pi,uparrow y3=y1*y2%（ 2 ）；x=linspace-2*pi,2*pi,100；y1=x.2；y2=cos2*xy3=y1.*y2；subplot1,3,1；%分区plotx,y1 titley1=x2；%设置标题subplot1,3,2；plotx,y2；titley2=cos2*x；subplot1,3,3；plotx,y3；titley3=x2*cos2*x；%（ 3 ）x=linspace-2*pi,2*pi
13、,20；y1=x.2；subplot2,2,1；%分区barx,y1；titley1=x2的条形图 ；%设置标题subplot2,2,2；stairsx,y1；titley1=x2的阶梯图 ；subplot2,2,3；stemx,y1；titley1=x2的杆图 ；subplot2,2,4；fillx,y1,rtitley1=x2；%假如少了 r的填充图 ；就会出错%其他的函数照样做；欢迎下载精品学习资源%第三题x=-5:0.01:5；y=；%起始设 y 为空向量for x0=xif x0=0 %不能写成 x0=0.5 n=n+1 ；endend p=n/30000%其次题%1欢迎下载精品学习
14、资源A=45+51*rand100,5Y,U=maxA；a,b=minA%2m=meanAs=stdA%3sumA,2Y,U=maxansa,b=minans%4zcj,xsxh=sortans%第三题h=6:2:18；x=6.5:2:17.5；t1=18,20,22,25,30,28,24；t2=15,19,24,28,34,32,30T1=splineh,t1,x；T2=splineh,t2,x%第四题x=1:0.1:101；y1=log10x；p=polyfitx,y1,5y2=polyvalp,x；plotx,y1,:,x,y2,-%第五题%1p1=1,2,4,0,5；p2=1,2p3
15、=1,2,3；p=p1+0,convp2,p3%2%为使两向量大小相同，所以补0A=rootsp%3A=-1,1.2,-1.4polyvalp,A；0.75,2,3.5；0,5,2.5；%4polyvalmp,A试验十程序：x=sym6；y=sym5；z=x+1/sqrt3+x-sqrty欢迎下载精品学习资源1 、 分解因式（1 ） 程序：syms x y；A=x4-y4；factorA（2 ） 程序：factorsym5135 3 、化简表达式（1 ）程序：syms beta1 beta2y=sinbeta1*cosbeta2-cosbeta1*sinbeta2 simpley（2 ）程序：
16、syms x y=4*x2+8*x+3/2*x+1 simpley5 、用符号方法求以下极限或导数（1 ） 程序： syms xf=x*expsinx+1-2*exptanx-1/sinx limitf（2 ）程序：syms xy=sqrtpi-sqrtacosx/sqrtx+1；limitf,x,-1,right（3 ） 程序： syms xy=1-cos2*x/x；y1=diffy y2=diffy,x,26 、用符号方法求以下积分（1 ） 程序： syms xf=1/1+x4+x8 intf（2 ）欢迎下载精品学习资源程序：syms x f=1/asinx2*sqrt1-x2 intf（
17、3 ）程序：syms x f=x2+1/x4+1 intf,x,0,inf（4 ）程序：syms x f=expx*1+expx2 y=intf,x,0,log2 doubley试验十一 级数与方程符号求解1. 级数符号求和；(1) 运算 ；(2) 求级数 的和函数，并求 之和；解：M 文件：clear all ；clc；n=symn ；x=symx ；S1=symsum1/2*n-1,n,1,10 S2=symsumn2*xn-1,n,1,infS3=symsumn2/5n,n,1,inf%vpaS3可以转化成小数运行结果：S1 = 31037876/14549535 S2 =piecewi
18、seabsx 1, -x2 + x/x*x - 13 S3 =15/322. 将 lnx 在 x=1 处按 5 次多项式开放为泰勒级数；解：M 文件： clear all ；clc；x=symx ；f=logx ；taylorf,x,6,1 运行结果： ans =x - x - 12/2 + x - 13/3 - x - 14/4 + x - 15/5 - 13. 求以下方程的符号解；欢迎下载精品学习资源解：M 文件：clear all ；clc；x1=solvelogx+1-5/1+sinx=2 x2=solvex2+9*sqrtx+1-1 x3=solve3*x*expx+5*sinx-7
19、8.5x4 y4=solvesqrtx2+y2-100,3*x+5*y-8运行结果：x1 = 521.67926389905839979437366649258x2 =-131/2*i*4/9*64651/2/2+2171/541/3-1/2*64651/2+2171/541/3/2- 64651/2/2 + 2171/541/3/2 - 2/9*64651/2/2 + 2171/541/3 + 1/31/3 - 64651/2/2+ 2171/541/3/2- 31/2*i*4/9*64651/2/2+ 2171/541/3- 1/2*64651/2 + 2171/541/3/2 - 2/9
20、*64651/2/2 + 2171/541/3x3 =2.3599419584772910151699327715486x4 =12/17 - 10*212461/2/17 10*212461/2/17 + 12/17y4 =6*212461/2/17 + 20/17 20/17 - 6*212461/2/174. 求微分方程初值问题的符号解，并与数值解进行比较；解：M 文件：clear all ；clc；dsolveD2y+4*Dy+29*y,y0=0,Dy0=15,x运行结果：ans = 3*sin5*x/exp2*x5. 求微分方程组的通解；解：M 文件：clear all ；clc；x
21、 y z=dsolveDx=2*x-3*y+3*z,. Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z,t运行结果：x =C1/expt + C2*exp2*t y =欢迎下载精品学习资源C1/expt + C2*exp2*t + C3/exp2*t z =C2*exp2*t + C3/exp2*t试验九 数值微积分与方程数值求解1. 求函数在指定点的数值导数；解： M 文件：clc；clear；x=1 ；i=1 ；f=inlinedetx x2 x3；1 2*x 3*x2 ；0 2 6*x ；while x0 % 非齐次方程组if rankA=rankA,b if rankA=n
22、disp 有唯独解 x；x=Ab ；elsedisp 有无穷个解，特解x，基础解系 y；x=Ab ；y=nullA,r ；end elsedisp无解 ；x= ；endelse %齐次方程组disp有零解 x；x=zerosn,1；if rankA In line_solution at 11 x =-2/11 10/11欢迎下载精品学习资源00y =1/11-9/11-5/111/111001所以原方程组的通解是：，其中 为任意常数；5. 求代数方程的数值解；(1) 3x+sinx-ex=0 在 x0=1.5 邻近的根；(2) 在给定的初值x0=1 ， y0=1， z0=1 下，求方程组的数
23、值解；解： M 文件： function g=fx g=3*x+sinx-expx；clc；clear；fzerof,1.5结果是：ans =1289/6822. M 文件： function F=funX x=X1 ；y=X2 ；z=X3 ；F1=sinx+y2+logz-7；F2=3*x+2-z3+1 ；F3=x+y+z-5 ；X=fsolvemyfun,1,1,1,optimsetDisplay,off 运行结果：6. 求函数在指定区间的极值；(1) 在0,1 内的最小值；(2) 在0,0 邻近的最小值点和最小值；解： M 文件：function f=gu x=u1 ； y=u2 ；f=
24、2*x.3+4*x.*y3-10*x.*y+y.2；clc；clear；format longf=inlinex3+cosx+x*logx/expx；x,fmin1=fminbndf,0,1 U,fmin2=fminsearchg,0,0运行结果欢迎下载精品学习资源7. 求微分方程的数值解；解： M 文件：function xdot= sys x,y xdot=y2 ；5*y2-y1/x；clc；clear；x0=1.0e-9 ；xf=20 ；x,y=ode45sys,x0,xf,0 0；x,y运行结果：8. 求微分方程组的数值解，并绘制解的曲线；解： 令 y1=x,y2=y,y3=z ； 这
25、样方程变为 :,自变量是 t M 文件：function xdot=sysx,yxdot=y2*y3；-y1*y3 ；-0.51*y1*y2；clc；clear；t0=0 ； tf=8 ；x,y=ode23sys,t0,tf,0,1,1 plotx,y版权申明本文部分内容，包括文字、图片、以及设计等在网上搜集整理；版权为潘宏亮个人全部This article includes some parts, including text,pictures, and design. Copyright is Pan Hongliangs personalownership.用户可将本文的内容或服务用于个
26、人学习、争论或观看，以及其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关法律的规定，不得侵害本网站及相关权益人的合法权益；除此以外，将本文任何内容或服务用于其他用途时，须征得本人及相关权利人的书面许可，并支付酬劳；Users may use the contents or services of this欢迎下载精品学习资源article for personal study, research or appreciation, and other non-commercial or non-profit purposes, but at the same time, they sh
27、all abide by the provisions of copyright law and other relevant laws, and shall not infringe upon the legitimate rights of this website and its relevant obligees. In addition, when any content or service of this article is used for other purposes, written permission and remuneration shall be obtained from the person concerned and the relevant