应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » 【求助】求积分程序中被积函数问题
161/2返回列表12下一页
查看: 1921  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

muqiuyan

木虫 (小有名气)

[交流] 【求助】求积分程序中被积函数问题已有5人参与

我用Simpson方法计算定积分。在程序中有两个地方用到积分，而被积分函数作为子程序调用，在计算积分的程序中该怎么定义被积函数呢？

这是 辛普森算法 的程序：
double fsimp(a,b,eps) /*辛普森算法*/  
double a,b,eps; /*a为积分下限，b为积分上限，eps是希望达到的精度*/  
{  
    int n,k;  
    double h,t1,t2,s1,s2,ep,p,x;  
    n=1; h=b-a;  
    t1=h*(fsimpf(a)+fsimpf(b))/2.0;  /*用梯形公式求出一个大概的估值*/  
    s1=t1;  
    printf("s1=%e\n",s1);
        ep=eps+1.0;  
    while (ep>=eps)  
{  
p=0.0;  
for (k=0;k<=n-1;k++)  
{  
x=a+(k+0.5)*h;  
p=p+fsimpf(x);  
}  
t2=(t1+h*p)/2.0;  
s2=(4.0*t2-t1)/3.0;  
printf("s2=%e\n",s2);
ep=fabs(s2-s1);  
t1=t2; s1=s2; n=n+n; h=h/2.0;  
}  
    return(s2);  
}



这是被积函数的程序：

double fsimpf(r)
double r; /* 被积分的函数*/

{  double y;
    y=被积函数f(x)的表达式；
  return(y);
}


两个被积函数不能都用fsimpf（）表示吧，可不用它，计算积分时调用被积函数又是用fsimpf（）表示的：


怎么区分啊
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-05-21 14:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

simthhuang

木虫 (小有名气)
★ ★
muqiuyan(金币+1):谢谢参与
resonant(金币+1):感谢积极参与：-） 2010-05-21 17:38:39
哪里有两个被积函数啊？明明就一个嘛！
赞一下(2人) 回复此楼
一贯知足二目远眺三餐有节四季不懒五谷皆食六欲不张七分忍让八方交往九薄烟断十分坦荡
2楼2010-05-21 17:04:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muqiuyan

木虫 (小有名气)
引用回帖:
Originally posted by simthhuang at 2010-05-21 17:04:29:
哪里有两个被积函数啊？明明就一个嘛！

主程序里有两个被积函数，我只列了两个子程序，一个是计算积分的，一个是被积函数的子程序，想问问第二个被积函数怎么弄，积分子程序里的f(a) f(b) f(x) 之类的需要调要被积函数子程序的地方要不要调整。
赞一下 回复此楼
3楼2010-05-21 21:46:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

simthhuang

木虫 (小有名气)
★
resonant(金币+1):感谢参与：-） 2010-05-22 21:01:25
第二个被积函数就是被积函数啊，就是积分号和dx之间那部分啊。f(a) f(b) f(x) 分别是用初值a,终值b和中间迭代的数x调用被积函数。MM，我的理解对吗？
赞一下(1人) 回复此楼
一贯知足二目远眺三餐有节四季不懒五谷皆食六欲不张七分忍让八方交往九薄烟断十分坦荡
4楼2010-05-22 19:41:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muqiuyan

木虫 (小有名气)
resonant:调用函数没有什么特别的吧？区分只要函数名不同不就可以区分么？ 2010-05-24 11:09:30
引用回帖:
Originally posted by simthhuang at 2010-05-22 19:41:04:
第二个被积函数就是被积函数啊，就是积分号和dx之间那部分啊。f(a) f(b) f(x) 分别是用初值a,终值b和中间迭代的数x调用被积函数。MM，我的理解对吗？

可能是我表达的不清楚。我要弄一个很复杂的程序，有很多公式，其中有两个计算积分的式子，所以我编了一个计算积分的子程序，就是fsimp(a,b,eps)；这其中就涉及到调用被积函数，那我的两个被积函数怎么调用呢，怎么区分开呢

PS：你是学计算机的吗，能认识下吗？
赞一下(1人) 回复此楼
5楼2010-05-23 10:37:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archun

铜虫 (初入文坛)
★ ★
muqiuyan(金币+1):谢谢参与
jjdg(金币+1):感谢热心应助 2010-05-26 04:51:38
muqiuyan(金币+2): 2010-05-26 09:21:19
我觉得你是这个意思吧。
你有两个函数要算积分，比如为 fsimpf1(x) 和 fsimpf2(x)
但是算积分的函数只有一个，就是 fsimp()
在 fsimp() 中有一句语句是 p=p+fsimpf(x)，但是在 fsimpf(x) 的定义中只定义了一个函数，比如是 fsimpf1(x)。这样的话，你算第二个积分时，仍然要调用 fsimp() 函数算。但是 fsimp() 是调用的 fsimpf1(x) 不是 fsimpf2(x)。你就不知道 fsimpf2(x) 该怎么算了。

你可以增加一个参数解决问题。
比如：
double fsimp(double a, double b, double eps, int func_index)
{
    .....
    .....
    p=p+fsimpf(x,func_index)
    .....
    .....
}

double fsimpf(double x, int func_index)
{
    swtich (func_index)
    {
    case 1:
        return (3*x+5); //第一个函数
        break;
    case 2:
        return (x*x+5); //第二个函数
        break;
    ........
    ........
    }
}

你积分时就写为
fsimp(a,b,eps,1) 积分第一个函数
fsimp(a,b,eps,2) 积分第二个函数
赞一下(2人) 回复此楼
6楼2010-05-25 22:04:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muqiuyan

木虫 (小有名气)
引用回帖:
Originally posted by archun at 2010-05-25 22:04:35:
我觉得你是这个意思吧。
你有两个函数要算积分，比如为 fsimpf1(x) 和 fsimpf2(x)
但是算积分的函数只有一个，就是 fsimp()
在 fsimp() 中有一句语句是 p=p+fsimpf(x)，但是在 fsimpf(x) 的定义中只定义了一个 ...

太感谢了，
那积分程序中的原来的fsimpf(a)就写成fsimpf(a,func_index)，fsimpf(b)写成fsimpf(b,func_index)，是吗？
另外我的两个被积函数是很复杂的式子，在哪里计算它们呢？下面的程序对吗
double fsimpf(double x, int func_index)
{


    swtich (func_index)
    {
    case 1:

double y1;
    y1=被积函数 fsimpf1(x) 的表达式；
      return (y1); //第一个函数
        break;

    case 2:
double y2;
    y2=被积函数 fsimpf1(x) 的表达式；
        return (y2); //第二个函数
        break;
    ........
    ........
    }
}
赞一下 回复此楼
7楼2010-05-26 09:21:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archun

铜虫 (初入文坛)
★
余泽成(金币+1):参与应助，欢迎常来程序语言版！ 2010-05-28 08:57:25
muqiuyan(金币+1): 2010-06-01 16:30:55
对的。
其实那个break语句可以不用写。只不过写也不错。
赞一下(1人) 回复此楼
8楼2010-05-28 06:04:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

muqiuyan

木虫 (小有名气)
引用回帖:
Originally posted by archun at 2010-05-28 06:04:01:
对的。
其实那个break语句可以不用写。只不过写也不错。

我的积分变量也是两个，怎么区分呢？一个是半径r；一个是温度T。哪些地方需要更改？
赞一下 回复此楼
9楼2010-06-01 16:32:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★ ★
muqiuyan(金币+1):谢谢参与
resonant(金币+1):感谢应助：-） 2010-06-02 19:03:48
muqiuyan(金币+1): 2010-06-03 10:15:13
最近比较忙，没有怎么回帖。

简单地说一下吧，建议你使用函数指针，不然程序照你这样写以后又要修改积分的函数那你还得改。
double fsimp(double (*f)(double ),double a, double b, double eps) /  /*a为积分下限，b为积分上限，eps是希望达到的精度,f为函数指针*/
在fsimp 的实现中，你就写f(a),f(b) 就是了。
前面申明double fsimpf(double)
然后在调用 fsimp时就直接
xx=fsimp(fsimpf, double a, double b, double eps) ;就是了
建议你看看the c programming language
赞一下(2人) 回复此楼
10楼2010-06-02 18:44:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 muqiuyan 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭