scilab中求积分
- CSDN博客
scilab中求积分
scilab中求积分
===========================================================
scilab中求一元积分函数：intg（a，b，f）
函数格式：[v,err]=intg(a,b,f [,ea [,er])
不过在scilab中求一元积分常用命令y＝intg（a，b，f）
这里的积分区间为（a，b）；
积分函数为 f；
deff('y=f(x);','y=x*x;'); //定义函数f（x）＝x*x
result=intg(1,2,f); //对函数f在（1，2）积分，
--&result //输入result，显示结果如下result =
求二元积分函数：int2d（）
函数格式：[I,err]=int2d(X,Y,f [,params])
I：积分误差估计值；
X：一个3&n的矩阵，包含n个三角形顶点的X坐标信息；
Y:一个3&n的矩阵，包含n个三角形顶点的x坐标信息；
f：被积函数；
params：几分参数为实数向量[tol，iclose，maxtri，mevals，iflag]，
X=[0,0;1,1;1,0];//X，Y为两个三角形的坐标Y=[0,0;0,1;1,1];deff('z=f(x,y)','z=cos(x+y)')[I,e]=int2d(X,Y,f)// computes the integrand over the square [0 1]x[0 1]
本文已收录于以下专栏：
相关文章推荐
Problem Description
Ignatius bought a land last week, but he didn't know the area of the land becaus...
运用Win32API实现
            题目地址：http://acm./showproblem.php?pid=1071
点击打开链接
            题目大意...
一.写一个用矩形法求定积分的函数，求sin(x)在（0,1）上的定积分。
#include 
#include 
float jifen(float a,float b)
{int i,l;
...
//复化梯形公式求定积分
//将[0,1]分成n等份，子区间长度为h=(b-a)/n
int main()
Time Limit: 1 Sec  Memory Limit:
Submit: 423  Solved: 269
[Submit][Status][Web
蒙特卡洛（Monte Carlo）法是一类随机算法的统称。随着二十世纪电子计算机的出现，蒙特卡洛法已经在诸多领域展现出了超强的能力。在机器学习和自然语言处理技术中，常常被用到的MCMC也是由此发展而来...
OpenMP 并行求数值积分， 临界区求数值积分圆周率（pi）方法
开启两个并行线程程序如下：
首先我们知道Matlab可以表示分段函数，一般都是是用if-else的结构，如下
x=-5:0.01:5;
for X=1:1:length(x);
他的最新文章
讲师：王禹华
讲师：宋宝华
您举报文章：
举报原因：
原文地址：
原因补充：
(最多只允许输入30个字)查看: 3149|回复: 14|关注: 0
已知积分值、积分上限求积分下限，大神看看我哪里错了
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
积分函数略复杂，为：(1-0.02*x)^4.176/(0.+0.02*x)^3.174*(0.886+0.0928*x-5.328e-3*x^2+1.1776e-4*x^3-8.96e-7*x^4)^3.174)
已知积分值为20492，积分上限为37.596，求积分下限。
我写的句子如下：
solve(int((1-0.02*x)^4.176/(0.+0.02*x)^3.174*(0.886+0.0928*x-5.328e-3*x^2+1.1776e-4*x^3-8.96e-7*x^4)^3.174),x,Q,37.596)-20492,Q)
即是以Q为变量作为下限，解出Q
也没报错，可为什么就是不出结果，程序一直busy。
想问问各位大神，我哪里不对，或者，有什么好的求下限的代码吗
<h1 style="color:# 麦片财富积分
本帖最后由 tzqgbt 于
11:43 编辑
可以换个算法试试比如说，直接用一个范围内的下限先按照较大的步长积分出结果，找到下限的粗略区间；然后在这个区间里按照较小的步长积分，确定一个更精确的区间，，，这样一次次地就可以找到足够精度的下限了
论坛优秀回答者
<h1 style="color:#17 麦片财富积分
关注者： 314
syms x
q=(1-0.02*x)^4.176/(0.+0.02*x)^3.174*(0.886+0.0928*x-5.328e-3*x^2+1.1776e-4*x^3-8.96e-7*x^4)^3.174);
q=matlabFunction(q);
k=@(Q)quadgk(q,Q,37.596)-20492
[x,fval]=fzero(k,[-6.6,-6])复制代码
<h1 style="color:# 麦片财富积分
第四行是有问题么，输出的k还是这串代码
论坛优秀回答者
<h1 style="color:#17 麦片财富积分
关注者： 314
k是匿名函数，输出k的话会把它的表达式输出来，也就还是那一串。。
结果应该没问题吧？
<h1 style="color:# 麦片财富积分
k是匿名函数，输出k的话会把它的表达式输出来，也就还是那一串。。
结果应该没问题吧？ ...
没太看懂。。想问我要求的下限值是哪个变量，输出的有x和fval。。。谢啦
论坛优秀回答者
<h1 style="color:#17 麦片财富积分
关注者： 314
是x, fval是求解误差
<h1 style="color:# 麦片财富积分
是x, fval是求解误差
为什么是在-6.6到6之间的，我预计这个值应该是个正数呢。想问一下大神是怎么预测到x会在-6.6到6之间的呢
论坛优秀回答者
<h1 style="color:#17 麦片财富积分
关注者： 314
我是作图找的，正数内也有解，你可以自己改变作图范围找：
syms x
q=(1-0.02*x)^4.176/(0.+0.02*x)^3.174*(0.886+0.0928*x-5.328e-3*x^2+1.1776e-4*x^3-8.96e-7*x^4)^3.174);
q=matlabFunction(q);
k=@(Q)quadgk(q,Q,37.596)-20492
x=linspace(77,78,0.5e2);& &%改变前两个参数作图
t=arrayfun(k,x);
plot(x,t);
% [x,fval]=fzero(k,[-6.6,-6])复制代码
<h1 style="color:# 麦片财富积分
我是作图找的，正数内也有解，你可以自己改变作图范围找：
我理解，目的不就是找出令t为0的x么，我自己找了个区间，但这里作出的图虽说是下降的趋势，但为什么纵坐标显示的一直是-2.0491e4呢
站长推荐 /1
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区}