- 1、本文档共4页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
第六单元
1求定积分的近似值
f=@(x)sqrt(cos(x.^2)+1);
%下面三种方法选择一种即可
%基于全局自适应计分方法
I1=integral(f,0,2*pi)
I1=6.0401
%辛普森方法
I2=quad(f,0,2*pi)
I2=6.0401
%Guss-Lobatto
I3=quadl(f,0,2*pi)
I3=6.0401
%梯形积分法
clear
f=@(x)sqrt(cos(x.^2)+1);
a=0;
b=2*pi;
n=1000;
h=(b-a)/n;
x=linspace(a,b,n+1);
y=f(x);
I_trapz=h*(sum(y)-0.5*(y(1)+y(end)))
I_trapz=6.0401
2在给定的初值x0=1,y0=1,z0=1下,求方程组的数值解
%方法一,定义独立的函数,下面这段代码需要一函数的形式存在单独的函数文件中
functionF=F0502(X)
x=X(1);
y=X(2);
z=X(3);
F(1)=sin(x)+y^2+log(z)-7;
1
F(2)=3*x+2^y-z^3+1;
F(3)=x+y+z-5;
end
xsolve=fsolve(@M0502,[1,1,1],optimset(Display,off))
xsolve=1×3
0.59912.39592.0050
%方法二
f=@(x)[sin(x(1))+x(2)^2+log(x(3))-7,3*x(1)+2^x(2)-x(3)^3+1,x(1)+x(2)+x(3)-5];
xsolve1=fsolve(f,[1,1,1],optimset(Display,off))
xsolve1=1×3
0.59912.39592.0050
%方法三
clear
equations=@(vars)[
sin(vars(1))+vars(2)^2+log(vars(3))-7;%第一个方程
3*vars(1)+2^vars(2)-vars(3)^3+1;%第二个方程
vars(1)+vars(2)+vars(3)-5%第三个方程
];
%初始猜测
x0=[1,1,1];%初值[x0,y0,z0]
%求解方程组
options=optimoptions(fsolve,Display,off);%不显示优化过程
solution=fsolve(equations,x0,options);
%输出结果
disp([x=,num2str(solution(1))]);
x=0.59905
disp([y=,num2str(solution(2))]);
y=2.3959
disp([z=,num2str(solution(3))]);
z=2.005
3求函数在(0,1)内的最小值点和最小值。
2
%方法一
f=inline((x^3+cos(x)+x*log(x))/exp(x))
f=
内联函数:
f(x)=(x^3+cos(x)+x*log(x))/exp(x)
[x,fmin]=fminbnd(f,0,1)
x=0.5223
fmin=0.3974
%方法二
f=@(x)(x^3+cos(x)+x*log(x))/exp(x)
f=包含以下值的function_handle:
@(x)(x^3+cos(x)+x*log(x))/exp(x)
[x,fmin]=fminbnd(f,0,1)
x=0.5223
fmin=0.3974
第四题
%方法1
symsxy;
y=dsolve(D2y-3*Dy+2*y-1=0,y(0)=1,Dy(0)=0)
警告:Suorto
文档评论(0)