一个通用数值积分的matlab实现(带附录)(3)

2019-08-03 14:14

输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,10,'mjx') 按回车键显示结果如下: f=0.7856 n =10

(2)100等份(n=100) 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'mjx') 按回车键显示结果: f =0.7854 n =100

(3)1000等份(n=1000) 输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,1000,'mjx') 按回车键显示结果如下: f =0.7854 n =1000

3.2.2梯形公式调用实例(计算不同等分份数的积分值)

(一)10等份 输入函数变量:

>>f=inline('1/(1+x*x)');

7

>>[f,n]=tyszjf(f,0,1,10,'tx') 按回车键显示结果如下: f =0.7850 n =10

(二)100等份 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'tx') 按回车键显示结果如下: f =0.7854 n =100

(三)1000等份 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,1000,'tx') 按回车键显示结果如下: f =0.7854 n=1000

3.2.3辛普森公式调用实例

(一)10等份 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,10,'sps') 按回车键显示结果如下:

8

f =0.7854 n =10

(二)100等份 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'sps') 按回车键显示结果如下: f =0.7854 n =100

(三)1000等份 输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,1000,'sps') 按回车键显示结果如下: f =0.7854 n =1000

(四)精度要求不同时

(1)wc=0.5*10^(-7) 输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,n,'sps',0.5*10^(-7)) 按回车键显示结果如下: f =0.7854 n =100

(2)wc=0.5*10^(-8)

9

输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,n,'sps',0.5*10^(-8)) 按回车键显示结果如下: f =0.7854 n =100

3.2.4牛顿-柯茨公式调用实例

(一)10等份 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,10,'nc') 按回车键显示结果如下: f =0.7854 n =10

(二)100等份 输入函数变量:

>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'nc') 按回车键显示结果如下: f =0.7854 n =100

(三)1000等份 输入函数变量:

>>f=inline('1/(1+x*x)');

10

>>[f,n]=tyszjf(f,0,1,1000,'nc') 按回车键显示结果如下: f =0.7854 n =1000

3.2.5自适应积分调用实例

(一)wc=0.5*10^(-6) 输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,n,'zsy',0.5*10^(-6)) 按回车键显示结果如下: f=0.7854 n=8

(二)wc=0.5*10^(-7) 输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,n,'zsy',0.5*10^(-7)) 按回车键显示结果如下: f =0.7854 n =8

(三)wc=0.5*10^(-5) 输入函数变量:

>>f=inline('1/(1+x*x)');

>>[f,n]=tyszjf(f,0,1,n,'zsy',0.5*10^(-5)) 按回车键显示结果如下:

11


一个通用数值积分的matlab实现(带附录)(3).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:PLC期末复习题

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: