实验报告(6) - 图文

怀化学院数学系 《数学试验》实验报告

实验序号： 实验项目名称：

学 号 实验地点 一、实验目的及要求 二、实验设备（环境）及要求 三、实验内容与步骤 (1.1) >> syms x1 x2; >> [x1,x2]=solve('x1^2-x2-1=0','(x1-2)^2+(x2-0.5)^2-1=0') x1 = -1.3068444845633173592407431426632-1.2136904451605911320167045558746*i -1.3068444845633173592407431426632+1.2136904451605911320167045558746*i 1.0673460858066897134085973128070 1.5463428833199450050728889725194 x2 = -.76520198984055124633885565606586+3.1722093284506318231178646481143*i -.76520198984055124633885565606586-3.1722093284506318231178646481143*i .13922766688686144048362498805141 1.3911763127942410521940863240803 >> err=[x1.^2-x2-1,(x1-2).^2+(x2-0.5).^2]; >> norm(double(eval(err))) ans = 2.0000 (1.2) >> syms x y z; 姓 名 指导教师 >> [x,y,z]=solve('x^2*y^2-z*x*y-4*x^2*y*z^2=x*z^2','x*y^3-2*y*z^2=3*x^3*z^2+4*x*z*y^2x = 0 0 x .31763330066708094899226744383312 -.57342208805666904699955413755011 -3.0414669315515393342891214854118 1.2406494592172202024946761753539+2.6726341184257446950051342796302*i 1.9395779085729820052691770221731-1.8460302354285923475099876321561*i -.12175724403625789114277654158124+.76608643264311922968528253940275*i -2.2119059377823750236571580666050-1.3098483298501442852199024916754*i .18091793647835958266865409918254-.43514167087309843831595260429073*i -.55869174189107447600021216575924+2.3297622385407573613585339888006*i -.68680307411564099103394996935634 2.2449864792628313778324119367757 .80149155267622824623322517104799 -.55869174189107447600021216575924-2.3297622385407573613585339888006*i .18091793647835958266865409918254+.43514167087309843831595260429073*i -2.2119059377823750236571580666050+1.3098483298501442852199024916754*i -.12175724403625789114277654158124-.76608643264311922968528253940275*i 1.9395779085729820052691770221731+1.8460302354285923475099876321561*i 1.2406494592172202024946761753539-2.6726341184257446950051342796302*i y = 0 y 0 -.10247418659393300993582067594236 .88186908582844666977541315922524 6.2206205278493911181266569155794 3.8883892104923440119866877297623+2.6775872796906346863164643744935*i -.19521984362257729205542902543114+2.1755659889967168917621095367677*i .19300043308429639136627772430946+.55085821682149293523314538512814*i -1.5345777137990494704652954417314+.42952982009383733582600047743750e-1*i .59285406023469299432404383306933e-1-.14560032849820344999561561093139*i -.13504742647411551604956518296015-1.3353390289175827351921401819303*i .23206681025858886578639216095528 1.2500352093099975887186243233379 -.40877757806122608090142625674869 -.13504742647411551604956518296015+1.3353390289175827351921401819303*i .59285406023469299432404383306933e-1+.14560032849820344999561561093139*i -1.5345777137990494704652954417314-.42952982009383733582600047743750e-1*i .19300043308429639136627772430946-.55085821682149293523314538512814*i -.19521984362257729205542902543114-2.1755659889967168917621095367677*i 3.8883892104923440119866877297623-2.6775872796906346863164643744935*i z = z 0 0 .14437409800864935567667171646044 1.2197465522537746902599296344473 1.2977198721910510602444950254432 1.0166830329716235095778013718903+.37011032131730244971789188681237*i .45233167057904744299654223927744+.51464117109451634939209195854989*i .31864673037377358785409663022781+.78415604801861425880017838557048*i .38457694363000428639271157033224e-1+.57967291207543468143911103709180*i -.83830187383058822380943737709146e-1+.20866770020547056525334588752981*i -.36349211548848392020797588683844+.46022200523351838974797619203293*i -.34096135567852262861343237947182 -.58923295037333858669309875602950 -.61423986723341834383214878861209 -.36349211548848392020797588683844-.46022200523351838974797619203293*i -.83830187383058822380943737709146e-1-.20866770020547056525334588752981*i .38457694363000428639271157033224e-1-.57967291207543468143911103709180*i .31864673037377358785409663022781-.78415604801861425880017838557048*i .45233167057904744299654223927744-.51464117109451634939209195854989*i 1.0166830329716235095778013718903-.37011032131730244971789188681237*i >> norm(double([x.^2.*y.^2-z.*x.*y-4*x.^2.*y.*z.^2-x.*z.^2, ans = 3.5192e-024 (2.1) x.*y.^3-2*y.*z.^2-3*x.^3.*z.^2-4*x.*z.*y.^2, y.^2.*x-7*x.*y.^2+3*x.*z.^2-x.^4.*z.*y]>> ezplot('exp(-(x+1)^2+pi/2)*sin(5*x+2)',[0 5]) >> hold on;line([0,5],[0,0]) (2.2) ezsurf('(x^2+y^2+x*y)*exp(-x^2-y^2-x*y)') (3.1) >> syms x; >> x4=fsolve(f,4) x4 = >> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x4) ans = -5.913350831018913e-013 >>x1=solve('exp(-(x+1)^2+pi/2)*sin(5*x+2)') x1 = -2/5 >> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x1) ans = 0 >> f=inline('exp(-(x+1).^2+pi/2).*sin(5*x+2)','x'); x2=fsolve(f,0) x2 = 0.22831852178755 >> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x2) ans = 4.750949292642762e-008 >> x3=fsolve(f,-1) x3 = -1.02831853071796 >> subs('exp(-(x+1)^2+pi/2)*sin(5*x+2)',x,x3) ans = -5.886413288211306e-016 (3.2) >> syms x; >> y1=solve('(x^2+y^2+x*y)*exp(-x^2-y^2-x*y)=0','y') y1 = (-1/2+1/2*i*3^(1/2))*x >> y2=simple(subs('(x^2+y^2+x*y)*exp(-x^2-y^2-x*y)','y',y1)) y2 = 0 (4) >> syms x c; >>y=int((exp(x)-c*x)^2,x,0,1) y = 1/2*exp(1)^2+1/3*c^2-1/2-2*c >>function y=exc6ff(c) y=1/2*exp(1)^2+1/3*c^2-1/2-2*c >> x=fminsearch('exc6ff',0) x = 3.00000000000000 >> ezplot(y,[0,5] (5) >> x=fminunc(f,ones(7,1)) x = 10.5464 111.2316 6.7823 -111.5042 1.0000 1.0000 1.0000 (6) >> y=inline('20+x(1)^2+x(2)^2-10*(cos(pi*x(1))+cos(pi*x(2)))','x'); >> x1=fminunc(y,[1;0]) x1 = 1.0e-008 * -0.3525 0.4427 >> y(x1) ans = 3.5527e-015 >> x2=fminunc(y,[0;2]) x2 = -0.0000 1.9602 >> y(x2) ans = 3.9205 >> [x1,x2]=meshgrid(-3:0.1:3); >> z=20+x1.^2+x2.^2-10*(cos(pi*x1)+cos(pi*x2)); >>f=inline(['100*(x(2)-x(1)^2)^2+(1-x(1))^2+','90*(x(4)-x(3)^2)+(1-x(3)^2)^2+','10.1*

本文来源：https://www.bwwdw.com/article/0e0d.html