噶米matlab16常用计算方法.docx

噶米matlab16常用计算方法.docx

ID:62504536

大小:144.83 KB

页数:13页

时间:2021-05-10

噶米matlab16常用计算方法.docx_第1页
噶米matlab16常用计算方法.docx_第2页
噶米matlab16常用计算方法.docx_第3页
噶米matlab16常用计算方法.docx_第4页
噶米matlab16常用计算方法.docx_第5页
资源描述:

《噶米matlab16常用计算方法.docx》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、常用计算方法1超越方程的求解一超越方程为x(2lnx-3)-100=0求超越方程的解。[算法]方法一:用迭代算法。将方程改为100x=2ln(冷)-3其中X0是一个初始值,由此计算终值值换成X0的值,重新进行计算;否则[程序]P1_1abs.m如下。%超越方程的迭代算法clearx0=30;xx=[];while1x=100/(2*log(x0)-3);xx=[xx,x];iflength(xx)>1000,break,endifabs(x0-x)<1e-4,break,endx0=x;endfigurex。取最大误差为e=10-4,当

2、x-x0

3、>e

4、时,就用x的

5、x-x°

6、

7、签text(length(xx),xx(end),num2str(xx(end)),[图示]用下标作为自变量画迭代的折线。要迭代19次才能达到精度,超越方程的解为'fontsize',fs)%显示结果如P0_20_1图所示,当最大误差为10-4时,需27.539。[算法]方法二:用求零函数和求解函数。将方程改为函数f(x)=21n(x)七-100xMATLAB求零函数为fzero,fzero函数的格式之一是x=fzero(f,x0)其中,f表示求解的函数文件,x0是估计值。fzero函数的格式之二是x=fzero(f,[x1,x2])其中,x1和x2表

8、示零点的范围。另外MATLAB还有求解函数solve,计算非线性方程和方程组的符号解。[程序]P1_2fzero.m如下。9超越方程的求法clearx=10:0.1:100;f=inline('2*log(x)-3-100./x')figureplot(x,f(x),'LineWidth',2)gridonx0=fzero(f,[20,30]);%x0=fzero(f,20);%holdon%清除变量%自变量向量%定义内线函数用的是字符窜%创建图形窗口%画曲线%加网格%求方程的零点求方程的零点%保持图像plot(x0,f(x0),'.')%画零点tit

9、le('超越方程的解','fontsize',16)%标题xlabel('itx','fontsize',16)%x标签ylabel('itf','fontsize',16)%y标签text(x0,0,num2str(x0),'fontsize',16)x0=solve('2*log(x)-3-100./x')plot(double(x0),0,'o')%标记零点%求超越方程的符号解%再画零点(double是把字符转化成数字)2.导数的计算正弦函数y=sinx的导数是余弦函数y'=cosx,余弦函数的导数是负的正弦函数,用MATLAB的数值导数和符

10、号导数求正弦函数的一阶和二阶导数,并与其解析解进行比较。[程序]P2diff.m如下。%正弦函数导数的计算方法cleardx=0.01*2*pi;x=0:dx:2*pi;y=sin(x);f1=diff(y)/dx;%清除变量%间隔%自变量向量%原函数%通过差分求导数f1=[f1(1),(f1(1:end-1)+f1(2:end))/2,f1(end)];%求平均值figure%创建图形窗口plot(x,cos(x),x,f1,'.')%plot(x,cos(x),x(1:end-1),f1,'.')%%plot(x,cos(x),x(2:end),f

11、1,'.')%symssxy=sin(sx);dy_dx=diff(y);df1=subs(dy_dx,sx,x);holdonplot(x,df1,'ro')gridon%画一阶导数和数值差分曲线数值导数(点)偏左数值导数(点)偏右%定义符号变量%建立符号函数%求符号导数%符号替换数值%保持图像%画符号导数曲线%加网格legend('解析导数’,‘数值差分’,‘符号导数’,4)%图例title('正弦函数的一阶导数','FontSize',16)%加标题f2=diff(f1)/dx;f2=[f2(1),(f2(1:end-1)+f2(2:end))/

12、2,f2(end)];d2y_dx2=diff(y,2);df2=subs(d2y_dx2,s

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

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

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