资源描述：
《自动化数学实验实验二》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、重庆大学学生实验报告实验课程名称数学实验开课实验室DS1404学院2011年级自动化专业班02班学生姓名朱梦磊学号20114983开课时间2012至2013学年第二学期总成绩教师签名数学与统计学院制开课学院、实验室:实验时间:2013年3月19日课程名称数学实验实验项目名称第二次实验实验项目类型验证演示综合设计其他指导教师温罗生成绩实验目的1)掌握方程求解的三种解法：解析法、数值解法以及图形表示解的方法；2)学会使用MATLAB软件求解析解、数值解和图形解；3)通过范例学习怎样建立方程模型和分析问题的思
2、想。基础实验一、实验内容—、(1),先建立一个函数M文件输入functionff=myf(x)ff=x*sin(x)+cos(x)-sin(x)-2*x;保存为myf.m。再建立M文件使用fzero函数，输入fzero(@myf,[・l,l])保存为ex_l_l,结果ans=0.3563；(2)先建立一个函数M文件输入functionff=myf2(x)ff=xA7+2*xA5+l;保存为myf2.m建立一个函数M文件输入functionff=mynewton(x)ff=(xA7+2*xA5+1)/(7
3、*xA64-10*xA4);保存为mynewton.m再建立M文件,输入p=[1,0,2,0,0,0,0,1];xl=roots(p)x2=fzero(@myf2,[-l,0])x3=fsolve(@myf2,-1)x0=-l;x4=x0-mynewton(x0);whileabs(x4-x0)>0.0001x0=x4;x4=x0-mynewton(x0);endx4保存为ex_l_2o结果xl=0.0577+1.4273i0.0577-1.4273i0.7063+0.4562i0.7063-0.456
4、2i-0.3533+0.8480i-0.3533・0.8480i-0.8214实根为-0.8214x2二-0.8214Optimizationterminated:first-orderoptimalityislessthanoptions.TolFun.x3二-0.8214x4二-0.8214二、建立M文件，输入A=[721-2;9153-2;-2-2115;13213];b=[4;7;・l;0];x=Ab保存为ex_2结果x二0.49790.14450.0629-0.0813三、(1)在m_fil
5、e中输入：functionff=myf3(x)ff(1)=2*x(1)-x(2)-exp(-x(1));ff(2)=-x(l)+2*x(2)・exp(・x(2));保存为myfBo在m_file中输入：x=fsolve('myf3',[l,1])保存为ex_3_lo结果x=0.56710.5671(2)在m_file中输入:functionff=myf4(x)ff(l)=x(l)A2-5*x(2)A2+7*x(3)A2+12;ff(2)=3*x(1)*x(2)+x(1)*x(3)-11*x(1);ff(
6、3)=2*x(2)*x(3)+40*x(l);保存为myf4・m在m_file中输入：x1=fsolve('myf4[1,1,1])x2=fsolve('myf4[-111])保存为ex_3_2o结果xl二0.00001.54920.0000x2=0.0000-1.5492-0.0000四、炮弹发射视为斜抛运动，已知初始速度为200m/s,问要击中水平距离360m、垂直距离160m的目标，当忽略空气阻力时，发射角应多大？进一步思考：如果要考虑水平方向的阻力，且设阻力与（水平方向）速度成正比，系数为0
7、」(1/skg),结果又如何？(本题要求：需要给出问题的分析和求解过程，最后给出问题的解。)解：设炮弹打中目标经过的时间为t,发射角为a,水平方向速度为vl,垂直方向速度为v2,则:不考虑水平方向阻力：v1=200cos(a),v2=200sin(a);考虑到竖直方向上的重力影响，有这2个方程：200tcos(a)=360(1)200tsin(a)-0.5gtA2=160(2),在此g取为9.8。两式消去(可得，360tan(a)-15.876/(cos(a)A2)=160在M.file中输入：a=so
8、lve(360*tan(x)-l5.876/(cos(x)A2)=160*)结果：-2.6783054372437130324777545077894-1.61585921356174772402503873810291.5257334400280455144376046451767.46328721634608020598488887549010考虑到实际情况：发射角为1.52或者0.46o考虑水平方向阻力：dv/dt=200*exp(-
