matlab实验六1数值微积分与常微分方程求解

matlab实验六1数值微积分与常微分方程求解

ID:31631404

大小:87.35 KB

页数:8页

时间:2019-01-16

matlab实验六1数值微积分与常微分方程求解_第1页
matlab实验六1数值微积分与常微分方程求解_第2页
matlab实验六1数值微积分与常微分方程求解_第3页
matlab实验六1数值微积分与常微分方程求解_第4页
matlab实验六1数值微积分与常微分方程求解_第5页
资源描述:

《matlab实验六1数值微积分与常微分方程求解》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实验六数值微积分与常微分方程求解、实验目的1.掌握微分与积分的数值计算方法2.掌握常微分方程的数值求解方法二、实验内容和步骤1.求函数指定点的数值导数X2Xx3f(x)=12x3x2,x=1,2,3026x2.求下列函数的导数(1)y=x,0+10v+logvl(),求y(2)y=ln(l+x),求yX=13.用数值方法求定积分⑵仕sinxf2fy+2⑷Q如歹⑶打:

2、cos(x+y'>jdxdy(1)£x2a/2x2+3dx4.求微分方程的数值解(1)y+(1.2+sinl0f)y=0,/()=0心=

3、5,y(/())=1(2)(1+t2)y+2ty+3y=2心=0占=5』仏)=0,y(/0)=1三、程序清单及运行结果1.程序清单:forx=l:1:3fx=[x,x2,x3;1,2*x,3*x2;0,2,6*x];diff(fx)end运行结果:1.(1)程序清单:x二0:0.05:1;y=x.°10+10."x+l・/loglO(x);dy=diff(y)/0.05运行结果:»Vtrninc:Locofrero.>InC:ttAILAB6^5to<>lb<>xT.atl»l.

4、oclO-«只liu«17InC:WAlLAK^5vorkVntitle4eatlifte2YornincDividebyzero.CTyp*offIUTLAB4)vi3«ByZ«ro"suppressthisyarning)>InC:«AILAKc5vorkVntitled.»吐line2dy=Coluvms1through6-12.S321-18894・1.2023-O8919-07380-06907ColuiJiz7through12・0.7468-0.9278-1.2787-1.87

5、57-2.8437・43920Coluroiz13throuch18・68872-110152-18.1724-315S42-597974-1334123Coluvms19through20•43&3387Inf(2)程序清单:x=0:0.01:1;y=log(l+x);dy=diff(y,2)/0.01;dy(99)运行结果:»ans=-0.00251.(1)程序清单:fx二inline('x."2.*sqrt(2*x.*x+3)');quad(fx,1,5)运行结果:»ans-232.8057»(2

6、)程序清单:fx=inline('x./sin(x).2');quad(fx,pi/4,pi/3)运行结果:(1)程序清单:fx二inline('abs(cos(x+y))');dblquad(fx,0,pi,0,pi)运行结果:ans-6.2832三»l_(2)程序清单:symsxy;fx=x*y;int(int(fx,y2,y+2),-1,2);运行结果:»ans=45/8»1.(1)程序清单:建立函数文件funst.mofunctionyp=funst(t,y)yp=-(l.2+sin(10*t)

7、)*y;求解微分方程。t0=0;tf=5;yO=l;[t,y]=ode23(©funst,[tO,tf],yO)运行结果:t二1.22462.55013.864201.30822.61403.92860.06671.37142.67793.99200.13751.43472.74194.06710.20031.50012.81934.14080.26951.58422.90464.21440.35281.65302.96854.28000.43621.72193.03234.34320.50331.78

8、583.09594.41590.56631.85013.17214.48120.63691.93193.23644.54650.69131.99533.30074.60940.74572.05873.36424.68120.80812.12363.43954.75670.87382.20533.53284.83220.95912.27453.59654.89901.02772.34383.66024.96201.09632.40803.72385.00001.16002.47193.7998y=0.51

9、820.29360.18191.00000.51060.25950.15850.90350.49760.24690.13870.78220.46560.24210.12250.68230.42800.23820.11640.59840.38440.22900.11400.54020.33550.20540.11220.10800.03610.01220.00540.09740.03160.01190.00500.08640.02780.01

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。