- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学实验三作业答案
实验三作业参考答案
1. 验证罗尔定理对函数在区间[,]上的正确性.
syms x
daoshu=diff(log(sin(x)))
daoshu=simple(daoshu)
g=inline(1/tan(x));
c=fzero(g,[pi/6,5*pi/6])
或
syms x
daoshu=simple(diff(log(sin(x))))
c=fzero(1/tan(x),[pi/6,5*pi/6])
输出:c = 1.5708
2. 验证拉格朗日定理对函数在区间[0.2, 1]上的正确性.
即验证存在 使得 ,
也即在区间[0.2, 1]找的零点.
syms x
f=inline(4*x^3-5*x^2+x-2)
k=(f(1)- f(0.2))/ 0.8
daoshu=diff(‘4*x^3-5*x^2+x-2’)
输出:k = -0.0400
daoshu =12*x^2-10*x+1
c=fzero(12*x^2-10*x+1+0.04 , [0.2, 1])
输出: c = 0.7115
3. 验证对函数应用拉格朗日中值定理时, 所求得的点总是位于区间的正中间.
即验证 .
解法1:
syms x p q r a b
daoshu=diff(p*x^2+q*x+r,x)
输出:
daoshu =2*p*x+q
x=(a+b)/2;
daoshuzhi=eval(2*p*x+q)
x=a;
ya=eval(p*x^2+q*x+r)
x=b;
yb=eval(p*x^2+q*x+r)
c=simple(daoshuzhi-(yb-ya)/(b-a))
解法2:
syms x p q r a b
daoshu=diff(p*x^2+q*x+r, x)
输出:daoshu =
2*p*x+q
在编辑窗口输入以下命令函数并保存为parabolic_fun.m
% M function parabolic_fun.m c=[p,q,r];d=[a,b];
function y=parabolic_fun(c,d, x)
y=2*c(1)*x+c(2)-(c(1)*d(2)^2+c(2)*d(2)+c(3)-(c(1)*d(1)^2+c(2)*d(1)+c(3)))/(d(2)-d(1));
在命令窗口输入
fzero(parabolic_fun([2,2,2], [2,4],x), [2,4]) % c=[p,q,r]= [2,2,2];d=[a,b]= [2,4]
ans =
3
再在命令窗口输入
fzero(parabolic_fun([4,2,5],[0,1],x), [0,1]) % c=[p,q,r]= [4,2,5];d=[a,b]= [0,1]
ans =
0.5000
注:c=[p,q,r]与d=[a,b]的值可以任意调整。
4. 验证柯西中值定理对函数及在区间[0, ]上的正确性.
即验证 , 而,故验证 .
syms x
dF=diff(x+cos(x))
df=diff(sin(x))
输出:dF =
1-sin(x)
df =
cos(x)
c=fzero((1-sin(x))/cos(x)-(pi/2+cos(pi/2)-1)/sin(pi/2),[0,pi/2])
输出:c =
0.5335
5. 求下列函数的一、二阶导数:
(1) ; syms x
diff(f(x^2),x)
ans =
2*D(f)(x^2)*x
(2) ;
syms x
diff(f(x)^2,x)
ans =
2*f(x)*diff(f(x),x)
(3);
diff(log(f(x)),x)
ans =
diff(f(x),x)/f(x)
(4)+.
diff(f(exp(x))+exp(f(x)),x)
ans =
D(f)(exp(x))*exp(x)+diff(f(x),x)*exp(f(x))
6. 求下列函数的高阶导数:
(1), 求;
diff(x*sinh(x),100)
ans =
100*cosh(x)+x*sinh(x)
(2) , 求;
diff(x^2*cos(x),10)
ans =
90*cos(x)-20*x*sin(x)-x^2*cos(x)
(3), 求.
diff(x^2*sin(2*x),50)
ans =
689613692941107200*sin(2*x)+56294995342131200*x*cos(2*x)-1125899906842624*x^2*sin(2*x)
好复杂,求吧.
diff(x^2
文档评论(0)