当前位置:文档之家› 数学建模matlab例题参考及练习

数学建模matlab例题参考及练习

数学实验与数学建模实验报告学院:专业班级:姓名:学号:完成时间:年月日承 诺 书本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成,所有结果都通过上机验证,无转载或抄袭他人,也未经他人转载或抄袭。

若承诺不实,本人愿意承担一切责任。

承诺人:年 月 日数学实验学习体会(每个人必须要写字数1200字以上,占总成绩的20%)练习1 一元函数的图形1. 画出x y arcsin =的图象.2. 画出x y sec =在],0[π之间的图象.3. 在同一坐标系中画出x y =,2x y =,3x y =,3x y =,x y =的图象. 4. 画出3232)1()1()(x x x f ++-=的图象,并根据图象特点指出函数)(x f 的奇偶性. 5. 画出)2ln(1++=x y 及其反函数的图象. 6. 画出321+=x y 及其反函数的图象.练习2 函数极限1.计算下列函数的极限.(1)xxx4cos12sin1lim4-+π→.程序:sym x;f=(1+sin(2*x))/(1-cos(4*x));limit(f,x,pi/4)运行结果:lx21ans =1(2).程序:sym x;f=(1+cos(x))^(3*sec(x));limit(f,x,pi/2)运行结果:lx22ans =exp(3)(3)22)2(sinlnlimxxx-ππ→.程序:sym x;f=log(sin(x))/(pi-2*x)^2;limit(f,x,pi/2)运行结果:lx23ans =-1/8(4)212lim xxex→.程序:x xx sec32)cos1(lim+π→sym x ;f=x^2*exp(1/x);limit(f,x,0)limit(f,x,0,'right')limit(f,x,0,'left')运行结果:lx24ans =NaNans =Infans =%左极限为零,存在,右极限为无穷大,在x 趋近于零时函数没有极限(5))215(lim 122x x x x +-∞→.程序:sym x ;f=5*x^2/(1-x^2)+2^(1/x);limit(f,x,inf)运行结果:>> lx25ans =-4(6)x x x x x -+-→32112lim .程序:sym x ;f=(x^2-2*x+1)/(x^3-x);limit(f,x,1)运行结果:>> lx26ans =0 (7)x x x 11lim 20-+→.程序:sym x ;f=(sqrt(1+x^2)-1)/x;limit(f,x,0)运行结果:>> lx27ans =0 (8))3sin(cos 21lim 3π--π→x xx .程序:sym x ;f=(1-2*cos(x))/sin(x-pi/3);limit(f,x,pi/3)运行结果:>> lx28ans =3^(1/2) (9)tgxx x )1(lim 0+→.程序:sym x ;f=(1/x)^tan(x);limit(f,x,0,'right')运行结果:>> lx29ans =(10)xx arctgx )2(lim π+∞→.程序:sym x ;f=(2/pi*atan(x))^x;limit(f,x,inf,'left')运行结果:>> lx210ans =Inf2. 解方程012=-⋅x x .程序:sym x ;X=solve(x*2^x-1)运行结果:>> lx202X =lambertw(0, log(2))/log(2)%方程有两个解3. 解方程1sin 3+=x x .程序:sym x ;X=solve(3*sin(x)+1-x)运行结果:>> lx203X =-0.53847936154. 解方程03=++q px x .(p 、q 为实数)程序:X=solve('x^3+p*x+q=0','x')运行结果:X =((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3) - p/(3*((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3))p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3)/2 - (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3)/2 + (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2练习 3 导数及偏导数计算1.求下列函数的导数. (1))11)(1(-+=x x y程序:sym x ;f=(sqrt(x)+1)*(1/sqrt(x)-1);diff(f)运行结果:>> lx31ans =(1/x^(1/2) - 1)/(2*x^(1/2)) - (x^(1/2) + 1)/(2*x^(3/2))(2)x x x y ln sin =程序:sym x ;f=x*sin(x)*log(x);diff(f)运行结果:>> lx32ans =sin(x) + log(x)*sin(x) + x*cos(x)*log(x)2.求下列参数方程所确定的函数的导数.(1)⎩⎨⎧==t y t x 44程序:sym t ;f1=t^4;f2=4*t;diff(f2)/diff(f1)运行结果:>> lx321ans =1/t^3(2)⎩⎨⎧-=+=arctgt t y t x )1ln(2程序:sym t ;f1=log(1+t^2);f2=t-atan(t);diff(f2)/diff(f1)运行结果:>> lx322ans =-((t^2 + 1)*(1/(t^2 + 1) - 1))/(2*t)3.求下列隐函数的导数. (1)22ln y x x y arctg +=程序:syms x y ;f=atan(y/x)-log(sqrt(x^2+y^2));yx=-diff(f,x)/diff(f,y)运行结果;>> lx331yx =(x/(x^2 + y^2) + y/(x^2*(y^2/x^2 + 1)))/(1/(x*(y^2/x^2 + 1)) - y/(x^2 + y^2))(2)xy y x =程序:syms x y ;f=x^y-y^xyx=-diff(f,x)/diff(f,y)运行结果:>> lx332f =x^y - y^xyx =(x^(y - 1)*y - y^x*log(y))/(x*y^(x - 1) - x^y*log(x))4.设x e y x cos =,求)4(y .程序:sym x ;f=exp(x)*sin(x);diff(f,x,4)运行结果:>> lx34ans =(-4)*exp(x)*sin(x)5.验证x e y x sin =满足关系式:22=+'-''y y y程序: sym x ;f=exp(x)*sin(x);y2=diff(f,x,2);y1=diff(f,x,1);y=f;y2-y1*2+2*y=='0'运行结果:>> lx35ans =1%运行结果为1表示y2-y1*2+2*y=='0'成立6.设)ln(y x x u +=,求22x u ∂∂,22y u∂∂,y x u ∂∂∂2. 程序:syms x y ;f=x*log(x+y);uxx=diff(f,x,2)uyy=diff(f,y,2)f1=diff(f,x);uxy=diff(f1,y)运行结果:>> lx36uxx =2/(x + y) - x/(x + y)^2uyy =-x/(x + y)^2uxy =1/(x + y) - x/(x + y)^27.求下列多元隐函数的偏导数y zx z ∂∂∂∂,.(1)1cos cos cos 222=++z y x程序:syms x y z ;f=(cos(x))^2+(cos(y))^2+(cos(z))^2-1;zx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果:>> lx371zx =-(cos(x)*sin(x))/(cos(z)*sin(z))zy =-(cos(y)*sin(y))/(cos(z)*sin(z))(2)xyz e z =程序:syms x y z ;f=exp(z)-x*y*zzx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果:>> lx372f =exp(z) - x*y*zzx =(y*z)/(exp(z) - x*y)zy =(x*z)/(exp(z) - x*y)8.证明函数22)()(ln b y a x u -+-=(b a ,为常数)满足拉普拉斯方程: 02222=∂∂+∂∂y u x u(提示:对结果用simplify 化简)练习4 积分计算1.计算下列不定积分. (1)⎰+dx x x 12 (2)⎰+x xdx 2sin 12sin2.计算下列定积分.(1)⎰e xdx x 1ln (2)⎰ππ342sin dx x x 3.求⎰+t dx x x x 12)ln (ln 1并用diff 对结果求导.4.求摆线)cos 1(),sin (t a y t t a x -=-=的一拱(π≤≤20t )与x 轴所围成的图形的面积.5.计算二重积分(1)⎰⎰≤++122)(y x dxdy y x (2)⎰⎰≤++x y x dxdy y x 22)(22 6.计算⎰+L ds y x 22 L 为圆周)0(22>=+a ax y x7.计算⎰++-L dy y x dx y x )()(2222,其中L 为抛物线2x y =上从点(0,0)到点(2,4)的一段弧.练习5 matlab 自定义函数与导数应用1.建立函数x x a a x f 3sin 31sin ),(+=,当a 为何值时,该函数在3π=x 处取得极值,它是极大值还是极小值,并求此极值.2.确定下列函数的单调区间.(1)7186223---=x x x y (2))0(82>+=x x x y3.求下列函数的最大值、最小值.(1)2332xx y -=41≤≤-x (2)312824≤≤-+-=x x x y练习6 matab 矩阵运算与数组运算1. 计算(1)⎪⎪⎪⎭⎫ ⎝⎛--521111204321+⎪⎪⎪⎭⎫ ⎝⎛--232002101041221 (2)⎪⎪⎭⎫ ⎝⎛-01301213⎪⎪⎪⎪⎪⎭⎫ ⎝⎛030101020501⎪⎪⎪⎭⎫ ⎝⎛-205101 (3)52422⎪⎪⎭⎫ ⎝⎛- 2.设⎪⎪⎪⎭⎫ ⎝⎛-=243121013A ,⎪⎪⎪⎭⎫ ⎝⎛-=112111201B ,求满足关系B X A =-23的X .练习7 矩阵与线性方程组1.求下列矩阵的秩.(1)⎪⎪⎪⎭⎫ ⎝⎛-321110021 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛4820322513454947513253947543173125 2.求下列矩阵的行列式,如可逆,试用不同的方法求其逆矩阵.(1)⎪⎪⎪⎭⎫ ⎝⎛--285421122 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---6201111121324321 3.设X ⎪⎪⎪⎭⎫ ⎝⎛-111012111=⎪⎪⎪⎭⎫ ⎝⎛521234311求X . 4.解下列线性方程组.(1)⎪⎪⎩⎪⎪⎨⎧=--+=+-+=+-+=+-+6223312433862344224221432143214321x x x x x x x x x x x x x x x x (2)⎪⎪⎩⎪⎪⎨⎧-=+--=+--=-+-212201432143214321x x x x x x x x x x x x练习8 常微分方程与级数 求1-6题微分方程的通解1.1222+='y y y x 2.x y x y dx dy -+= 3.x x x y y +='cos 4.1)2sin cos (='+y y y x 5.x e y y y x 2cos 3=-'+'' 6.x x y y sin 14++=+'' 求7、8题初值问题的解7.⎪⎩⎪⎨⎧==-++-+=10)2(212222x y dx dy x xy y y xy x8.⎪⎪⎩⎪⎪⎨⎧===++==0000222,02V dt dxx x x a dt dx n dt x d t t9.给出函数x x e x f x x cos 2sin )(+=在点0=x 的7阶taylor 展开式以及在x=1处的 5阶taylor 展开式.10.判别下列级数的敛散性,若收敛求其和.(1)++++7151311(2)∑∞=+112n nntgπ11.求幂级数∑∞=--22)1(nnnnnx的和函数.12.求函数项级数∑∞=-1)2sin)1(nnnn xπ的和函数.。

相关主题