科学计算与MATLAB语言-第六单元.pdf

  1. 1、本文档共4页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 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)

***** + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档