| 查看: 1036 | 回复: 5 | ||
破碎东瓜金虫 (初入文坛)
|
[求助]
谁能帮我把一段matlab程序转为fortran啊！！！（很短的一段）
|
|
while err>errlim s=find(P>0); sn=find(P<=0); dd=real(ifft2(fA.*fft2(P,n*2,n*2))); u=dd(1:n,1:n)*2.0*dx/pi^2; rk=h+u; do=mean(rk(s)); rk=rk-do; G=sum(sum(rk(s).*rk(s))); pk(s)=rk(s)+G/Gold*pk(s); pk(sn)=0; Gold=G; dd=real(ifft2(fA.*fft2(pk,2*n,2*n))); qk=dd(1:n,1:n)*2.0*dx/pi^2; rb=mean(qk(s)); qk=qk-rb; dp=sum(sum(rk(s).*pk(s)))/sum(sum(qk(s).*pk(s))); P(s)=P(s)-dp*pk(s); s=find(P<=0);P(s)=0; sol=find((P==0)&(rk<0)); P(sol)=P(sol)-dp*rk(sol); W=sum(sum(P))*dx*dx; P=charge/W*P; err=sqrt(Gold*dx*dx); end |
回复此楼» 猜你喜欢
基金申报
已经有5人回复
-
基金委咋了？2026年的指南还没有出来？
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
-
常年博士招收(双一流，工科)
已经有4人回复
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
2楼2015-12-12 09:57:31
破碎东瓜
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1520.2
- 散金: 53
- 帖子: 46
- 在线: 32.1小时
- 虫号: 2231948
- 注册: 2013-01-09
- 性别: GG
- 专业: 机械摩擦学与表面技术
3楼2015-12-12 15:36:58
wxws.2008
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 5015.5
- 散金: 256
- 红花: 9
- 帖子: 354
- 在线: 147.9小时
- 虫号: 639701
- 注册: 2008-10-29
- 性别: GG
- 专业: 荒漠化与水土保持
4楼2015-12-12 18:07:54
lishouyin
铁杆木虫 (著名写手)
- 应助: 25 (小学生)
- 金币: 12773.3
- 散金: 614
- 红花: 10
- 沙发: 1
- 帖子: 2824
- 在线: 793.6小时
- 虫号: 1983857
- 注册: 2012-09-07
- 性别: GG
- 专业: 电路与系统
5楼2015-12-13 10:32:35
破碎东瓜
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1520.2
- 散金: 53
- 帖子: 46
- 在线: 32.1小时
- 虫号: 2231948
- 注册: 2013-01-09
- 性别: GG
- 专业: 机械摩擦学与表面技术
6楼2015-12-13 14:49:55