| 查看: 521 | 回复: 8 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请教,C语言中pow(x,y)函数对幂指数的类型有限制吗?
|
|||
|
如题,最近算个程序,当调用pow()函数时,当幂指数取0.0时输出正常,但是当幂指数取0.6时输出全是0,不知道为什么? 请问会不会是与C语言中pow()函数原代码中幂指数的类型有限制呢? 希望大家指点迷津,或者是因为别的问题,谢谢!!! [ Last edited by cadick on 2009-12-11 at 02:49 ] |
» 猜你喜欢
求国际会议网站
已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有278人回复
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
反铁磁体中的磁性切换:两种不同的机制已成功可视化
已经有0人回复
26申博推荐:南京航空航天大学国际前沿院光学方向招收博士生!
已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复

wgpig
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2504.6
- 帖子: 191
- 在线: 273.7小时
- 虫号: 35477
- 注册: 2004-02-03
- 性别: GG
- 专业: 生物大分子结构与功能
★ ★
mingdong(金币+1,VIP+0):感谢交流! 10-10 11:01
cll3380(金币+1,VIP+0):谢谢,不过我就是找不到代码哦! 10-12 10:41
mingdong(金币+1,VIP+0):感谢交流! 10-10 11:01
cll3380(金币+1,VIP+0):谢谢,不过我就是找不到代码哦! 10-12 10:41
| 无。贴代码看看。 |
2楼2009-10-09 22:32:41

3楼2009-10-12 10:48:51
wgpig
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2504.6
- 帖子: 191
- 在线: 273.7小时
- 虫号: 35477
- 注册: 2004-02-03
- 性别: GG
- 专业: 生物大分子结构与功能
4楼2009-10-12 16:18:33

5楼2009-10-12 22:10:49
6楼2009-10-13 12:10:34
wgpig
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2504.6
- 帖子: 191
- 在线: 273.7小时
- 虫号: 35477
- 注册: 2004-02-03
- 性别: GG
- 专业: 生物大分子结构与功能
★ ★ ★ ★
mingdong(金币+2,VIP+0):感谢交流! 10-13 17:31
cll3380(金币+2,VIP+0):问题解决,非常非常感谢!!! 10-14 09:32
mingdong(金币+2,VIP+0):感谢交流! 10-13 17:31
cll3380(金币+2,VIP+0):问题解决,非常非常感谢!!! 10-14 09:32
|
$ cat test.cpp #include #include #include #include #include #include "fstream.h" //变量及函数的预定义-------------------------------------------------------------------------- const int m=500,n=70000; const double pi=3.1415926; double xx,x[n+2],noise1[n+2]; double stren1=0.1,tau1=1.0; double a=0.6,h=0.001; double g1(double xx) { return pow(fabs(xx),a); } //产生0-1均匀分布的随机数--------------------------------------------------------------------- static double r[98]; static int iff,ix1,ix2,ix3; double ran(int &idum) { const int ia=16807,im=2147483647,iq=127773,ir=2836,ntab=32; const int ndiv=(1+(im-1)/ntab); const double eps=3.0e-16,am=1.0/im,rnmx=(1.0-eps); static int iy=0; static int iv[ntab]; int j,k; double temp; if (idum<=0||!iy) { if (-idum<1) idum=1; else idum=-idum; for(j=ntab+7;j>=0;j--) { k=idum/iq; idum=ia*(idum-k*iq)-ir*k; if(idum<0) idum+=im; if(j iy=iv[0]; } k=idum/iq; idum=ia*(idum-k*iq)-ir*k; if(idum<0) idum+=im; j=iy/ndiv; iy=iv[j]; iv[j]=idum; if ((temp=am*iy)>rnmx) return rnmx; else return temp; } //主程序-------------------------------------------------------------------------------------- ofstream xzuob("x1.txt" ,yzuob("y1.txt" ;int main() { double gw_noise1,gw_noisea,cum[n+2]; int i,j; int idum=-13; for(i=1;i<=n;i++) { cum=0.0; } for(j=1;j<=m;j++) { x[1]=0.1,noise1[1]=0.0; for(i=1;i<=n;i++) { double y1,y2; y1=ran(idum); y2=ran(idum); gw_noise1=sqrt(-4.0*log(y1))*cos(2.0*pi*y2); gw_noisea=sqrt(stren1*h)*gw_noise1; x[i+1]=x+g1(x)*noise1*h; noise1[i+1]=noise1-noise1/tau1*h+gw_noisea/tau1; cum=cum+x*x/m; } } for(i=1;i<=n;i++) { xzuob< yzuob< return 0; } |
7楼2009-10-13 13:50:11
wgpig
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2504.6
- 帖子: 191
- 在线: 273.7小时
- 虫号: 35477
- 注册: 2004-02-03
- 性别: GG
- 专业: 生物大分子结构与功能
8楼2009-10-13 13:53:09
wgpig
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2504.6
- 帖子: 191
- 在线: 273.7小时
- 虫号: 35477
- 注册: 2004-02-03
- 性别: GG
- 专业: 生物大分子结构与功能
9楼2009-10-13 13:55:53













回复此楼
,yzuob("y1.txt"