matlab代码--数值微分.doc

下载文档 降价啦

42
0
约7.78千字
约 9页
2017-12-21 发布于江西
举报
版权申诉
保障服务

matlab代码--数值微分.doc

1、本文档共9页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

matlab代码--数值微分

数值微分 MidPoint中点公式求取导数 function df=MidPoint(func,x0,h) if nargin == 2 h = 0.1; else if (nargin == 3 h == 0.0) disp(h不能为0！); return; end end y1 = subs(sym(func), findsym(sym(func)),x0+h); y2 = subs(sym(func), findsym(sym(func)),x0-h); df = (y1-y2)/(2*h); ThreePoint三点法求函数的导数 function df=ThreePoint(func,x0,type,h) if nargin == 3 h = 0.1; else if (nargin == 4 h == 0.0) disp(h不能为0！); return; end end y0 = subs(sym(func), findsym(sym(func)),x0); y1 = subs(sym(func), findsym(sym(func)),x0+h); y2 = subs(sym(func), findsym(sym(func)),x0+2*h); y_1 = subs(sym(func), findsym(sym(func)),x0-h); y_2 = subs(sym(func), findsym(sym(func)),x0-2*h); switch type case 1, df = (-3*y0+4*y1-y2)/(2*h); %用第一个公式求导数 case 2, df = (3*y0-4*y_1+y_2)/(2*h); %用第二个公式求导数 case 3, df = (y1-y_1)/(2*h); %用第三个公式求导数 end FivePoint五点法求函数的导数 function df=FivePoint(func,x0,type,h) if nargin == 3 h = 0.1; else if (nargin == 4 h == 0.0) disp(h不能为0！); return; end end y0 = subs(sym(func), findsym(sym(func)),x0); y1 = subs(sym(func), findsym(sym(func)),x0+h); y2 = subs(sym(func), findsym(sym(func)),x0+2*h); y3 = subs(sym(func), findsym(sym(func)),x0+3*h); y4 = subs(sym(func), findsym(sym(func)),x0+4*h); y_1 = subs(sym(func), findsym(sym(func)),x0-h); y_2 = subs(sym(func), findsym(sym(func)),x0-2*h); y_3 = subs(sym(func), findsym(sym(func)),x0-3*h); y_4 = subs(sym(func), findsym(sym(func)),x0-4*h); switch type case 1, df = (-25*y0+48*y1-36*y2+16*y3-3*y4)/(12*h); %用第一个公式求导数 case 2, df = (-3*y_1-10*y0+18*y1-6*y2+y3)/(12*h); %用第二个公式求导数 case 3, df = (y_2-8*y_1+8*y1-y2)/(12*h); %用第三个公式求导数 case 4, df = (3*y1+10*y0-18*y_1+6*y_2-y_3)/(12*h); %用第四个公式求导数 case 5, df = (25*y0-48*y_1+36*y_2-16*y_3+3*y_4)/(12*h); %用第五个公式求导数 end Di